Microsoft Word - numero_51_art_22_2683 R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 275 Focussed on IGF25 – Fracture and Structural Integrity International Conference 2019 Local zigzag effects and brittle delamination fracture of n-layered beams using a structural theory with three displacement variables Roberta Massabò, Ilaria Monetto Department of Civil, Chemical and Environmental Engineering, University of Genova, Italy Roberta.massabo@unige.it, http://orcid.org/0000-0002-9746-1801 ABSTRACT. Equivalent single layer theories for layered beams effectively and accurately predict global displacements and internal force and moment resultants using a limited number of displacement variables. However, they cannot reproduce local effects due to material architecture or weak/imperfect bonding of the layers, such as zigzag displacement fields, displacement jumps at the layer interfaces and complex transverse stress fields, nor can they simulate delamination damage growth. In this work we will present some applications and discuss advantages and limitations of a recently formulated zigzag model. The model, through a modification of the equilibrium equations of an equivalent single layer theory, which maintains the same number of variables, reproduces local effects and delamination fracture under mode II dominant conditions. The approach is based on a local enrichment of the displacement field of first order shear deformation theory, the introduction of cohesive interfaces and homogenization. KEYWORDS. Structural theory; Zigzag; Homogenization; Cohesive zone model; Brittle fracture; Laminate. Citation: Massabò, R., Monetto, I., Local zigzag effects and brittle delamination fracture of n-layered beams using a structural theory with three displacement variables, Frattura ed Integrità Strutturale, 51 (2020) 275-287. Received: 13.11.2019 Accepted: 27.11.2019 Published: 01.01.2020 Copyright: © 2020 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION he mechanical response of layered structures, such as laminates, sandwich composites, laminated glass, and laminated wood, is controlled by local effects due to the elastic mismatch of the layers, the presence of thin interlayers (adhesives, resin rich regions, polymeric layers), and interfacial imperfections, defects and delaminations, which are typical consequences of manufacturing processes. Zigzag displacement fields and complex transverse stress distributions, which are due to the multilayer architecture, and displacement discontinuities, due to the presence of interlayers or interfacial damage/delaminations, are not captured by the widely used equivalent single layer theories (ESLT), such as first order shear deformation theory (FSDT). These theories assume a C1 continuous displacement field through the thickness which introduces discontinuities in transverse stresses and the impossibility to satisfy local equilibrium. Modeling local effects requires the use of layerwise, discrete-layer or 2D/3D finite-element models, where layers and interfacial regions between layers are distinctly modelled and continuity of displacements, but not of their derivatives (slopes), is imposed at the interfaces [1]. Modeling delamination evolution using these techniques requires in addition a T http://www.gruppofrattura.it/VA/51/2683.mp4 R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 276 cohesive zone approach to capture the displacement discontinuities and an in-plane discretization to follow the evolution of the cracks. The main drawback of these approaches is that the number of displacement variables is related to the number of mathematical layers and to the kinematics assumed to describe the layers. To model a n-layered beam with perfectly bonded layers and transverse inextensibility using Timoshenko beam theory in each layer leads to a total of (n+2) displacement unknowns (namely (2n+1) initial displacements minus (n-1) displacement continuity conditions at the interfaces). Delamination damage evolution in the beam can be studied by releasing the displacement continuity at the layer interfaces and introducing cohesive interfacial tractions, normal and tangential, and cohesive traction laws which relate the tractions to the relative sliding and opening displacements [2–4]. This approach requires 3n unknowns for each beam region (or (2n+1) unknowns for mode II dominant problems, where crack opening is neglected and a single displacement variable describes the transverse displacements of the layers). An effective approach for modeling multilayered structures was proposed in the early works in [5,6] and is known in the community as “zigzag approach”. The basic idea of the global-local zigzag approach is to enrich the displacement field of an equivalent single layer theory (ELST) by local zigzag functions [5,6, and subsequent works] in order to account for the effects of the multilayer structure (zigzag). This adds a large number of local kinematic variables to the global variables, which are then removed from the problem by imposing continuities conditions on interfacial tractions at the layer interfaces [6]. Variational techniques are then used to derive the equilibrium equations, which depend on the global variables only. Other refined approaches have been proposed, e.g. [7,8], which introduce additional global variables to try and overcome some drawbacks of the original theories. The presence of thin interlayers and interfacial damage has been modelled following two different techniques. The first, known as “compliant layer concept”, is based on a description of thin interlayers or interfacial damaged regions as regular layers of the stack, with elastic constants apt to describe their response; the method then essentially follows the procedure described above for fully bonded structures [9,10]. The second technique, known as “spring-layer approach”, introduces zero thickness interfaces and interfacial traction laws which describe the response of the interlayer/damaged regions [11– 14]. In addition to the zigzag functions, the displacement field of the global equivalent single layer theory is then enriched by interfacial jumps. Considering the n-layered plate in Fig. 1a, the displacements ( )k v (=1,2,3) in an arbitrary layer k , with upper and lower coordinates 13 kx  and 3 kx , are described by the following equations: 1 1 ( ) 1 1 2 3 01 1 2 3 1 1 2 1 1 2 3 3 1 1 2 1 1 ˆ( , , , ) ( , , ) ( , , ) ( , )( ) ( , , ) k k k i i i i i v x x x t v x x t x x x t x x x x v x x t            1 1 ( ) 2 1 2 3 02 1 2 3 2 1 2 2 1 2 3 3 2 1 2 1 1 ˆ( , , , ) ( , , ) ( , , ) ( , )( ) ( , , ) k k k i i i i i v x x x t v x x t x x x t x x x x v x x t            (1) 1 1 ( ) 3 1 2 3 0 1 2 3 3 1 2 3 1 2 3 3 3 1 2 1 1 ˆ( , , , ) ( , , ) ( , , ) ( , )( ) ( , , ) k k k i i i i i v x x x t w x x t x x x t x x x x v x x t            where 01 02 0 1 2 3, , , , ,v v w    are the primary variables of the equivalent single layer model (here a first order shear and normal deformation plate theory), 1 2 3 3( , )( ) i ix x x x  for 1, 2, 3  and 1,..., 1i n  are piecewise linear zigzag functions and ( 1) ( ) 1 2 3 3 3 3ˆ ( , , ) ( ) ( ) i k k k kv x x t v x x v x x       for 1, 2, 3  and 1,..., 1i n  are displacement jumps at the layer interfaces [14]. Higher order equivalent single layer models can be used which increase the number of the global variables and better describe the transverse stresses in fully bonded problems [1]. Zigzag functions and displacement jumps are then defined by imposing continuity of transverse shear and normal tractions at the layer interfaces and the interfacial constitutive laws to relate tractions and jumps. This yields the macro-scale displacement field and homogenized equilibrium equations which depend on the global variables only. Using this second approach care must be put on properly accounting for the strain energy absorbed during deformation within the interface [15] in the variational equations used for the derivation of the equations of motion. This term was neglected in the early models in [11–13] and corrected only later in [14]. Zigzag models have been formulated to study interfacial damage evolution using the compliant-layer concept and continuum damage approaches in [9,16,17]. Recently, the model in [14], which uses a cohesive zone approach and the zigzag kinematic approximation in Eqn. (1), has been applied to simulate single and multiple delamination growth under mode II dominant conditions in wide plates within the framework of fracture mechanics [18]. R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 277 In this paper the model in [14] is briefly reviewed, the equations which describe the response of n-layered beams with imperfect interfaces and delaminations under mode II dominant conditions, as shown in Figs. 1b,c, are recalled and some applications of the model, novel and from the literature, are presented and discussed. The applications are chosen in order to highlight the applicability of the homogenized approach to predict local stress and displacement fields in layered beams with imperfect interfaces and delaminations and to analyze single and multiple delamination crack propagation in homogeneous and layered beams with different boundary conditions. The results are used to discuss advantages and limitations of the homogenized approach. Figure 1: (a) n-layered plate with weakly bonded layers modelled as imperfect interfaces and delaminations; (b1,2) n-layered beam with imperfect interfaces and related homogenized description; (c1,2) n-layered beam with imperfect interfaces and delaminations and related homogenized description. The homogenized description eliminates the need for through-thickness and longitudinal discretization in the presence of continuous imperfect interfaces (b2), and requires only longitudinal discretization in the presence of delaminations, (c2). (a) (b) Figure 2: (a) Exemplary interfacial traction laws used to describe a n-layered beam with continuous fully bonded, imperfectly bonded and fully detached interfaces, as shown in Fig. 1b (after [14]); (b) interfacial traction law used to analyze perfectly brittle fracture in the n-layered beam in Fig. 1c; the crack tip is defined at the coordinate where the sliding displacement becomes critical, . 2 2ˆ ˆ cv v ., and the area under the first branch up to 2ˆ cv is the critical energy release rate or mode II fracture toughness of the interface (after [18]). MULTISCALE STRUCTURAL MODEL FOR N-LAYERED BEAMS WITH INTERFACIAL IMPERFECTIONS AND DELAMINATIONS n this section the equations of the model formulated in [14] for plates with weakly bonded layers, oriented along the geometrical axes, and delaminations shown in Fig. 1a are particularized to a beam (wide plate) with imperfect sliding interfaces and delaminations and used to describe the two schematics shown in Figs. 1b and 1c. Assuming transverse inextensibility and the absence of normal separation between the layers, in order to describe mode II dominant problems with layers in constrained contact, the displacement field in Eqn. (1) modifies since the dependence on 1x is removed and 3 3 3ˆ 0 i iv     , for 1,..., 1i n  . The displacement unknowns of the problem are then the global displacements 02 0 2, ,v w  , the zigzag functions 2 3 3( ) i ix x  for 1,..., 1i n  and the jumps 2ˆ iv for 1,..., 1i n  , for a total of I R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 278 3 2( 1)n  . The ( 1)n  zigzag functions are defined in terms of the global variables by imposing continuity of the shear tractions at the ( 1)n  layer interfaces (the shear tractions are calculated using the elastic constitutive laws in Eqn. (7) in the Appendix). The remaining local unknowns, which describe the interfacial jumps, depend on the mechanisms acting at the interfaces which are described through piece-wise linear interfacial traction laws relating interfacial shear tractions and jumps. To describe the beam in Fig. 1b, with continuous imperfect interfaces, such as those due to thin elastic interlayers (adhesives, resin rich layers), the linear traction laws shown in Fig. 2a, 2ˆ ˆ k k k S SK v  , are used, which relates k SK , the interfacial stiffness of the interface k, and the jump 2ˆ kv ; the interfacial stiffness may be different in the different interfaces and 0kSK  is used to describe fully debonded layers, so that 2ˆ ˆ0, 0 k k S v   , while 1/ 0 k SK  describes perfect bonding, where 2ˆ 0 kv  . In the beam in Fig. 1c, with imperfect interfaces and delaminations, the interfacial stiffness also varies between the different longitudinal regions and a longitudinal discretization of the problem is required; different interfacial traction laws are used to define the traction free delaminations and the bonded or partially bonded ligaments ahead of the delamination tips. To simulate perfectly brittle fracture of the layers, the law in Fig. 2b, with 1/ 0kSK  in the intact regions and 0 k SK  along the traction free delaminations, is used and the critical energy release rate is defined by the area under the first branch. The crack tips are located where the crack sliding displacement equals the critical value, 2 2ˆ ˆ k cv v . The interfacial jumps are defined in terms of the global variables by imposing equilibrium conditions at the interface boundaries, which relate interfacial tractions and shear stresses in the layers. Using Eqns. (1) and (7), the macroscale displacement field then takes the form:  ( ) 2 2 3 02 2 3 2 2 0 2 2 2 2 22 3 ( ) 3 2 0 2 ( , ) ( ) ( ) , ( ) ( ) ( ) ( ) ( ) k k S k v x x v x x x w x x R x v x w x       (2) with   ( 1)1 (1) (1)23 22 3 23 3 3 23( 1) ( ) ( 1) ( ) 1 123 23 23 23 1 1 1 1 ( ) 1 ik i k i S i i i j j i jS G R x G x x G G G K G G                               (3) where ( ) 23 i G is the shear modulus of the layer i. In a beam with layers having the same elastic properties, as a unidirectionally reinforced laminate, ( 1) ( )23 23 k kG G  and Eqn. (3) simplifies in ( 1)22 3 231.. 1( ) / k i i S Si k R x G K     ; in the intact portions, where 1/ 0kSK  , 22 3( ) 0 k SR x  and the displacement field coincides with that of first order shear deformation theory. If the elastic constants of the layers differ, in the intact portions of the plate where 1/ 0kSK  the displacement field coincides with that assumed by the original zigzag theory in [6]. The low order of the global model and Eqns. (2),(7) imply that the shear strains obtained from (2) through compatibility are constants in and between the layers in the intact regions and vanish in the delaminated regions. Following the approach which is commonly used for the structural low order theories, accurate shear stresses and strains can be obtained a posteriori through the imposition of local equilibrium: ( ) ( )22 2 23 3, , 0 k k post   and ( ) ( ) ( )23 23 23 2 / k post k post k G  . On the other hand, the absence of shear strains in the delaminated regions induces under-predictions of the transverse displacements in shear deformable plates and the need for corrections or adjustments [19,20]. The problem will be discussed later. Variationally consistent equilibrium equations and boundary conditions are derived through the principle of Virtual Works [14]. They are shown below for a beam subjected to distributed transverse load, 3 2( )f x : 22 2 2 22 2 2 2 2 2 2 2 3 2 22 2 2 02 02 22 2 2 2 2 2 2 3 0 0 22 2 2 0 2 0 2 , ( ) 0; , ( ) ( ) 0; , ( ) ( ) 0 or ; or ; or ; or , , b g g b b zS zS g N x M x Q x Q x f x N n N v v M n M Q n N w w M n M w w                    (4) R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 279 where 22 22 2, , b gN M Q are axial force, bending moment and generalized shear force, given by (8) in the Appendix, where the tilde defines prescribed values at the beam edges and 2 1n   is the second component of the outward unit normal. The equations have the same form of those of first order shear deformation theory and the effects of the local enrichments (zigzag and interfacial jumps) enter the generalized shear force, Eqn. (8), which is variationally consistent and coincides with the resultant of the a posteriori calculated shear stresses, 3 1 3 ( ) 2 23 31.. k k x k post g k n x Q dx     , as in classical structural theories. The constitutive equations of the homogenized beam have been derived in [18]:   0 0 22 22 22 02 2 02 22 1 1 0 1 2 2 22 22 22 2 2 2 44 2 0 2 2 22 0 1 2 2 0 22 0 22222 , , , ; , , , , S S b S S S S S S g zS S S S S N A B C C v v M B D C C Q A w C C C C w wM C C C C                                               (5) with 3 1 3 ( ) 2 22 22 22 22 3 3 3 44 44 44 1 , , D ) (1, , ) ; ( k k n xk P S x k C x x dx A kA B C C       where  ( )( ) 22 22 23 32 33/ kk C C C C C  with ( )k ijC the coefficients of the 6×6 stiffness matrix in the layer k ; the remaining coefficients are given in (10) in the Appendix; 44k is a shear correction factor which can be set equal to 5/6, as for a homogeneous beam, since the effects of the layered structure are already accounted for through the zigzag functions. The equilibrium equations in terms of global displacement variables (for a beam subjected to transverse surface tractions, 3 2( )f x ) are [19]: 0 0 22 02 22 22 2 22 0 222 0 1 2 1 2 22 02 22 22 2 22 0 222 44 0 2 2 0 1 2 2 02 222 2 222 0 2222 44 0 22 2 2 3 , ( ) , , 0 ( ) , ( 2 ) , ( ) , ( , ) 0 , ( ) , , ( , , ) 0 S S S S S S S S S S S A v B C C w B C v D C C C C w A w C v C C C w A w f                          (6) They return the equilibrium equations of first order shear deformation theory in a n-layered beam with equal fully bonded layers, where all coefficients with upper index S vanish. And they return the equilibrium equations of the zigzag theory in [6] for fully bonded unequal layers. LOCAL STRESS AND DISPLACEMENT FIELDS AND DELAMINATION FRACTURE OF N-LAYERED BEAMS n this section various applications of the model are presented in order to highlight its main features. Stress and displacement fields in n-layered beams subjected to mechanical loadings The equilibrium Eqns. (6) with boundary conditions in (4) are easily solved for n-layered beams and wide plates with continuous imperfect interfaces, shown in Fig. 1(b); the model is applicable to beams with various boundary conditions (simply supported, clamped) and loadings (distributed, concentrated loads). Exemplary closed form solutions have been presented in [19] for simply supported plates under thermo-mechanical loadings; in [21] closed form solutions have been derived for the wave propagation problem. The model accurately captures the transverse displacements in fully bonded beams while requires a correction in shear deformable beams with fully debonded or weakly bonded interfaces where the I R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 280 effects of the shear strains on the global transverse displacement are not fully captured. The local fields are accurately captured also in very thick beams and wide plates. Predictions are less accurate near clamped edges, due to the presence of boundary regions which will be discussed later [19,20,22]. The diagrams in Fig. 3 show transverse shear stresses through the thickness in a unidirectionally reinforced orthotropic wide plate with fibers (L) along 2x , two layers (n=2) and one interface. Results in the figure refer to a simply supported plate subjected to a sinusoidal transverse load. The solutions of the homogenized zigzag model are compared with exact elasticity solutions for perfectly bonded, partially bonded and fully debonded layers, after [14]. Geometry and elastic constants are defined in the caption. Fig. 4 shows axial displacements and transverse shear stresses in a thick wide-plate with three symmetrically oriented orthotropic layers (0/90/0) and two equally spaced weak interfaces. The boundary conditions are the same of the previous case. The homogenized model accurately captures the local effects due to the material architecture, e.g. the zigzag displacements, the jumps at the layers interfaces and their effects on transverse shear stresses. Figure 3: Dimensionless transverse shear stresses shown through thickness at 2x L in the thick unidirectionally reinforced wide plate with two layers and one imperfect interface shown on the right. The plate is subjected to a sinusoidal transverse load, 0 2sin( / )p q x L  . Materials and geometry: / 5L h  , / 25T LE E  , / 50LT LG E  , , 0.25LT TT   (modified after [14] ). Figure 4: Axial displacements and transverse shear stresses at 2x L shown through thickness in a three-layer laminate (0/90/0), thick simply supported wide plate subjected to sinusoidal transverse load, 0 2sin( / )p q x L  . Materials and geometry: / 4L h  , / 25T LE E  , / 50LT LG E  , / 125TT LG E  , 0.25LT TT   (modified after [19]). R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 281 Single and multiple delamination fracture in n-layered beams The equilibrium Eqns. (6) can be applied directly to study the response of beams (wide plates) with stationary delaminations, which are regions where the layers are fully debonded, as shown in Fig. 1c. For these problems a discretization in the axial direction is required, as shown in Fig. 1c2 since the coefficients in Eqns. (6) and (4) differ between regions. Continuity conditions between the different regions are then applied and the problem is solved analytically or semi-analytically. Displacement and stress fields accurately reproduce the 2D elasticity solutions but for boundary regions, which occur near clamped edges or crack tip cross sections, and will be discussed later in the paper. The energy release rate for the collinear propagation of the delaminations can be calculated using different techniques, namely the compliance method and the J- integral along different paths [18]. Different fracture problems have been examined in [18]: bend beams with two and three layers (sandwiches), with single and multiple delaminations and various boundary conditions (simply supported and clamped). Here two relevant applications after [18] are recalled and novel applications are presented in order to highlight the capability of the model, which uses only three displacement variables as a classical single layer theory, to describe the discrete event of brittle delamination fracture, to reproduce the interaction effects of multiple delaminations and to follow the evolution of multiple cracks, in homogeneous and layered beams. (a) (b) Figure 5: (a) Critical load for crack propagation versus load point displacement in a ENF specimen (modified after [18] ). Geometry and materials: L = 50 mm, 2h = 3.4 mm, 0 25a  mm, b = 25 mm (width). Graphite/epoxy AS-4/828 with 139 16.7LE   GPa and 6LTG  GPa and Mode II fracture energy 1.04 0.17IIC  G N/mm [24]. (b) Dimensionless critical load for crack propagation versus crack length in an End Load Split specimen. Unidirectional E-glass-epoxy with 41LE  GPa, 10.4TE  GPa, 4.3LTG  GPa , 0.28, 0.5LT TT   . Exact solution in [27]. Crack propagation in unidirectionally reinforced End-Notched Flexure and End-Load Split specimens Fig. 5a shows the macro-structural response of the unidirectional End-Notched Flexure (ENF) specimen shown in the inset, through the critical load for crack propagation versus load-point deflection (after [18]). The layup and elastic constants are described in the caption. The energy release rate has been obtained using the compliance method and the transverse displacement calculated through Eqns. (6). Local effects in the layers due to the presence of the delamination are well captured by the model and the distribution of the shear stresses is qualitatively similar to that presented in Fig. 3. This is a pure mode II problem and the dimensionless critical load is obtained in closed form and is 2/cr IICP hE G  4 / 3 /h a with IICG the critical fracture toughness [18]. The critical load coincides with predictions using classical structural mechanics assumptions and a discrete-layer approach. The exact solution of the problem [23] has two additional terms, which account for the effects of crack tip shear on the near tip deformations and are important when the crack is short. The response diagram is compared with the experimental results on a graphite/epoxy [0]24 laminate tested in [24]. The two theoretical curves shown in the diagram have been obtained using the average values of the elastic constants and energy release rate, R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 282 red dashed lines Model (a), and the maximum values, black solid lines Model (b). The experimental results, obtained under displacement, control show a load drop in the critical load at the onset of propagation; this is due to an unstable propagation of the crack which grows catastrophically and arrests near the mid-span. The homogenized model, which is under crack- length control [18,25], is able to capture the snap-back instability and follow the virtual branch where crack growth is associated to a reduction of the load-point displacement. Crack propagation is modelled also in the region beyond the mid- span to show that the curve stably approaches the limiting solution (dotted line) corresponding to two fully delaminated layers. Fig. 5b shows the critical load for crack propagation in an End Loaded Split (ELS) specimen. The material is a unidirectional E-glass-epoxy laminate with elastic constants given in the caption of Fig. 5b. This is a pure mode II problem. Local effects are similar to those of the ENF specimen and described in Fig. 3. The critical load for crack propagation obtained using the homogenized model is /cr IIC LP hE G  2 / 3 /h a . The solution coincides with predictions using discrete models and elementary beam theory. An accurate solution of the ELS specimen based on dimensional analysis, orthotropy rescaling and 2D finite element analyses has been obtained in [26] and has been recently confirmed by the analytical solution in [27]. For the E-glass-epoxy examined here, the exact solution is /cr IIC LP hE G   1/42 / 3 / (1 ( ) / )IIh a Y h a   , with 1/40.33, =2.56 and ( ) 0.4IIY      . The 2D solution accounts for the near tip deformations which are not accounted for in the homogenized solution since continuity conditions at the crack tip are imposed on global quantities. The relative error on the critical load for crack propagation is 4%, 2%,1% for / 10, 20, 40a h  , respectively. (a) (b) Figure 6: (a) Critical load for crack propagation versus load point displacement in the ELS specimen tested in [28]. Geometry and materials: L = 130 mm, 2h = 5.98 mm, 0 60a  mm, b = 20 mm. Cross-ply carbon/epoxy laminate with layup [(02/90)7/02//02/(90/02)7] and 144LE  GPa, 10TE  GPa, 4.2LT TTG G  GPa, 0.25LT  and 0.3TT  . Mode II fracture energy 0.5, 0.45IIC G N/mm. (b) Comparison of theoretical results using the homogenized model and 2D finite element models. Crack propagation in End-Load Split cross-ply laminate [(02/90)7/02//02/(90/02)7] The diagram in Fig. 6 shows the critical load for crack propagation versus load point displacement in the End Load Split (ELS) specimen tested in [28]. The material is a cross ply laminate with 46 unidirectional carbon plies of thickness 0.13 mm, layup [(02/90)7/02//02/(90/02)7] and elastic constants given in the caption of Fig. 6. The specimen has total length 160TL  mm , width 20B  mm and initial crack length 0 60a  mm. The theoretical results have been obtained by assuming that a part of the specimen, of length 30 mm, has been clamped so that the free length is 130L  mm (this is necessary to match 2D finite element results to the slope of the experimental linear branch; using TL would produce unacceptably large R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 283 displacements). The diagrams show the critical load for crack propagation versus load-point deflection obtained using the homogenized model, the compliance method, the crack propagation criterion II IICG =G and crack length control [18,25]. The results are compared with the experimental results in [28] and 2D finite element results obtained here. The finite element results have been obtained by using ANSYS, the VCCT technique, the crack propagation criterion II IICG =G and a crack length control. The finite element model uses plane strain isoparametric and quadrilateral 8-noded elements to mesh each layer of the beam with elements of equal size 0.075 0.13 mm in order to ensure convergence on both load point displacement and energy release rate. The results of the homogenized model agree well with the finite element predictions and the relative error on the critical load is always less than 3%; this error is due to an underestimation of the compliance in the homogenized model which neglects shear deformations along the delaminated portion of the specimen. The theoretical curves are presented for two values of the mode II fracture energy 0.45IIC G N/mm and 0.5IIC G N/mm. The first value better matches the experimental results at crack initiation and the second the post-critical branches; this difference is probably due to the presence of nonlinear cohesive mechanisms which are not reproduced by the linear elastic fracture mechanics solutions. The theoretical curves, obtained under crack length control, show snap-back instabilities, with virtual branches having negative slope, which are also evident in the almost vertical drops of some of the experimental curves obtained using displacement control. In the simulations in [28] the fracture toughness was assumed as 0.7IIC G N/mm; according to the LEFM analyses, this large value of IICG could be explained only if crack propagation occured under large scale bridging, with a large process zone, during the entire loading process. Further analyses are required to clarify this point; however propagation under extensive large scale bridging appears to be supported by the experimental measurements of the traction free crack lengths during propagation, which are shorter than those obtained using linear elastic fracture mechanics (not shown here). (a) (b) Figure 7: Critical load for crack propagation versus load point displacement in the homogeneous and isotropic cantilever beams with multiple delaminations shown in the insets (results for homogenized model in [18] and discrete layer model in [4]): (a) symmetrically positioned cracks; (b) asymmetrical positioned cracks (lengths and positions in the main text). Interaction effects of multiple delaminations Fig. 7 (after [18]) is used to highlight the capability of the homogenized model to analyze multiple delamination fracture and reproduce the effects of the interaction between delaminations. The dimensionless diagrams depict the critical load for the propagation of the cracks in the homogeneous and isotropic cantilever beams in the insets in Fig. 8 versus load-point displacement. The diagram (a) refers to a beam where ( 2 ) (1) 2h h  , ( 3 ) (1) 1h h  and the initial upper and lower crack lengths are 0 / 5.5Ua h  and 0 / 6La h  . The diagram (b) refers to a homogeneous isotropic beam where ( 2 ) (1) ( 3 ) (1) 1/ 3h h h h  and the initial upper and lower crack lengths are 0 / 5Ua h  and 0 / 4La h  . The delaminations R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 284 are assumed to propagate when the energy release rate, calculated using the compliance method, approaches the critical value, II IICG G . The results of the homogenized model in [18] are compared with the results of the discrete-layer cohesive- crack model with spring-contact in [4]. A local snap-through instability is observed in the diagram (a) when the upper crack starts to propagate in A and approaches the lower crack tip, in B. Then the load to propagate the crack must be increased, due to a shielding phenomenon, up to point C where the two cracks propagate together unstably. In the diagram (b) the lower crack, which is shorter, starts to propagate at the maximum load, point A; crack propagation is unstable and characterized by a snap-back instability up to point B. Then there is a sudden drop in the load, to point C, caused by a sudden amplification discontinuity. After point C the lower crack continues to propagate alone. DISCUSSION he homogenized model uses a zigzag kinematic approximation based on first order shear deformation theory (FSDT) which is enriched by local zigzag functions and interfacial jumps, as shown in Eqn. (1). The low order of the global model and the small number of unknowns, which are those of FSDT, implies that shear strains are constants in and between the layers in the intact regions of the beam and vanish in the delaminated regions. Following the approach which is commonly used for the structural low order theories, accurate shear stresses and strains can be obtained a posteriori through the imposition of local equilibrium. The diagrams in Figs. 3 and 4 show that the complex local fields which arise in multilayered wide plates with imperfect or fully debonded interfaces are accurately reproduced, also in very thick and highly anisotropic cases. However, stress recovery using the equation above is not very accurate in a finite element framework [8] and this limits the applicability of the model. The introduction of a shear correction factor, 44 5 / 6k  , in the equilibrium Eqns. (6) allows to account for the shear strains in the solution of the problem and the displacements of fully bonded beams are accurately predicted. On the other hand, due to the vanishing of the shear strains in the delaminated regions of the beams, the transverse displacements in Eqn. (3) and the equilibrium Eqns. (6) will only account for the bending contributions and this may affect the load re-distribution between different regions in statically indeterminate systems. This problem cannot be solved by increasing the order of the ESL theory used as global model and, in our opinion, could be solved only by modifying the zigzag formulation and leaving the displacement jumps as part of the unknowns of the problem, as it has been done for instance in [29]. Another approach to solve this problem could be through the use of the refined zigzag theory [7]; however, preliminary analyses in [20,30] highlight other still unresolved difficulties of the approach in dealing with in-plane discontinuities. The effects of neglecting shear deformations along the delaminated portions of the beam, however, are not relevant in all problems dominated by bending defomations (as in the fracture analyses of the cross ply laminated ELS specimen in Fig. 6). Another limitation of the homogenized approach presented here is related to the imposition of continuity and boundary conditions in terms of global quantities, either displacement or force and moment resultants. This implies that continuity is not satisfied at the local level and displacements and force sub-resultants within the single layers may differ. This typically occurs near the crack tip and clamped ends, within boundary regions whose size depends on the mismatch of the elastic constants of the layers and on the stiffness of the interfaces and is very small when the mismatch is small and the interfacial stiffness is very large or very small [20,22]. The model and analyses presented here are limited to mode II dominated problems, were transverse extensibility of the layers and interfacial openings may be neglected. The homogenized approach is currently being extended (work in progress [31]) to general mixed-mode problems using the extended formulation in [14], which uses first order shear, first order normal deformation theory as global model. CONCLUSIONS homogenized structural model based on a zigzag approach has been presented for the analysis of beams with imperfect interfaces and delaminations and layers oriented along the geometrical axes. The model enriches the displacement field of first order shear deformation theory, in order to introduce zigzag effects due to the layered architecture and interfacial sliding jumps due to the presence of soft interlayers, imperfections or delaminations. Piece-wise linear cohesive traction laws describe the constitutive behavior of the interfaces. Homogenization and variational techniques are used to define the local variables in terms of the global variables and derive the equilibrium equations of the problem. The model has been validated through various novel applications and using experimental results and 2D elasticity solutions. T A R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 285 Advantages and limitations of the approach have been discussed. Local stress and displacement fields are accurately predicted also in very thick and highly anisotropic beams and wide plates with continuous imperfect interfaces and delaminations. Crack propagation under mode II dominant conditions has been successfully reproduced in various fracture specimens and bend beams: unidirectionally reinforced End-Notched Flexure and End-Load Split specimens; End-Load Split specimen made of a 46-layers cross-ply laminate; isotropic cantilever beams with multiple delaminations interacting during propagation. The solution is analytical or semi-analytical, requires only in-plane discretization of the problem and the displacement unknowns are only three as in first order shear deformation theory. ACKNOWLEDGEMENTS upport by the U.S. Navy Office of Naval Research, grant N00014-17-1-2914 and by the Italian Dept. for University and Scientific and Technological Research, MIUR Prin15 project 2015LYYXA8. REFERENCES [1] Abrate, S., Di Sciuva, M. (2018).Multilayer models for composite and sandwich structures. Comprehensive Composite Materials II, (editors Zweben, C.H., Beaumont, P.W.R.), Elsevier, USA, pp. 399–425. [2] Williams, T.O., Addessio, F.L. (1997). A general theory for laminated plates with delaminations, Int. J. Solids Struct., 34(16), pp. 2003–2024, DOI: 10.1016/s0020-7683(96)00131. [3] Andrews, M.G., Massabò, R., Cavicchi, A., Cox, B.N. (2009). Dynamic interaction effects of multiple delaminations in plates subject to cylindrical bending, Int. J. Solids Struct., 46(9), pp. 1815–1833, DOI: 10.1016/j.ijsolstr.2008.11.027. [4] Andrews, M.G., Massabò, R. (2008). Delamination in flat sheet geometries with material imperfections and thickness variations, Compos. Part B Eng., 39(1), pp. 139–150, DOI: 10.1016/j.compositesb.2007.02.017. [5] Murakami, H. (1986). Laminated composite plate theory with improved in-plane responses, J. Appl. Mech. Trans. ASME, 53(3), pp. 661–666, DOI: 10.1115/1.3171828. [6] Di Sciuva, M. (1986). Bending, vibration and buckling of simply supported thick multilayered orthotropic plates: An evaluation of a new displacement model, J. Sound Vib., 105(3), pp. 425–442, DOI: 10.1016/0022-460X(86)90169-0. [7] Tessler, A., Di Sciuva, M., Gherlone, M. (2009). A refined zigzag beam theory for composite and sandwich beams, J. Compos. Mater., 43(9), pp. 1051–1081, DOI: 10.1177/0021998308097730. [8] Tessler, A. (2015). Refined zigzag theory for homogeneous, laminated composite, and sandwich beams derived from Reissner’s mixed variational principle, Meccanica, 50(10), pp. 2621–2648, DOI: 10.1007/s11012-015-0222-0. [9] Averill, R.C. (1994). Static and dynamic response of moderately thick laminated beams with damage, Compos. Eng., 4(4), pp. 381–395, DOI: 10.1016/S0961-9526(09)80013-0. [10] Groh, R.M.J., Tessler, A. (2017). Computationally efficient beam elements for accurate stresses in sandwich laminates and laminated composites with delaminations, Comput. Methods Appl. Mech. Eng., 320, pp. 369–395, DOI: 10.1016/j.cma.2017.03.035. [11] Di Sciuva, M. (1997). Geometrically nonlinear theory of multilayered plates with interlayer slips, AIAA J., 35(11), pp. 1753–1759, DOI: 10.2514/2.23. [12] Cheng, Z.Q., Jemah, A.K., Williams, F.W. (1996). Theory for multilayered anisotropic plates with weakened interfaces, J. Appl. Mech. Trans. ASME, 63, pp. 1019–1026, DOI: 10.1115/1.2787221. [13] Schmidt, R., Librescu, L. (1996). Geometrically nonlinear theory of laminated anisotropic composite plates featuring interlayer slips, N J Math Game Theory Algebr., 5, pp. 131–137. [14] Massabò, R., Campi, F. (2015). Assessment and correction of theories for multilayered plates with imperfect interfaces, Meccanica, 50(4), pp. 1045–1071, DOI: 10.1007/s11012-014-9994-x. [15] Shu, X., Soldatos, K.P. (2001). An accurate stress analysis model for angle-ply laminates with weakly bonded layers, Acta Mech., 150(3-4), pp. 161–178, DOI: 10.1007/BF01181809. [16] Eijo, A., Oñate, E., Oller, S. (2013). A numerical model of delamination in composite laminated beams using the LRZ beam element based on the refined zigzag theory, Compos. Struct., 104, pp. 270–280, DOI: 10.1016/j.compstruct.2013.04.035. [17] Groh, R.M.J., Weaver, P.M., Tessler, A. (2015). Application of the refined zigzag theory to the modeling of delaminations in laminated composites, NASA/TM-2015-218808, pp. 1–22, DOI: 10.13140/RG.2.1.3147.0804. S R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 286 [18] Massabo, R, Darban, H. (2019). Mode II dominant fracture of layered composite beams and wide-plates: a homogenized structural approach, Eng. Fract. Mech., 213, pp. 280–301, DOI: 10.1016/j.engfracmech.2019.03.002. [19] Pelassa, M., Massabò, R. (2015). Explicit solutions for multi-layered wide plates and beams with perfect and imperfect bonding and delaminations under thermo-mechanical loading, Meccanica, 50(10), pp. 2497–2524, DOI: 10.1007/s11012-015-0147-7. [20] Darban, H., Massabò, R. (2018). A homogenized structural model for shear deformable composites with compliant interlayers, Multiscale Multidiscip. Model. Exp. Des., 1, pp. 269–290, DOI: 10.1007/s41939-018-0032-x. [21] Massabò, R. (2017). Propagation of Rayleigh-Lamb waves in multilayered plates through a multiscale structural model, Int. J. Solids Struct., 124, pp. 108–124, DOI: 10.1016/j.ijsolstr.2017.06.020. [22] Massabò, R. (2014). Influence of boundary conditions on the response of multilayered plates with cohesive interfaces and delaminations using a homogenized approach, Frat. Ed Integrita Strutt., 8(29), pp. 230–240, DOI: 10.3221/IGF-ESIS.29.20. [23] Andrews, M.G., Massabò, R. (2007). The effects of shear and near tip deformations on energy release rate and mode mixity of edge-cracked orthotropic layers, Eng. Fract. Mech., 74(17), pp. 2700–2720, DOI: 10.1016/j.engfracmech.2007.01.013. [24] Madhukar, M.S., Drzal, L.T. (1992). Fiber-matrix adhesion and its effect on composite mechanical properties: IV. mode I and mode II fracture toughness of graphite/epoxy composites, J. Compos. Mater., 26(7), pp. 936–968, DOI: 10.1177/002199839202600701. [25] Monetto, I., Campi, F. (2017). Numerical analysis of two-layer beams with interlayer slip and step-wise linear interface law, Eng. Struct., 144, pp. 201–209, DOI: 10.1016/j.engstruct.2017.04.010. [26] Suo, Z., Bao, G., Fan, B., Wang, T.C. (1991). Orthotropy rescaling and implications for fracture in composites, Int. J. Solids Struct., 28(2), pp. 235–248, DOI: 10.1016/0020-7683(91)90208-W. [27] Ustinov, K.B., Massabò, R., Lisovenko, D.S. (2019). Orthotropic strip with central semi-infinite crack under arbitrary loads applied far apart from the crack tip. Analytical solution, Submitted. [28] Tawk, I., Navarro, P., Ferrero, J.F., Barrau, J.J., Abdullah, E. (2010). Composite delamination modelling using a multi- layered solid element, Compos. Sci. Technol., 70(2), pp. 207–214, DOI: 10.1016/j.compscitech.2009.10.008. [29] Kim, H.S., Chattopadhyay, A., Ghoshal, A. (2003). Characterization of delamination effect on composite laminates using a new generalized layerwise approach, Comput. Struct., 81(5), pp. 1555–1566, DOI: 10.1016/S0045-7949(03)00150-0. [30] Flores, F.G., Oller, S., Nallim, L.G. (2018). On the analysis of non-homogeneous laminates using the refined zigzag theory, Compos. Struct., 204, pp. 791–802, DOI: 10.1016/j.compstruct.2018.08.018. [31] Massabò, R. (2019). Mixed-mode delamination in layered beams using a homogenized structural approach, internal report. University of Genova, Italy. APPENDIX Relevant constitutive, equilibrium and compatibility equations in layer k: ( ) ( ) ( ) 22 22 22 k k kC  ; ( ) ( ) ( )23 23 232 k k kG  ; ( ) ( ) 22 2 23 3, , 0 k k   ; (7)     22 2 2, k k v  ;      23 2 3 3 22 , , k k kv v   Force and moment resultants:    3 1 3 ( ) 22 22 22 3 3 2 2 2 2 2 2 22 2 2 2 2 1 ˆ, 1, ; ( ) ( ) ( ) , ( ) ( ) k k n xb k b z zS gx k N M x dx Q x Q x Q x M x x         (8) with: R. Massabò et alii, Frattura ed Integrità Strutturale, 51 (2020) 275-287; DOI: 10.3221/IGF-ESIS.51.22 287 3 3 1 1 3 3 3 1 3 ( ) ( ) 2 23 3 22 22 22 3 1 1 1 ( ) (1) ( 1) (1) 2 23 23 3 2 23 23( 1) ( ) ( 1) ( ) 1 1 123 23 23 23 ; 1 1 1 1 ˆ ; 1 k k k k k k n nx xb k zS k k Sx x k k n k kxz k k i i j jx k i j Q dx M R dx Q G dx G G G G G G                                          1 2 1 ˆ n k k v        (9) Coefficients of the equilibrium equations:   3 3 1 1 3 3 3 1 3 ( ) 2 ( ) 2 22 3 22 3 22 22 3 1 1 2 2( ) ( 1) (1) 44 23 22 3 3 23 23 ( 1) ( ) 1 1 23 23 ( ) ( ) 1 1 1 1 , 1 k k k k k k n nx xrS k r k S k k S Sx x k k n kxP k k S k S k j jx k jS C C x R dx C C R dx C G R dx C G G K G G                                  1 1 n k    (10) << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice