Review of Force, Stress, and Strain Tensors

Size: px
Start display at page:

Download "Review of Force, Stress, and Strain Tensors"

Transcription

1 Review of Fore, Stress, and Strain Tensors.1 The Fore Vetor Fores an be grouped into two broad ategories: surfae fores and body fores. Surfae fores are those that at over a surfae (as the name implies), and result from diret physial ontat between two bodies. In ontrast, body fores are those that at at a distane, and do not result from diret physial ontat of one body with another. The fore of gravity is the most ommon type of body fore. In this hapter we are primarily onerned with surfae fores, the effets of body fores (suh as the weight of a struture) will be ignored. A fore is a three-dimensional (3D) vetor. A fore is defined by a magnitude and a line of ation. In SI units, the magnitude of a fore is expressed in Newtons, abbreviated N, whereas in English units the magnitude of a fore is expressed in pounds-fore, abbreviated lbf. A fore vetor F ating at a point P and referened to a right-handed x y oordinate system is shown in Figure.1. Components of F ating parallel to the x y oordinate axes, F x, F y, and F, respetively, are also shown in the figure. The algebrai sign of eah fore omponent is defined in aordane with the algebraially positive diretion of the orresponding oordinate axis. All fore omponents shown in Figure.1 are algebraially positive, as eah omponent points in the orresponding positive oordinate diretion. The reader is likely to have enountered several different ways of expressing fore vetors in a mathematial sense. Three methods will be desribed here. The first is alled vetor notation, and involves the use of unit vetors. Unit vetors parallel to the x-, y-, and -oordinate axes are typially labeled î, ĵ, and ˆk, respetively, and by definition have a magnitude equal to unity. A fore vetor F is written in vetor notation as follows: F = F iˆ + F ˆj + F kˆ (.1) x y The magnitude of the fore is given by F = F + F + F x y (.) 41 K13483_C00.indd 41 9/13/01 11:04:14 AM

2 4 Strutural Analysis of Polymeri Composite Materials +x +y F F x F y P F + Figure.1 A fore vetor F ating at point P. Fore omponents F x, F y, and F ating parallel to the x y oordinate axes, respetively, are also shown. A seond method of defining a fore vetor is through the use of indiial notation. In this ase, a subsript is used to denote individual omponents of the vetoral quantity: F = ( Fx, Fy, F ) The subsript denotes the oordinate diretion of eah fore omponent. One of the advantages of indiial notation is that it allows a shorthand notation to be used, as follows: F = F, where i = x, y, or (.3) i Note that a range has been expliitly speified for the subsript i in Equation.3. That is, it is expliitly stated that the subsript i may take on values of x, y, or. Usually, however, the range of a subsript(s) is not stated expliitly but rather is implied. For example, Equation.3 is normally written simply as F = where it is understood that the subsript i takes on values of x, y, and. The third approah is alled matrix notation. In this ase, individual omponents of the fore vetor are listed within braes in the form of a olumn array: F i Fx F = Fy F (.4) K13483_C00.indd 4 9/13/01 11:04:16 AM

3 Review of Fore, Stress, and Strain Tensors 43 Indiial notation is sometimes ombined with matrix notation as follows: F = { F i } (.5). Transformation of a Fore Vetor One of the most ommon requirements in the study of mehanis is the need to desribe a vetor in more than one oordinate system. For example, suppose all omponents of a fore vetor F i are known in one oordinate system (the x y oordinate system, say) and it is desired to express this fore vetor in a seond oordinate system (the x y oordinate system, say). To desribe the fore vetor in the new oordinate system, we must alulate the omponents of the fore parallel to the x, y, and axes that is, we must alulate F x, F y, and F. The proess of relating fore omponents in one oordinate system to those in another oordinate system is alled transformation of the fore vetor. This terminology is perhaps unfortunate, in the sense that the fore vetor itself is not transformed but rather our desription of the fore vetor transforms as we hange from one oordinate system to another. It an be shown [1,] that the fore omponents in the x y oordinate system (F x, F y, and F ) are related to the omponents in x y oordinate system (F x, F y, F ) aording to: F = F + F + F x x x x x y y x F = F + F + F y y x x y y y y (.6a) F = F + F + F x x y y The terms i j that appear in Equation.6a are alled diretion osines, and are to equal the osine of the angle between the axes of the new and original oordinate systems. An angle of rotation is defined from the original x y oordinate system to the new x y oordinate system. The algebrai sign of the angle of rotation is defined in aordane with the right-hand rule. Equation.6a an be suintly written using the summation onvention as follows: F = F (.6b) i i j j Alternatively, these three equations an be written using matrix notation as Fx Fy = F x x x y x y x y y y x y Fx Fy F (.6) K13483_C00.indd 43 9/13/01 11:04:16 AM

4 44 Strutural Analysis of Polymeri Composite Materials Note that although values of individual fore omponents vary as we hange from one oordinate system to another, the magnitude of the fore vetor (given by Equation.) does not. The magnitude is independent of the oordinate system used, and is alled an invariant of the fore tensor. Diretion osines relate unit vetors in the new and old oordinate systems. For example, a unit vetor direted along the x -axis (i.e., unit vetor i ˆ ) is related to unit vetors in the x y oordinate system as follows: iˆ = iˆ + ˆj + kˆ (.7) x x x y x As i ˆ is a unit vetor, then in aordane with Equation.: ( ) + ( ) + ( ) = 1 (.8) x x x y x To this point we have referred to a fore as a vetor. A fore vetor an also be alled a fore tensor. The term tensor refers to any quantity that transforms in a physially meaningful way from one Cartesian oordinate system to another. The rank of a tensor equals the number of subsripts that must be used to desribe the tensor. A fore an be desribed using a single subsript, F i, and therefore a fore is said to be a tensor of rank one, or equivalently, a firstorder tensor. Equation.6 is alled the transformation law for a first-order tensor. It is likely that the reader is already familiar with two other tensors: the stress tensor, σ ij, and the strain tensor, ε ij. The stress and strain tensors will be reviewed later in this hapter, but at this point it an be noted that two subsripts are used to desribe stress and strain tensors. Hene, stress and strain tensors are said to be tensors of rank two, or equivalently, seond-order tensors. Example Problem.1 Given: All omponents of a fore vetor F are known in a given x y oordinate system. It is desired to express this fore in a new x y oordinate system, where the x y system is generated from the original x y system by the following two rotations (see Figure.): A rotation of θ-degrees about the original -axis (whih defines an intermediate x y oordinate system), followed by A rotation of β-degrees about the x -axis (whih defines the final x y oordinate system) Problem a. Determine the diretion osines i j relating the original x y oordinate system to the new x y oordinate system. b. Obtain a general expression for the fore vetor F in the x y oordinate system. K13483_C00.indd 44 9/13/01 11:04:18 AM

5 Review of Fore, Stress, and Strain Tensors 45 (a) +x +θ +x +y +y +θ +, + +θ (b) +x, +x +β +y +y +β + +β + Figure. Generation of the x y oordinate system from the x y oordinate system. (a) Rotation of θ-degrees about the original -axis (whih defines an intermediate x y oordinate system); (b) rotation of β-degrees about the x -axis (whih defines the final x y oordinate system).. Calulate numerial values of the fore vetor F in the x y oordinate system if θ = 0, β = 60, and F x = 1000 N, F y = 00 N, F = 600 N. Solution a. One way to determine diretion osines i j is to rotate unit vetors. In this approah unit vetors are first rotated from the original x y oordinate system to the intermediate x y oordinate system, and then from the x y system to the final x y oordinate system. Define a unit vetor I that is aligned with the x-axis: I ( 1) i That is, vetor I is a vetor for whih I x = 1, I y = 0, and I = 0. The vetor I an be rotated to the intermediate x y oordinate system using Equation.6: K13483_C00.indd 45 9/13/01 11:04:19 AM

6 46 Strutural Analysis of Polymeri Composite Materials I = I + I + I x x x x x y y x I = I + I + I y y x x y y y y I = I + I + I x x y y The diretion osines assoiated with a transformation from the x y oordinate system to the intermediate x y oordinate system an be determined by inspetion (see Figure.a), and are given by = osine(angle between x - and x-axes) = os θ x x x y x = osine( angle between x - and y-axes) = os( 90 θ) = sin θ = osine( angle between x - and -axes) = os( 90 ) = 0 y x y y = osine(angle between the y - and x-axes) = os( 90 + θ) = sin θ = osine( angle between the y - and y-axes) = os θ = osine(angle between the y - and -axes) = os( 90 ) = 0 y x y = osine(angle between the - and x-axes) = os( 90 ) = 0 = osine(angle between the - and y-axes) = os( 90 ) = 0 = osine(angle between the - and -axes) = os( 0 ) = 1 Using these diretion osines: Ix = x xix + x yiy + x I = (os θ)( 1) + (sin θ)( 0) + ( 0)( 0) = os θ Iy = y xix + y yiy + y I = ( sin θ)( 1) + (os θ)( 0) + ( 0)( 0) = sinθ I = xix + yiy + I = ( 0)( 1) + ( 0)( 0) + ( 1)( 0) = 0 Therefore, in the x y oordinate system, the vetor I is written: I = (os θ) i + ( sin θ) j Now define two additional unit vetors, one aligned with the original y-axis (vetor J ) and one aligned with the original -axis (vetor K ); that is, let J = ( 1 ) j and K = ( 1 ) k. Transforming these vetors to the x y oordinate system, again using the diretion osines listed above, results in J = (sin θ) i + ( os θ) j K13483_C00.indd 46 9/13/01 11:04:19 AM

