Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userWeak Boundary Condition Enforcement for Linear Kirchhoff-Love Shells: Formulation, Error Analysis, and Verification
68
0
0.0
(
0
)
Added by
Rasmus Tamstorf
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Stable and accurate modeling of thin shells requires proper enforcement of all types of boundary conditions. Unfortunately, for Kirchhoff-Love shells, strong enforcement of Dirichlet boundary conditions is difficult because both functional and derivative boundary conditions must be applied. A popular alternative is to employ Nitsches method to weakly enforce all boundary conditions. However, while many Nitsche-based formulations have been proposed in the literature, they lack comprehensive error analyses and verifications. In fact, existing formulations are variationally inconsistent and yield sub-optimal convergence rates when used with common boundary condition specifications. In this paper, we present a novel Nitsche-based formulation for the linear Kirchhoff-Love shell that is provably stable and optimally convergent for general sets of admissible boundary conditions. To arrive at our formulation, we first present a framework for constructing Nitsches method for any abstract variational constrained minimization problem. We then apply this framework to the linear Kirchhoff-Love shell and, for the particular case of NURBS-based isogeometric analysis, we prove that the resulting formulation yields optimal convergence rates in both the shell energy norm and the standard $L^2$-norm. In the process, we derive the Euler-Lagrange equations for general sets of admissible boundary conditions and show that the Euler-Lagrange boundary conditions typically presented in the literature is incorrect. We verify our formulation by manufacturing solutions for a new shell obstacle course that encompasses flat, parabolic, hyperbolic, and elliptic geometric configurations. These manufactured solutions allow us to robustly measure the error across the entire shell in contrast with current best practices where displacement and stress errors are only measured at specific locations.
rate research
Read More
We present a comprehensive rotation-free Kirchhoff-Love (KL) shell formulation for peridynamics (PD) that is capable of modeling large elasto-plastic deformations and fracture in thin-walled structures. To remove the need for a predefined global parametric domain, Principal Component Analysis is employed in a meshfree setting to develop a local parameterization of the shell midsurface. The KL shell kinematics is utilized to develop a correspondence-based PD formulation. A bond-stabilization technique is employed to naturally achieve stability of the discrete solution. Only the mid-surface velocity degrees of freedom are used in the governing thin-shell equations. 3D rate-form material models are employed to enable simulating a wide range of material behavior. A bond-associative damage correspondence modeling approach is adopted to use classical failure criteria at the bond level, which readily enables the simulation of brittle and ductile fracture. NAT{Discretizing the model with asymptotically compatible meshfree approximation provides a scheme which converges to the classical KL shell model while providing an accurate and flexible framework for treating fracture.} A wide range of numerical examples, ranging from elastostatics to problems involving plasticity, fracture, and fragmentation, are conducted to validate the accuracy, convergence, and robustness of the developed PD thin-shell formulation. It is also worth noting that the present method naturally enables the discretization of a shell theory requiring higher-order smoothness on a completely unstructured surface mesh.
An isogeometric Galerkin approach for analysing the free vibrations of piezoelectric shells is presented. The shell kinematics is specialised to infinitesimal deformations and follow the Kirchhoff-Love hypothesis. Both the geometry and physical fields are discretised using Catmull-Clark subdivision bases. It provides the required C1 continuous discretisation for the Kirchhoff-Love theory. The crystalline structure of piezoelectric materials is described using an anisotropic constitutive relation. Hamiltons variational principle is applied to the dynamic analysis to derive the weak form of the governing equations. The coupled eigenvalue problem is formulated by considering the problem of harmonic vibration in the absence of external load. The formulation for the purely elastic case is verified using a spherical thin shell benchmark. Thereafter, the piezoelectric effect and vibration modes of a transverse isotropic curved plate are analysed and evaluated for the Scordelis-Lo roof problem. Finally, the eigenvalue analysis of a CAD model of a piezoelectric speaker shell structure showcases the ability of the proposed method to handle complex geometries.
We develop and analyze an ultraweak variational formulation for a variant of the Kirchhoff-Love plate bending model. Based on this formulation, we introduce a discretization of the discontinuous Petrov-Galerkin type with optimal test functions (DPG). We prove well-posedness of the ultraweak formulation and quasi-optimal convergence of the DPG scheme. The variational formulation and its analysis require tools that control traces and jumps in $H^2$ (standard Sobolev space of scalar functions) and $H(mathrm{div,Div})$ (symmetric tensor functions with $L_2$-components whose twice iterated divergence is in $L_2$), and their dualities. These tools are developed in two and three spatial dimensions. One specific result concerns localized traces in a dense subspace of $H(mathrm{div,Div})$. They are essential to construct basis functions for an approximation of $H(mathrm{div,Div})$. To illustrate the theory we construct basis functions of the lowest order and perform numerical experiments for a smooth and a singular model solution. They confirm the expected convergence behavior of the DPG method both for uniform and adaptively refined meshes.
A hybrid surface integral equation partial differential equation (SIE-PDE) formulation without the boundary condition requirement is proposed to solve the electromagnetic problems. In the proposed formulation, the computational domain is decomposed into two emph{overlapping} domains: the SIE and PDE domains. In the SIE domain, complex structures with piecewise homogeneous media, e.g., highly conductive media, are included. An equivalent model for those structures is constructed through replacing them by the background medium and introducing a surface equivalent electric current density on an enclosed boundary to represent their electromagnetic effects. The remaining computational domain and homogeneous background medium replaced domain consist of the PDE domain, in which inhomogeneous or non-isotropic media are included. Through combining the surface equivalent electric current density and the inhomogeneous Helmholtz equation, a hybrid SIE-PDE formulation is derived. Unlike other hybrid formulations, where the transmission condition is usually used, no boundary conditions are required in the proposed SIE-PDE formulation, and it is mathematically equivalent to the original physical model. Through careful construction of basis functions to expand electric fields and the equivalent current density, the discretized formulation is compatible on the interface of the SIE and PDE domain. Finally, its accuracy and efficiency are validated through two numerical examples. Results show that the proposed SIE-PDE formulation can obtain accurate results including both near and far fields, and significant performance improvements in terms of CPU time and memory consumption compared with the FEM are achieved.
We investigate the piecewise linear nonconforming Crouzeix-Raviar and the lowest order Raviart-Thomas finite-element methods for the Poisson problem on three-dimensional anisotropic meshes. We first give error estimates of the Crouzeix-Raviart and the Raviart-Thomas finite-element approximate problems. We next present the equivalence between the Raviart-Thomas finite-element method and the enriched Crouzeix-Raviart finite-element method. We emphasise that we do not impose either shape-regular or maximum-angle condition during mesh partitioning. Numerical results confirm the results that we obtained.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
