Analysis of Intact/Delaminated Composite and Sandwich Beams Using a Higher-Order Modeling Technique

: A simple higher-order model (HOM) is presented in this study for the bending analysis of an intact or delaminated composite and sandwich beam. This model adopts the concept of sub-laminates to simulate multilayered structures, and each sub-laminate takes cubic variation for axial displacement and linear variation for transverse displacement through the thickness. A sub-laminate possesses displacement components at its surfaces (bottom and top) that provide a straightforward way to improve the accuracy of prediction by stacking several sub-laminates. Thus, analysts will have the flexibility to balance the computational cost and the accuracy by selecting an appropriate sub-lamination scheme. The proposed model was implemented by developing a C 0 beam element that has only displacement unknowns. The model was used to solve numerical examples of composite and sandwich beams to demonstrate


Introduction
Composite laminates are assembled by several layers consisting of fiber-reinforced polymer, which are stacked one over the other and bonded to behave as a single monolithic structure.These structures are widely used in mechanical, aerospace, marine, and civil engineering applications owing to their outstanding mechanical features (e.g., high stiffness and strength, and being lightweight) [1].One of the specific types of composite laminates is sandwich structures, where a thick core layer (extremely lightweight with comparatively low stiffness and strength) is stacked in the middle of two thin face sheets (high strength and stiffness) for a more lightweight construction.
The precise modeling of these structures is challenging due to the complex material configurations of the structural system.Generally, the transverse shear deformations of composites are more due to lower shear stiffness compared to extensional rigidity.This has drawn a significant research interest in the recent past and different laminate theories [2][3][4][5] have been proposed.Initially, single-layer-based theories have evolved where a single formulation is employed to describe the through-thickness deformation kinematics so that the multilayered structures are simulated by a single-layer-based model [2].These belong to first-order shear-deformation theory (FSDT) [6][7][8][9] or higher-order shear-deformation theory (HSDT) [10][11][12][13][14][15].The FSDT needs a shear correction factor because it assumes a constant through-thickness shear strain, but an accurate calculation of this factor is cumbersome for a multilayered composite structure [16].
The performance of HSDT is better than FSDT, since it includes cross-sectional warping, and the shear correction factor is no longer needed.However, these equivalent single layers models give a continuous through-thickness variation for displacement fields.This fails to model the interfacial discontinuity of transverse strain which is caused by an abrupt change in material properties between layers having a different fiber orientation.Layer-wise (LW) models [17][18][19][20][21] are developed for accurately modeling the interfacial strain/stress of multilayered structures but they require enormous computational efforts, especially for laminates having a high number of layers, since variables are taken at all interfaces.In order to improve computational efficiency, zigzag theories have been proposed.The main difference between the zigzag theory and LW theory is that the stress continuities of interfacial transverse shear stress are utilized to replace interfacial variables in terms of variables at the reference plane [22][23][24][25][26][27].
The existing zigzag theories can have different forms, with higher-order or linear through-thickness displacement variation for every single layer, but they need a C 1 continuous finite element formulation, which is the main issue unless developing a particular element [28,29].To avoid the requirement of C 1 continuity, different strategies have been adopted for modifying zigzag models with different forms to have a C 0 continuous formulation [30][31][32].The other limitation of the conventional/improved zigzag model is the absence of transverse deformability that must be included for sandwich structures.Due to the relatively high deformability of the thick, soft-core layer, the deviation of transverse displacements between the top and bottom surface cannot be ignored, especially when a localized loading is applied.Layer-wise modes can address this issue, although they are computationally expensive.In addition, a dedicated sandwich model [33,34] has also been developed, where a quadratic variation of transverse displacement is applied only to the core.Moreover, some studies tried to unify various theories with different ranges of generality [35][36][37][38][39].For this purpose, a thorough study was conducted by Carrera and his group using mixed formulations involving stress/strain fields in addition to displacements [36,38,39].
Delamination occurs as a common form of damage in composite laminates that would significantly reduce the strength stiffness of the laminates.It was observed that the number of existing studies on the modeling of delaminated sandwich structures and composite laminates is limited.Moreover, the suitability of existing models for multilayered structures with any arbitrary number/location of delaminations, except the computationally expensive LW [40][41][42][43][44] model, is an issue.Thus, this study attempted to develop a pure displacement-based higher-order model for simulating laminated composite and sandwich beams with any possible combination of delaminations.The model is simple but at the same time it is sufficiently powerful due to its sub-laminate-based modeling capability that can accommodate several bonded layers with a sub-laminate [45][46][47][48] to achieve computational efficiency with adequate accuracy.The model provides cubic variation for axial displacement and linear variation for transverse displacement through the thickness within a sub-laminate.This was achieved by taking unknowns at two internal reference planes in addition to the surfaces (top and bottom) of a sub-laminate.This helps to connect sub-laminates easily due to nodes at the top and bottom surfaces.
The finite element (FE) implementation of the model was performed through the development of a three-nodes beam element with C 0 continuity.The present model was validated by solving static responses of intact laminated composite and sandwich beams and comparing the results to Pagano's exact solution [49] and Reddy's conventional higherorder-equivalent single-layer model [11].As results available in the literature are not adequate for the validation of beams containing delamination, detailed 2D FE models developed by a commercial FE code (ABAQUS) were used for analyzing the delaminated beam.The validation of predictions showed an encouraging performance of the model.Some new results are also reported for future benchmarking, especially for composite beams with delamination.