7 Review of Fore, Stress, and Strain Tensors 47 K = ( 1 ) k We now rotate vetors I, J and K from the intermediate x y oordinate system to the final x y oordinate system. The diretion osines assoiated with a transformation from the x y oordinate system to the final x y oordinate system are easily determined by inspetion (see Figure.b) and are given by = 1 = 0 = 0 x x x y x = 0 = os β = sin β y x y y y = 0 = sin β = os β x y These diretion osines together with Equation.6 an be used to rotate the vetor I from the intermediate x y oordinate system to the final x y oordinate system: I = I + I + I = ( 1)(os θ) + ( 0)( sin θ) + ( 0)( 0) x x x x x y y x I x = os θ I = I + I + I = ( 0)(os θ) + (os β)( sin θ) + (sin β)( 0) y y x x y y y y I y = os β sin θ I = I + I + I = ( 0)(os θ) + ( sin β)( sin θ) + (os β)( 0) x x y y I = sin β sin θ Therefore, in the final x y oordinate system, the vetor I is written: I = (os θ) i + ( os βsin θ) j + (sin β sin θ) k (a) Reall that in the original x y oordinate system I is simply a unit vetor aligned with the original x-axis: I ( 1 ) i. Therefore, result (a) defines the diretion osines assoiated with the angle between the original x-axis and the final x -, y -, and -axes. That is x x = osθ y x x = os β sin θ = sin β sin θ K13483_C00.indd 47 9/13/01 11:04:19 AM

8 48 Strutural Analysis of Polymeri Composite Materials A similar proedure is used to rotate the unit vetors J and K from the intermediate x y oordinate system to the final x y oordinate system. These rotations result in J = (sin θ) i + (os βos θ) j + ( sin β os θ) k (b) K = ( 0) i + (sin β) j + (os β ) k () As vetor J is a unit vetor aligned with the original y-axis, J = ( 1) j, result (b) defines the diretion osines assoiated with the angle between the original y-axis and the final x -, y -, and -axes: = sin θ x y y y y = os β os θ = sin β os θ Finally, result () defines the diretion osines assoiated with the angle between the original -axis and the final x -, y -, and -axes: = 0 x y = sin β = os β Assembling the preeding results, the set of diretion osines relating the original x y oordinate system to the final x y oordinate system an be written: x x x y x y x y y y x y osθ sinθ 0 = osβ sinθ osβ osθ sinβ sinβ sinθ sinβ osθ osβ b. As diretion osines have been determined, transformation of fore vetor F an be aomplished using any version of Equation.6. For example, using matrix notation, Equation.6: Fx Fy = F x x x y x y x y y y x y Fx Fy F os θ sin θ 0 Fx = os β sin θ os β os θ sin β Fy sin β sin θ sin β os θ os β F K13483_C00.indd 48 9/13/01 11:04:0 AM

9 Review of Fore, Stress, and Strain Tensors 49 Fx (os θ) Fx + (sin θ) Fy Fy = ( osβsin θ) Fx + (osβos θ) Fy + (sin β) F F (sinβsin θ) Fx + ( sinβos θ) Fy + (os β) F. Using the speified numerial values and the results of part (b): [os( 0 )]( 1000 N) + [sin( 0 )] 00 N Fx [ os( 60 )(sin 0 )]( 1000 N) + [os( 60 )os( 0 )]( 00 N) Fy = + [sin( 60 )]( 600 N) F [sin( 60 )sin( 0 )]( 1000 N) + [ sin( 60 )os( 0 )]( 00 N) + (os 60 )( 600 N) Fx N Fy = N F N Using vetor notation, F an now be expressed in the two different oordinate systems as or, equivalently F = (1000 N) iˆ + (00 N) ˆj + ( 600 N) kˆ F = (873.1 N) i + (784.6 N) j + ( N) k Where iˆ, ˆ, j k ˆ and i ˆ, ˆ j, k ˆ are unit vetors in the x y and x y oordinate systems, respetively. Fore vetor F drawn in the x y and x y oordinate systems is shown in Figure.3a and b, respetively. The two desriptions of F are entirely equivalent. A onvenient way of (partially) verifying this equivalene is to alulate the magnitude of the original and transformed fore vetors. As the magnitude is an invariant, it is independent of the oordinate system used to desribe the fore vetor. Using Equation., the magnitude of the fore vetor in the x y oordinate system is F = F + F + F = ( 1000 N) + ( 00 N) + ( 600 N) = 1183 N x y The magnitude of the fore vetor in the x y oordinate system is x y F = ( F ) + ( F ) + ( F ) = ( N) + ( N) + ( N) = 1183 N (agrees) K13483_C00.indd 49 9/13/01 11:04:0 AM

10 50 Strutural Analysis of Polymeri Composite Materials (a) +x +y F F x = 1000 N F = 600 N F y = 00 N + (b) F +x F x = N +y F y = N F = N + Figure.3 Fore vetor F drawn in two different oordinate systems. (a) Fore vetor F in the original x y oordinate system; (b) fore vetor F in a new x y oordinate system..3 Normal Fores, Shear Fores, and Free-Body Diagrams Fore F ating at an angle to a planar surfae is shown in Figure.4. As fore is a vetor, it an always be deomposed into two fore omponents, a normal fore omponent and a shear fore omponent. The line-of-ation of the normal fore omponent is orthogonal to the surfae, whereas the line-of-ation of the shear fore omponent is tangent to the surfae. Internal fores indued within a solid body by externally applied fores an be investigated with the aid of free-body diagrams. A simple example is shown in Figure.5, whih shows a straight irular rod with onstant diameter subjeted to two external fores of equal magnitude (R) but opposite diretion. The internal fore (FI, say) indued at any ross-setion of the rod an be investigated by making an imaginary ut along the plane of interest. Suppose an imaginary ut is made along plane a a a a, whih is perpendiular to the axis of the rod. The resulting free-body diagram for the lower half of the rod is shown in Figure.5a, where a x y oordinate system has been assigned suh that the x-axis is parallel to the rod axis, as shown. On the basis of this free-body diagram, it is onluded that an internal fore FI = ( R) i + ( 0) j + ( 0) k is indued at ross-setion a a. That K13483_C00.indd 50 9/13/01 11:04:0 AM

11 Review of Fore, Stress, and Strain Tensors 51 F Normal fore omponent Shear fore omponent Planar surfae Figure.4 A fore F ating at an angle to a planar surfae. (a) R +x R a a a a +y + R R (b) R +x +θ +x a b b a a b b a +θ R os θ +y +θ +y +,+ R os θ R R Figure.5 The use of free-body diagrams to determine internal fores ating on planes a a a a and b b b b. (a) Free-body diagram based on plane a a a a, perpendiular to rod axis; (b) free-body diagram based on plane b b b b, inlined at angle +θ to the rod axis. is, only a normal fore of magnitude R is indued at ross-setion a a a a, whih has been defined to be perpendiular to the axis of the rod. In ontrast, the imaginary ut need not be made perpendiular to the axis of the rod. Suppose the imaginary ut is made along plane b b b b, whih is inlined at an angle of +θ with respet to the axis of the rod. The resulting K13483_C00.indd 51 9/13/01 11:04:0 AM

12 5 Strutural Analysis of Polymeri Composite Materials free-body diagram for the lower half of the rod is shown in Figure.5b. A new x y oordinate system has been assigned so that the x -axis is perpendiular to plane b b b b and the -axis is oinident with the -axis, that is, the x y oordinate system is generated from the x y oordinate system by a rotation of +θ about the original -axis. The internal fore FI an be expressed with respet to the x y oordinate system by transforming FI from the x y oordinate system to the x y oordinate system. This oordinate transformation is a speial ase of the transformation onsidered in Example Problem.1. The diretion osines now beome (with β = 0 ): = os θ = sin θ = 0 = sin θ = os θ = 0 x x x y x y x y y y = 0 x y Applying Equation.6, we have Fx Fy = F x x x y x y x y y y x y = 0 = 1 Fx os θ sin θ 0 R (os θ) R Fy = sin θ os θ 0 0 = ( sin θ) R F In the x y oordinate system the internal fore is FI = ( Ros θ) i ( Rsin θ ) j + ( 0 ) k. Hene, by defining a oordinate system whih is inlined to the axis of the rod, we onlude that both a normal fore (R os θ) and a shear fore ( R sin θ) are indued in the rod. Although the preeding disussion may seem simplisti, it has been inluded to demonstrate the following: A speifi oordinate system must be speified before a fore vetor an be defined in a mathematial sense. In general, the oordinate system is defined by the imaginary ut(s) used to form the free-body diagram. All omponents of a fore must be speified to fully define the fore vetor. Further, the individual omponents of a fore hange as the vetor is transformed from one oordinate system to another. These two observations are valid for all tensors, not just for fore vetors. In partiular, these observations hold in the ase of stress and strain tensors, whih will be reviewed in the following setions..4 Definition of Stress There are two fundamental types of stress: normal stress and shear stress. Both types of stress are defined as a fore divided by the area over whih it ats. K13483_C00.indd 5 9/13/01 11:04:1 AM

