Interfacial partial debonding and its in uence on the elasticity of a two-phase composite

Size: px
Start display at page:

Download "Interfacial partial debonding and its in uence on the elasticity of a two-phase composite"

Transcription

1 Mechanics of Materials 32 (2000) 695±709 Interfacial partial debonding and its in uence on the elasticity of a two-phase composite S.F. Zheng, M. Denda, G.J. Weng * Department of Mechanical and Aerospace Engineering, Rutgers University, New Brunswick, NJ 08903, USA Received 8 July 1999; received in revised form 4 May 2000 Abstract A dual e ective-medium and nite-element study is carried out to examine the in uence of interfacial partial debonding on the elastic sti ness of a two-phase composite containing aligned elliptic bers. In the e ective-medium approach double debonding on the top and bottom of the elliptic interface is considered, but in the nite-element approach both double and single debonding are examined. The e ective-medium approach makes use of the concept of a ctitious ber whose load-carrying capacity is taken to be lost in the debonding direction but remains intact in the transverse direction. The nite-element analysis allows one to examine further the in uence of debonding angle on the sti ness of the composite and it also provides the needed magnitude for this angle at which the e ective-medium approach applies. It is found that the angle has to be su ciently wide, but not so wide as to lead to a potential complete debonding. Such an angle also increases when the bers become more ribbon-like, but it decreases with increasing volume concentration. The basic assumption that the ctitious ber has zero, or very low, load-carrying capacity along the debonding direction is also veri ed by the nite-element results. For both double and single debonding the Young's modulus of the composite along the debonding direction can be signi cantly reduced as the debonding angle increases. The local stress distributions inside the ber and the matrix are also illustrated for both types of debonding. Ó 2000 Elsevier Science Ltd. All rights reserved. Keywords: Interfacial partial debonding; Elastic sti ness; Two-phase composites 1. Introduction In a recent paper Zhao and Weng (1997) developed a simple concept of ctitious inclusions to mimic the load-transfer capacity of two kinds of partially debonded inclusions in a composite. The rst kind involves partial debonding on the top and bottom of the oblate inclusions while the * Corresponding author. Tel.: ; fax: address: weng@jove.rutgers.edu (G.J. Weng). second kind involves partial debonding on the lateral surface of the prolate ones. A ctitious inclusion was used to represent the smeared e ect of the original inclusion and the surrounding interface cracks. When debonding is su ciently wide, the load-carrying capacity of the ctitious inclusion was taken to be zero in the debonding direction, but its ability to carry the load in the bonded direction was taken to remain present. As such, the isotropic inclusion and the interfacial cracks together were replaced by the ctitious inclusion which now becomes perfectly bonded to the matrix but with transversely isotropic moduli. Following /00/$ - see front matter Ó 2000 Elsevier Science Ltd. All rights reserved. PII: S ( 0 0 )

2 696 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 this concept, the e ective moduli of the transversely isotropic composite with such partially debonded inclusions were estimated using a well established micromechanical model for a twophase composite containing aligned, transversely isotropic spheroidal inclusions. The advantage of such an approach is that it delivers simple ± albeit approximate ± explicit results for the e ective moduli of an otherwise rather complicated composite. What remains unclear is how wide the partial debonding has to be in order for the concept to apply. Apparently, the loadcarrying capacity of the ctitious inclusion can not be zero in the debonding direction if the debonding occurs only over a small portion of the interface, and that its load-transfer ability in the perpendicular direction can not be retained if the debonding extends over the entire interface. It is thus highly desirable to have this question answered. Moreover, the extent of the partial debonding for this concept to hold is also likely to depend on the inclusion shape as a platelet inclusion will need a wider debonding than a spherical one, and the overall sti ness of the debonded composite will also depend on the extent of debonding. The e ective-medium method as usual does not deliver this sort of information and it also gives no local stress and strain elds which can serve to detect the subsequent damage process. These questions can perhaps be best answered by a companion nite-element analysis. It is with this perspective that the present study was undertaken. Here we are concerned with a two-dimensional problem involving partially debonded, aligned elliptic bers. The unit cells are shown in Figs. 1(a) and (b) for double debonding on the top and bottom, and for single debonding on the top, respectively. Fibers are taken to be aligned along x 1, and the debonded direction along x 3. The cross-sectional aspect ratio or the thickness to the width ratio, t=w, is represented by a, so that by varying a the e ect of the ber shape can also be taken into account. In both cases the composite as a whole is orthotropic. The e ective-medium method based on the concept of ctitious bers can only be applied to con guration (a), for in (b) it is not realistic to expect the smeared region to lose its load-carrying capacity. Moreover, for the e ective-medium method to hold the cross-sectional aspect ratio should be limited to the range a < 1, as for a > 1 ± especially when it is very large, the shear stress (shear lag) may still transfer the load from the matrix to the ber. For the - nite-element method both con gurations will be examined so that a comparison between these two types of partial debonding can be made. In both gures the extent of partial debonding is speci ed by the debonding angle /, with / ˆ 0 corresponding to perfect bonding and / ˆ 90 to complete debonding. It is intuitively intriguing to see how the sti ness of the composite is weakened as the debonding angle increases. The outer cell of the unit cell is chosen to have an identical shape as the inclusion for the reason that the e ective-medium theory to be used corresponds to such a microgeometry. A rectangular outer cell would lead to a composite with a periodic structure and this is not the case for the model to be adopted. The radii of the inner and outer cells are chosen in accordance with the volume concentration of the bers, and when the aspect ratio of inclusions becomes one, the unit cell re- Fig. 1. Unit cells for: (a) double interfacial debonding; (b) single interfacial debonding.

