Solid-shell element for modeling consolidation-induced wrinkle defect formation
Maarten M. Hodzelmans, Stephen R. Hallett, Joris J.C. Remmers, Jonathan P.-H. Belnoue
Abstract. When consolidating thick fiber-reinforced composite laminates through compaction, the reduction in thickness may cause wrinkle defects to form. These defects considerably reduce the structural properties of components. The lack of robustness of this manufacturing process is one of the reasons why, in the aerospace industry, composite structures are designed using large safety factors. A numerical method to predict these defects was proposed in the literature [1], wherein a ply is modeled as having a homogeneous constitutive response. Multiple plies and the resin rich regions between them are then homogenized through kinematic enrichment. The authors observed that an overly coarse mesh in the through-thickness direction did not converge. In this work, their method is combined with a solid-shell element formulation [2] to allow the use of a single element through the thickness, aiming to increase the speed of the simulations. The initial approach makes use of an implicit finite element solver. This works well for compaction, as the timescale involved is typically in the order of hours. However, it is impractical for forming, as this process occurs in second(s). Meanwhile, the aerospace industry is striving for higher throughput, which could be enabled by multi-ply forming. This scheme could aid in finding adequate deposition conditions. For this reason, as well as to further reduce the computational cost of simulations, the element and constitutive behavior were implemented in the explicit in-house finite element solver SimTex with three alterations. The adaptive hourglassing-stabilization stiffness was altered to work with the anisotropic material, an iterative step was removed from the kinematic enrichment scheme, and the compaction limit was modelled as an elastic- rather than viscous phenomenon to bound the apparent stiffness. Although the predictive capabilities of the resulting scheme are limited by its dependence on the timestep, wrinkling can be predicted correctly with the solid-shell element.
Keywords
Solid Shell Element, Kinematic Enrichment, Hyperviscoelasticity, Consolidation
Published online 5/7/2025, 10 pages
Copyright © 2025 by the author(s)
Published under license by Materials Research Forum LLC., Millersville PA, USA
Citation: Maarten M. Hodzelmans, Stephen R. Hallett, Joris J.C. Remmers, Jonathan P.-H. Belnoue, Solid-shell element for modeling consolidation-induced wrinkle defect formation, Materials Research Proceedings, Vol. 54, pp 1935-1944, 2025
DOI: https://doi.org/10.21741/9781644903599-208
The article was published as article 208 of the book Material Forming
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
References
[1] J.P.H. Belnoue and S.R. Hallett, A rapid multi-scale design tool for the prediction of wrinkle defect formation in composite components, Materials & Design 187 (2020) 108388. https://doi.org/10.1016/j.matdes.2019.108388
[2] M. Schwarze and S. Reese, A reduced integration solid-shell finite element based on the EAS and the ANS concept—Large deformation problems, Int. J. Numer. Meth. Engng 85 (2011) 289-329. https://doi.org/10.1002/nme.2966
[3] J. P.-H. Belnoue, O. J. Nixon-Pearson, A. J. Thompson, D. S. Ivanov, K. D. Potter, and S. R. Hallett, Consolidation-Driven Defect Generation in Thick Composite Parts, J. Manuf. Sci. Eng 140 (2018) 071006. https://doi.org/10.1115/1.4039555
[4] J.P.-H. Belnoue, O.J. Nixon-Pearson, Dymitr Ivanov, and S. R. Hallett, A novel hyper-viscoelastic model for consolidation of toughened prepregs under processing conditions. Mechanics of Materials 97 (2016), 118–134. https://doi.org/10.1016/J.MECHMAT.2016.02.019
[5] R. Hauptmann and K. Schweizerhof, A systematic development of ‘solid‐shell’element formulations for linear and non‐linear analyses employing only displacement degrees of freedom, Int. J. Numer. Meth. Engng 42 (1998), 49-69. https://doi.org/10.1002/(SICI)1097-0207(19980515)42:1%3C49::AID-NME349%3E3.0.CO;2-2
[6] E.N. Dvorkin, K.J. Bathe, A continuum mechanics based four‐node shell element for general non‐linear analysis, Engineering computations, 1.1 (1984), 77-88. https://doi.org/10.1108/eb023562
[7] M. Pagani, S. Reese, and U. Perego, Computationally efficient explicit nonlinear analyses using reduced integration-based solid-shell finite elements. Comput. Methods in Appl. Mech. Engng 268 (2014), 141–159. https://doi.org/10.1016/J.CMA.2013.09.005
[8] M.A. Puso, A highly efficient enhanced assumed strain physically stabilized hexahedral element. Int. J. Numer. Meth. Engng 49 (2000), 1029–1064 https://doi.org/10.1002/1097-0207(20001120)49:8%3C1029::AID-NME990%3E3.0.CO;2-3
[9] J. P.H. Belnoue, M.A. Valverde, M. Onoufriou, X. (Ric) Sun, D.S. Ivanov, and S.R. Hallett, On the physical relevance of power law-based equations to describe the compaction behaviour of resin infused fibrous materials, International Journal of Mechanical Sciences 199 (2021), 106425. https://doi.org/10.1016/j.ijmecsci.2021.106425
[10] M.A. Valverde, J.P.H. Belnoue, R. Kupfer, L.F. Kawashita, M. Gude, and S.R. Hallett, Compaction behaviour of continuous fibre-reinforced thermoplastic composites under rapid processing conditions, Composites Part A: Applied Science and Manufacturing 149 (2021), 106549. https://doi.org/10.1016/J.COMPOSITESA.2021.106549
[11] J. Bonet and A.J. Burton, A simple orthotropic, transversely isotropic hyperelastic constitutive equation for large strain computations, Computer Methods in Applied Mechanics and Engineering 162 (1998), 151–164. https://doi.org/10.1016/S0045-7825(97)00339-3
[12] M.M. Hodzelmans, Numerical prediction of wrinkle formation in laminar composites, Student thesis: Master (2024). https://research.tue.nl/en/studentTheses/numerical-prediction-of-wrinkle-formation-in-laminar-composites