13 Review of Fore, Stress, and Strain Tensors 53 A general 3D solid body subjeted to a system of external fores is shown in Figure.6a. It is assumed that the body is in stati equilibrium, that is, it is assumed that the sum of all external fores is ero, ΣFi = 0. These external fores indue internal fores ating within the body. In general, the internal fores will vary in both magnitude and diretion throughout the body. An illustration of the variation of internal fores along a single line within an internal plane is shown in Figure.6b. A small area (ΔA) isolated from this plane is shown in Figure.6. Area ΔA is assumed to be infinitesimally small. That (a) F 1 F F 3 F 5 F 4 (b) Internal fores along a single line F 3 F 5 F 4 () V N ΔA F 3 F 5 F 4 Figure.6 A solid 3-D body in equilibrium. (a) A solid 3-D body subjet to external fores F1 F5; (b) variation of internal fores along an internal line; () internal fore ating over infinitesimal area ΔA. K13483_C00.indd 53 9/13/01 11:04:3 AM

14 54 Strutural Analysis of Polymeri Composite Materials is, the area ΔA is small enough suh that the internal fores ating over ΔA an be assumed to be of onstant magnitude and diretion. Therefore, the internal fores ating over ΔA an be represented by a fore vetor whih an be deomposed into a normal fore, N, and a shear fore, V, as shown in Figure.6. Normal stress (usually denoted σ ), and shear stress (usually denoted τ ) are defined as the fore per unit area ating perpendiular and tangent to the area ΔA, respetively. That is, N σ lim A V τ lim A A 0 A 0 (.9) Note that by definition the area ΔA shrinks to ero: ΔA 0. Stresses σ and τ are therefore said to exist at a point. As internal fores generally vary from point-to-point (as shown in Figure.6), stresses also vary from point-to-point. Stress has units of fore per unit area. In SI units stress is reported in terms of Pasals (abbreviated Pa ), where 1 Pa = 1 N/m. In English units stress is reported in terms of pounds-fore per square inh (abbreviated psi ), that is, 1 psi = 1 lbf/in. Conversion fators between the two systems of measurement are 1 psi = 6895 Pa, or equivalently, 1 Pa = psi. Common abbreviations used throughout this hapter are as follows: Pa = 1 kilo-pasals = 1 kpa psi = 1 kilo-psi = 1 ksi Pa = 1 Mega-Pasals = 1 MPa psi = 1 mega-psi = 1 Msi Pa = 1 Giga-Pasals = 1 GPa..5 The Stress Tensor A general 3D solid body subjeted to a system of external fores is shown in Figure.7a. It is assumed that the body is in stati equilibrium and that body fores are negligible, that is, is it assumed that the sum of all external fores is ero, ΣFi = 0. A free-body diagram of an infinitesimally small ube removed from the body is shown in Figure.7b. The ube is referened to a x y oordinate system, and the ube edges are aligned with these axes. The lengths of the ube edges are denoted dx, dy, and d. Although (in general) internal fores are indued over all six faes of the ube, for larity the fores ating on only three faes have been shown. The fore ating over eah ube fae an be deomposed into a normal fore omponent and two shear fore omponents, as illustrated in Figure.7. Although eah fore omponent ould be identified with a single subsript (as fore is a first-order tensor), for onveniene two subsripts have been used. The first subsript identifies the fae over whih the fore is distributed, whereas the seond subsript identifies the diretion in whih the K13483_C00.indd 54 9/13/01 11:04:3 AM

15 Review of Fore, Stress, and Strain Tensors 55 (a) (b) F 1 +x F +y F 3 F 5 F 4 + () V xy N xx +x (d) τ xy σ xx +x V x V +y yx +y τ x τ yx V y τ y V y τ y N N yy σ σ yy + V x + τ x Figure.7 Free-body diagrams used to define stress indued in a solid body. (a) 3-D solid body in equilibrium; (b) infinitesimal ube removed from the solid body (internal fores ating on three faes shown); () normal fore and two shear fores at over eah fae of the ube; (d) normal stress and two shear stresses at over eah fae of the ube. fore is oriented. For example, N xx refers to a normal fore omponent whih is distributed over the x-fae and whih is ating parallel to the x-diretion. Similarly, V y refers to a shear fore distributed over the -fae whih is ating parallel to the y-diretion. Three stress omponents an now be defined for eah ube fae, in aordane with Equation.9. For example, for the three faes of the infinitesimal element shown in Figure.7: Stresses ating on the -x-fae: Nxx Vxy σxx = lim τxy τ y y = lim d, d 0 d d dy, d 0 dy d x V = lim x dy, d 0 dy d Stresses ating on the y-fae: Nyy Vyx σ yy = lim τyx τ x x = lim d, d 0 d d dx, d 0 dx d y V y = lim d x, d 0 dx d K13483_C00.indd 55 9/13/01 11:04:4 AM

16 56 Strutural Analysis of Polymeri Composite Materials Stresses ating on the +-fae: N Vx V y σ = lim τx τy x y x y = lim x y x y = lim d, d 0 d d d, d 0 d d d x, d y 0 dx dy As three fore omponents (and therefore three stress omponents) exist on eah of the six faes of the ube, it would initially appear that there are 18 independent fore (stress) omponents. However, it is easily shown that for stati equilibrium to be maintained (assuming body fores are negligible): Normal fores ating on opposite faes of the infinitesimal element must be of equal magnitude and opposite diretion. Shear fores ating within a plane of the element must be orientated either tip-to-tip (e.g., fores V x and V x in Figure.7) or tail-totail (e.g., fores V xy and V yx ), and be of equal magnitude. That is V xy = V yx, V x = V x, V y = V y. These restritions redue the number of independent fore (stress) omponents from 18 to 6, as follows: Independent Fore Component(s) Corresponding Stress Component(s) N xx N yy σ xx σ yy N σ V xy (=V yx ) τ xy (= τ yx ) V x (=V x ) τ x (= τ x ) V y (=V y ) τ y (= τ y ) We must next define the algebrai sign onvention we will use to desribe individual stress omponents. The omponents ating on three faes of an infinitesimal element are shown in Figure.8. We first assoiate an algebrai sign with eah fae of the infinitesimal element. A ube fae is positive if the outward unit normal of the fae (i.e., the unit normal pointing away from the interior of the element) points in a positive oordinate diretion; otherwise, the fae is negative. For example, faes (ABCD) and (ADEG) in Figure.8 are a positive faes, whereas fae (DCFE) is a negative fae. Having identified the positive and negative faes of the element, a stress omponent is positive if The stress omponent ats on a positive fae and points in a positive oordinate diretion, or if The stress omponent ats on a negative fae and points in a negative oordinate diretion. If neither of these onditions are met, then the stress omponent is negative. This onvention an be used to onfirm that all stress omponents K13483_C00.indd 56 9/13/01 11:04:4 AM

17 Review of Fore, Stress, and Strain Tensors 57 +x σ xx B +y A τ x C τ y τ xy D τ yx τ y σ G τ x F σ yy E + Figure.8 An infinitesimal stress element (all stress omponents shown in a positive sense). shown in Figure.8 are algebraially positive. For example, to determine the algebrai sign of the normal stress σ xx whih ats on fae ABCD in Figure.8, note that (a) fae ABCD is positive, and (b) the normal stress σ xx whih ats on this fae points in the positive x-diretion. Therefore, σ xx is positive. As a seond example, the shear stress τ y whih ats on ube fae DCFE is positive beause (a) fae DCFE is a negative fae, and (b) τ y points in the negative -diretion. The preeding disussion shows that the state of stress at a point is defined by six omponents of stress: three normal stress omponents and three shear stress omponents. The state of stress is written using matrix notation as follows: σ τ τ τ σ τ τ τ σ xx xy x yx yy y x y σ τ τ = τ σ τ τx τ σ xx xy x xy yy y y (.10) To express the state of stress using indiial notation we must first make the following hange in notation: τ σ xy xy τ τ τ τ τ x yx y x y σ σ σ σ σ x yx y x y K13483_C00.indd 57 9/13/01 11:04:5 AM

18 58 Strutural Analysis of Polymeri Composite Materials With this hange the matrix on the left side of the equality sign in Equation.10 beomes: σ τ τ τ σ τ τ τ σ xx xy x yx yy y x y σ σ σ σ σ σ σ σ σ xx xy x yx yy y x y Whih an be suintly written using indiial notation as σ ij, i, j = x, y, or (.11) In Setion.1 it was noted that a fore vetor is a first-order tensor, as only one subsript is required to desribe a fore tensor, F i. From Equation.11 is it lear that stress is a seond-order tensor (or equivalently, a tensor of rank two), as two subsripts are required to desribe a state of stress. Example Problem. Given: The stress element referened to a x y oordinate system and subjet to the stress omponents shown in Figure.9. Determine: Label all stress omponents, inluding algebrai sign. Solution The magnitude and algebrai sign of eah stress omponent is determined using the sign onvention defined above. The proedure will +x B 100 +y A D C 5 G 30 F 50 E + Figure.9 Stress omponents ating on an infinitesimal element (all stresses in MPa). K13483_C00.indd 58 9/13/01 11:04:6 AM