Mathematical Formulations
This study presents a sub-lamination-based higher-order model (HOM), which is capable of simulating a composite beam having multiple layers with different approaches: (1) treat the whole multilayered beam as a single HOM based sub-laminate (Figure 1a); (2) represent every individual physical layer of the beam by HOM, and the entire multilay-ered beam is simulated by stacking multiple HOMs over the beam's thickness (Figure 1b); or (3) represent several (less than the total number of layers) physical layers in the form of sub-laminates by the HOM and stack them to simulate the beam (Figure 1c).For every sub-laminate applying formulations of HOM (shown in Figure 2), the displacement fields are taken in terms of pure displacements components through the thickness as follows: where L i is the Lagrange polynomial in terms of z, which can be expressed as Taking displacement variables at sub-laminate bottom and top (see u 1 , u 4 , w 1 , w 4 in Figure 2) and two internal axial displacements (u 2 , u 3 ), the proposed model provides cubic variation for axial displacement and linear vibration for transverse displacement through the thickness.It should be noted that the proposed model has benefit in staking multiple sub-laminates to model a composite beam having many layers by incorporating surface unknowns u 4 , u 1 , w 2 , and w 1 , where u 1 and w 1 represent axial and transverse displacements at the sub-laminate bottom, while w 2 and u 4 are transverse and axial displacement at the sub-laminate top.These surface unknowns (u 4 , u 1 , w 2 , w 1 ) help the model directly obtain the displacements at surfaces of a sub-laminate without any additional transformation, which is required in a conventional model only having unknowns at mid-plane.The benefits of this include that several sub-laminates can be stacked by easily connecting the corresponding surface unknowns to simulate a whole multilayered beam, and a delamination is explicitly modeled by debonding them.Moreover, stacking several HOM-based sub-laminates in the thickness direction also improves the accuracy of the modeling in addition to the higher order approximation used for displacement fields of a sub-laminate.

