16.21 Techniques of Structural Analysis and sig Spring 2003 Unit #2 - Stress and momentum balance Stress at a point We are going to consider the forces exerted on a material. These can be external or internal. External forces come in two flavors: body forces(given per unit mass or volume) and surface forces(given per unit area). If we cut a body of material in equilibrium under a set of external forces along a plane as shown in fig. 1. and consider one side of it we draw two conclusions: 1 the equilibrium provided by the loads from the side taken out is provided by a set of forces that are distributed among the material particles adjacent to the cut plane and that should provide an equivalent set of forces to the ones loading the part taken out, 2)these forces can now be considered as external surface forces acting on the part of material under consideration The stress vector at a point on As is defined as s-0△S If the cut had gone through the same point under consideration but along a plane with a different normal, the stress vector t would have been different Let's consider the three stress vectors t () acting on the planes normal to the coordinate axes. Let's also decompose each to) in its three components in the coordinate system e:(this can be done for any vector)as( see Fig. 2) )=可ej
16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #2 - Stress and Momentum balance Stress at a point We are going to consider the forces exerted on a material. These can be external or internal. External forces come in two flavors: body forces (given per unit mass or volume) and surface forces (given per unit area). If we cut a body of material in equilibrium under a set of external forces along a plane, as shown in Fig.1, and consider one side of it, we draw two conclusions: 1) the equilibrium provided by the loads from the side taken out is provided by a set of forces that are distributed among the material particles adjacent to the cut plane and that should provide an equivalent set of forces to the ones loading the part taken out, 2) these forces can now be considered as external surface forces acting on the part of material under consideration. The stress vector at a point on ΔS is defined as: t = lim ΔS→0 f (1) ΔS If the cut had gone through the same point under consideration but along a plane with a different normal, the stress vector t would have been different. Let’s consider the three stress vectors t(i) acting on the planes normal to the coordinate axes. Let’s also decompose each t(i) in its three components in the coordinate system ei (this can be done for any vector) as (see Fig.2): t(i) = σijej (2) 1
body forces Figure 1: Surface force f on area AS of the cross section by plane whose normal is n Jii is the component of the stress vector th along the e, direction x1 Figure 2: Stress components Stress tensor We could keep analyzing different planes passing through the point with different normals and therefore. different stress vectors t(n)and one might wonder if there is any relation among them or if they are all independent. The answer to this question is given by invoking equilibrium on the(shrinking)
n body forces ∆ s f surface forces body forces n Figure 1: Surface force f on area ΔS of the cross section by plane whose normal is n σij is the component of the stress vector t(i) along the ej direction. σ11 σ12 σ13 σ21 σ22 σ23 σ31 σ32 σ33 t t t (1) (2) (3) x1 x2 x3 e e e1 2 3 Figure 2: Stress components Stress tensor We could keep analyzing different planes passing through the point with different normals and, therefore, different stress vectors t(n) and one might wonder if there is any relation among them or if they are all independent. The answer to this question is given by invoking equilibrium on the (shrinking) 2
tetrahedron of material of Fig 3. The area of the faces of the tetrahedron are △S1,△S2,△S3and△S. The stress vectors on planes with reversed normals t(ei have been replaced with -t() using Newton's third law of action and reaction(which is in fact derivable from equilibrium): t(-n)=-t(n Enforcing equilibrium we have t)△S-t)△S1-t2)△S2-t(△s3=0 (3) n) Figure 3: Cauchy 's tetrahedron representing the equilibrium of a tetrahedron shrinking to a point where AV is the volume of the tetrahedron and f is the body force per unit volume. The following relation: ASni= ASi derived in the following mathematical aside y virtue of Green's Theorem 3
� � tetrahedron of material of Fig.3. The area of the faces of the tetrahedron are ΔS1, ΔS2, ΔS3 and ΔS. The stress vectors on planes with reversed normals t(−ei) have been replaced with −t(i) using Newton’s third law of action and reaction (which is in fact derivable from equilibrium): t(−n) = −t(n) . Enforcing equilibrium we have: t(n) ΔS − t(1)ΔS1 − t(2)ΔS2 − t(3)ΔS3 = 0 (3) t (3) − t (2) − t (1) − −e −e −e 1 2 3 n t(n) Figure 3: Cauchy’s tetrahedron representing the equilibrium of a tetrahedron shrinking to a point where ΔV is the volume of the tetrahedron and f is the body force per unit volume. The following relation: ΔSni = ΔSi derived in the following mathematical aside: By virtue of Green’s Theorem: �φdV = nφdS V S 3
applied to the function g= l, we get which applied to our tetrahedron gives 0=△Sn-△S1 If we take the scalar product of this equation with ei, we obtain AS(n:e)=△S S;=△S can then be replaced in equation 3 to obtain t2)+(n:e3)t③) The factor in parenthesis is the definition of the Cauchy stress tensor o elt/+e2t+est Note it is a tensorial expression (independent of the vector and tensor com- ponenents in a particular coordinate system). To obtain the tensorial com- ponenents in our rectangular system we replace the expressions of to from Eqn. 2 (6) Replacing in Eqn. 4: n0ie,e;=0i(n-;=(0igniej ti okit
� � � � � applied to the function φ = 1, we get 0 = S ndS which applied to our tetrahedron gives: 0 = ΔSn − ΔS1e1 − ΔS2e2 − ΔS3e3 If we take the scalar product of this equation with ei, we obtain: ΔS(n · ei) = ΔSi or ΔSi = ΔSni can then be replaced in equation 3 to obtain: ΔS � t(n) − (n · e1)t(1) + (n · e2)t(2) + (n · e3)t(3)� = 0 or � � t(n) = n · e1t(1) + e2t(2) + e3t(3) (4) The factor in parenthesis is the definition of the Cauchy stress tensor σ: σ = e1t(1) + e2t(2) + e3t(3) = eit(i) (5) Note it is a tensorial expression (independent of the vector and tensor componenents in a particular coordinate system). To obtain the tensorial componenents in our rectangular system we replace the expressions of t(i) from Eqn.2 σ = eiσijej (6) Replacing in Eqn.4: t(n) = n · σ (7) or: t(n) = n · σijeiej = σij n · ei ej = σijni ej (8) ti = σkink (9) 4
Transformation of stress components Consider a different system of cartesian coordinates e. We can express our tensor in either one We would like to relate the stress components in one the two systems. To this end, we take the scalar product of (10) with e and e e;oe=ok(e; ek(el e=omn(e em)(e OI Ok(e;ek)(ere The factors in parenthesis are the cosine directors of the angles between the original and primed coordinate axes Principal stresses and directions Given the components of the stress tensor in a given coordinate system, the determination of the maximum normal and shear stresses is critical for the design of structures. The normal and shear stress components on a plane with normal n are given by n -Kinin It is obvious from these equations that the normal component achieves its maximum tN=lt( n)ll when the shear components are zero. In this case =An= al or in components Kink= asking (oki-Xoki nk=0 5
� � � � � �� � � � � Transformation of stress components Consider a different system of cartesian coordinates e� i. We can express our tensor in either one: σ = σklekel = σ� mnemen (10) We would like to relate the stress components in one the two systems. To this end, we take the scalar product of (10) with e� i and e� j : � �� � � �� � mn ei · em en · ej = σ� ei ij � · σ · ej � = σkl ei � · ek el · ej � = σ� � � � � mnδimδnj = σ� or σ� ij = σkl ei · ek el · ej (11) The factors in parenthesis are the cosine directors of the angles between the original and primed coordinate axes. Principal stresses and directions Given the components of the stress tensor in a given coordinate system, the determination of the maximum normal and shear stresses is critical for the design of structures. The normal and shear stress components on a plane with normal n are given by: tN = t(n) · n = σkinkni tS = �t(n)�2 − t2 N It is obvious from these equations that the normal component achieves its maximum tN = �t(n) � when the shear components are zero. In this case: t(n) = n · σ = λn = λIn or in components: σkink = λni σkink = λδkink (12) σki − λδki nk = 0 5
which means that the principal stresses are obtained by solving the previous eigenvalue problem, the principal directions are the eigenvectors of the prob- lem. The eigenvalues A are obtained by noticing that the last identity can b satisfied for non-trivial n only if the factor is singular, i. e, if its determinant vanishes 01-入12013 22 23 which leads to the characteristic equation 3+12-I2入+I3=0 where 1=0i=011+022+a33 h2=5(10y-oy1)=0102+0203+0301-(v12+2+)(14) I 3= det[o]=loiill are called the stress invariants because they do not depend on the coordinate ystem of choic Linear and angular momentum balance We are going to derive the equations of momentum balance in integral form, since this is the formulation that is more aligned with our integral"approach in this course. We start from the definition of linear and angular momentum For an element of material at position x of volume dv, density p, mass pdv which remains constant, moving at a velocity v, the linear momentum is pvdv and the angular momentum xx(pvdv). The total momenta of the body are obtained by integration over the volume as pvdv and/x×pvd respectively. The principle of conservation of linear momentum states that the rate of change of linear momentum is equal to the sum of all the external forces acting on the body Dt/pvdv fdV+tds 16) 6
� � � � � � � which means that the principal stresses are obtained by solving the previous eigenvalue problem, the principal directions are the eigenvectors of the problem. The eigenvalues λ are obtained by noticing that the last identity can be satisfied for non-trivial n only if the factor is singular, i.e., if its determinant vanishes: � � �σ11 − λ σ12 σ13 � � � σ21 σ22 − λ σ23 � � = 0 � σ31 σ32 σ33 − λ� which leads to the characteristic equation: −λ3 + I1λ2 − I2λ + I3 = 0 where: I1 = σii = σ11 + σ22 + σ33 (13) 1� � � � I2 = σiiσjj − σijσji = σ11σ22 + σ22σ33 + σ33σ11 − σ2 23 + σ2 12 + σ2 31 (14) 2 I3 = det[σ] = �σij� (15) are called the stress invariants because they do not depend on the coordinate system of choice. Linear and angular momentum balance We are going to derive the equations of momentum balance in integral form, since this is the formulation that is more aligned with our “integral” approach in this course. We start from the definition of linear and angular momentum. For an element of material at position x of volume dV , density ρ, mass ρdV which remains constant, moving at a velocity v, the linear momentum is ρvdV and the angular momentum x × (ρvdV ). The total momenta of the body are obtained by integration over the volume as: ρvdV and x × ρvdV V V respectively. The principle of conservation of linear momentum states that the rate of change of linear momentum is equal to the sum of all the external forces acting on the body: D ρvdV = fdV + tdS (16) Dt V V S 6
ere B is the total derivative. The Ths can be expanded as D (pdv)v+ but di (pdv)=0 from conservation of mass, so the principle reads fdv+tdS Now, using what we've learned about the tractions and their relation to the stress tensor. Padv=/fdv+/n ads This is the linear momentum balance equation in integral form. We can replace the surface integral with a volume integral with the aid of the diver- gence theorem and then(18)becomes f-V·adv=0 Since this principle applies to an arbitrary volume of material, the integrand must vanish f = 0 at This is the linear momentum balance equation in differential form. In com- 可x+f=p at Angular momentum balance and the symmetry of the stress tensor The principle of conservation of angular momentum states that the rate of change of angular momentum is equal to the sum of the moment of all the xternal forces acting on the body Dt/ ex xvdv x×fdV+)x×tdS
� � � � � � � � � � � � where D is the total derivative. The lhs can be expanded as: Dt D D ∂v ρvdV = Dt (ρdV )v + ρ dV Dt V V V ∂t D but Dt(ρdV ) = 0 from conservation of mass, so the principle reads: ∂v ρ dV = fdV + tdS (17) V ∂t V S Now, using what we’ve learned about the tractions and their relation to the stress tensor: � V ρ ∂v ∂t dV = � V fdV + � S n · σdS (18) This is the linear momentum balance equation in integral form. We can replace the surface integral with a volume integral with the aid of the divergence theorem: � � n · σdS = � · σdV S V and then (18) becomes: ∂v ρ − f − � · σ dV = 0 V ∂t Since this principle applies to an arbitrary volume of material, the integrand must vanish: ∂v ρ − f − � · σ = 0 (19) ∂t This is the linear momentum balance equation in differential form. In components: ∂vi σji,j + fi = ρ ∂t Angular momentum balance and the symmetry of the stress tensor The principle of conservation of angular momentum states that the rate of change of angular momentum is equal to the sum of the moment of all the external forces acting on the body: D ρx × vdV = x × fdV + x × tdS (20) Dt V V S 7
It can be conveniently written as (4-s+/(b=21=-/(m-mm Using t;= Kink, the divergence theorem and(19), this expression leads to ee homework pro (oij-dji)dv=0 which applies to an arbitrary volume V if the integrand vanishes. This implies. > and therefore, can only be satisfied
� � � � � � � It can be conveniently written as xitj − xj ti dS + xifj − xjfi � dV = � xi ∂vj − xj ∂vi � dV S V V ∂t ∂t Using ti = σkink, the divergence theorem and (19), this expression leads to (see homework problem): (σij − σji)dV = 0 V which applies to an arbitrary volume V , and therefore, can only be satisfied if the integrand vanishes. This implies: σij = σji (21) 8