3 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± turns to the composite cylinder assemblage. The unit cell implies that every ber has an identical damage. This is a highly idealized situation and the results must be viewed in this light. We also note in passing that the elastic eld of an interfacial arc crack has been derived by England (1966) and Toya (1974) for a circular ber (a ˆ 1) in an in nitely extended matrix, and the result has been applied by Ju (1991) to develop a damage model for a ber-reinforced composite. For an elliptic interfacial crack in a nite matrix, however, the problem remains unsolved. In what follows the concept of zero load-transfer capacity in the debonding direction will be used to establish the orthotropic property of the ctitious ber. Then based on an e ective-medium method the nine e ective moduli of the partially debonded composite will be calculated as a function of the cross-sectional aspect ratio a and volume concentration c 1 of the elliptic bers. The nite-element study will subsequently be presented. Comparison between the two for double debonding, and the additional results by the nite element for single debonding, will be given at the end. 2. Properties of the ctitious bers and the e ectivemedium method In the two-phase composite the elliptic bers will be referred to as phase 1 and the matrix as phase 0. For simplicity both phases will be taken to be isotropic, with the bulk and shear moduli denoted by the Greek j r and l r, respectively, for phase r. Its plane-strain bulk modulus will be denoted by the Roman k r, and volume concentration by c r. The orthogonal symmetry of the debonding con guration as depicted in Fig. 1(a) necessitates that the ctitious ber be orthotropic. For an orthotropic material its elastic constitutive equations carry the form r 11 ˆ c c c 13 ; r 22 ˆ c c c 23 ; r ˆ c c c ; r 32 ˆ 2c ; r 31 ˆ 2c ; r 12 ˆ 2c ; 1 in terms of stress r ij, strain ij, and elastic constants c ij ( ˆ c ji ). We shall denote the nine elastic constants of the ctitious ber by c 1 ij. When interfacial debonding or the debonding angle / is su ciently wide, we may assume that the load-carrying capacity of the ctitious ber be zero in the x 3 -direction so that c 1 ˆ c1 32 ˆ c1 31 ˆ 0: 2 No shear stress on the x 3 -surface of the ctitious ber is taken to endure either, with c 1 44 ˆ c1 55 ˆ 0: 3 The r 12 -loading is on the bonded edge of the elliptic cylinder and is taken to be una ected by the partial debonding c 1 66 ˆ l 1: 4 There remain c 1 11, c1 22, and c1 12 to be determined, all in the x 1 ±x 2 plane. These three constants could still be somewhat weakened by the partial loss of load-transfer ability along x 1 and x 2. We expect their values to be less than those of the original ber, which carries c 11 ˆ c 22 ˆ c ˆ k 1 l 1 and c 12 ˆ c 13 ˆ c 23 ˆ k 1 l 1. These values will be determined by arti cially imposing r ˆ 0 in the constitutive equations of the isotropic ber so that ˆ k 1 l 1 k 1 l ; and consequently r 11 ˆ 4k 1l 1 k 1 l l 1 k 1 l 1 k 1 l 1 22 ; r 22 ˆ 2l 1 k 1 l 1 k 1 l k 1l 1 k 1 l 1 22 : The three constants are then identi ed as c 1 11 ˆ c1 22 ˆ 4k 1l 1 k 1 l 1 ; c 1 12 ˆ 2l 1 k 1 l 1 k 1 l 1 : The reductions in c 1 11, c1 22, and c1 12 from the original k 1 l 1 and k 1 l 1 are found to be all k 1 l 1 2 = k 1 l 1.

4 698 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 With this scheme, the partially debonded isotropic bers are replaced by the perfectly bonded orthotropic bers of the same cross-section. The original problem ± which in essence is a three-phase one involving bers, matrix, and interface cracks ± has been replaced by a two-phase problem with only the ctitious bers and the matrix. In this context it is important to recognize that the bers in the threephase problem are not the same as the ctitious bers in the two-phase problem. Rather, each ctitious ber now represents the smeared medium of the original ber and the two interface cracks. We shall represent these orthotropic constants by the tensor L 1. It is then straightforward to nd the e ective moduli of such a two-phase composite by some well established e ective-medium method. For the present problem the simplest one is likely to be Willis' (1977) approach, by taking the elastic moduli of the comparison material as those of the matrix, denoted by the tensor L 0. It follows that the e ective moduli tensor is given by L ˆ c 1 L 1 A 1 c 0 L 0 I c 1 A 1 c 0 I 1 ; 8 where I is the fourth-order identity tensor, and A 1 is the strain concentration tensor of a single ctitious ber embedded in the in nitely extended matrix. It is given by A 1 ˆ I SL 1 0 L 1 L 0 Š 1 ; 9 in terms of Eshelby's (1957) S-tensor for an elliptic inclusion. Its components can be found in Mura (1987). As proved in Weng (1992), these moduli are identical to those derived by the Mori and Tanaka (1973) method. Since Willis' model involves a distribution function that is identical to the inclusion shape, and the outer cell of the unit cell also represents the distribution function, its microgeometry is exactly what is depicted in Fig. 1. It then makes sense to compare the nite-element calculations based on such a unit cell with the theoretical results. In order to show how the double debonding on the top and bottom of the elliptic cylinders would a ect the overall sti ness of the composite, we have applied this approach to calculate the nine independent elastic constants of the partially debonded composite, using the Young's modulus and Poisson's ratio of silicon carbide and aluminum (Arsenault, 1984; Nieh and Chellman, 1984): Silicon carbide : E 1 ˆ 490 GPa; m 1 ˆ 0:17; Aluminum : E 0 ˆ 68:3 GPa; m 0 ˆ 0:: 10 The results are shown in Fig. 2, normalized by the properties of the matrix. Here the dependence of the three Young's moduli, three shear moduli, and three Poisson's ratios on the volume concentration c 1, of the debonded bers are illustrated. In each gure the in uence of the ber cross-sectional aspect ratio a, ranging from 1.0 to 0.01, are also given. For comparison, the corresponding moduli with perfectly bonded bers are plotted in Fig. 3. A quick glance over these gures indicates that E 11, E 22, l 12, and m 12 are almost una ected by the partial debonding, but E, l 13, and l 23 are signi cantly reduced. Partial debonding has also changed the Poisson's ratios m 13 and m 32, which now can go out of the range of the constituent phases. A closer examination on the in uence of the cross-sectional shape further reveals that, for the Young's modulus E and shear moduli l 13 and l 23, the reduction of the sti ness is far more signi cant when bers take the shape of thin ribbons (a ˆ 0.01). The moduli displayed here must be viewed only under the condition of a crack opening mode. Consequently E is a tensile Young's modulus, not a compressive one. 3. Description of the nite-element analysis As the most important property change is the transverse Young's modulus E, our nite-element investigation will focus on the tensile loading along x 3. We consider a plane-strain problem for the unit cell subjected to a boundary traction giving rise to a uniform stress r (an overbar signi es the volume, or rather the area, average over the unit cell here). The traction boundary condition gives rise to a lower bound value for the e ective moduli; it was chosen over the displacement boundary condition for the e ective-