Mathematical Formulations
This study presents a sub-lamination-based higher-order model (HOM), which is capable of simulating a composite beam having multiple layers with different approaches: (1) treat the whole multilayered beam as a single HOM based sub-laminate (Figure 1a); (2) represent every individual physical layer of the beam by HOM, and the entire multilayered beam is simulated by stacking multiple HOMs over the beam's thickness (Figure 1b); or (3) represent several (less than the total number of layers) physical layers in the form of sub-laminates by the HOM and stack them to simulate the beam (Figure 1c).For every sub-laminate applying formulations of HOM (shown in Figure 2), the displacement fields are taken in terms of pure displacements components through the thickness as follows: where Li is the Lagrange polynomial in terms of z, which can be expressed as  () = ∏ , 1  4,   ( = − ,  = − ,  = ,  = ) .Taking displacement variables at sub-laminate bottom and top (see u1, u4, w1, w4 in Figure 2) and two internal axial displacements (u2, u3), the proposed model provides cubic variation for axial displacement and linear vibration for transverse displacement through the thickness.It should be noted that the proposed model has benefit in staking multiple sub-laminates to model a composite beam having many layers by incorporating surface unknowns u4, u1, w2, and w1, where u1 and w1 represent axial and transverse displacements at the sub-laminate bottom, while w2 and u4 are transverse and axial displacement at the sub-laminate top.These surface unknowns (u4, u1, w2, w1) help the model directly obtain the displacements at surfaces of a sub-laminate without any additional transformation, which is required in a conventional model only having unknowns at mid-plane.The benefits of this include that several sublaminates can be stacked by easily connecting the corresponding surface unknowns to simulate a whole multilayered beam, and a delamination is explicitly modeled by debonding them.Moreover, stacking several HOM-based sub-laminates in the thickness direction also improves the accuracy of the modeling in addition to the higher order approximation used for displacement fields of a sub-laminate.Due to the fact that the dimension of a beam is small in its width direction ( compared to other directions, the plane stress assumption is implemented in the xof the beam (Figure 2).The strain components at an arbitrary location are expres the derivative of displacements at this point as Due to the fact that the dimension of a beam is small in its width direction (y-axis) compared to other directions, the plane stress assumption is implemented in the x-z plane of the beam (Figure 2).The strain components at an arbitrary location are expressed by the derivative of displacements at this point as Substituting Equations ( 1) and (2) into Equation (3), and decoupling the two-dimensional strain vector into a cross-sectional matrix [H] and a one-dimensional strain vector {ε} as where [H] is expressed in terms of functions of z: where {ε} is expressed in terms of function of x: In order to retail a good degree of generality for analysis, standard FE procedures were used to implement this model.The quadratic Lagrangian interpolation function was employed as the shape function to approximate the six displacement components (u 1 , u 2, u 3, u 4 , w 1 , and w 2 ) that appeared in Equation ( 6).As a result, a quadratic C 0 beam element was adopted to approximate all six displacement variables along the sub-laminate length.
where N i represents the shape function corresponding to node j [50] and the nodal displacement of variable c at node i is represented by χ c j .Substituting Equation (7) into Equation ( 6) to express {ε} in terms of the strain-displacement matrix [B] and nodal displacement vector {χ}, and as {ε} = [B]{χ} where [B] comprises 1D shape functions and their derivatives.Layers/plies made of composites are idealized as a homogeneous orthogonal material, and their 3D constitutive relationship is expressed as follows [51]: The material stiffness matrix [C] used in the above equation can be written in terms of different direction-dependent elastic modulus (E 1 , E 2 , E3), shear modulus (G 12 , G 23 , G 31 ) and Poisson's ratio (ν 12 , ν 23 , ν 31 ) as where A composite laminate is usually a combination of number of layers having different fiber orientations.Thus, it is necessary to express the stress-strain relationship of the individual layer in terms of the same axis system of the structure so as to combine their individual contributions.This needs transformation of the constitutive relationship of a physical layer expressed in terms of its principal material axis system (1-2-3) to the global laminates axis (x-y-z): The transformed material stiffness matrix C in the above equation can be obtained utilising the fibre orientation (θ) [51]: where [T] is the transform matrix, which is dependent on fibre orientation (θ) as [51] follows: To achieve the plane stress material condition at plane x-z (σ y = τ xy = τ yz = 0), the strainstress relationship of a sub-laminate is expressed with a compliance matrix The 2D plane stress condition can now be easily obtained by eliminating the corresponding rows and columns.(Reduce system of equations from 6 to 3.) This will reduce the 3D problem to plane stress problem, and the reduced strain-stress relationship is reversed again to form the final constitutive relationship of the beam as Now the virtual work principle (δU = δW) is applied to obtain the governing equation within an element.And δU represents the variation of strain energy, which is expressed as Equations ( 4), (8), and (15) are substituted in the above equation to express it in the following form: The above equation can now be expressed in the form of a 1D integral by operating integration along z direction for the cross-section considering a unity width as where The beam can be subjected to any transversely distributed load (q = q(x)) along the x-axis, and the variation of work performed by it (δW) can be written as where w q can be w 1 or w 2 depending on the position of the load.Equation ( 12) can be substituted in the above equation to express it in terms of nodal displacement vector as where w q = [N] q {∆}.Substituting Equations ( 18) and ( 21) into the equation for virtual work principle (δU = δW), the governing equation can finally be obtained and expressed as where the stiffness matrix [K] e and the nodal load vector {P} e of an element can be written as follows: The numerical integration technique Gauss quadrature rule was used to undertake the integration in the above equations, where four integration points were used for integrations within each physical layer along the z-axis, while three integration points were used for integrations in terms of x for a single element.With these equations, [K] e and {P} e of all elements along the sub-laminate length were calculated.The global stiffness matrix [K] and global load vector {P} were then formed by assembling all [K] e and {P} e , and the global beam governing equation in terms of global nodal displacements {χ} g as [K]{χ} g = {P}.The boundary condition was incorporated by the penalty method to restrain the corresponding degrees of freedom, and then {χ} g can be solved.Stresses at an arbitrary location within an element can then be calculated using nodal displacements by substituting Equations ( 5) and ( 9) into ( 16): The utilization of the surface displacements in HOM provides ease in modeling delamination, which is shown in Figure 3.For the contact region, the adjunct sublaminate/element is directly connected in the thickness direction by taking the nodes at the bottom surface of the upper sub-laminate and nodes at the top surface of the lower sub-laminate as common nodes.As for the delaminated region, these common nodes are separated into two sets of nodes without physical gaps so that the continuity of displacements is relaxed in this region.The figure also presents a typical beam element, where unknowns are taken at four reference planes; thus, four nodes are taken along the thickness direction.This model was implemented by developing a beam element where 3 nodes are taken along the axial direction.It should be noted that the 2D strain was decoupled into x (axial) and z (thickness direction) separately (see Equations ( 4)-( 6)).Thus, the finite element interpolation was only applied to the one-dimension strain {ε} in terms of x (see Equations ( 6)-( 8)).The integration along z direction was separately performed by using Gauss quadrature integration (see Equations ( 18) and ( 19)).
is relaxed in this region.The figure also presents a typical beam element, where unknowns are taken at four reference planes; thus, four nodes are taken along the thickness direction.This model was implemented by developing a beam element where 3 nodes are taken along the axial direction.It should be noted that the 2D strain was decoupled into x (axial) and z (thickness direction) separately (see Equations ( 4)-( 6)).Thus, the finite element interpolation was only applied to the one-dimension strain ̅ in terms of x (see Equations ( 6)-( 8)).The integration along z direction was separately performed by using Gauss quadrature integration (see Equations ( 18) and ( 19)).
The above formulation was implemented by developing a code using MATLAB, which was employed to analyze intact/delaminated composite and sandwich beams having various configurations.