19 Review of Fore, Stress, and Strain Tensors 59 be illustrated using the stress omponents ating on fae DCFE. First, note that fae DCFE is a negative fae, as an outward unit normal for this fae points in the negative y-diretion. The normal stress whih ats on fae DCFE has a magnitude of 50 MPa, and points in the positive y-diretion. Hene, this stress omponent is negative and is labeled σ yy = 50 MPa. One of the shear stress omponents ating on fae DCFE has a magnitude of 75 MPa, and points in the positive x-diretion. Hene, this stress omponent is also negative and is labeled τ yx = 75 MPa (or equivalently, τ xy = 75 MPa). Finally, the seond shear fore omponent ating on fae DCFE has a magnitude of 50 MPa, and points in the positive -diretion. Hene, this omponent is labeled τ y = 50 MPa (or equivalently, τ y = 50 MPa). Following this proess for all faes of the element, the state of stress represented by the element shown in Figure.9 an be written: σ τ τ τ σ τ τ τ σ xx xy x yx yy y x y 100 MPa 75MPa 30 MPa = 75MPa 50 MPa 50 MPa 30 MPa 50 MPa 5 MPa.6 Transformation of the Stress Tensor In Setion.5 the stress tensor was defined using a free-body diagram of an infinitesimal element removed from a 3D body in stati equilibrium. This onept is again illustrated in Figure.10a, whih shows the stress element referened to an x y oordinate system. Now, the infinitesimal element need not be removed in the orientation shown in Figure.10a. An infinitesimal element removed from preisely the same point within the body but at a different orientation is shown in Figure.10b. This stress element is referened to a new x y oordinate system. The state of stress at the point of interest is ditated by the external loads applied to the body, and is independent of the oordinate system used to desribe it. Hene, the stress tensor referened to the x y oordinate system is equivalent to the stress tensor referened to the x y oordinate system, although the diretion and magnitude of individual stress omponents will differ. The proess of relating stress omponents in one oordinate system to those in another is alled transformation of the stress tensor. This terminology is perhaps unfortunate, in the sense that the state of stress itself is not transformed but rather our desription of the state of stress transforms as we hange from one oordinate system to another. K13483_C00.indd 59 9/13/01 11:04:6 AM

20 60 Strutural Analysis of Polymeri Composite Materials (a) F 1 F +x τ xy σ xx τ x +y τ yx τ y F 3 τ y σ σ yy F 5 F 4 + τ x (b) F 1 σ τ F x x x +x +y τ xy τ y x F 3 τ y τ y σ y y σ τ x F 5 + F 4 Figure.10 Infinitesimal elements removed from the same point within a 3D solid but in two different orientations. (a) Infinitesimal element referened to the x y oordinate system; (b) infinitesimal element referened to the x y oordinate system. It an be shown [1,] that the stress omponents in the new x y oordinate system (σ i j ) are related to the omponents in the original x y oordinate system (σ ij ) aording to σi j = i k j l σkl where i, j, k, l = x, y, (.1a) Or, equivalently (using matrix notation): [ σ ] = [ ][ σ ][ ] T i j i j ij i j K13483_C00.indd 60 9/13/01 11:04:7 AM

21 Review of Fore, Stress, and Strain Tensors 61 where [ i j ] T is the transpose of the diretion osine array. Writing in full matrix form: σ σ σ σ σ σ σ σ σ x x x y x y x y y y x y = x x x y x y x y y y x y σ σ σ σ yx σ yy σ σ σ σ xx xy x y x y x x y x x x y y y y x y (.1b) As disussed in Setion., the terms i j whih appear in Equation.1a,b are diretion osines and equal the osine of the angle between the axes of the x y and x y oordinate systems. Reall that the algebrai sign of an angle of rotation is defined in aordane with the right-hand rule, and that angles are defined from the x y oordinate system to x y oordinate system. Equation.1a,b is alled the transformation law for a seondorder tensor. If an analysis is being performed with the aid of a digital omputer, whih nowadays is almost always the ase, then matrix notation (Equation.1b) will most likely be used to transform a stress tensor from one oordinate system to another. Conversely, if a stress transformation is to be aomplished using hand alulations, then indiial notation (Equation.1a) may be the preferred hoie. To apply Equation.1a, the stress omponent of interest is speified by seleting the appropriate values for subsripts i and j, and then the terms on the right side of the equality are summed over the entire range of the remaining two subsripts, k and l. For example, suppose we wish write the relationship between σ x and the stress omponents in the x y oordinate system in expanded form. We first speify that i = x and j =, and Equation.1a beomes: σx = x k l σkl where k, l = x, y, We then sum all terms on the right side of the equality formed by yling through the entire range of k and l. In expanded form, we have σ = σ + σ + σ x x x x xx x x y xy x x x + σ + σ + σ + x y x yx x y y yy x y y σ + σ + σ (.13) x x x x y y x K13483_C00.indd 61 9/13/01 11:04:8 AM

22 6 Strutural Analysis of Polymeri Composite Materials Equations.1a and.1b show that the value of any individual stress omponent σ i j varies as the stress tensor is transformed from one oordinate system to another. However, it an be shown [1,] that there are features of the total stress tensor that do not vary when the tensor is transformed from one oordinate system to another. These features are alled the stress invariants. For a seond-order tensor three independent stress invariants exist, and are defined as follows: First stress invariant = Θ = σ ii (.14a) Seond stress invariant = Φ = 1 ( σ σ σ σ ) ii jj ij ij (.14b) Third stress invariant = Ψ = 1 (σ σ σ 3σ σ σ + σ σ σ 6 ii jj kk ii jk jk ij jk ki) (.14) Alternatively, by expanding these equations over the range i,j,k = x,y, and simplifying, the stress invariants an be written: First stress invariant = Θ = σxx + σ yy + σ (.15a) Seond stress invariant = Φ = σ σ + σ σ + σ σ ( σ + σ + xx yy xx yy xy x Thirdstressinvariant = Ψ = σ σ σ σ σ σ σ σ σ + xx yy xx y yy x xy σ y ) (.15b) σ σ σ xy x y (.15) The three stress invariants are oneptually similar to the magnitude of a fore tensor. That is, the value of the three stress invariants is independent of the oordinate used to desribe the stress tensor, just as the magnitude of a fore vetor is independent of the oordinate system used to desribe the fore. This invariane will be illustrated in the following Example Problem. Example Problem.3 Given: A state of stress referened to a x y oordinate is known to be σ σ σ σ σ σ σ σ σ xx xy x yx yy y x y = ( ksi) K13483_C00.indd 6 9/13/01 11:04:9 AM

23 Review of Fore, Stress, and Strain Tensors 63 It is desired to express this state of stress in an x y oordinate system, generated by the following two sequential rotations: i. Rotation of θ = 0 about the original -axis (whih defines an intermediate x y oordinate system), followed by ii. Rotation of β = 35 about the x -axis (whih defines the final x y oordinate system) Problem a. Rotate the stress tensor to the x y oordinate system. b. Calulate the first, seond, and third invariants of the stress tensor using (i) elements of the stress tensor referened to the x y oordinate system, σ ij, and (ii) elements of the stress tensor referened to the x y oordinate system, σ i j. Solution a. General expressions for diretion osines relating the x y and x y oordinate systems were determined as a part of Example Problem.1. The diretion osines were found to be x x x y x y x y y y x y = osθ = sin θ = 0 = os β sin θ = os β os θ = sin β = sin β sin θ = sin β os θ = os β As in this problem θ = 0 and β = 35, the numerial values of the diretion osines are x x x y x y x = os( 0 ) = = sin( 0 ) = = 0 = os( 35 )sin( 0 ) = y y = os( 35 )os( 0 ) = y = sin( 35 ) = x = sin( 35 )sin( 0 ) = y = sin( 35 )os( 0 ) = = os( 35 ) = K13483_C00.indd 63 9/13/01 11:04:30 AM

24 64 Strutural Analysis of Polymeri Composite Materials Eah omponent of the transformed stress tensor is now found through appliation of either Equation.1a or.1b. For example, if indiial notation is used stress omponent σ x an be found using Equation.13: σ = σ + σ + σ x x x x xx x x y xy x x x + σ + σ + σ x y x yx x y y yy x y y + σ + σ + σ x x x x y y x σ x = ( )( ) ( 50ksi) + ( )( ) ( 10ksi) + ( )( ) ( 15ksi) + ( ) ( ) ( 10ksi) + ( )( ) ( 5ksi) + ( )( ) ( 30ksi) + ( 0) ( ) ( 15ksi) + ( 0)( )( 30ksi) + ( 0)( )( 5ksi) σ x = 8. 95ksi Alternatively, if matrix notation is used, then Equation.1b beomes: σ σ σ σ σ σ σ σ σ x x x y x y x y y y x y = Completing the matrix multipliation indiated, there results: σ σ σ σ σ σ σ σ σ x x x y x y x y y y x y = ( ksi) Notie that the value of σ x determined through matrix multipliation is idential to that obtained using indiial notation, as previously desribed. The stress element is shown in the original and final oordinate systems in Figure.11. b. The first, seond, and third stress invariants will now be alulated using omponents of both σ ij and σ i j. It is expeted K13483_C00.indd 64 9/13/01 11:04:31 AM

