Determination of Transverse Shear Stresses and Delamination in Composite Laminates Using Finite Elements

Author(s):  
George Z. Voyiadjis ◽  
Pawel Woelke
2012 ◽  
Vol 19 (3) ◽  
pp. 459-475 ◽  
Author(s):  
Pawel Woelke ◽  
Ka-Kin Chan ◽  
Raymond Daddazio ◽  
Najib Abboud ◽  
George Z. Voyiadjis

We present a simple and efficient method for the analysis of shear flexible isotropic and orthotropic composite shells. Classical thin shell constitutive equations used in the explicit finite element code EPSA to model homogenous isotropic shells using "through-the-thickness-integration" and layered orthotropic composite shells [1–3,5] are modified to account for transverse shear deformation. This effect is important in the analysis of thick plates and shells as well as composite laminates, where interlaminar effects matter. Transverse shear stresses are calculated using a linear normal strain distribution, where first the shear forces are calculated and then the stresses are calculated by means of the generalized section properties, i.e., first and second moments of area. The formulation is a generalization of the analytical method of analyzing composite beams. It is simple and computationally inexpensive, and it yields accurate results without employing higher order displacement interpolations. In the case of isotropic shells, the transverse shear stresses are distributed parabolically, based on the assumption of linear normal strain distribution through the thickness and on application of the quadratic shape function to transverse shear strains. The transverse shear stresses are included in the elastic-perfectly plastic yield function of the Huber-Mises-Hencky type.


2020 ◽  
Vol 86 (2) ◽  
pp. 44-53
Author(s):  
Yu. I. Dudarkov ◽  
M. V. Limonin

An engineering approach to estimation of the transverse shear stresses in layered composites is developed. The technique is based on the well-known D. I. Zhuravsky equation for shear stresses in an isotropic beam upon transverse bending. In general, application of this equation to a composite beam is incorrect due to the heterogeneity of the composite structure. According to the proposed method, at the first stage of its implementation, a transition to the equivalent model of a homogeneous beam is made, for which the Zhuravsky formula is valid. The transition is carried out by changing the shape of the cross section of the beam, provided that the bending stiffness and generalized elastic modulus remain the same. The calculated shear stresses in the equivalent beam are then converted to the stress values in the original composite beam from the equilibrium condition. The main equations and definitions of the method as well as the analytical equation for estimation of the transverse shear stress in a composite beam are presented. The method is verified by comparing the analytical solution and the results of the numerical solution of the problem by finite element method (FEM). It is shown that laminate stacking sequence has a significant impact both on the character and on the value of the transverse shear stress distribution. The limits of the applicability of the developed technique attributed to the conditions of the validity of the hypothesis of straight normal are considered. It is noted that under this hypothesis the shear stresses do not depend on the layer shear modulus, which explains the absence of this parameter in the obtained equation. The classical theory of laminate composites is based on the similar assumptions, which gives ground to use this equation for an approximate estimation of the transverse shear stresses in in a layered composite package.


1997 ◽  
Vol 64 (3) ◽  
pp. 525-531 ◽  
Author(s):  
Wen-Hwa Chen ◽  
Tain-Fu Huang

By utilizing the general solutions derived for the plies with arbitrary fiber orientations under uniform axial strain (Huang and Chen, 1994), the explicit solutions of the edge-delamination stress singularities for the angle-ply and cross-ply laminates are obtained. The dominant edge-delamination stress singularities for the angle-ply laminates are found to be a real constant, −1/2, and a pair of complex conjugates, −1/2±i/2πln{(b+b2−a2)/a}. For the cross-ply laminates, the significant effect of transverse shear stresses of the laminate is considered and the dominant edge-delamination stress singularities are shown as −1/2 and −1/2±i/2πln{(c2+c22−4c1c3)/2c1}. a, b, cl, c2, and c3 are the corresponding combined complex constants. In addition, two elementary forms of edge-delamination stress singularity, say, r−1/2 and rδr(lnr)n(δr=n−3/2,n=1,2...) exist for both the angle-ply and cross-ply laminates. Excellent correlations between the present results and available solutions show the validity of the approach. The deficiencies of the solutions available in the literature are compensated. New results for other angle-ply and cross-ply laminates are also provided.


2012 ◽  
Vol 29 (2) ◽  
pp. 241-252 ◽  
Author(s):  
A. S. Sayyad ◽  
Y. M. Ghugal

AbstractThis paper deals with the problem of stress distribution in orthotropic and laminated plates subjected to central concentrated load. An equivalent single layer trigonometric shear deformation theory taking into account transverse shear deformation effect as well as transverse normal strain effect is used to obtain in-plane normal and transverse shear stresses through the thickness of plate. Governing equations and boundary conditions of the theory are obtained using the principle of virtual work. A simply supported plate with central concentrated load is considered for the numerical analysis. Anomalous behavior of inplane normal and transverse shear stresses is observed due to effect of stress concentration compared to classical plate theory and first order shear deformation theory.


2017 ◽  
Vol 21 (6) ◽  
pp. 1820-1842
Author(s):  
Wu Zhen ◽  
Ma Rui ◽  
Chen Wanji

This paper will try to overcome two difficulties encountered by the C0 three-node triangular element based on the displacement-based higher-order models. They are (i) transverse shear stresses computed from constitutive equations vanish at the clamped edges, and (ii) it is difficult to accurately produce the transverse shear stresses even using the integration of the three-dimensional equilibrium equation. Invalidation of the equilibrium equation approach ought to attribute to the higher-order derivations of displacement parameters involved in transverse shear stress components after integrating three-dimensional equilibrium equation. Thus, the higher-order derivatives of displacement parameters will be taken out from transverse shear stress field by using the three-field Hu–Washizu variational principle before the finite element procedure is implemented. Therefore, such method is named as the preprocessing method for transverse shear stresses in present work. Because the higher-order derivatives of displacement parameters have been eliminated, a C0 three-node triangular element based on the higher-order zig-zag theory can be presented by using the linear interpolation function. Performance of the proposed element is numerically evaluated by analyzing multilayered sandwich plates with different loading conditions, lamination sequences, material constants and boundary conditions, and it can be found that the present model works well in the finite element framework.


2011 ◽  
Vol 110-116 ◽  
pp. 1458-1465 ◽  
Author(s):  
M. Khadem ◽  
M. M. Kheirikhah

Nowadays Shape Memory Alloys (SMAs) are used as actuators in many applications such as aerospace structures. In sandwich structures, the SMA wires or plates are used in the skins for shape control of the structure or vibration damping. In this paper, bending behavior of sandwich plates with embedded SMA wires in their skins is studied. 3D finite element method is used for construction and analysis of the sandwich plate with a flexible core and two stiff skins. Some important points such as continuity conditions of the displacements, satisfaction of interlaminar transverse shear stresses, the conditions of zero transverse shear stresses on the upper and lower surfaces and in-plane and transverse flexibility of soft core are considered for accurate modeling and analysis of sandwich structures. Solution for bending analysis of sandwich plates under various transverse loads are presented and the effect of many parameters such as plate dimensions, loading conditions, material properties of core, skins and SMA wires are studied. Comparison of the present results in special case with those of the three-dimensional theory of elasticity and some plate theories confirms the accuracy of the proposed model.


2018 ◽  
Vol 204 ◽  
pp. 896-911 ◽  
Author(s):  
J. Awrejcewicz ◽  
V.A. Krysko ◽  
M.V. Zhigalov ◽  
I.V. Papkova ◽  
V.A. Krysko

Sign in / Sign up

Export Citation Format

Share Document