5 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± Fig. 2. The nine e ective moduli of the doubly debonded composite as calculated by the e ective-medium method. medium approach also provides a lower bound estimate (even though not always the same type of lower bound) when the matrix is softer than the inclusion in the perfectly bounded case. Since the direction of the loading is perpendicular to the debonding con guration, contact of the debonded surfaces will not occur when the crosssectional aspect ratio is small and the debonding

6 700 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 Fig. 3. The nine e ective moduli of the perfectly bonded composite as calculated by the e ective-medium method. angle is far from 90. But in every step of the way we have checked for the contact condition and whenever such a contact was detected the calculation for the speci c debonding angle / was halted. Indeed such a contact tends to occur under the in-plane shear r 23. The treatment of such a contact problem is essentially nonlinear and is not included here.

7 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± The transverse Young's modulus is calculated by E ˆ r ; 11 where is the overall average strain in the x 3 - direction. In a cracked body, the contributions to the overall strain ± as pointed out by Benveniste (1986) and others ± must include the displacement jump across the crack surfaces. For the present problem as depicted in Figs. 1(a) and (b), it comes from the average strains of matrix (phase 0) and the ber (phase 1), with their weighted sum now denoted by 0 1, and the displacement jump, u 3 Š, across the debonded surfaces in the x 3 -direction ˆ c 0 0 c X Z u 3 Šn 3 ds k V k S k ˆ 0 1 uš ; 12 where n 3 is the component of the unit normal to the debonding surface in the x 3 -direction, S k the kth debonded surface and V is the volume of the unit cell. The summation is taken on all the debonded surfaces in the unit cell. The symbol uš represents the last term from the displacement jump across the interfaces. It turns out that, when the extent of debonding is small, uš contributes little to the overall strain, but when it is wide this term will dominate. Typical nite-element meshes used for both double and single debonding are shown in Figs. 4(a) and (b), respectively. Here the illustrations are given for the aspect ratio a ˆ 1/4, ber concentration c 1 ˆ 30%, and debonding angle / ˆ 67 (for reasons to become clear later). The debonding along the ber±matrix interface is simulated by a sharp crack-like notch introduced by removing a notched region from the matrix adjacent to the interface. The maximum opening of the notch d max is determined to be 1% of the half debonding length. Since our primary concern is with the overall modulus E ± not with the singular behavior of the stress eld at the crack tip ± the debonding area is represented by a sharp cracklike notch, rather than by a crack. To justify this approximation we have introduced an alternative representation of the debonding by a slit of constant width d max with rounded tips. This model eliminates the singular behavior at all. Our numerical results have shown that these two models di er only in the local stress distribution near the tips and practically have no e ects on the average strain of the unit cell. Even the local notch tip stress distributions do not di er much if the radius of the rounded notch is small enough, such as less than 1% of the half debonding length. We have used ANSYS nite-element program. Eight-node isoparametric quadratic elements have been used except at the debonding tip, where triangular quarter-node singular elements have been employed. Each of the singular element surrounding the debonding tip is a six-node triangular element with its radial length less than one hundredth of the half debonding length. The next layer of elements consists of quadratic elements, edge lengths of which are magni ed by 120% of those for the inner singular elements. The mesh, re ned in the debonding tip region, gradually becomes coarse as we move away from it. Strictly speaking, the quarter-node p singular element that numerically generates the 1= r singular stress behavior at the tip, is valid only for the crack tip in a homogeneous elastic material. Since the region of the oscillatory singularity is con ned to the near interface crack tip region and the envelope p of this singularity is still of the order of 1= r, use of the quarter-node is not a bad approximation even if we replace the crack by a notch whose singularity is slightly weaker than that of the crack. Exploiting the symmetry of the geometry and loading, we have analyzed one quarter of the unit cell for the double debonding case and the half for the single debonding. 4. Dependence of the tensile Young's modulus E on the debonding angle / The variations of tensile Young's modulus E as a function of the debonding angle / are shown in Figs. 5(a)±(c) for three cross-sectional aspect ratios a ˆ 1=10, 1=4, and 1, respectively. In each gure two volume concentrations of the bers c 1 ˆ 30% and 10% have been calculated. Fig. 5(a) indicates that, when the aspect ratio is low ± namely when

8 702 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 Fig. 4. Typical nite-element mesh for the unit cell under: (a) double debonding (one quadrant); (b) single debonding (one half). the bers are ribbon-like ± the tensile modulus remains una ected when the debonding angle is small, but the reduction becomes very sharp after / reaches a certain critical angle. For instance with c 1 ˆ 30%, E remains at until the debonding angle is close to 30, and with c 1 ˆ 10% it remains almost una ected until it reaches 60. As the crosssectional aspect ratio increases, such a sustained region also shortens, and when a ˆ 1 the variation of E is almost instant as / increases. Such a dependence is easy to conceive as, when a ˆ 1=10, most portion of the ber would be between the perfectly bonded interfaces if the debonding angle is small, but such is not the case with circular bers. On each of these curves a heavy dot has also been marked; theses points correspond to the E calculated by the e ective-medium method, as shown in Fig. 2. The e ective-medium method is based on the assumption that the debonding angle / is su ciently wide, but it does not provide the information on how wide this angle has to be. These calculations indicate that, with a ˆ 1=10, this angle is about 72 if c 1 ˆ 30%, and it is close to 78 if c 1 ˆ 10%. With the more rounded a ˆ 1=4, these angles are somewhat lower, about 67 and 74, respectively (the choice of a ˆ 1=4, c 1 ˆ 30%, and / ˆ 67 in Fig. 4 is now explained). The greater angle needed for the more ribbon-like a ˆ 1=10 is understandable as, again, when the

9 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± Fig. 5. Dependence of the transverse Young's modulus E on the debonding angle /.