25 Review of Fore, Stress, and Strain Tensors 65 (a) 50 +x (b) x y y Figure.11 Stress tensor of Example Problem.3, referened to two different oordinate systems (magnitude of all stress omponents in ksi). (a) Referened to x y oordinate system; (b) referened to x y oordinate system. that idential values will be obtained, as the stress invariants are independent of oordinate system. First Stress Invariant: x y oordinate system: Θ = σ = σ + σ + σ ii xx yy Θ = ( ) ksi Θ = 70 ksi x y oordinate system: Θ = σ = σ + σ + σ i i x x y y Θ = ( ) Θ = 70ksi As expeted, the first stress invariant is independent of oordinate system. Seond Stress Invariant: x y oordinate system: Φ = ( ) 1 ( σ σ σ σ ) = σ σ + σ σ + σ σ σ + σ + σ ii jj ij ij xx yy xx yy xy x y K13483_C00.indd 65 9/13/01 11:04:3 AM

26 66 Strutural Analysis of Polymeri Composite Materials { + + } Φ = ( 50)( 5) + ( 50)( 5) + ( 5)( 5) [( 10) ( 15) ( 30) ] ( ksi) Φ = 350( ksi) x y oordinate system: Φ = 1 ( ) σ i i σ j j σ i j σ i j ( ) Φ = σ σ + σ σ + σ σ σ + σ + σ x x y y x x y y x y x y { Φ = ( )( ) + ( )( 13. 7) + ( )( 13. 7) [( ) + ( 8. 95) + ( ) ]}( ksi) Φ = 350( ksi) As expeted, the seond stress invariant is independent of oordinate system. Third Stress Invariant: x y oordinate system: 1 Ψ = σ σ σ σ σ σ σ σ σ ( ii jj kk ii jk jk ij jk ki ) Ψ = σ σ σ σ σ σ σ σ σ + σ σ σ xx yy xx y yy x xy xy x y Ψ = [( 50)( 5)( 5) ( 50)( 30) ( 5)( 15) ( 5)( 10) + ( 10)( 15)( 30)]( ksi) 3 Ψ = 65375( ksi) 3 x y oordinate system: 1 Ψ = σ σ σ σ σ σ σ σ σ ( i i j j k k i i j k j k i j j k k i ) Ψ = σ σ σ σ σ σ σ σ σ + σ x x y y x x y y y x x y x y σx σ y K13483_C00.indd 66 9/13/01 11:04:36 AM

27 Review of Fore, Stress, and Strain Tensors 67 Ψ = [( )( )( 13. 7) ( )( ) ( )( 8. 95) ( 13. 7)( ) + ( )( 8. 95)( )]( ksi) 3 Ψ = 65375( ksi) 3 As expeted, the third stress invariant is independent of oordinate system..7 Prinipal Stresses It is always possible to rotate the stress tensor to a speial oordinate system in whih no shear stresses exist. This oordinate system is alled the prinipal stress oordinate system, and the normal stresses that exist in this oordinate system are alled prinipal stresses. Prinipal stresses an be used to predit failure of isotropi materials. Therefore, knowledge of the prinipal stresses and orientation of the prinipal stress oordinate system is of vital importane during design and analysis of isotropi metal strutures. This is not the ase for anisotropi omposite materials. Failure of omposite material is not governed by prinipal stresses. Prinipal stresses are only of oasional interest to the omposite engineer and are reviewed here only in the interests of ompleteness. Prinipal stresses are usually denoted as σ 1, σ, and σ 3. However, in the study of omposites, the labels 1,, and 3 are used to label a speial oordinate system alled the prinipal material oordinate system. Therefore, in this hapter the axes assoiated with the prinipal stress oordinate system will be labeled the p 1, p, and p 3 axes, and the prinipal stresses will be denoted as σ p1, σ p, and σ p3. Prinipal stresses may be related to stress omponents in an x y oordinate system using the free-body diagram shown in Figure.1. It is assumed that plane ABC is one of the three prinipal planes (i.e., n = 1,, or 3), and therefore no shear stress exists on this plane. The line-of-ation of prinipal stress σ pn defines one axis of the prinipal stress oordinate system. The diretion osines between this prinipal axis and the x-, y-, and -axes are pnx, pny, and pn, respetively. The surfae area of triangle ABC is denoted A ABC. The normal fore ating over triangle ABC therefore equals (σ pn )(A ABC ). The omponents of this normal fore ating in the x-, y-, and -diretions equal ( pnx )(σ pn )(A ABC ), ( pny )(σ pn )(A ABC ) and ( pn )(σ pn )(A ABC ), respetively. K13483_C00.indd 67 9/13/01 11:04:36 AM

28 68 Strutural Analysis of Polymeri Composite Materials +x σ pn τ y σ τ x C +y B τ xy τ yx D σ yy τ y τ x + A σ xx Figure.1 Free-body diagram used to relate strress omponents in the x y oordinate system to a prinipal stress. The area of the other triangular faes are given by Area of triangle ABD = ( pnx )(A ABC ) Area of triangle ACD = ( pny )(A ABC ) Area of triangle BCD = ( pn )(A ABC ). Summing fores in the x-diretion and equating to ero, we obtain σ A σ A τ A τ A = 0 pnx pn ABC xx pnx ABC xy pny ABC x pn ABC whih an be redued and simplified to ( σ pn σxx) pnx τxypny τxpn = 0 (.16a) Similarly, summing fores in the y- and -diretions results in τxypnx + ( σ pn σ yy) pny τypn = 0 (.16b) τxpnx τypny + ( σ pn σ) pn = 0 (.16) Equation.16 represent three linear homogeneous equations whih must be satisfied simultaneously. As diretion osines pnx, pny, and pn must also satisfy Equation.8, and therefore annot all equal ero, the solution an be obtained by requiring that the determinant of the oeffiients of pnx, pny, and pn equal ero: ( σ σ ) τ τ pn xx xy x τ ( σ σ ) τ xy pn yy y τ τ ( σ σ ) x y pn = 0 K13483_C00.indd 68 9/13/01 11:04:38 AM

29 Review of Fore, Stress, and Strain Tensors 69 Equating the determinant to ero results in the following ubi equation: 3 σ Θσ + Φσ Ψ = 0 (.17) pn pn pn where Θ, Φ, and Ψ are the first, seond, and third stress invariants, respetively, and have been previously listed as Equations.14 and.15. The three roots of this ubi equation represent the three prinipal stresses and may be found by appliation of the standard approah [3]. By onvention, the prinipal stresses are numbered suh that σ p1 is the algebraially greatest prinipal stress, whereas σ p3 is the algebraially least. That is, σ p1 > σ p > σ p3. One the prinipal stresses are determined the three sets of diretion osines (whih define the prinipal oordinate diretions) are found by substituting the three prinipal stresses given by Equation.17 into Equation.16 in turn. As only two of Equation.16 are independent, Equation.8 is used as a third independent equation involving the three unknown onstants, pnx, pny }, and pn. The proess of finding prinipal stresses and diretion osines will be demonstrated in the following Example Problem. Example Problem.4 Given: A state of stress referened to an x y oordinate is known to be: σ σ σ σ σ σ σ σ σ xx xy x yx yy y x y = ( ksi) Problem Find (a) the prinipal stresses and (b) the diretion osines that define the prinipal stress oordinate system. Solution This is the same stress tensor onsidered in Example Problem.3. As a part of that problem the first, seond, and third stress invariants were found to be Θ = 70 ksi Φ = 350 ( ksi) Ψ = ( ksi) 3 K13483_C00.indd 69 9/13/01 11:04:40 AM

30 70 Strutural Analysis of Polymeri Composite Materials a. Determining the Prinipal Stresses: In aordane with Equation.17, the three prinipal stresses are the roots of the following ubi equation: 3 σ 70σ 350σ = 0 The three roots of this equation represent the three prinipal stresses, and are given by σ = ksi, σ = ksi, and σ = 7. 7 ksi p1 p p3 b. Determining the Diretion Cosines: The first two of Equations.16 and.8 are used to form three independent equations in three unknowns. We have ( σ σ ) τ τ = 0 pn xx pnx xy pny x pn τ + ( σ σ ) τ = 0 xy pnx pn yy pny y pn ( ) + ( ) + ( ) = 1 pnx pny pn Diretion osines for σ p1 : The three independent equations beome: ( ) = 0 p1x p1y p ( ) 30 = 0 p1x p1y p1 ( ) + ( ) + ( ) = 1 p1x p1y p1 Solving simultaneously, we obtain: p1x = p1y = p1 = Diretion osines for σ p : The three independent equations beome: ( ) = 0 px py p 10 + ( ) 30 = 0 px py p ( ) + ( ) + ( ) = 1 px py p K13483_C00.indd 70 9/13/01 11:04:43 AM