Intact Laminated Composite Beams
Simply supported intact laminated composite beams under a sinusoidal distributed load ( =  sin (/) where l is the beam length, and q0 is the load intensity) were first analyzed to show the performance of the proposed model in analyzing both symmetric (0/90/90/0) and asymmetric (0/90/0/90) composite beams.There were three aspect ratios ( = /ℎ where h represents beam depth).Different sub-laminated schemes were used to model those beams, which were represented in terms of the total number of sub-laminates (n) by present (NSn).A convergence study showed the model requires 20 elements in the beam axis to produce stable converged predictions of displacements and stresses.Therefore, the results in this section and other examples were calculated by using at least 20 elements along beam length.These beams were made by four equal-thickness composite plies whose material properties along the fiber orientation were as follows: E1/E2 = 25.0,E2 = E3, G12 = G13 = 0.5E2, G23 = 0.2E2, and ν12 = ν23 = ν13 = 0.25.In addition, separate MATLAB (R2020a, MathWorks, New York, USA) codes were developed based on Pagano's exact The above formulation was implemented by developing a code using MATLAB, which was employed to analyze intact/delaminated composite and sandwich beams having various configurations.

Intact Laminated Composite Beams
Simply supported intact laminated composite beams under a sinusoidal distributed load (q = q 0 sin(πx/l) where l is the beam length, and q 0 is the load intensity) were first analyzed to show the performance of the proposed model in analyzing both symmetric (0/90/90/0) and asymmetric (0/90/0/90) composite beams.There were three aspect ratios (S = l/h where h represents beam depth).Different sub-laminated schemes were used to model those beams, which were represented in terms of the total number of sub-laminates (n) by present (NS n ).A convergence study showed the model requires 20 elements in the beam axis to produce stable converged predictions of displacements and stresses.Therefore, the results in this section and other examples were calculated by using at least 20 elements along beam length.These beams were made by four equal-thickness composite plies whose material properties along the fiber orientation were as follows: E 1 /E 2 = 25.0,E 2 = E 3 , G 12 = G 13 = 0.5E 2 , G 23 = 0.2E 2 , and ν 12 = ν 23 = ν 13 = 0.25.In addition, separate MAT-LAB (R2020a, MathWorks, New York, USA) codes were developed based on Pagano's exact solution [49] and Reddy's higher-order analytical model [11], which were used for validation.
Table 1 presents the mid-span transverse displacements (w) predicted at the midplane by the present mode with three sub-lamination schemes and the analytical models [11,49].All the results in this table are expressed in their non-dimensional form ( w = 100E 2 h 3 w/ l 4 q 0 .They show that the deflection was significantly underestimated by the single-layer-based models for beams with S = 5 (i.e., thick beams).However, the present model NS 1 (using one sub-laminate for modeling) had a relatively better performance compared to Reddy's equivalent single-layer model [11].In addition, the capability of sub-lamination modeling of the proposed model provided good flexibility to users, with a tradeoff between solution accuracy and computational efficiency.To study the performance of the present model in predicting stresses of laminated composites, the mid-span non-dimensional bending stresses (σ xx = σ xx /q 0 ) at the beam's top surface and the non- dimensional mid-plane transverse shear stresses (τ xz = τ xz /q 0 ) at the beam's end (left) are shown in Table 2.The stresses were predicted in these locations because they were expected to be high there.In addition, the elastic exact solution [48] and results produced by Reddy's model [11] are also included in the table.Moreover, the through-thickness variations of σ xx , and τ xz predicted by different models at the mid-span of the symmetrical composite beam are plotted in Figure 4.It should be noted that all transverse shear stresses were calculated directly from the constitutive relationship.A similar performance found in calculating displacements could also be observed for the prediction of stresses from Table 2 and Figure 4.Although modeling the four-layer composite beam with two sub-laminates had a large deviation for calculating the mid-plane transverse shear stresses through a constitutive relationship, it still provided good predictions for bending stresses.Figure 4 also shows that the zero transverse shear stresses conditions and interlaminar continuity of transverse shear stresses were fulfilled when four sub-laminates were used.Thus, the sub-lamination scheme S2 is a good option for balancing the accuracy and computational cost unless an accurate prediction of transverse shear stress is needed.

Intact Laminated Sandwich Beams
A simply supported five-layer sandwich beam (0/90/C/90/0) was analyzed by the proposed model.Two different load cases applied on the beam's top surface were considered: (1) localized load applied at beam mid-span (i.e., point load q = P), and (2) a load that was uniformly distributed ( =  ) along the x-axis.The stiff composite layers used as face z/h z/h

Intact Laminated Sandwich Beams
A simply supported five-layer sandwich beam (0/90/C/90/0) was analyzed by the proposed model.Two different load cases applied on the beam's top surface were considered: (1) localized load applied at beam mid-span (i.e., point load q = P), and (2) a load that was uniformly distributed (q = q 0 ) along the x-axis.The stiff composite layers used as face sheets, and the properties along their fiber direction were as follows: and ν 12 = ν 13 = ν 23 = 0.25.The material properties of the soft-core layer were E 1 =4.0E,E 2 =E 3 =0.05E,G 13 =0.06E, and ν 12 = ν 23 = ν 13 = 0.25.The composite face sheets had a small distribution in material properties since the thickness of each of them was 0.1h, while that of the core layer was 0.8 h.Similarly, Pagano's exact solution [49] and Reddy's equivalent single-layer higherorder model [11] were also employed to analyze the same sandwich beam for the validation of the proposed HOM.
Table 3 shows the mid-span deflection calculated by the above models at the beam's top surface (z = h/2), bottom surface (z = −h/2), and mid-plane (z = 0) for both loads.Different aspect ratios (S) were taken (S = 5, 20, and 50) for this simulation.The deflection (w) of the beam is expressed in non-dimensional form as w = 100E 2 h 3 w/ l 3 P for the point load, while this is the same expression as used in the previous example (Section 3.1) for the uniformly distributed load.It can be observed from Table 3 that the mid-span transverse displacement under the point load was significantly underestimated by Reddy's model [11], especially for the thickest beam (S = 5).Due to the absence of through-thickness variation of transverse displacement in Reddy's model [11], this deviation was most significant (more than 38%) at the surface where the load was applied (top surface).And it can be expected that with the increasing aspect ratio and the loading of a less localized force, this difference would have gradually diminished.A similar problem can also be observed in the present model when the whole beam was modeled with one sub-laminate (NS 1 ).However, the introduction of variation for deflection over the thickness helped the present model (NS 1 ), providing an improved prediction compared to the higher-order model of Reddy [11].For the transverse displacement predicted at the beam's top surface with S = 5 (the worst case), the error of the present model (NS 1 ) was approximately 19.7% (corresponding to Pagano's solution [49]), which was 38.7% for Reddy's results [11].In order to improve the performance of the present model, two different sub-lamination schemes of the present model were also employed to model the sandwich beam: (a) NS 4 , where two sublaminates were used for the thick core layer and every two composite plies were modeled by one sub-laminate; and (b) NS 6 , where each physical composite ply was simulated by one sub-laminate, and the thick-core was simulated by two sub-laminates.Table 3 presents that NS 4 and NS 6 have evident benefits in improving the accuracy of the present model for most cases.Moreover, the results calculated by using these sub-laminate schemes showed good agreement with the exact solution [48], with errors of 0.33% (NS 4 ) and 0.84% (NS 6 ) for the worst case mentioned before.Moreover, the predictions of NS 6 showed a minor improvement compared to the results of NS 4 , and this was expected since the soft-core layer has a much more prominent contribution of material prosperity than the composite sheets along the thickness direction.Thus, the sub-laminate S4 should be an optimum option in this scenario considering its better computational efficiency.4.This table presents transverse displacements (w) and axial displacements (u) (u = 100E 2 h 3 u/ l 3 P ) predicted at four different points: 1, beam top; 2, surface above the delamination; 3, surface below the delamination; and 4, beam bottom at the beam free end.Equations w = 100E 2 h 3 w/ l 3 P and u = 100E 2 h 3 u/ l 3 P were used to nondimensionalize these results.The value in this table shows that the surface above the delamination had a larger deflection than that below the delamination, which is an indicator of the opening of the delamination.For better visualization, the through-beam length variation of deflection at two sides of the delamination was predicted and plotted in Figure 6.Since suitable results available in the literature are not adequate to validate the results of delaminated beams, a reliable commercial finite element package ABAQUS(Version 6.14, DASSAULT SYSTEMS, Paris, French) was employed to develop a detailed FE model for the validation.To accommodate the delamination by using the proposed model, at least two sub-laminates were needed.NS 8 , which provides the highest accuracy, was also used for modeling in addition to an intermediate option, NS 4 .
It can be observed from Table 4 that NS 2 was adequate to simulate the delaminated composite beam in most cases.

Conclusions
This study presented a higher-order model that incorporates the concept of sublamination for the static analysis of intact and delaminated composite and sandwich beams.By taking displacement components at the beam top, bottom, and two internal planes, the proposed model achieved a cubic variation for axial displacement and linear variation for transverse displacement through the thickness.The presence of surface displacement parameters showed a significant benefit in modeling multilayered structures owing to the ease in connection of several sub-laminates in the thickness direction of the structure.The model was implemented by developing a pure displacement-based C 0 beam element with three nodes.It was validated by analyzing different intact or delaminated composite and sandwich beams and comparing the results with existing solutions and results produced by a commercial finite element program.The following conclusions were made based on the numerical results.
• The present model shows encouraging performance in analyzing multilayered structures subjected to different types of static loading.The incorporation of transverse displacement variation was beneficial in modeling the sandwich beams, especially subjected to localized load (e.g., point load).The proposed model also showed a good estimation for transverse shear stresses, which were calculated directly from the constitutive relationship.Moreover, modeling the whole structure with a single sub-laminate showed a better performance than the typical higher-order model with the equivalent single-layer hypothesis.
• The model provides ease in the simulation of delamination by inserting it between two sub-laminates.The numerical results also showed that the model has adequate accuracy in analyzing delaminated composite beams.

•
The HOM-based model has high flexibility to choose the level of detail in the variation of displacement through the thickness for the modeling of multilayer composite and sandwich beams with or without delamination by choosing an appropriate sub-lamination scheme.Therefore, users can balance computational efficiency and prediction accuracy.

Figure 4 .
Figure 4. Through-thickness variation of (a)   , and (b) ̅  for the symmetrical composite beams predicted at the beam midspan by different models [11,49].

Figure 4 .
Figure 4. Through-thickness variation of (a) σ xx , and (b) τ xz for the symmetrical composite beams predicted at the beam midspan by different models [11,49].
Top surface of the delaminationBottom surface of the delamination

Figure 5 .
Figure 5.A cantilever laminated beam [0/90/0/90] s with single delamination under a point load.The present model with three sub-lamination schemes (NS 2 ,NS 4 ,NS 8 ) was used to simulate the beam and the results are shown in Table4.This table presents transverse displacements (w) and axial displacements (u) (u = 100E 2 h 3 u/ l 3 P ) predicted at four different points: 1, beam top; 2, surface above the delamination; 3, surface below the delamination; and 4, beam bottom at the beam free end.Equations w = 100E 2 h 3 w/ l 3 P and u = 100E 2 h 3 u/ l 3 P were used to nondimensionalize these results.The value in this table shows that the surface above the delamination had a larger deflection than that below the delamination, which is an indicator of the opening of the delamination.For better visualization, the through-beam length variation of deflection at two sides of the delamination was predicted and plotted in Figure6.

Figure 6 .Figure 6 .
Figure 6.Variation of transverse displacement captured at surfaces above and below the delamination.

Table 1 .
Non-dimensional middle-plane deflection (w) at mid-span for simply supported symmetrical and antisymmetric intact composite beams under sinusoidal load.
[49]rror was calculated corresponding to the exact solution of Pagano[49].

Table 2 .
Non-dimensional bending stresses and transverse shear stresses of simply supported symmetrical composite beams under distributed transverse load.

Table 3 .
(a): Nondimensionalized deflection (w) at the mid-span of a simply supported sandwich beam under localized loading.(b) Nondimensionalized deflection (w) at the mid-span of a simply supported sandwich beam under uniformly distributed load.

Table 4 .
Non -dimensional displacements obtained from different depths (z) of the free end of the cantilever laminated beam [0/90/0/90]s with single delamination.a surface above the delamination; b surface below the delamination.

Table 4 .
Non-dimensional displacements obtained from different depths (z) of the free end of the cantilever laminated beam [0/90/0/90]s with single delamination.
a surface above the delamination; b surface below the delamination.

Table 4 .
Non-dimensional displacements obtained from different depths (z) of the free end of the cantilever laminated beam [0/90/0/90] s with single delamination.
a surface above the delamination; b surface below the delamination.