10 704 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 Fig. 6. The variation of the debonding angle / in the e ectivemedium theory as the ber concentration c 1 increases. angle is small the inclusion can still transfer substantial amount of stress to the matrix. The heavy dots in Fig. 5(c) lie on the extended dashed lines, which in essence can not be attained as our niteelement calculations indicated that, when / is so large, the surfaces of the circular cylinder and the matrix would come to contact under the tensile loading. Such a possibility has also been noted by Toya (1974). This implies that, for a ˆ 1, the concept of ctitious ber is not valid and the E -curve in Fig. 2 for a ˆ 1 should be discarded. This all the more reiterates that the concept of ctitious bers should be restricted to the range a < 1, as stated earlier. The fact that the two dotted points in (a) and (b) do not lie vertically implies that the debonding angle / associated with the e ective-medium theory is somewhat dependent upon the ber concentration. Such a dependence is shown in Fig. 6 over a continuous range of c 1. These results display an almost linear correlation, and further indicate that, for the e ective-medium theory to hold, the debonding angle indeed has to be su ciently wide, typically ranging around 60±80. The increasing signi cance of the displacement jump uš as the debonding angle increases is illustrated in Fig. 7, for a ˆ 1=10 and c 1 ˆ 30%. It is evident from Fig. 7(a) that when / > 57 the contribution to the overall strain by the displacement jump ± represented by the uš -curve ± has exceeded that by the weighted sum of the average strains of the bers and the matrix, and that at / ˆ 72 ± which corresponds to the e ective-medium theory ± the displacement jump accounts for 90% of the total strain. To verify the validity of the assumption that the load-carrying capacity of the ctitious ber is zero in the debonding direction, we have plotted in Fig. 7(b) the ratio of the jump strain uš to the inclusion strain 1. Here uš, written without an overbar, is calculated as uš in Eq. (12) but divided by the volume V 1 of the ber ( uš uš ˆ / c 1 ) and, together with 1, form the strain of the ctitious ber. Again, as / increases the contribution of the displacement jump to the strain of the smeared region also increases, and at / ˆ 72 where the e ective-medium theory is intended for application, the strain from the cracks is as high as 158-times of that of the ber. As the crack regions have zero moduli, this means that the sti ness of the ctitious ber in the debonding direction is only 1=158th, or less than 1%, of the sti ness of the real ber. As such, the basic assumption that the ctitious ber has no load-carrying capacity in the debonding direction is justi ed. With single debonding as sketched in Fig. 1(b), the reduction of the sti ness as a function of debonding angle / is shown in Fig. 8, for the same three aspect ratios. For comparison the corresponding results for double debonding are also depicted as dashed lines. Both classes of curves show similar characteristics, such as an incubation period for / when a is small. Even though the overall appearance of these two groups of curves does not seem to di er signi cantly, the magnitudes of E at a given / within the sharply a ected region can be quite di erent. For instance at / ˆ 60 and when a ˆ 1=10 and c 1 ˆ 30%, the overall sti ness under single debonding is 38% higher than that under double debonding. The local stress distributions for r x 2 ; x 3 under double debonding at c 1 ˆ 30% are shown for a ˆ 1=10 and 1=4 in Figs. 9(a) and (b), respectively. The debonding angles are deliberately chosen at 72 and 67 here so that the plotted results would also shed some light for the elds based on the e ective-medium method. The x 2 -x 3 plane at

11 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± Fig. 7. The relative contribution of the displacement jump over the interface cracks to: (a) the overall strain of the composite; (b) the overall strain of the smeared ber and interface cracks region. the bottom of each gure re ects the geometry of the rst quadrant of Fig. 1(a), with the ber±matrix interface marked by the sign ``interface'' on both axes. Looking along the x 3 -axis, it is apparent that the stress in the ber is zero. Continuous increase along the x 3 -axis takes one into the matrix region and the stress begins to increase from zero towards a greater quantity. The stress along the x 2 - axis shows a sustained region of zero stress also inside the ber. Its value, however, begins to increase as one moves below the notch tip. Similar features are also present in Fig. 8(b), with a ˆ 1=4. It is evident that stress distribution in the ber and the matrix is quite heterogeneous due to the presence of the interface cracks and elastic inhomogeneity. It may be reminded that, while the stress in the ber is not zero, this is not to be interpreted as a direct violation to the basic concept of the ctitious ber. Instead as we did in Fig. 7(b) such a validity should be examined from the strain consideration. The local stress elds under single debonding are shown in Figs. 10(a) and (b), chosen for the same ber aspect ratios and debonding angles as in Figs. 9(a) and (b) so as to render a comparison. These illustrations are displayed for the right half of the con guration in Fig. 1(b) as the stress eld is no longer symmetric with respect to the x 2 -axis. The stress r x 2 ; x 3 along the x 3 -axis is seen to be zero on the top but nite at the bottom of the interface. The lower portion of the ber ± which remains bonded to the matrix ± now carries a noticeable magnitude of tensile stress everywhere, leading to a higher sti ness for the overall composite as re ected in Fig Concluding remarks In this paper we have presented a dual e ectivemedium and nite-element study on the overall elastic behavior of a partially debonded ber-re-

12 706 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 Fig. 8. Comparison between the e ects of single debonding and double debonding.

13 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± Fig. 9. Local stress distribution in the unit cell for r x 2 ; x 3 under double debonding. inforced composite. The e ective-medium method is aimed at double debonding, while the nite-element analysis is concerned with both double and single debonding and it has been further used to examine how the debonding angle / a ects the reduction of the overall sti ness of the composite.

14 708 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695±709 Fig. 10. Local stress distribution in the unit cell for r x 2 ; x 3 under single debonding. The e ective-medium method is based on the concept of ctitious bers whose load-carrying capacity is taken to be lost in the debonding direction and whose orthotropic properties are then derived from the isotropic ones of the original - bers. The debonded composite is then treated as a perfectly bonded two-phase composite. The niteelement analysis on the other hand treats the problem as a three-phase material including the bers, matrix, and the interface cracks, the extent of the latter being represented by the debonding angle /. A detailed examination on the nite-element results indicates that, for ribbon-like bers, the debonding angle has very little e ect on the

