The Multiscale Damage Mechanics in Objected-oriented Fortran Framework

We develop a dual-purpose damage model (DPDM) that can simultaneously model intralayer damage (ply failure) and interlayer damage (delamination) as an alternative to conventional practices that models ply failure by continuum damage mechanics (CDM) and delamination by cohesive elements. From purely computational point of view, if successful, the proposed approach will significantly reduce computational cost by eliminating the need for having double nodes at ply interfaces. At the core, DPDM is based on the regularized continuum damage mechanics approach with vectorial representation of damage and ellipsoidal damage surface. Shear correction factors are introduced to match the mixed mode fracture toughness of an analytical cohesive zone model. A predictor-corrector local-nonlocal regularization scheme, which treats intralayer portion of damage as nonlocal and interlayer damage as local, is developed and verified. Two variants of the DPDM are studied: a single- and two- scale DPDM. For the two-scale DPDM, reduced-order-homogenization (ROH) framework is employed with matrix phase modeled by the DPDM while the inclusion phase modeled by the CDM. The proposed DPDM is verified on several multi-layer laminates with various ply orientations including double-cantilever beam (DCB), end-notch-flexure (ENF), mixed-mode-bending (MMB), and three-point-bending (TPB). The simulation is executed in the platform of FOOF (Finite element solver based on Object-Oriented Fortran). The objective of FOOF is to develop a new architecture of the nonlinear multiphysics finite element code in object oriented Fortran environment. The salient features of FOOF are reusability, extensibility, and performance. Computational efficiency stems from the intrinsic optimization of numerical computing intrinsic to Fortran, while reusability and extensibility is inherited from the support of object-oriented programming style in Fortran 2003 and its later versions. The shortcomings of the object oriented style in Fortran 2003 (in comparison to C++) are alleviated by introducing the class hierarchy and by utilizing a multilevel programming style.