31 Review of Fore, Stress, and Strain Tensors 71 Solving simultaneously, we obtain px = py = p = Diretion osines for σ p3 : The three independent equations beome: ( ) = 0 p3x p3y p ( ) 30 = 0 p3x p3y p3 ( ) + ( ) + ( ) = 1 p3x p3y p3 Solving simultaneously, we obtain: p3x = p3y = p3 = Plane Stress A stress tensor is always defined by six omponents of stress: three normal stress omponents and three shear stress omponents. However, in pratie a state of stress is often enountered in whih the magnitudes of three stress omponents in one oordinate diretion are known to be ero a priori. For example, suppose σ = τ x = τ y = 0, as shown in Figure.13a. As the three remaining nonero stress omponents (σ xx, σ yy, and τ xy ), all lie within the x y plane, suh a ondition is alled a state of plane stress. Plane stress onditions our most often beause of the geometry of the struture of interest. Speifially, the plane stress ondition usually exists in relatively thin, platelike strutures. Examples inlude the web of an I-beam, the body panel of an automobile, or the skin of an airplane fuselage. In these instanes the stresses indued normal to the plane of the struture are very small ompared with those indued within the plane of the struture. Hene, the small out-of-plane stresses are assumed to be ero, and attention is foused on the relatively high stress omponents ating within the plane of the struture. K13483_C00.indd 71 9/13/01 11:04:46 AM

32 7 Strutural Analysis of Polymeri Composite Materials (a) σ xx +x +y τ xy τ yx σ yy + (b) +x σ xx τ xy σ yy σ yy +y τ xy σ xx Figure.13 Stress elements subjeted to a state of plane stress. (a) 3-D stress element subjeted to a plane stress state (all stress omponents shown in a positive sense); (b) plane stress element drawn as a square rather than a ube (positive -axis out of the plane of the figure; all stress omponents shown in a positive sense). As laminated omposites are often used in the form of thin plates or shells, the plane stress assumption is widely appliable in omposite strutures and will be used throughout most of the analyses disussed in this book. As the out-of-plane stresses are negligibly small, for onveniene an infinitesimal stress element subjeted to plane stress will usually be drawn as a square rather than a ube, as shown in Figure.13b. Results disussed in earlier setions for general 3D state of stress will now be speialied for the plane stress ondition. It will be assumed that the nonero stresses lie in the x y plane (i.e., σ = τ x = τ y = 0). This allows the remaining omponents of stress to be written in the form of a olumn array, rather than a 3 3 array: σxx τxy 0 σ τxy σ yy 0 σ τ xx yy xy K13483_C00.indd 7 9/13/01 11:04:46 AM

33 Review of Fore, Stress, and Strain Tensors 73 Note that when a plane stress state is desribed, stress appears to be a first-order tensor, as (apparently) only three omponents of stress (σ xx, σ yy, and τ xy ) need be speified to desribe the state of stress. This is, of ourse, not the ase. Stress is a seond-order tensor in all instanes, and six omponents of stress must always be speified to define a state of stress. When we invoke the plane stress assumption, we have simply assumed a priori that the magnitude of three stress omponents (σ, τ x, and τ y ) are ero. Reall that either Equation.1a or Equation.1b governs the transformation of a stress tensor from one oordinate system to another. Equation.1b is repeated here for onveniene: σ σ σ σ σ σ σ σ σ x x x y x y x y y y x y = x x x y x y x y y y x y σxx σxy σx σ yx σ yy σ y σx σy σ x x y x x x y y y y x y (.1b) (repeated) When transformation of a plane stress tensor is onsidered, it will be assumed that the x y oordinate system is generated from the x y system by a rotation θ about the -axis. That is, the - and -axes are oinident, as shown in Figure.14. In this ase, the diretion osines are x x x y x y x y y y x y = os( θ) = os( 90 θ) = sin( θ) = os( 90 ) = 0 = os( 90 + θ) = sin( θ) = os( θ) = os( 90 ) = 0 = os( 90 ) = 0 = os( 90 ) = 0 = os( 0 ) = 1 If we now (a) substitute these diretion osines into Equation.1b, (b) label the shear stresses using the symbol using τ rather σ, K13483_C00.indd 73 9/13/01 11:04:47 AM

34 74 Strutural Analysis of Polymeri Composite Materials (a) +x σ xx τ xy σ yy σ yy +y τ xy σ xx (b) +x +θ +x σ x x σ y y τ x y +y τ x y +θ σ y y σ x x +y Figure.14 Transformation of a plane stress element from one oordinate system to another. (a) Plane stress element referened to the x y oordinate system; (b) plane stress element referened to the x y oordinate system, oriented θ-degrees ounter-lokwise from the x y oordinate system. and () note that σ = τ x = τ y = 0 by assumption, then Equation.1b beomes: σ τ τ τ σ τ τ τ σ x x x y x y x y y y x y osθ sinθ 0 σxx τxy 0 = sinθ osθ 0 τyx σ yy osθ sinθ 0 sinθ osθ K13483_C00.indd 74 9/13/01 11:04:48 AM

35 Review of Fore, Stress, and Strain Tensors 75 Completing the matrix multipliation indiated results in: σ τ τ τ σ τ τ τ σ x x x y x y x y y y x y os θσxx + sin θσ yy + osθsinθτxy osθsinθσ xx + osθsin θσ yy + (os θ sin θ) τxy 0 = osθsinθσ xx + osθsin θσ yy + (os θ sin θ) τxy sin θσxx + os θσ yy osθsinθτ xy (.18) As would be expeted, the out-of-plane stresses are ero: σ = τ x =τ y =0. The remaining stress omponents are σ = os ( θ) σ + sin ( θ) σ + os( θ)sin( θ) τ x x xx yy xy σ y y = sin ( θ) σ + os ( θ) σ os( θ)sin( θ) τ xx yy xy (.19) τ = os( θ)sin( θ) σ + os( θ)sin( θ) σ + [os ( θ) sin ( θ)] τ x y xx yy xy Equation.19 an be written using matrix notation as σ σ τ x x y y x y os ( θ) sin ( θ) os( θ)sin( θ) σ = sin ( θ) os ( θ) os( θ)sin( θ) σ os( θ)sin( θ) os( θ)sin( θ) os ( θ) sin ( θ) τ xx yy xy (.0) It should be kept in mind that these results are valid only for a state of plane stress. More preisely, Equations.19 and.0 represent stress transformations within the x y plane, and are only valid if the -axis is a prinipal stress axis. The 3 3 array that appears in Equation.0 is alled the transformation matrix, and is abbreviated as [T]: os ( θ) sin ( θ) os( θ)sin( θ) [ T ] = sin ( θ) os ( θ) os( θ)sin( θ) os( θ)sin( θ) os( θ)sin( θ) os ( θ) sin ( θ) (.1) The stress invariants (given by Equation.14 or.15) are onsiderably simplified in the ase of plane stress. Sine by definition σ = τ x = τ y = 0, the stress invariants beome: K13483_C00.indd 75 9/13/01 11:04:49 AM

Maximum Entropy and Exponential Families

Maximum Entropy and Exponential Families Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It

More information

4 Puck s action plane fracture criteria

4 Puck s action plane fracture criteria 4 Puk s ation plane frature riteria 4. Fiber frature riteria Fiber frature is primarily aused by a stressing σ whih ats parallel to the fibers. For (σ, σ, τ )-ombinations the use of a simple maximum stress

More information

LECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes.

LECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes. V. DEMENKO MECHNCS OF MTERLS 05 LECTURE Geometrial Properties of Rod Cross Setions (Part ) Moments of nertia Transformation with Parallel Transfer of xes. Parallel-xes Theorems S Given: a b = S = 0. z

More information

Beams on Elastic Foundation

Beams on Elastic Foundation Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating

More information

BEAMS: SHEARING STRESS

BEAMS: SHEARING STRESS LECTURE Third Edition BEAMS: SHEARNG STRESS A. J. Clark Shool of Engineering Department of Civil and Environmental Engineering 14 Chapter 6.1 6.4 b Dr. brahim A. Assakkaf SPRNG 200 ENES 220 Mehanis of

More information

The Hanging Chain. John McCuan. January 19, 2006

The Hanging Chain. John McCuan. January 19, 2006 The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

Modeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers

Modeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation

More information

The casing is subjected to the following:

The casing is subjected to the following: 16.50 Leture 13 Subjet: Roket asing design; Strutural modeling Thus far all our modeling has dealt with the fluid mehanis and thermodynamis of rokets. This is appropriate beause it is these features that

More information

Torsion. Torsion is a moment that twists/deforms a member about its longitudinal axis

Torsion. Torsion is a moment that twists/deforms a member about its longitudinal axis Mehanis of Solids I Torsion Torsional loads on Cirular Shafts Torsion is a moment that twists/deforms a member about its longitudinal axis 1 Shearing Stresses due to Torque o Net of the internal shearing