15 S.F. Zheng et al. / Mechanics of Materials 32 (2000) 695± reduction of the transverse Young's modulus E until it reaches a su ciently large value, but that for circular bers its e ect is immediate. As the debonding angle increases, the interfacial displacement jump is found to contribute the major portion of the overall strain. The nite-element results also provide the precise / at which the effective-medium method can be best applied. For debonding angles signi cantly di erent from this one, the ctitious- ber approach is not recommended. The validity of the basic assumption that the ctitious ber has zero, or very low, load-carrying capacity in the debonding direction has also been veri ed. The results for single debonding suggest that, while the general trend of the E variation remains similar to the double debonding case, the di erence in their magnitudes can be signi cant for certain critical range of debonding angle. The local stress elds for both debonding types have also been illustrated. Acknowledgements This work was supported by the National Science Foundation, under Grant CMS , and by the O ce of Naval Research, under Grant N J References Arsenault, R.J., The strengthening of aluminum alloy 6061 by ber and platelet silicon carbide. Mater. Sci. Eng. 64, 171±181. Benveniste, Y., On the Mori±Tanaka's method in cracked bodies. Mech. Res. Comm. 13, 193±201. England, A.H., An arc crack around a circular elastic inclusion. ASME J. Appl. Mech., 637±640. Eshelby, J.D., The determination of the elastic eld of an ellipsoidal inclusion and related problems. Proc. Roy. Soc. London. Ser. A 241, 376±396. Ju, J.W., A micromechanical damage model for uniaxially reinforced composites weakened by interfacial arc microcracks. ASME J. Appl. Mech. 58, 923±930. Mori, T., Tanaka, K., Average stress in the matrix and average elastic energy of materials with mis tting inclusions. Acta Metall. 21, 571±574. Mura, T., Micromechanics of Defects in Solids, second ed. Martinus Niho, Dordrecht. Nieh, T.G., Chellman, D.J., Modulus measurements in discontinuous reinforced aluminum composites. Scr. Metall. 8, 925±928. Toya, M., A crack along the interface of a circular inclusion embedded in an in nite solid. J. Mech. Phys. Solids 22, 325±348. Weng, G.J., Explicit evaluation of Willis' bounds with ellipsoidal inclusions. Int. J. Eng. Sci. 30, 83±92. Willis, J.R., Bounds and self-consistent estimates for the overall properties of anisotropic composites. J. Mech. Phys. Solids 25, 185±202. Zhao, Y.H., Weng, G.J., Transversely isotropic moduli of two partially debonded composites. Int. J. Solids Struct. 34, 493±507.

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without

More information

The overall elastic moduli of orthotropic composite and description of orthotropic damage of materials

The overall elastic moduli of orthotropic composite and description of orthotropic damage of materials International Journal of Solids and Structures 37 (2000) 6755±6771 www.elsevier.com/locate/ijsolstr The overall elastic moduli of orthotropic composite and description of orthotropic damage of materials

More information

MODELING OF THE ANISOTROPIC ELASTIC PROPERTIES OF PLASMA-SPRAYED COATINGS IN RELATION TO THEIR MICROSTRUCTURE

MODELING OF THE ANISOTROPIC ELASTIC PROPERTIES OF PLASMA-SPRAYED COATINGS IN RELATION TO THEIR MICROSTRUCTURE Acta mater. 48 (2000) 1361±1370 www.elsevier.com/locate/actamat MODELING OF THE ANISOTROPIC ELASTIC PROPERTIES OF PLASMA-SPRAYED COATINGS IN RELATION TO THEIR MICROSTRUCTURE I. SEVOSTIANOV{ and M. KACHANOV

More information

Finite Element of Composite Materials with Partial Debonding along Boundary. Kohei HIGUCHI, Tetsuo IWAKUMA

Finite Element of Composite Materials with Partial Debonding along Boundary. Kohei HIGUCHI, Tetsuo IWAKUMA Vol. 8 (2005 8 ) Finite Element of Composite aterials with Partial Debonding along Boundary Kohei HIGUCHI, Tetsuo IWAKUA JFE PhD 980-8579 6-6-06 Damages in composite materials are mostly governed by interfacial

More information

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model

More information

Surface crack subject to mixed mode loading

Surface crack subject to mixed mode loading Engineering Fracture Mechanics 65 (2000) 1±14 www.elsevier.com/locate/engfracmech Surface crack subject to mixed mode loading M.Y. He a, J.W. Hutchinson b, * a Materials Engineering Department, University

More information

A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds

A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds International Journal of Solids and Structures 42 (2005) 4179 4191 www.elsevier.com/locate/ijsolstr A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial

More information

A simple plane-strain solution for functionally graded multilayered isotropic cylinders

A simple plane-strain solution for functionally graded multilayered isotropic cylinders Structural Engineering and Mechanics, Vol. 24, o. 6 (2006) 000-000 1 A simple plane-strain solution for functionally graded multilayered isotropic cylinders E. Pan Department of Civil Engineering, The

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

SPLIT SINGULARITIES: STRESS FIELD NEAR THE EDGE OF A SILICON DIE ON A POLYMER SUBSTRATE

SPLIT SINGULARITIES: STRESS FIELD NEAR THE EDGE OF A SILICON DIE ON A POLYMER SUBSTRATE Acta mater. Vol. 47, No. 1, pp. 67±76, 1999 # 1999 Acta Metallurgica Inc. Published by Elsevier Science Ltd. All rights reserved Printed in Great Britain PII: S1359-6454(98)00345-0 1359-6454/99 $19.00

More information

Three-dimensional thermoelastic deformations of a functionally graded elliptic plate

Three-dimensional thermoelastic deformations of a functionally graded elliptic plate JCOMB Composites: Part B () 9 6 www.elsevier.com/locate/compositesb Three-dimensional thermoelastic deformations of a functionally graded elliptic plate Z.-Q. Cheng a, R.C. Batra b, * a Department of Modern

More information

A modified quarter point element for fracture analysis of cracks

A modified quarter point element for fracture analysis of cracks ndian Journal of Engineering & Materials Sciences Vol. 14, February 007, pp. 31-38 A modified quarter point element for fracture analysis of cracks Sayantan Paul & B N Rao* Structural Engineering Division,

More information

Micromechanical simulation of the failure of ber reinforced composites

Micromechanical simulation of the failure of ber reinforced composites Journal of the Mechanics and Physics of Solids 48 (2000) 621±648 www.elsevier.com/locate/jmps Micromechanical simulation of the failure of ber reinforced composites Chad M. Landis a, *, Irene J. Beyerlein

More information

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,

More information

Prediction of Elastic Constants on 3D Four-directional Braided

Prediction of Elastic Constants on 3D Four-directional Braided Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

More information

On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar

On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar NDT&E International 33 (2000) 401 407 www.elsevier.com/locate/ndteint On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar T.-T. Wu*, J.-H. Sun, J.-H.

More information

Buckling of cylindrical sandwich shells with metal foam cores

Buckling of cylindrical sandwich shells with metal foam cores International Journal of Solids and Structures 37 (2000) 6777±6794 www.elsevier.com/locate/ijsolstr Buckling of cylindrical sandwich shells with metal foam cores J.W. Hutchinson a, *, M.Y. He b a Division

More information

Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites

Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites Computational Materials Science 39 (2007) 198 204 www.elsevier.com/locate/commatsci Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites S. Kari

More information

Ball bearing skidding under radial and axial loads

Ball bearing skidding under radial and axial loads Mechanism and Machine Theory 37 2002) 91±113 www.elsevier.com/locate/mechmt Ball bearing skidding under radial and axial loads Neng Tung Liao a, Jen Fin Lin b,* a Department of Mechanical Engineering,

