Module-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
|
|
- Bethanie Ward
- 5 years ago
- Views:
Transcription
1 Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination theory it is clear that if we are given the strains and curvatures at a point (x, y) on the reference surface ε o, ε o, γ o, K o, K o and x y xy x y K o then we can compute the strain distributions through the thickness of the laminate. By xy using the stress-strain relations, we can compute the stress distributions. By using the transformation relations, we can determine the stresses and strains in the principal material system, and by using the definitions of the stress resultants, we can compute the force and moment resultants acting at that point. If we specify that the given strains and curvatures are the same at every point on the reference surface of a laminate of a given length and width, then we can determine the total forces and moments acting on the edges of the laminate. These steps are all a result of the plane stress assumption and the Kirchhoff s hypothesis. Figure M6.1.1 illustrates the connection between these steps. What remains in the development of classical lamination theory is to be able to specify the force and moment resultants acting at a point (x, y) on the reference surface, and then to be able to compute the strains and curvatures of that point on the reference surface that these resultants cause. We want to fill in the missing link in the upper left portion of Figure M6.1.1, namely, compute the reference surface strains and curvatures, knowing the stress resultants. With these computed reference surface strains and curvatures, we can then, as in the previous examples, compute the strain and stress distributions through the thickness of the laminate. Relating the stress resultants to the reference surface strains and curvatures is an important step. In the application of composite materials to structures, we are often given the forces and moments that act on a laminate, and we want to know the stresses and strains that are caused by these loads. Having a relation between the force and moment resultants and the strains and curvatures at a point (x, y) on the reference surface will allow us to do this. In what follows we shall develop this key relationship.
2 Figure M6.1.1 Interrelations among concepts developed so far M6.1.1 Formal Definition of Force and Moment Resultants Though we have introduced and computed them in the examples, we have not formally defined in one location the three force resultants and three moment resultants that are a result of integrating the stresses through the thickness of the laminate. The stress resultants in the x direction, N in the y direction, N and in shear, N are defined as, x y xy
3 (M6.1.1) The resultants N and x N referred to as the shear force resultant. y will be referred to as the normal force resultants, and N will be xy Figure M6.1.2 illustrates a small element of laminate surrounding a point (x, y) on the geometric mid-plane, and the figure indicates the directions of these three stress resultants. The units of the stress resultants are force per unit length; the unit of length is a unit of length in the x direction or in the y direction. If the small element of laminate has dimensions dx by dy, then the force in the x direction due to N, the force in the x direction due to N, the force in the y direction due x to N y, and the force in the y direction due to xy N are, respectively, xy (M6.1.2) The moment resultants, M, M and M are defined as x y xy (M6.1.3)
4 The resultants M and M will be referred to as bending moment resultants and xy will be x y referred to as the twisting moment resultant. In Figure M6.1.3 the sense of these three moment resultants is illustrated, and the senses shown are important. Figure M6.1.2 Definition of force resultants N, N and x y N xy Figure: M6.1.3 Definition of moment resultants M, M, and M x y xy Shortly we shall be relating a positive M, for example, with a positive curvature in the x x direction, K o x. Knowing the correct sense of positive quantities is critical. The sense illustrated in Figure M6.1.3 is consistent with the cross-product definition of moments learned in elementary mechanics courses (i.e.( M = r F ). The units of the moment resultants are moment per unit length, again the unit of length being a unit of length in either the x direction or y direction. If a small element of laminate has dimensions dx by dy, then the moments about the positive y and the positive x axis due to M, M and M are, respectively, x y xy (M6.1.4)
5 We will use these formal definitions of N,..., M in computing the laminate stiffness x matrix. Specifically, we will use the force and moment resultants when studying the equilibrium of a small element of laminated plate. As the definitions of the moment resultants involve integrals with respect to the coordinate z, and because z is measured relative to the reference surface, the moment resultants have to be thought of as producing moments on the reference surface. The force resultants can be considered to act anywhere through the thickness of the laminate. But when force resultants are used in conjunction with moment resultants, it is most consistent to think of the force resultants as producing forces on the reference surface. This will also be consistent when we begin discussing the relationship between the resultants and reference surface response, namely, ε o,..., K o. From this point forward, then, the force x xy and moment resultants will be considered to act on the reference surface. M6.1.2 Laminate Stiffness: The ABD Matrix Equation (M6.1.1) and Equation (M6.1.3) give the definitions of the force and moment resultants. These definitions involve the stresses σ, σ, andτ x y xy. From the stress-strain relations, the stresses can be written in terms of the strains by the now-familiar relation, xy (M6.1.5) and the strains can, in turn, be written in terms of the strains and curvatures of the reference surface to produce another familiar relation, namely, (M6.1.6) We have emphasized repeatedly that in a laminate the stresses are functions of z because the transformed reduced stiffness Q ; change from layer to layer, and also because of the linear dependence on z - The reference surface strains and curvatures, by definition, do not depend on z - let us examine the expression for the normal force resultant N by substituting for σ from x x Equation (M6.1.6) into the integrand of the first equation of Equation (M6.1.1), specifically: (M6.1.7)
6 and expanding the integrand leads to, (M6.1.8) The integration can be distributed over the six terms to give, (M6.1.9) (M6.1.9) Because the reference surface deformations are not functions of z, they can be taken outside of the integrals, resulting in, (M6.1.10) The integrals through the thickness now involve just the material properties, specifically, the reduced stiffnesses. We need to look at these integrals in detail. The first integral in Equation (M6.1.10) is (M6.1.11)
7 Because the reduced stiffnesses are materials properties, they are constant within a layer. However, from layer to layer the values of the reduced stiffnesses change. Because of this piecewise-constant property of the reduced stiffnesses, the integral through the thickness can be expanded to give, (M6.1.12) In the above expansion the locations of the layer interfaces have been used and the subscript on Q 11 indicates layer number. The systematic scheme introduced earlier to identify the layer numbers and the z location of the layer interfaces will be used to advantage in what follows. Because within each layer Q 11 is a constant; in the above series of integrals it can be removed from each integral and the series of integrals written as (M6.1.13) But the integrations with respect to z are quite simple, Equation (M6.1.13) becoming, (M6.1.14) or (M6.1.15)
8 th The difference in the z values is recognized as simply the thickness of the layer, h so k Equation (M6.1.15) becomes, k (M6.1.16) The integral of Q 11 through the thickness is denoted by A 11 ; that is: (M6.1.17) In a similar manner, we find, referring to the third and fifth terms in Equation (M6.1.10), that (M6.1.18) and (M6.1.19) Thus the first, third, and fifth terms on the right side of Equation (M6.1.10) become, (M6.1.20) The second, fourth, and sixth terms on the right side of Equation (M6.1.10) involve integrals of Q11 z, Q12 zand Q16 z rather than justq11, Q12 and Q 16. Specifically, the second term is, (M6.1.21)
9 Again using the fact that Q 11 is constant within each layer but changes from layer to layer, we find that, (M6.1.22) Here the integrals on z are also quite simple. Equation (M6.1.22) becomes, (M6.1.23) or (M6.1.24) This integral is denoted as B 11 ; that is: (M6.1.25) In like manner, referring to the fourth and sixth terms in Equation (M6.1.10), we find, (M6.1.26) and (M6.1.27)
10 The expression for N in Equation (M6.1.10) becomes, after slight rearrangement, x (M6.1.28) Because they are summations of material and geometric properties, A, A, A, B, B,andB are easy to compute and are inherent properties of the laminate Because the A an db ; are integrals, they represent smeared or integrated properties of the laminate. Equation (M6.1.28) represents a relation between the normal force resultant N and x the six reference surface deformations. This is part of the missing link in Figure M We will develop five other equations that will lead to a total of six relations between the six stress resultants and the six reference surface deformations. The missing link will then be complete. For the four-layer [0/90] 5 laminate, as each layer is the same thickness h, by Equation (M6.1.17), (M6.1.29) or (M6.1.30) By knowing the values ofq, we can get value of, and, A11 A12 A16 A11 = some numerical value1( = C 1) (M6.1.31) Likewise, A11 = some numerical value1( = C 2) A16 = 0 (For cross ply la min ate) (M6.1.32) For all cross-ply laminates A 16 will be zero because for each layerq 16 = 0. The B can be computed from Equations (M6.1.25)-(M6.1.27) as,
11 (M6.1.33) or (M6.1.34) We find, using the numerical values of the interface locations for the [ 0/90 S laminate, Equation (M6.1.21) leads to, ] (M6.1.35) The quantity is exactly zero because the contributions to B from the layers at negative z- B11 11 B 11 locations exactly cancel the contributions to from the layers at the positive z locations. This is a characteristic of a symmetric laminate. The values of all the B will be zero for all symmetric laminates because of this cancellation effect. For symmetric laminates there is no need to go through the algebra to compute the values of the B. They are all exactly zero; thus, (M6.1.36) B16 is zero not only because it is a symmetric laminate, but also because theq 16, are zero for every layer. For the six-layer [ ± 30/ 0] S laminate, (M6.1.37) Using layer orientation, we find that, (M6.1.38) For j = 1, substituting numerical values for the appropriateq, we find
12 and for j = 2, and for j = 6, A11 = some numerical value1( = C 3) (M6.1.39) A12 = some numerical value1( = C 4) (M6.1.40) A 0 16 = (M6.1.41) The term A is exactly zero because the contributions from the +30 layers exactly cancel the 16 contributions from the 30 layers. This is characteristic of the 16-component of the A, for what we will come to call a balanced laminate. Balanced laminates are characterized by ± θ pairs of layers. Because the laminate is symmetric, the components, B, B, and B are exactly zero That is, B = 0 B B = 0 16 = 0 for Symmetric Laminates. ±30/ 0] T exactly one-half the number of layers in any particular direction, the Ay- for the [ ± ] laminate are one-half the values of the A for the [ 30/ 0] S the[ ], is not a symmetric laminate, the When turning to the [ laminate, because relative to the six-layer laminate there are 30/ 0 T ± laminate. However, as ± 30/ 0 T B will not be zero. Specifically, referring to equations (M6.1.39), (M6.1.40) and (M6.1.41), we find that, For the B, A 0 N/ 16 = m (M6.1.42)
13 (M6.1.43) For j = 1 and using the interface locations, we find, B11 = some numerical value1( = C 5) (M6.1.44) For j = 2and6, B12 = some numerical value1( = C 6) B16 = some numerical value1( = C 7) (M6.1.45) To continue with the development of the laminate stiffness matrix, we can substitute the expressions for σ y and τ xy into the definitions in Equation (M6.1.1) for N and xy, y respectively, to yield, (M6.1.46) and (M6.1.47) As in going from Equation (M6.1.8) to (6.1.9), Equation (M6.1.9) to (6.1.10), and finally to Equation (M6.1.28), and continuing to define quantities A and B rearranging, we find that (M6.1.48) and (M6.1.49)
14 If we combine these definitions with the definitions associated with N, we can write the results x in matrix notation as, (M6.1.50) The above matrix equation constitutes three of the six relations between the six stress resultants and the six reference surface deformations. Three equations have yet to be derived, but with the above notation, it is clear that the A, and B are components of three-by-three matrices. To derive the remaining three equations, we use the moment expressions, Equation (M6.1.3), by substituting for ax from Equation (M6.1.6) into the first equation of Equation (M6.1.3) to yield, (M6.1.51) where, as before, the integration can be distributed over the six terms to give, (M6.1.52) Removing the reference surface deformations from the integrals, we find that, (M6.1.53) The first integral on the right hand side of Equation (M6.1.56) is,
15 (M6.1.54) and this is immediately recognizable, from Equation (M6.1.25), as and fifth terms on the right side of Equation (M6.1.53) become, B 11. In fact, the first, third, (M6.1.55) The appearance of B important to remember. in both the force resultant and the moment resultant expressions; this is The second, fourth, and sixth terms on the right hand side of Equation (M6.1.53) are new; the second term is, (M6.1.56) which can be expanded to, (M6.1.57) Because of the constancy of the reduced stiffnesses within layers, (M6.1.58) or
16 (M6.1.59) This integral is denoted as D 11 ; that is: (M6.1.60) Referring to the fourth and sixth terms in Equation (M6.1.53), we find that, (M6.1.61) And (M6.1.62) With the definitions of the D, and the earlier B, the expression for the bending moment resultant M in Equation (M6.1.53) becomes, x (M6.1.63) The D, are similar to the interfaces. B in that they involve powers of the z locations of the layer Although the A also involved z location of the layer interfaces, it was really only the thickness of each layer that was involved. The actual location of the layers through the thickness is not important. However, with the locations B, and through the thickness of the laminate, as well as the thickness of the layers, are important. To continue with the numerical examples, the values of D, D and D, for the example laminates will be computed for illustration. For the four-layer 0/90 laminate, D [ ] S
17 (M6.1.64) or (M6.1.65) The ensuing algebra leads to, For D 12, D11 = some numerical value1( = C 8) (M6.1.66) (M6.1.67) or, (M6.1.68) Numerically, D12 = some numerical value1( = C 9) (M6.1.69) Because Q 16 is zero for each layer, D 0 16 = (M6.1.70) For the [ ± 30/ 0] S laminate,
18 (M6.1.71) or (M6.1.72) We find, referring to the interface locations for this laminate that the algebra leads to, D11 = some numerical value1( = C 8) D12 = some numericalvalue1( = C 9) D16 = some numericalvalue1( = C 10) (M6.1.73) Note that the value of D is not zero. Unlike the situation with A, where the contribution of the layers to the A 16, cancels the contribution of the +30 layers, the contributions of the layers with opposite orientation do not cancel in the calculation of D 16. This is because the locations of the layers, in addition to the layer thicknesses, are involved in the computation of the. We will see later in this module that the existence of a nonzero, term has important consequences. Also, we will see the dependence of the A on the ± θ layer arrangement is 16 identical to the dependence of A, and the dependence of D is identical to the dependence of. D 16 Finally, for the three-layer [ ±30/ 0] T laminate, D 16 D (M6.1.74) or,
19 (M6.1.75) Using numerical values leads to, D11 = some numerical value1( = C 8) N m D12 = some numericalvalue1( = C 9) N m D16 = some numericalvalue1( = C 10) N m (M6.1.76) To finish the development of the laminate stiffness matrix, we can substitute the expressions for σ and τ into the definitions in Equation (M6.1.3) for M and M. The result is, y xy y xy (M6.1.77) and, (M6.1.78) As in going from Equation (M6.1.51) to (6.1.52), Equation (M6.1.52) to (6.1.53), and finally to Equation (M6.1.63), and continuing to define additional, we find that, D (M6.1.79) and, (M6.1.80) We can write the results for the three moments M, M andm in matrix notation as, y y xy
20 (M6.1.81) The above analysis constitutes the remaining three relations of the six relations between the six stress resultants and the six reference surface deformations. We can combine the six relations just derived to form one matrix relation between the six stress resultants and the six reference surface deformations. The combined matrix relation is, (M6.1.82) Of course, in the above, (M6.1.83) The six-by-six matrix in Equation (M6.1.82) consisting of the components D, i= 1,2,6; j = 1,2,6, is called the laminate stiffness matrix. For obvious reasons it is also known as the ABD matrix. The ABD matrix defines a relationship between the stress resultants (i.e., loads) applied to a laminate, and the reference surface strains and curvatures (i.e., deformations). This form is a direct result of the Kirchhoff hypothesis, the planestress assumption, and the definition of the stress resultants. The laminate stiffness matrix involves everything that is used to define the laminate layer material properties, fiber orientation, thickness, and location.
21 The relation of Equation (M6.1.82) can be inverted to give, (M6.1.84a) where, (M6.1.84b) Equation (M6.1.82) and its inverse Equation (M6.1.84) are very important in the analysis of composite structures, particularly the inverse form. Knowing the loads, we can determine the reference surface strains and curvatures uniquely from Equation (M6.1.84). If we know the reference surface strains and curvatures, as we have seen in the previous examples, then we can determine the strains and stresses in each layer. Figure M6.1.4 illustrates the importance of the ABD matrix, namely the fact that it closes the loop on the analysis at composite laminates (see Figure M6.1.1). It is appropriate at this point to present the numerical values of the remaining components of the ABD matrix for each of the three laminates we have been discussing. So far only limited numerical values have been illustrated. It will be left as an exercise to verify the numerical values presented for all the components of the ABD matrix for the three laminates. It is important to note the units on the various components of the ABD matrix. The units are consistent with the units of the reference surface strains and curvatures and the units of the stress resultants.
22 Figure M6.1.4 Closing the loop on concepts developed Writing the numerical values in full form is quite revealing; the coupling, or lack thereof, between the various stress resultants and reference deformations becomes quite obvious. The coupling is important in determining how much interaction there will be between the specific loads and specific reference surface strains and curvatures. Essentially, the more zeros there are in the ABD matrix, or its inverse, the less coupling there is between the stress resultants and reference surface deformations. For example, for the [0/90] S and [0/30] S laminates, because all components of the B matrix are zero, the reference surface curvatures are decoupled from the force resultants, and the reference surface strains are decoupled from the moment resultants. The six-by-six relation of Equation (M6.1.82) reduces to two three-by-three relations, namely, (M6.1.85) And
23 (M6.1.86) Such a reduction is not possible when there are elements of the B matrix present (i.e., when the laminate is unsymmetric). We will have more to say regarding the couplings within the ABD matrix in the next section, where the degree of coupling, and its effects on laminate response, are categorized in a systematic way. Moreover, we can do the categorization a priori, before any analysis is started, and we can use the advantages of knowing the level of coupling to simplify the analysis procedure. Learning Unit-2: M6.2 M6.2 Special Classification of Laminates We have seen in the numerical results presented in the previous section that the form of the ABD matrix depends strongly on whether the laminate is what we referred to as symmetric, whether it consists of just 0 and 90 layers, or whether for every layer at orientation + θ there is a layer at θ. In this section we would like to formally classify laminates as to their stacking arrangement, and indicate the effects of the various classifications on the ABD matrix. In a subsequent section, physical interpretations will be given to particular components of the ABD matrix. As a sidenote, in connection with stress-strain equations and their inverse relations, it was stated that the presence of the zeros in the compliance and stiffness matrices in the principal material coordinate system was associated with a fiber-reinforced material being orthotropic. Similarly, with the lack of any zero terms in Equation (M6.1.82), the elastic behavior of a general laminate is sometimes termed nonorthotropic, or anisotropic. M6.2.1 Symmetric Laminates A laminate is said to be symmetric if for every layer to one side of the laminate reference surface with a specific thickness, specific material properties, and specific fiber orientation, there is another layer the identical distance on the opposite side of the reference surface with the identical thickness, material properties, and fiber orientation. If a laminate is not symmetric, then it is referred to as an unsymmetric laminate. Note that this pairing on opposite sides of the reference surface must occur for every layer. Although we have restricted our examples to laminates with layers of the same material and same thickness, this does not have to be the case for a laminate to be symmetric. A four-layer laminate with the two outer layers made of glassreinforced material oriented at +45 and the two inner layers made of graphite-reinforced material oriented at 30 is a symmetric laminate. For a symmetric laminate all the components of the B matrix are identically zero. Consequently, the full six-by-six set of equations in Equation (M6.1.82) decouples into two three-by-three sets of equations, namely,
24 (M6.2.1a) And (M6.2.1b) With symmetric laminates, the inverse relation, Equation (M6.1.84), also degenerates into two three-by-three relations, specifically, (M6.2.2a) (M6.2.2b) In the above the a matrix is the inverse of the A matrix and has components,
25 (M6.2.3a) and the d matrix is the inverse of the D matrix and has components, (M6.2.3b) where the definitions,
26 (M6.2.3c) and, (M6.2.3d) have been used. Because neither the A nor D matrices have any zero terms, a general symmetric laminate is anisotropic in both inplane and bending behavior. M6.2.2 Balanced Laminates A laminate is said to be balanced if for every layer with a specified thickness, specific material properties, and specific fiber orientation, there is another layer with the identical thickness, material properties, but opposite fiber orientation somewhere in the laminate. The layer with opposite fiber orientation does not have to be on the opposite side of the reference surface, nor immediately adjacent to the other layer, nor anywhere in particular. The other layer can be anywhere within the thickness. A laminate does not have to be symmetric to be balanced. The symmetric [ ± 30/ 0] S laminate and the unsymmetric [ ± 30/ 0] T laminate are both balanced laminates. If a laminate is balanced, then the stiffness matrix components A, and A are always zero because the Q 16 and Q 26 from the layer pairs with opposite orientation are of opposite sign, and the net contribution to A 16, and A from these layer pairs is then zero. To 26 classify as balanced, all off-axis layers must occur in pairs. As a special case, all laminates consisting entirely of layers with their fibers oriented at either 0 or 90, to be discussed shortly, are balanced, as Q 16 and Q26 are zero for every layer, resulting in A, and A being zero. The ABD matrix of a balanced but otherwise general laminate is not that much simpler than the ABD matrix of Classification of Laminates and Their Effect on a general unsymmetric, unbalanced laminate. The full six-by-six form, Equation (M6.1.85), applies but with the A 16, and A 26 components set to zero. To obtain the inverse relation, Equation (M6.1.87a), the full six-by-six with the zero entries for A, and A must be inverted M6.2.3 Symmetric Balanced Laminates A laminate is said to be a symmetric balanced laminate if it meets both the criterion for being symmetric and the criterion for being balanced. If this is the case, then Equation (M6.1.82) takes the decoupled form,
27 (M6.2.4a) with, and, (M6.2.4b) The inverse form involving the force resultants is, (M6.2.5a) with, The inverse form involving the moments is no different than for symmetric but unbalanced laminates, namely, (M6.2.5b) For the symmetric balanced case, the third of Equation (M6.1.82) decouples from all other equations. The components of the a matrix in Equation (M6.2.5a) are,
28 (M6.2.6a) and the components of the d matrix are the same as the symmetric but balanced case, Equation (M6.2.3b), namely, (M6.2.6b) with det[d] being given by Equation (M6.2.3d). Because A, and A are zero, a symmetric balanced laminate is orthotropic with respect to inplane behavior. M6.2.4 Cross-Ply Laminates A laminate is said to be a cross-ply laminate if every layer has its fibers oriented at either 0 or 90. If this is the case, because Q 16 and Q26 are zero for every layer, then A, A, B, B, D andd are zero, and there is some decoupling of the six equations. In
29 particular, the six equations decouple to a set of four equations and a set of two equations, namely, (M6.2.7a) (M6.2.7b) M6.2.5 Symmetric Cross-Ply Laminates A laminate is said to be a symmetric cross-ply laminate if it meets the criterion for being symmetric and the criterion for being cross-ply. This results in the simplest form of the ABD matrix. All the B, are zero, and A, A, D andd are zero The six equations of Equation (M6.1.82) decouple significantly. Specifically, (M6.2.8a) with, and, (M6.2.8b) with
30 The inverted form is, (M6.2.9a) and with (M6.2.9b) with, The components of the a matrix are the same as the symmetric balanced case, (M6.2.10a) and the components of the d matrix are,
31 (M6.2.10b) Compared to the most general case, the symmetric cross-ply laminate is trivial. Because A, A, D and D are zero, a symmetric cross-ply laminate is orthotropic with respect to both inplane and bending behavior. M6.2.6 Single Isotropic Layer As a comparison, it is of interest to record the ABD matrix for the case of a single isotropic layer, the subject of studies in classical plate theory. What we have presented in the development of classical lamination theory is a generalization of classical plate theory. The isotropic case is a special case of this generalization, and for a single isotropic layer of thickness H, (M6.2.11) In the above equation, E is the extensional modulus, referred to as Young's modulus when dealing with isotropic materials, and ν is Poisson's ratio of the single layer of material. Because of Equation (M6.2.11), the ABD matrix for a single isotropic layer decouples and greatly simplifies, and Equation (M6.1.82) becomes, (M6.2.12a)
32 with, and, (M6.2.12b) with, The inverse relation is, (M6.2.13a) with and (M6.2.13b) with,
33 Finally, as we have implied, there are no simplifications in the ABD matrix if the laminate is unsymmetric, if it is unbalanced, and if it is not a cross-ply. One must work with the full six-bysix matrix of Equation (M6.1.82). In Equation (M6.1.84), one must invert the full six-by-six relation. M6.2.7 Comments Whether we invert the full six-by-six ABD matrix or take advantage of the simplifications just discussed, it is clear we are now in a position to start at any point in the loop of Figure M6.4 and calculate our way to any other point. Despite this fact, it is the loads that are usually known, and these are converted to force and moment resultants by dividing by the appropriate lengths. We use the laminate material properties (materials, fiber orientations, etc.) to compute the ABD matrix. We can then compute the reference surface strains and curvatures and we can determine the stress and strain response, in the x-y-z global coordinate system, or the principal material coordinate system. An important alternative to this process is determining the laminate material properties that will produce a specific response to a known set of loads. This process is called laminate design. Optimal laminate design goes one step further and seeks the laminate material properties that minimize, for example, laminate stresses in the presence of a given set of loads. Design and optimization are important facets of composite analysis, and in many instances they are the primary activities. However, they both depend on having a thorough understanding of the mechanics of composite materials, and all of the steps that go into closing the loop of Figure M Certain subtleties are involved in some of the steps of Figure M6.1.4, and one subtlety, the coupling effects that occur in laminates, is discussed next. These couplings are reflected in the various components of the ABD matrix and its inverse. Because the ABD matrix is a stiffness matrix, these couplings are properly referred to as elastic couplings. We have discussed two of these important couplings. One was the coupling that occurs in an unsymmetric laminate due to the presence of the B matrix. These couplings do not occur with metallic materials and hence are not fully appreciated by many who do not work with composites. M6.2.8 Effective Engineering Properties of a Laminate It is often convenient to have what can be referred to as effective engineering properties of a laminate. These consist of the effective extensional modulus in the x direction, the effective extensional modulus in the y direction, the effective Poisson's ratios, and the effective shear modulus in the x-y plane. These particular effective properties can be defined for general laminates but they make the most sense when considering the inplane loading of symmetric balanced laminates. Let us restrict our discussion to the inplane loading of symmetric balanced or symmetric cross-ply laminates. In this case the identical equations (6.2.5a) and (6.2.9a) are applicable. To introduce effective engineering properties, define the average laminate stress in the x direction, the average laminate stress in the y direction, and the average laminate shear stress to be, respectively,
34 (M6.2.14) The average stresses, like all average quantities, do not really exist. They are simply a definition. In the definitions of the average stresses, the integrals are actually the definitions of the normal and shear force resultants, Equation (M6.1.1); that is: (M6.2.15) or (M6.2.16) with H being the thickness of the laminate. Knowing the force resultants and the laminate thickness, we can compute the average stresses easily. Substituting the above form into Equation (M a), we find
35 (M6.2.17) If we compare the form of this equation with the form of the compliance for a single layer in a state of plane stress, from the compliances equations, it is possible to define the above threeby-three matrix as the laminate compliance matrix. Accordingly, by analogy with compliances equations, the laminate's effective extensional modulus in the x direction, effective extensional modulus in the y direction, effective shear modulus, and two Poisson's ratios are given by, (M6.2.18) In terms of the elements of the A matrix, from Equation (M6.2.6a), (M6.2.19)
36 Of course, the reciprocity relation, (M6.2.20) is valid and so ν yx andν xy are not independent. As an example, for the four-layer[0/90] S laminate for which thickness H = Numerical value=c mm, substituting into Equation (M6.2.19) from 12 Equation (M6.1.88a), we find that the effective engineering properties are, E x = Numerical value =C GPa 13 E y = Numerical value =C GPa 14 G xy = Numerical value =C GPa 15 ν ν xy yx = Numerical value=c = Numerical value=c (M6.2.21) For this particular laminate, the x and y directions have identical effective extensional moduli, and the effective shear modulus is the same as the shear modulus for a single layer G 12. This is not surprising, considering that there are two layers with their fibers in the x direction and two with their fibers in the y direction, and that the two directions respond the same to an applied load. Furthermore, because both the 0 and the 90 layers have their fibers aligned with the global coordinate directions, the shear modulus for the laminate is identical to the shear modulus for a single layer.
Classical Lamination Theory: The Kirchhoff Hypothesis
212 Chapter 6 Classical Lamination Theory: The Kirchhoff Hypothesis In the preceding chapters we developed the tools needed to understand the elastic response of a small volume of fiber-reinforced material
More informationModule 5: Laminate Theory Lecture 17: Laminate Constitutive Relations. The Lecture Contains: Laminate Constitutive Relations
Lecture 17: Laminate Constitutive Relations The Lecture Contains: Laminate Constitutive Relations Classification of Laminates Cross-Ply Laminates Specially Orthotropic Laminates Examples Homework References
More informationChapter 5 Plane-Stress Stress-Strain Relations in a Global Coordinate System
Chapter 5 Plane-Stress Stress-Strain Relations in a Global Coordinate System One of the most important characteristics of structures made of fiber-reinforced materials, and one which dictates the manner
More informationDESIGN OF LAMINATES FOR IN-PLANE LOADING
DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationComposite Structures. Indian Institute of Technology Kanpur
Mechanics of Laminated Composite Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 23 Analysis of an Orthotropic Ply Lecture Overview Introduction Engineering constants for an 2
More informationComparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost
Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost 1 Mihir A. Mehta, 2 Satyen D. Ramani 1 PG Student, Department
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationChapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.
Chapter - Macromechanical Analysis of a Lamina Exercise Set. The number of independent elastic constants in three dimensions are: Anisotropic Monoclinic 3 Orthotropic 9 Transversely Orthotropic 5 Isotropic.
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationVIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS
VIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS Mechanical Engineering Department, Indian Institute of Technology, New Delhi 110 016, India (Received 22 January 1992,
More informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationTABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA TABLE OF CONTENTS 1. INTRODUCTION TO COMPOSITE MATERIALS 1.1 Introduction... 1.2 Classification... 1.2.1
More information3. BEAMS: STRAIN, STRESS, DEFLECTIONS
3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets
More informationFig. 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 informationLAMINATED COMPOSITE PLATES
LAMINATED COMPOSITE PLATES David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 February 10, 2000 Introduction This document is intended
More informationModule III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina
Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the
More informationON THE NUMERICAL ANALYSIS OF COMPOSITE MATERIAL
The 4th International Conference Advanced Composite Materials Engineering COMAT 8- October, Brasov, Romania O THE UMERICAL AALYSIS OF COMPOSITE MATERIAL D.D. icoara Transilvania University, City Brasov,
More informationLAMINATED COMPOSITE PLATES
LAMINATED COMPOSITE PLATES David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 February 10, 2000 Introduction This document is intended
More informationThe stiffness of plates
The stiffness of plates 1. Introduction The word plate is a collective term for elements in which forces can be transferred in two directions. Floors, walls, bridge slabs and laminates are all plates.
More informationMechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA What programs are in PROMAL? Master Menu The master menu screen with five separate applications from
More informationComposite Plate Documentation. John Robinson
Composite Plate Documentation Release v ersion John Robinson February 18, 2016 Contents 1 Indices and tables 3 2 How to Use Composite Plate 5 3 Modules 7 Python Module Index 15 i ii Contents: Contents
More informationHygrothermal stresses in laminates
Hygrothermal stresses in laminates Changing environment conditions (temperature and moisture) have an important effect on the properties which are matrix dominated. Change in temperaturet and moisture
More informationComposite Structural Mechanics using MATLAB
Session 2520 Composite Structural Mechanics using MATLAB Oscar Barton, Jr., Jacob B. Wallace United States Naval Academy Annapolis, Md 21402 Abstract In this paper MATLAB is adopted as the programming
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
More informationLAMINATION 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 informationPIEZOELECTRIC TECHNOLOGY PRIMER
PIEZOELECTRIC TECHNOLOGY PRIMER James R. Phillips Sr. Member of Technical Staff CTS Wireless Components 4800 Alameda Blvd. N.E. Albuquerque, New Mexico 87113 Piezoelectricity The piezoelectric effect is
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationMechanical Behavior of Composite Tapered Lamina
International Journal of Engineering Research and Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 10, Issue 8 (August 2014), PP.19-27 Mechanical Behavior of Composite Tapered Lamina
More informationFinite element analysis of FRP pipelines dynamic stability
Boundary Elements and Other esh Reduction ethods XXXVIII 39 Finite element analysis of FRP pipelines dynamic stability D. G. Pavlou & A. I. Nergaard Department of echanical and Structural Engineering and
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationAERSYS KNOWLEDGE UNIT
-7016 1. INTRODUCTION The scope of this document is to provide a clarification and a deeper understanding of the two different ways to move the mid plane of the element out of the nodal plane. Although
More informationMaterials and Structures. Indian Institute of Technology Kanpur
Introduction to Composite Materials and Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 34 Thermal Stresses in Plates Lecture Overview Introduction Mechanical and Thermal Strains
More informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationLinear elastic analysis of thin laminated beams with uniform and symmetric cross-section
Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,
More information(Refer Slide Time: 01:00 01:01)
Strength of Materials Prof: S.K.Bhattacharya Department of Civil Engineering Indian institute of Technology Kharagpur Lecture no 27 Lecture Title: Stresses in Beams- II Welcome to the second lesson of
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationMechanical Design in Optical Engineering
Torsion Torsion: Torsion refers to the twisting of a structural member that is loaded by couples (torque) that produce rotation about the member s longitudinal axis. In other words, the member is loaded
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationStacking sequences for Extensionally Isotropic, Fully Isotropic and Quasi-Homogeneous Orthotropic Laminates.
49th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials 7-10 April 2008, Schaumburg, IL AIAA 2008-1940 Stacking sequences for Extensionally Isotropic, Fully Isotropic and Quasi-Homogeneous
More informationTorsion/Axial Illustration: 1 (3/30/00)
Torsion/Axial Illustration: 1 (3/30/00) Table of Contents Intro / General Strategy Axial: Different Materia The Displacement Method 1 2 Calculate the Stresses General Strategy The same structure is loaded
More informationStrength of Materials Prof: S.K.Bhattacharya Dept of Civil Engineering, IIT, Kharagpur Lecture no 28 Stresses in Beams- III
Strength of Materials Prof: S.K.Bhattacharya Dept of Civil Engineering, IIT, Kharagpur Lecture no 28 Stresses in Beams- III Welcome to the 3 rd lesson of the 6 th module which is on Stresses in Beams part
More informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
More informationNomenclature. Length of the panel between the supports. Width of the panel between the supports/ width of the beam
omenclature a b c f h Length of the panel between the supports Width of the panel between the supports/ width of the beam Sandwich beam/ panel core thickness Thickness of the panel face sheet Sandwich
More informationME 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 informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
More informationExercise: 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 informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationXI. 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 informationMATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS
. MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain
More informationFinite Element Analysis of Composite Laminate By Using ABDH Matrix(Stiffness Matrix)
Finite Element Analysis of Composite Laminate By Using ABDH Matrix(Stiffness Matrix) Nikhil J. Chaudhari 1 Post Graduate Student Department of Mechanical Engineering Veermata Jijabai Technological Institute
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More information3 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 informationQUESTION BANK Composite Materials
QUESTION BANK Composite Materials 1. Define composite material. 2. What is the need for composite material? 3. Mention important characterits of composite material 4. Give examples for fiber material 5.
More informationComputational 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 informationCLASSICAL TORSION AND AIST TORSION THEORY
CLASSICAL TORSION AND AIST TORSION THEORY Background The design of a crane runway girder has not been an easy task for most structural engineers. Many difficult issues must be addressed if these members
More informationName (Print) ME Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM
Name (Print) (Last) (First) Instructions: ME 323 - Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM Circle your lecturer s name and your class meeting time. Gonzalez Krousgrill
More informationPROMAL2012 SOFTWARE PACKAGE A USER GUIDE
PROMAL2012 SOFTWARE PACKAGE A USER GUIDE 1. This manual is only for VISTA, WINDOWS 7 and WINDOWS 8 users. The PROMAL2012 software and manual are available at http://www.eng.usf.edu/~kaw/promal2012/ 2.
More informationPLAT DAN CANGKANG (TKS 4219)
PLAT DAN CANGKANG (TKS 4219) SESI I: PLATES Dr.Eng. Achfas Zacoeb Dept. of Civil Engineering Brawijaya University INTRODUCTION Plates are straight, plane, two-dimensional structural components of which
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationEFFECT OF LAMINATION ANGLE AND THICKNESS ON ANALYSIS OF COMPOSITE PLATE UNDER THERMO MECHANICAL LOADING
Journal of MECHANICAL ENGINEERING Strojnícky časopis, VOL 67 (217), NO 1, 5-22 EFFECT OF LAMINATION ANGLE AND THICKNESS ON ANALYSIS OF COMPOSITE PLATE UNDER THERMO MECHANICAL LOADING Arnab Choudhury 1,
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationCOMPOSITE PLATE THEORIES
CHAPTER2 COMPOSITE PLATE THEORIES 2.1 GENERAL Analysis of composite plates is usually done based on one of the following the ries. 1. Equivalent single-layer theories a. Classical laminate theory b. Shear
More informationAs an example we consider the stacking and analysis of a 2-ply symmetric laminate. First we clear the Matlab space and close figures.
Program Lam 1 Introduction The collection of Matlab command and function files in the package lam allows the stacking and analysis of a laminate. In fact only a representative piece of a laminate is considered,
More information8 Properties of Lamina
8 Properties of Lamina 8- ORTHOTROPIC LAMINA An orthotropic lamina is a sheet with unique and predictable properties and consists of an assemblage of fibers ling in the plane of the sheet and held in place
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More information2. Mechanics of Materials: Strain. 3. Hookes's Law
Mechanics of Materials Course: WB3413, Dredging Processes 1 Fundamental Theory Required for Sand, Clay and Rock Cutting 1. Mechanics of Materials: Stress 1. Introduction 2. Plane Stress and Coordinate
More informationINTERNATIONAL 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 informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationStresses in Curved Beam
Stresses in Curved Beam Consider a curved beam subjected to bending moment M b as shown in the figure. The distribution of stress in curved flexural member is determined by using the following assumptions:
More informationDESIGN OF COMPOSITE LAMINATED STRUCTURES BY POLAR METHOD AND TOPOLOGY OPTIMISATION
DESIGN OF COMPOSIE LAMINAED SUCUES BY POLA MEOD AND OPOLOGY OPIMISAION A. Jibawy,,3, C. Julien,,3, B. Desmorat,,4, (*), A. Vincenti UPMC Univ Paris 6, UM 79, Institut Jean Le ond d Alembert B.P. 6 4, place
More informationMECHANICS OF COMPOSITE STRUCTURES
MECHANICS OF COMPOSITE STRUCTURES LÁSZLÓ P. KOLLÁR Budapest University of Technology and Economics GEORGE S. SPRINGER Stanford University PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE
More informationModule 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 informationPassive Damping Characteristics of Carbon Epoxy Composite Plates
Journal of Materials Science and Engineering A 6 (-) 35-4 doi:.765/6-63/6.-.5 D DAVID PUBLISHING Passive Damping Characteristics of Carbon Epoxy Composite Plates Dileep Kumar K * and V V Subba Rao Faculty
More informationANALYSIS OF STRAINS CONCEPT OF STRAIN
ANALYSIS OF STRAINS CONCEPT OF STRAIN Concept of strain : if a bar is subjected to a direct load, and hence a stress the bar will change in length. If the bar has an original length L and changes by an
More informationChapter 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 informationChapter 2: Elasticity
OHP 1 Mechanical Properties of Materials Chapter 2: lasticity Prof. Wenjea J. Tseng ( 曾文甲 ) Department of Materials ngineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W.F.
More informationOn elastic lines of double curvature
Sur les lignes élastiques à double coubure Correspondance sur l École Poly. 3 (1816) 355-360. On elastic lines of double curvature By Siméon Denis Poisson Translated by D. H. Delphenich In this article
More informationEffect of Thermal Stresses on the Failure Criteria of Fiber Composites
Effect of Thermal Stresses on the Failure Criteria of Fiber Composites Martin Leong * Institute of Mechanical Engineering Aalborg University, Aalborg, Denmark Bhavani V. Sankar Department of Mechanical
More informationMATERIAL ELASTIC ANISOTROPIC command
MATERIAL ELASTIC ANISOTROPIC command.. Synopsis The MATERIAL ELASTIC ANISOTROPIC command is used to specify the parameters associated with an anisotropic linear elastic material idealization. Syntax The
More informationBeams III -- Shear Stress: 1
Beams III -- Shear Stress: 1 Internal Shear Force Shear Stress Formula for Beams The First Area Moment, Q Shear Stresses in Beam Flanges Shear Distribution on an I Beam 1 2 In this stack we will derive
More information4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL
4. BEMS: CURVED, COMPOSITE, UNSYMMETRICL Discussions of beams in bending are usually limited to beams with at least one longitudinal plane of symmetry with the load applied in the plane of symmetry or
More informationStress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina
Chapter 2 Macromechanical Analysis of a Lamina Stress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina Islamic
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More information