More information

Advanced Computational Fluid Dynamics AA215A Lecture 4

Advanced Computational Fluid Dynamics AA215A Lecture 4 Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas

More information

). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become

). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become Relativity and quantum mehanis: Jorgensen 1 revisited 1. Introdution Bernhard Rothenstein, Politehnia University of Timisoara, Physis Department, Timisoara, Romania. brothenstein@gmail.om Abstrat. We first

More information

Chapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3

Chapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3 hapter 3 eture 7 Drag polar Topis 3..3 Summary of lift oeffiient, drag oeffiient, pithing moment oeffiient, entre of pressure and aerodynami entre of an airfoil 3..4 Examples of pressure oeffiient distributions

More information

MOLECULAR ORBITAL THEORY- PART I

MOLECULAR ORBITAL THEORY- PART I 5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.

More information

Slenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method

Slenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory

More information

Aharonov-Bohm effect. Dan Solomon.

Aharonov-Bohm effect. Dan Solomon. Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that

More information

Slenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method

Slenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') 22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),

More information

Heat exchangers: Heat exchanger types:

Heat exchangers: Heat exchanger types: Heat exhangers: he proess of heat exhange between two fluids that are at different temperatures and separated by a solid wall ours in many engineering appliations. he devie used to implement this exhange

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

Relativity in Classical Physics

Relativity in Classical Physics Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of

More information

TORSION By Prof. Ahmed Amer

TORSION By Prof. Ahmed Amer ORSION By Prof. Ahmed Amer orque wisting moments or torques are fores ating through distane so as to promote rotation. Example Using a wrenh to tighten a nut in a bolt. If the bolt, wrenh and fore are

More information

The Lorenz Transform

The Lorenz Transform The Lorenz Transform Flameno Chuk Keyser Part I The Einstein/Bergmann deriation of the Lorentz Transform I follow the deriation of the Lorentz Transform, following Peter S Bergmann in Introdution to the

More information

Green s function for the wave equation

Green s function for the wave equation Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0

More information

A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM.

A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. S. Kanagaraj Eulidean Relativity s.kana.raj@gmail.om (1 August 009) Abstrat By re-interpreting the speial relativity (SR) postulates based on Eulidean

More information

Velocity Addition in Space/Time David Barwacz 4/23/

Velocity Addition in Space/Time David Barwacz 4/23/ Veloity Addition in Spae/Time 003 David arwaz 4/3/003 daveb@triton.net http://members.triton.net/daveb Abstrat Using the spae/time geometry developed in the previous paper ( Non-orthogonal Spae- Time geometry,

More information

Fig Review of Granta-gravel

Fig Review of Granta-gravel 0 Conlusion 0. Sope We have introdued the new ritial state onept among older onepts of lassial soil mehanis, but it would be wrong to leave any impression at the end of this book that the new onept merely

More information

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia 3 r d International Conferene on New Developments in Soil Mehanis and Geotehnial Engineering, 8-30 June 01, Near East University, Niosia, North Cyprus Values of of partial fators for for EC EC 7 7 slope

More information

ON DYNAMICALLY EQUIVALENT FORCE SYSTEMS AND THEIR APPLICATION TO THE BALANCING OF A BROOM OR THE STABILITY OF A SHOE BOX

ON DYNAMICALLY EQUIVALENT FORCE SYSTEMS AND THEIR APPLICATION TO THE BALANCING OF A BROOM OR THE STABILITY OF A SHOE BOX Proeedings of DEC 04 ASME 004 Design Engineering ehnial Conferenes and Computers and Information in Engineering Conferene September 8-Otober, 004, Salt Lake City, Utah, USA DE C0 04-5 7 188 ON DYNAMICALLY

More information

SURFACE WAVES OF NON-RAYLEIGH TYPE

SURFACE WAVES OF NON-RAYLEIGH TYPE SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh

More information

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential

More information

1 sin 2 r = 1 n 2 sin 2 i

1 sin 2 r = 1 n 2 sin 2 i Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with

More information

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California

More information

Chapter 2 Linear Elastic Fracture Mechanics