More information

Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip

Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip Giovanni Formica 1, Michela Talò, Walter Lacarbonara 1 Dipartimento di Architettura, Università degli Studi Roma Tre

More information

3D Compression Molding

3D Compression Molding Autodesk Simulation Moldflow Insight 2014 3D Compression Molding Executive summary In this work, the simulation results from a program developed for the three-dimensional analysis of compression molding

More information

Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics

Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Neelam Mishra 1, and Kaushik Das *1 1 School of Minerals Metallurgical and Materials Engineering, Indian

More information

On the uniformity of stresses inside an inhomogeneity of arbitrary shape

On the uniformity of stresses inside an inhomogeneity of arbitrary shape IMA Journal of Applied Mathematics 23) 68, 299 311 On the uniformity of stresses inside an inhomogeneity of arbitrary shape Y. A. ANTIPOV Department of Mathematics, Louisiana State University, Baton ouge,

More information

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading Indian Journal of Engineering & Materials Sciences Vol. 15, October 2008, pp. 382-390 Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading K Sivaji Babu a *, K Mohana

More information

Fatigue Damage Development in a Steel Based MMC

Fatigue Damage Development in a Steel Based MMC Fatigue Damage Development in a Steel Based MMC V. Tvergaard 1,T.O/ rts Pedersen 1 Abstract: The development of fatigue damage in a toolsteel metal matrix discontinuously reinforced with TiC particulates

More information

A fracture parameter for welded structures with residual stresses

A fracture parameter for welded structures with residual stresses Computational Mechanics 22 (1998) 281±288 Ó Springer-Verlag 1998 A fracture parameter for welded structures with residual stresses Y.-C. Hou, J. Pan Abstract In this paper, a nite element procedure to

More information

Enhancing Prediction Accuracy In Sift Theory

Enhancing Prediction Accuracy In Sift Theory 18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department

More information

Soufiane Belhouideg. Keywords Composite materials, homogenization method, Eshelby s inclusion problem.

Soufiane Belhouideg. Keywords Composite materials, homogenization method, Eshelby s inclusion problem. Prediction of effective mechanical properties of composite materials using homogenization approach: Application to tungsten fiber reinforced bulk metallic glass matrix composite Soufiane Belhouideg Abstract

More information

Analysis of intersonic crack growth in unidirectional ber-reinforced composites

Analysis of intersonic crack growth in unidirectional ber-reinforced composites Journal of the Mechanics and Physics of Solids 47 (999) 893±96 Analysis of intersonic crack growth in unidirectional ber-reinforced composites Y. Huang a, *, W. Wang b, C. Liu c, A.J. Rosakis d a Department

More information

Effect of fibre shape on transverse thermal conductivity of unidirectional composites

Effect of fibre shape on transverse thermal conductivity of unidirectional composites Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and

More information

Authors: Correspondence: ABSTRACT:

Authors: Correspondence: ABSTRACT: Multi-Scale Modeling of the Impact and Failure of Fiber Reinforced Polymer Structures using DIGIMAT to LS-DYNA Interface Authors: L. Adam, R. Assaker & R. Ramaya e-xstream engineering S.A. 7, Rue du Bosquet.

More information

Mode selection of guided waves for ultrasonic inspection of gas pipelines with thick coating

Mode selection of guided waves for ultrasonic inspection of gas pipelines with thick coating Mechanics of Materials 31 (1999) 165±174 Mode selection of guided waves for ultrasonic inspection of gas pipelines with thick coating E. Pan a, J. Rogers b, S.K. Datta a, *, A.H. Shah b a Department of

More information

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE 1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for

More information

strain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip

strain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this

More information

STEEL FIBER CURVATURE IN CONCRETE COMPOSITES: MODULUS PREDICTIONS USING EFFECTIVE STEEL FIBER PROPERTIES

