P Kerfriden, S Claus, I Mihai, A mixed-dimensional CutFEM methodology for the simulation of fibre-reinforced composites, Advanced Modeling and Simulation in Engineering Science, 2020
We develop a novel unfitted finite element solver for composite materials with quasi-1D fibrous reinforcements. The method belongs to the class of mixed-dimensional non-conforming finite element solvers. The fibres are treated as 1D structural elements that may intersect the mesh of the embedding structure in an arbitrary manner. No meshing of the unidimensional elements is required. Instead, fibre solution fields are described using the trace of the background mesh. A regularised “cut” finite element formulation is carefully designed to ensure that analyses using such non-conforming finite element descriptions are stable. We also design a dedicated primal/dual operator splitting scheme to resolve the coupling between structure and fibrous reinforcements efficiently. The novel computational strategy is applied to the solution of stiff computational models whereby fibrous reinforcements may lose their bond to the embedding material above a certain level of stress. It is shown that the primal-dual 1D/3D CutFEM scheme is convergent and well-behaved in variety of scenarios involving such highly nonlinear structural computations.
Simulations by S. Claus, S. Bigot and P. Kerfriden in the FEniCS library CutFEM.
S. Claus, S. Bigot and P. Kerfriden, CutFEM Method for Stefan--Signorini Problems with Application in Pulsed Laser Ablation, SIAM J. Sci. Comput., 40(5), 2018
Funding: Sêr Cymru National Research Network
Eddy is held by the tails and gravity acts in direction [1 -1] in the plane of the picture. The dinosaur will either slip or stumble forward depending on the roughness of the contact between its feet and the support.![]() |
| Harmonic Dirichlet boundary conditions |