Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Auto-Generated Geometry-Based Discrete Finite Element Model for Damage Evolution in Composite Laminates with Arbitrary Stacking Sequence
204
0
0.0
(
0
)
Added by
Jiakun Liu
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Stiffness degradation and progressive failure of composite laminates are complex processes involving evolution and multi-mode interactions among fiber fractures, intra-ply matrix cracks and inter-ply delaminations. This paper presents a novel finite element model capable of explicitly treating such discrete failures in laminates of random layup. Matching of nodes is guaranteed at potential crack bifurcations to ensure correct displacement jumps near crack tips and explicit load transfer among cracks. The model is entirely geometry-based (no mesh prerequisite) with distinct segments assembled together using surface-based tie constraints, and thus requires no element partitioning or enrichment. Several numerical examples are included to demonstrate the models ability to generate results that are in qualitative and quantitative agreement with experimental observations on both damage evolution and tensile strength of specimens. The present model is believed unique in realizing simultaneous and accurate coupling of all three types of failures in laminates having arbitrary ply angles and layup.
rate research
Read More
Notched components are commonly used in engineering structures, where stress concentration may easily lead to crack initiation and development. The main goal of this work is to develop a simple numerical method to predict the structural strength and crack-growth-path of U-notched specimens made of brittle materials. For this purpose, the Fragile Points Method (FPM), as previously proposed by the authors, has been augmented by an interface damage model at the interfaces of the FPM domains, to simulate crack initiation and development. The formulations of FPM are based on a discontinuous Galerkin weak form where point-based piece-wise-continuous polynomial test and trial functions are used instead of element-based basis functions. In this work, the numerical fluxes introduced across interior interfaces between subdomains are postulated as the tractions acting on the interface derived from an interface damage model. The interface damage is triggered when the numerical flux reaches the interface strength, and the process of crack-surface separation is governed by the fracture energy. In this way, arbitrary crack initiation and propagation can be naturally simulated without the need for knowing the fracture-patch before-hand. Additionally, a small penalty parameter is sufficient to enforce the weak-form continuity condition before damage initiation, without causing problems such as artificial compliance and numerical ill-conditioning. As validations, the proposed FPM method with the interface damage model is used to predict the structural strength and crack-development from U-notched structures made of brittle materials, which is useful but challenging in engineering structural design practices.
The Lorentz equations describe the motion of electrically charged particles in electric and magnetic fields and are used widely in plasma physics. The most popular numerical algorithm for solving them is the Boris method, a variant of the Stormer-Verlet algorithm. Boris method is phase space volume conserving and simulated particles typically remain near the correct trajectory. However, it is only second order accurate. Therefore, in scenarios where it is not enough to know that a particle stays on the right trajectory but one needs to know where on the trajectory the particle is at a given time, Boris method requires very small time steps to deliver accurate phase information, making it computationally expensive. We derive an improved version of the high-order Boris spectral deferred correction algorithm (Boris-SDC) by adopting a convergence acceleration strategy for second order problems based on the Generalised Minimum Residual (GMRES) method. Our new algorithm is easy to implement as it still relies on the standard Boris method. Like Boris-SDC it can deliver arbitrary order of accuracy through simple changes of runtime parameter but possesses better long-term energy stability. We demonstrate for two examples, a magnetic mirror trap and the Solevev equilibrium, that the new method can deliver better accuracy at lower computational cost compared to the standard Boris method. While our examples are motivated by tracking ions in the magnetic field of a nuclear fusion reactor, the introduced algorithm can potentially deliver similar improvements in efficiency for other applications.
We present a 3D hybrid method which combines the Finite Element Method (FEM) and the Spectral Boundary Integral method (SBIM) to model nonlinear problems in unbounded domains. The flexibility of FEM is used to model the complex, heterogeneous, and nonlinear part -- such as the dynamic rupture along a fault with near fault plasticity -- and the high accuracy and computational efficiency of SBIM is used to simulate the exterior half spaces perfectly truncating all incident waves. The exact truncation allows us to greatly reduce the domain of spatial discretization compared to a traditional FEM approach, leading to considerable savings in computational cost and memory requirements. The coupling of FEM and SBIM is achieved by the exchange of traction and displacement boundary conditions at the computationally defined boundary. The method is suited to implementation on massively parallel computers. We validate the developed method by means of a benchmark problem. Three more complex examples with a low velocity fault zone, low velocity off-fault inclusion, and interaction of multiple faults, respectively, demonstrate the capability of the hybrid scheme in solving problems of very large sizes. Finally, we discuss potential applications of the hybrid method for problems in geophysics and engineering.
In this paper, we propose a local-global multiscale method for highly heterogeneous stochastic groundwater flow problems under the framework of reduced basis method and the generalized multiscale finite element method (GMsFEM). Due to incomplete characterization of the medium properties of the groundwater flow problems, random variables are used to parameterize the uncertainty. As a result, solving the problem repeatedly is required to obtain statistical quantities. Besides, the medium properties are usually highly heterogeneous, which will result in a large linear system that needs to be solved. Therefore, it is intrinsically inevitable to seek a computational-efficient model reduction method to overcome the difficulty. We will explore the combination of the reduced basis method and the GMsFEM. In particular, we will use residual-driven basis functions, which are key ingredients in GMsFEM. This local-global multiscale method is more efficient than applying the GMsFEM or reduced basis method individually. We first construct parameter-independent multiscale basis functions that include both local and global information of the permeability fields, and then use these basis functions to construct several global snapshots and global basis functions for fast online computation with different parameter inputs. We provide rigorous analysis of the proposed method and extensive numerical examples to demonstrate the accuracy and efficiency of the local-global multiscale method.
In this work, we develop a discretisation method for the mixed formulation of the magnetostatic problem supporting arbitrary orders and polyhedral meshes. The method is based on a global discrete de Rham (DDR) sequence, obtained by patching the local spaces constructed in [Di Pietro, Droniou, Rapetti, Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra, arXiv:1911.03616] by enforcing the single-valuedness of the components attached to the boundary of each element. The first main contribution of this paper is a proof of exactness relations for this global DDR sequence, obtained leveraging the exactness of the corresponding local sequence and a topological assembly of the mesh valid for domains that do not enclose any void. The second main contribution is the formulation and well-posedness analysis of the method, which includes the proof of uniform Poincare inequalities for the discrete divergence and curl operators. The convergence rate in the natural energy norm is numerically evaluated on standard and polyhedral meshes. When the DDR sequence of degree $kge 0$ is used, the error converges as $h^{k+1}$, with $h$ denoting the meshsize.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