Chapter 2 Linear Elastic Fracture Mechanics Chapter 2 Linear Elasti Frature Mehanis 2.1 Introdution Beginning with the fabriation of stone-age axes, instint and experiene about the strength of various materials (as well as appearane, ost, availability

More information

Determination of the reaction order

Determination of the reaction order 5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order

More information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

Uniaxial Concrete Material Behavior

Uniaxial Concrete Material Behavior COMPUTERS AND STRUCTURES, INC., JULY 215 TECHNICAL NOTE MODIFIED DARWIN-PECKNOLD 2-D REINFORCED CONCRETE MATERIAL MODEL Overview This tehnial note desribes the Modified Darwin-Peknold reinfored onrete

More information

Physics 523, General Relativity Homework 4 Due Wednesday, 25 th October 2006

Physics 523, General Relativity Homework 4 Due Wednesday, 25 th October 2006 Physis 523, General Relativity Homework 4 Due Wednesday, 25 th Otober 2006 Jaob Lewis Bourjaily Problem Reall that the worldline of a ontinuously aelerated observer in flat spae relative to some inertial

More information

TENSOR FORM OF SPECIAL RELATIVITY

TENSOR FORM OF SPECIAL RELATIVITY TENSOR FORM OF SPECIAL RELATIVITY We begin by realling that the fundamental priniple of Speial Relativity is that all physial laws must look the same to all inertial observers. This is easiest done by

More information

THE REFRACTION OF LIGHT IN STATIONARY AND MOVING REFRACTIVE MEDIA

THE REFRACTION OF LIGHT IN STATIONARY AND MOVING REFRACTIVE MEDIA HDRONIC JOURNL 24, 113-129 (2001) THE REFRCTION OF LIGHT IN STTIONRY ND MOVING REFRCTIVE MEDI C. K. Thornhill 39 Crofton Road Orpington, Kent, BR6 8E United Kingdom Reeived Deember 10, 2000 Revised: Marh

More information

Modes are solutions, of Maxwell s equation applied to a specific device.

Modes are solutions, of Maxwell s equation applied to a specific device. Mirowave Integrated Ciruits Prof. Jayanta Mukherjee Department of Eletrial Engineering Indian Institute of Tehnology, Bombay Mod 01, Le 06 Mirowave omponents Welome to another module of this NPTEL mok

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

Complexity of Regularization RBF Networks

Complexity of Regularization RBF Networks Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw

More information

The homopolar generator: an analytical example

The homopolar generator: an analytical example The homopolar generator: an analytial example Hendrik van Hees August 7, 214 1 Introdution It is surprising that the homopolar generator, invented in one of Faraday s ingenious experiments in 1831, still

More information

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u Strauss PDEs e: Setion 3.4 - Exerise 3 Page 1 of 13 Exerise 3 Solve u tt = u xx + os x, u(x, ) = sin x, u t (x, ) = 1 + x. Solution Solution by Operator Fatorization Bring u xx to the other side. Write

More information

Simplified Buckling Analysis of Skeletal Structures

Simplified Buckling Analysis of Skeletal Structures Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation

More information

An Effective Photon Momentum in a Dielectric Medium: A Relativistic Approach. Abstract

An Effective Photon Momentum in a Dielectric Medium: A Relativistic Approach. Abstract An Effetive Photon Momentum in a Dieletri Medium: A Relativisti Approah Bradley W. Carroll, Farhang Amiri, and J. Ronald Galli Department of Physis, Weber State University, Ogden, UT 84408 Dated: August

More information

After the completion of this section the student should recall

After the completion of this section the student should recall Chapter I MTH FUNDMENTLS I. Sets, Numbers, Coordinates, Funtions ugust 30, 08 3 I. SETS, NUMERS, COORDINTES, FUNCTIONS Objetives: fter the ompletion of this setion the student should reall - the definition

More information

WRAP-AROUND GUSSET PLATES

WRAP-AROUND GUSSET PLATES WRAP-AROUND GUSSET PLATES Where a horizontal brae is loated at a beam-to-olumn intersetion, the gusset plate must be ut out around the olumn as shown in Figure. These are alled wrap-around gusset plates.

More information

The Unified Geometrical Theory of Fields and Particles

The Unified Geometrical Theory of Fields and Particles Applied Mathematis, 014, 5, 347-351 Published Online February 014 (http://www.sirp.org/journal/am) http://dx.doi.org/10.436/am.014.53036 The Unified Geometrial Theory of Fields and Partiles Amagh Nduka

More information

Q2. [40 points] Bishop-Hill Model: Calculation of Taylor Factors for Multiple Slip

Q2. [40 points] Bishop-Hill Model: Calculation of Taylor Factors for Multiple Slip 27-750, A.D. Rollett Due: 20 th Ot., 2011. Homework 5, Volume Frations, Single and Multiple Slip Crystal Plastiity Note the 2 extra redit questions (at the end). Q1. [40 points] Single Slip: Calulating

More information

n n=1 (air) n 1 sin 2 r =

n n=1 (air) n 1 sin 2 r = Physis 55 Fall 7 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.4, 7.6, 7.8 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with index

More information

Control Theory association of mathematics and engineering

Control Theory association of mathematics and engineering Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology

More information

BUCKLING OF LONG COMPRESSION-LOADED ANISOTROPIC PLATES RESTRAINED AGAINST INPLANE LATERAL AND SHEAR DEFORMATIONS

BUCKLING OF LONG COMPRESSION-LOADED ANISOTROPIC PLATES RESTRAINED AGAINST INPLANE LATERAL AND SHEAR DEFORMATIONS AIAA-3-789 BUCKLING OF LONG COMPRESSION-LOADED ANISOTROPIC PLATES RESTRAINED AGAINST INPLANE LATERAL AND SHEAR DEFORMATIONS Mihael P. Nemeth * Mehanis and Durability Branh, NASA Langley Researh Center

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS

More information

max min z i i=1 x j k s.t. j=1 x j j:i T j

max min z i i=1 x j k s.t. j=1 x j j:i T j AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be

More information

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J.

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J. A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE W. J. Haering* Senior Projet Engineer General Motors Corporation Warren, Mihigan R.

More information

Lecture 15 (Nov. 1, 2017)

Lecture 15 (Nov. 1, 2017) Leture 5 8.3 Quantum Theor I, Fall 07 74 Leture 5 (Nov., 07 5. Charged Partile in a Uniform Magneti Field Last time, we disussed the quantum mehanis of a harged partile moving in a uniform magneti field

More information

arxiv:gr-qc/ v2 6 Feb 2004

arxiv:gr-qc/ v2 6 Feb 2004 Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this

More information

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading 2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.1 Tuning controllers, in the large view Quoting from Section 16.7: Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output

More information

Practice Exam 2 Solutions

Practice Exam 2 Solutions MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department o Physis Physis 801T Fall Term 004 Problem 1: stati equilibrium Pratie Exam Solutions You are able to hold out your arm in an outstrethed horizontal position

More information

On the Quantum Theory of Radiation.

On the Quantum Theory of Radiation. Physikalishe Zeitshrift, Band 18, Seite 121-128 1917) On the Quantum Theory of Radiation. Albert Einstein The formal similarity between the hromati distribution urve for thermal radiation and the Maxwell

More information

Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames

Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames. Virtual Work for Frames IL 32 /9 ppling the virtual work equations to a frame struture is as simple as separating the frame into a series of beams and summing the virtual work for eah setion. In addition, when evaluating the

More information

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW P. М. Меdnis Novosibirs State Pedagogial University, Chair of the General and Theoretial Physis, Russia, 636, Novosibirs,Viljujsy, 8 e-mail: pmednis@inbox.ru

More information

Spinning Charged Bodies and the Linearized Kerr Metric. Abstract

Spinning Charged Bodies and the Linearized Kerr Metric. Abstract Spinning Charged Bodies and the Linearized Kerr Metri J. Franklin Department of Physis, Reed College, Portland, OR 97202, USA. Abstrat The physis of the Kerr metri of general relativity (GR) an be understood

More information

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013 Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it

More information

Chapter 2: One-dimensional Steady State Conduction

Chapter 2: One-dimensional Steady State Conduction 1 Chapter : One-imensional Steay State Conution.1 Eamples of One-imensional Conution Eample.1: Plate with Energy Generation an Variable Conutivity Sine k is variable it must remain insie the ifferentiation

More information

arxiv: v2 [math.pr] 9 Dec 2016

arxiv: v2 [math.pr] 9 Dec 2016 Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First

More information

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures Downloaded from orbit.dtu.dk on: Jan 19, 2019 Development of a user element in ABAQUS for modelling of ohesive laws in omposite strutures Feih, Stefanie Publiation date: 2006 Doument Version Publisher's

More information

A Characterization of Wavelet Convergence in Sobolev Spaces

A Characterization of Wavelet Convergence in Sobolev Spaces A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat

More information

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum

More information

STRUCTURAL BEHAVIOR OF R/C DEEP BEAM WITH HEADED LONGITUDINAL REINFORCEMENTS

STRUCTURAL BEHAVIOR OF R/C DEEP BEAM WITH HEADED LONGITUDINAL REINFORCEMENTS 13 th World Conferene on Earthquake Engineering anouver, B.C., Canada August 1-6, 24 Paper No. 58 STRUCTURAL BEHAIOR OF R/C DEEP BEAM WITH HEADED LONGITUDINAL REINFORCEMENTS Soo-Yeon SEO 1, Seung-Joe YOON

More information

Part G-4: Sample Exams

Part G-4: Sample Exams Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink

More information

The coefficients a and b are expressed in terms of three other parameters. b = exp

The coefficients a and b are expressed in terms of three other parameters. b = exp T73S04 Session 34: elaxation & Elasti Follow-Up Last Update: 5/4/2015 elates to Knowledge & Skills items 1.22, 1.28, 1.29, 1.30, 1.31 Evaluation of relaxation: integration of forward reep and limitations

More information

Collinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b

Collinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti

More information

Relativistic Dynamics

Relativistic Dynamics Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable

More information

Lecture 3 - Lorentz Transformations

Lecture 3 - Lorentz Transformations Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the

More information

Berry s phase for coherent states of Landau levels

Berry s phase for coherent states of Landau levels Berry s phase for oherent states of Landau levels Wen-Long Yang 1 and Jing-Ling Chen 1, 1 Theoretial Physis Division, Chern Institute of Mathematis, Nankai University, Tianjin 300071, P.R.China Adiabati

More information

A simple expression for radial distribution functions of pure fluids and mixtures

A simple expression for radial distribution functions of pure fluids and mixtures A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.

More information

Likelihood-confidence intervals for quantiles in Extreme Value Distributions

Likelihood-confidence intervals for quantiles in Extreme Value Distributions Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio

More information

Analysis of discretization in the direct simulation Monte Carlo

Analysis of discretization in the direct simulation Monte Carlo PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,

More information

Properties of Quarks

Properties of Quarks PHY04 Partile Physis 9 Dr C N Booth Properties of Quarks In the earlier part of this ourse, we have disussed three families of leptons but prinipally onentrated on one doublet of quarks, the u and d. We

More information

City, University of London Institutional Repository

City, University of London Institutional Repository City Researh Online City, University of London Institutional Repository Citation: Labib, M., Moslehy, Y. & Ayoub, A. (07). Softening Coeffiient of Reinfored Conrete Elements Subjeted to Three-Dimensional

More information

Frequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach

Frequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach Frequeny Domain Analysis of Conrete Gravity Dam-Reservoir Systems by Wavenumber Approah V. Lotfi & A. Samii Department of Civil and Environmental Engineering, Amirkabir University of Tehnology, Tehran,

More information

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution. arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat

More information

Illustrating the relativity of simultaneity Bernhard Rothenstein 1), Stefan Popescu 2) and George J. Spix 3)

Illustrating the relativity of simultaneity Bernhard Rothenstein 1), Stefan Popescu 2) and George J. Spix 3) Illustrating the relativity of simultaneity ernhard Rothenstein 1), Stefan Popesu ) and George J. Spix 3) 1) Politehnia University of Timisoara, Physis Department, Timisoara, Romania, bernhard_rothenstein@yahoo.om

More information

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines General Closed-form Analytial Expressions of Air-gap Indutanes for Surfaemounted Permanent Magnet and Indution Mahines Ronghai Qu, Member, IEEE Eletroni & Photoni Systems Tehnologies General Eletri Company

More information

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. All systems with interation of some type have normal modes. One may desribe them as solutions in absene of soures; they are exitations of the system

More information

Bäcklund Transformations: Some Old and New Perspectives

Bäcklund Transformations: Some Old and New Perspectives Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **

More information

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG

More information

Case Study in Reinforced Concrete adapted from Simplified Design of Concrete Structures, James Ambrose, 7 th ed.

Case Study in Reinforced Concrete adapted from Simplified Design of Concrete Structures, James Ambrose, 7 th ed. ARCH 631 Note Set 11 F015abn Case Study in Reinfored Conrete adapted from Simplified Design of Conrete Strutures, James Ambrose, 7 th ed. Building desription The building is a three-story offie building

More information

Chapter 9. The excitation process

Chapter 9. The excitation process Chapter 9 The exitation proess qualitative explanation of the formation of negative ion states Ne and He in He-Ne ollisions an be given by using a state orrelation diagram. state orrelation diagram is

More information

On Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations

On Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations On Certain Singular Integral Equations Arising in the Analysis of Wellbore Reharge in Anisotropi Formations C. Atkinson a, E. Sarris b, E. Gravanis b, P. Papanastasiou a Department of Mathematis, Imperial

More information

Geometry of Transformations of Random Variables

Geometry of Transformations of Random Variables Geometry of Transformations of Random Variables Univariate distributions We are interested in the problem of finding the distribution of Y = h(x) when the transformation h is one-to-one so that there is

More information

General Equilibrium. What happens to cause a reaction to come to equilibrium?

General Equilibrium. What happens to cause a reaction to come to equilibrium? General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember

More information

3 Tidal systems modelling: ASMITA model

3 Tidal systems modelling: ASMITA model 3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level

More information