Multilayer Modeling of Delaminated Bilayer Beams and Its Application to Fracture and Buckling Analysis

A new beam element only having the translational displacement degrees of freedom is developed to analyze fracture and buckling behaviors of bilayer beams or columns with delamination using the multilayer modeling method. By comparing their degrees of freedom, it is found that the present beam element is equivalent to the conventional Timoshenko beam element. The upper and lower layers are both divided into one or several sublayers, which are modeled by the developed beam elements. The element stiffness matrix, load vector, and geometric stiffness matrix using the linear interpolation polynomial are derived by the principle of minimum potential energy. To illustrate the effectiveness, accuracy, and efficiency, both the fracture and buckling analyses are conducted. Four typical fracture specimens (i.e., homogeneous and asymmetric double cantilever beams, single leg bending, and asymmetric end-notched flexure specimens) are analyzed using the developed beam element, and the obtained energy release rate and its components using the virtual crack closure technique are compared with those of analytical and conventional two-dimensional (2D) finite element solutions. Further, buckling analysis of bilayer beam-columns with delamination is conducted, focusing on the effect of the number of sublayers. The results show that the developed beam element is capable of effectively and efficiently calculating the convergent stress intensity factor ratio, total energy release rate and its components, and critical buckling load of delaminated bilayer beams or columns with relatively few elements, but without introducing the interface continuity condition in comparison to the conventional modeling method.