STEEL FIBER CURVATURE IN CONCRETE COMPOSITES: MODULUS PREDICTIONS USING EFFECTIVE STEEL FIBER PROPERTIES American Journal of Applied Sciences 11 (1): 145-151, 2014 ISSN: 1546-9239 2014 Science Publication doi:10.3844/ajassp.2014.145.151 Published Online 11 (1) 2014 (http://www.thescipub.com/ajas.toc) STEEL

More information

ME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites

ME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-

More information

Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact

Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Philip P. Garland 1 and Robert J. Rogers 2 1 School of Biomedical Engineering, Dalhousie University, Canada 2 Department

More information

Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation

Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation Charles L. Tucker III Department of Mechanical and Industrial Engineering University of Illinois at Urbana-Champaign

More information

Effective elastic moduli of two dimensional solids with distributed cohesive microcracks

Effective elastic moduli of two dimensional solids with distributed cohesive microcracks European Journal of Mechanics A/Solids 23 (2004) 925 933 Effective elastic moduli of two dimensional solids with distributed cohesive microcracks Shaofan Li a,, Gang Wang a, Elise Morgan b a Department

More information

Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear

Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear Xi an 2-25 th August 217 Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension compression or shear Heyin Qi 1 Mingming Chen 2 Yonghong Duan 3 Daxu

More information

Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip

Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 1 1 JULY 2003 Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip Dejun Ma Department

More information

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and

More information

Dynamic response of tubular joints with an annular void subjected to a harmonic torsional load

Dynamic response of tubular joints with an annular void subjected to a harmonic torsional load SPECIAL ISSUE PAPER 361 Dynamic response of tubular joints with an annular void subjected to a harmonic torsional load A Vaziri and H Nayeb-Hashemi* Department of Mechanical, Industrial and Manufacturing

More information

On the Path-Dependence of the J-Integral Near a Stationary Crack in an Elastic-Plastic Material

On the Path-Dependence of the J-Integral Near a Stationary Crack in an Elastic-Plastic Material Dorinamaria Carka Chad M. Landis e-mail: landis@mail.utexas.edu Department of Aerospace Engineering and Engineering Mechanics, University of Texas at Austin, 10 East 4th Street, C0600 Austin, TX 7871-035

More information

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A

More information

Apparent thermal conductivity of periodic two-dimensional composites

Apparent thermal conductivity of periodic two-dimensional composites Computational Materials Science 25 (2002) 329 338 www.elsevier.com/locate/commatsci Apparent thermal conductivity of periodic two-dimensional composites M. Jiang a, I. Jasiuk b, *, M. Ostoja-Starzewski

More information

International Journal of Engineering Science

International Journal of Engineering Science International Journal of Engineering Science 63 (2013) 1 9 Contents lists available at SciVerse ScienceDirect International Journal of Engineering Science journal homepage: www.elsevier.com/locate/ijengsci

More information

Received 5 April 2003; received in revised form 29 May 2003

Received 5 April 2003; received in revised form 29 May 2003 Materials Science and Engineering A360 (2003) 339/344 www.elsevier.com/locate/msea Connection between elastic moduli and thermal conductivities of anisotropic short fiber reinforced thermoplastics: theory

More information

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Biplab Chatterjee, Prasanta Sahoo 1 Department of Mechanical Engineering, Jadavpur University

More information

Practice Final Examination. Please initial the statement below to show that you have read it

Practice Final Examination. Please initial the statement below to show that you have read it EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use

More information

Thermal load-induced notch stress intensity factors derived from averaged strain energy density

Thermal load-induced notch stress intensity factors derived from averaged strain energy density Available online at www.sciencedirect.com Draft ScienceDirect Draft Draft Structural Integrity Procedia 00 (2016) 000 000 www.elsevier.com/locate/procedia 21st European Conference on Fracture, ECF21, 20-24

More information

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds

More information

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011 Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.

More information

Effect of Interlayers on Mechanical Properties and Interfacial Stress Transfer of 2D Layered Graphene- Polymer Nanocompsites

Effect of Interlayers on Mechanical Properties and Interfacial Stress Transfer of 2D Layered Graphene- Polymer Nanocompsites University of Kentucky UKnowledge Theses and Dissertations--Mechanical Engineering Mechanical Engineering 2017 Effect of Interlayers on Mechanical Properties and Interfacial Stress Transfer of 2D Layered

More information

CRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th

CRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA

More information

Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model

Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model J.B. Sang*, L.F. Sun, S.F. Xing,

More information

Computational Analysis for Composites

Computational Analysis for Composites Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics

More information

Finite element analysis of longitudinal debonding between fibre and matrix interface

Finite element analysis of longitudinal debonding between fibre and matrix interface Indian Journal of Engineering & Materials Sciences Vol. 11, February 2004, pp. 43-48 Finite element analysis of longitudinal debonding between fibre and matrix interface K Aslantaş & S Taşgetiren Department

More information

Finite solid circular cylinders subjected to arbitrary surface load. Part II Ð Application to double-punch test

Finite solid circular cylinders subjected to arbitrary surface load. Part II Ð Application to double-punch test International Journal of Solids and Structures 37 (2000) 5733±5744 www.elsevier.com/locate/ijsolstr Finite solid circular cylinders subjected to arbitrary surface load. Part II Ð Application to double-punch

More information

Material instability criterion near a notch-tip under locally adiabatic deformations of thermoviscoplastic materials

Material instability criterion near a notch-tip under locally adiabatic deformations of thermoviscoplastic materials Theoretical and Applied Fracture Mechanics 30 (1998) 153±158 Material instability criterion near a notch-tip under locally adiabatic deformations of thermoviscoplastic materials L. Chen, R.C. Batra * Department

More information

Analysis of microstructural elds in heterogeneous piezoelectric solids

Analysis of microstructural elds in heterogeneous piezoelectric solids PERGAMON International Journal of Engineering Science 37 (1999) 665±685 Analysis of microstructural elds in heterogeneous piezoelectric solids J.Y. Li, M.L. Dunn * Center for Acoustics, Mechanics, and

More information

SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION

SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION 8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas

More information

A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core

A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core Commun. Theor. Phys. 56 774 778 Vol. 56, No. 4, October 5, A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core JIANG

More information

Composite Structures

Composite Structures Composite Structures 96 (2013) 601 605 Contents lists available at SciVerse ScienceDirect Composite Structures journal homepage: www.elsevier.com/locate/compstruct Meso-scale modelling of shock wave propagation

More information

Most of the material in this package is based on a recently published book. This is:

Most of the material in this package is based on a recently published book. This is: Mechanics of Composite Materials Version 2.1 Bill Clyne, University of Cambridge Boban Tanovic, MATTER Assumed Pre-knowledge It is assumed that the student is familiar with simple concepts of mechanical

More information

Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites

Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites International OPEN ACCESS Journal Of Modern Engineering Research (IJMER) Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites Sonali Gholap 1, Dr. Dhananjay R. Panchagade

More information

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Sandhyarani Biswas Department of Mechanical Engineering, N.I.T Rourkela, INDIA Abstract: Now-a-days,

More information

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES L.V. Smith 1 *, M. Salavatian 1 1 School of Mechanical and Materials

More information

Viscosity of magmas containing highly deformable bubbles

Viscosity of magmas containing highly deformable bubbles Journal of Volcanology and Geothermal Research 105 (2001) 19±24 www.elsevier.nl/locate/jvolgeores Viscosity of magmas containing highly deformable bubbles M. Manga a, *, M. Loewenberg b a Department of

More information

Surface stress and relaxation in metals

Surface stress and relaxation in metals J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S0953-8984(00)11386-4 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown

More information

On the Ballooning Motion of Hyperelastic Strings

On the Ballooning Motion of Hyperelastic Strings On the Ballooning Motion of Hyperelastic Strings S. Sarangi 1, R. Bhattacharyya, A. K. Samantaray Indian Institute of Technology, 7130 Kharagpur, India. Abstract The e ect of material nonlinearity on the

More information

XI. NANOMECHANICS OF GRAPHENE

XI. NANOMECHANICS OF GRAPHENE XI. NANOMECHANICS OF GRAPHENE Carbon is an element of extraordinary properties. The carbon-carbon bond possesses large magnitude cohesive strength through its covalent bonds. Elemental carbon appears in

More information

CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES

CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES 4. Introduction Fillers added to polymer matrices as additives are generally intended for decreasing the cost (by increase in bulk) of

More information

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS Djoković, J. M., et.al.: Influence of Temperature on Behavior of the Interfacial THERMAL SCIENCE: Year 010, Vol. 14, Suppl., pp. S59-S68 S59 INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK

More information

Contemporary Research in Engineering Science

Contemporary Research in Engineering Science Bazant, Z.P., and Prat, P.C. (1995). 'Elastic material with systems of growing or closing cracks: Tangential Stiffness.' in Contemporary Research in Engineering Science (Proc., Eringen Medal Symp. honoring

More information

Chapter 3. Load and Stress Analysis

Chapter 3. Load and Stress Analysis Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3

More information

Module 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites

Module 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites Module 4: Behaviour of a Laminae-II Learning Unit 1: M1 M4.1 Mechanics of Composites M4.1.1 Introduction to Mechanics of Composites The relation between ply uniaxial strengths and constituent properties

More information

Finite Element Analysis of FRP Debonding Failure at the Tip of Flexural/Shear Crack in Concrete Beam

Finite Element Analysis of FRP Debonding Failure at the Tip of Flexural/Shear Crack in Concrete Beam Marquette University e-publications@marquette Civil and Environmental Engineering Faculty Research and Publications Civil and Environmental Engineering, Department of 12-1-2013 Finite Element Analysis

More information

Open-hole compressive strength prediction of CFRP composite laminates

Open-hole compressive strength prediction of CFRP composite laminates Open-hole compressive strength prediction of CFRP composite laminates O. İnal 1, A. Ataş 2,* 1 Department of Mechanical Engineering, Balikesir University, Balikesir, 10145, Turkey, inal@balikesir.edu.tr

More information

A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests

A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests Bolton, M. D. & Whittle, R. W. (999). GeÂotechnique 49, No., 33±4 TECHNICAL NOTE A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests M. D. BOLTON and R. W. WHITTLE{

More information

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING N. Krishna Vihari 1, P. Phani Prasanthi 1, V. Bala Krishna Murthy 2* and A. Srihari Prasad 3 1 Mech. Engg. Dept., P. V. P. Siddhartha

More information

An empirical correlation for natural convection heat transfer from thin isothermal circular disks at arbitrary angles of inclination

An empirical correlation for natural convection heat transfer from thin isothermal circular disks at arbitrary angles of inclination International Journal of Heat and Mass Transfer 45 2002) 1159±1163 Technical Note An empirical correlation for natural convection heat transfer from thin isothermal circular disks at arbitrary angles of

More information

The Effects of Transverse Shear on the Delamination of Edge-Notch Flexure and 3-Point Bend Geometries

The Effects of Transverse Shear on the Delamination of Edge-Notch Flexure and 3-Point Bend Geometries The Effects of Transverse Shear on the Delamination of Edge-Notch Flexure and 3-Point Bend Geometries M. D. Thouless Department of Mechanical Engineering Department of Materials Science & Engineering University

More information

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES

EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC

More information

CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles

CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles Mark Hughes 14 th March 2017 Today s learning outcomes To understand the role of reinforcement, matrix and

More information

Analytical formulation of Modified Upper Bound theorem

Analytical formulation of Modified Upper Bound theorem CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are

More information

ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS

ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS Dieter H. Pahr and Helmut J. Böhm Institute of Lightweight Design and Structural

More information

Finite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements

Finite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements Proceedings o the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. Finite Element Modeling o Residual Thermal Stresses in Fiber-Reinorced Composites Using Dierent Representative

More information

CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement

CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement Mark Hughes 19 th September 2016 Outline Fibre architecture Volume fraction and the rule of mixtures Principle

More information

Debonding process in composites using BEM

Debonding process in composites using BEM Boundary Elements XXVII 331 Debonding process in composites using BEM P. Prochazka & M. Valek Czech Technical University, Prague, Czech Republic Abstract The paper deals with the debonding fiber-matrix

More information

Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load

Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Sriram Chintapalli 1, S.Srilakshmi 1 1 Dept. of Mech. Engg., P. V. P. Siddhartha Institute of Technology.

More information

THERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH CYLINDERS OF FUNCTIONALLY GRADED MATERIALS

THERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH CYLINDERS OF FUNCTIONALLY GRADED MATERIALS Journal of Thermal Stresses, 28: 391 408, 2005 Copyright # Taylor & Francis Inc. ISSN: 0149-5739 print/1521-074x online DOI: 10.1080/01495730590916623 THERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH

More information

Magnetoelectroelastic multi-inclusion and inhomogeneity problems and their applications in composite materials

Magnetoelectroelastic multi-inclusion and inhomogeneity problems and their applications in composite materials International Journal of Engineering Science 38 (2000) 1993±2011 www.elsevier.com/locate/ijengsci Magnetoelectroelastic multi-inclusion and inhomogeneity problems and their applications in composite materials

More information

J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia

J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia Evaluation of fracture parameters for functionally gradient materials J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, 842 20 Bratislava, Slovakia

More information

Finite element modeling of the resin transfer molding process based on homogenization techniques

Finite element modeling of the resin transfer molding process based on homogenization techniques Computers and Structures 76 (2000) 379±397 www.elsevier.com/locate/compstruc Finite element modeling of the resin transfer molding process based on homogenization techniques Z. Dimitrovova a, *, L. Faria

More information

Shock wave speed and stress-strain relation of aluminium honeycombs under dynamic compression

Shock wave speed and stress-strain relation of aluminium honeycombs under dynamic compression EPJ Web of Conferences 83, 7 (8) DYMAT 8 https://doi.org/.5/epjconf/8837 Shock wave speed and stress-strain relation of aluminium honeycombs under dynamic compression Peng Wang,*, Jun Zhang, Haiying Huang,

More information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

Failure modes of glass panels subjected to soft missile impact

Failure modes of glass panels subjected to soft missile impact Failure modes of glass panels subjected to soft missile impact L. R. Dharani & J. Yu Dept. of Mech. and Aerospace Engineering and Engineering Mechanics, University of Missouri-Rolla, U.S.A. Abstract Damage

More information

L e c t u r e. D r. S a s s a n M o h a s s e b

L e c t u r e. D r. S a s s a n M o h a s s e b The Scaled smteam@gmx.ch Boundary Finite www.erdbebenschutz.ch Element Method Lecture A L e c t u r e A1 D r. S a s s a n M o h a s s e b V i s i t i n g P r o f e s s o r M. I. T. C a m b r i d g e December

More information