16.21 Techniques of Structural Analysis and Design pring 2003 Unit#8-Principle of Virtual Displacements Principle of Virtual Displacements Consider a body in equilibrium. We know that the stress field must satisfy the differential equations of equilibrium. Multiply the differential equations of equilibrium by an "arbitrary"displacement field T (0+f1)1=0 Note that the field ii is noT the actual displacement field u; corresponding to the solution of the problem but a virtual displacement field. Therefore equation(1)can be interpreted as the local expression of virtual work done by the actual stresses and the body forces on the virtual displacement ii and that it must be zero. The total virtual work done on the body is obtained by integration over the volume 0jij+fiiidv=0
� � � � � 16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #8 - Principle of Virtual Displacements Principle of Virtual Displacements Consider a body in equilibrium. We know that the stress field must satisfy the differential equations of equilibrium. Multiply the differential equations of equilibrium by an “arbitrary” displacement field u¯i: σji,j + fi u¯i = 0 (1) Note that the field u¯i is NOT the actual displacement field ui corresponding to the solution of the problem but a virtual displacement field. Therefore, equation (1) can be interpreted as the local expression of virtual work done by the actual stresses and the body forces on the virtual displacement u¯i and that it must be zero. The total virtual work done on the body is obtained by integration over the volume: σji,j + fi u¯idV = 0 (2) V 1
and it must also be zero since the integrand is zero every where in the domain oji,ju di fiuid=0 (a),-0nin,)d+/fa=0 o lin:ds dV+/fi dV=0 The integral over the surface can be decomposed into two: an integral over the portion of the boundary where the actual external surface loads(tractions Qre specified St and an integral over the portion of the boundary where the placements are specified(supports)Su. This assumes that these sets are disjoint and complementary, i.e S=SUS,S∩St= (6) Oije We will require that the virtual displacements i; vanish on Su, i.e., that the virtual displacement field satisfy the homogeneous essential boundary condi- a1(x)=0,ifx;∈S Then, the second integral vanishes. The resulting expression is a statement of the Principle of virtual Displacements(PVD) aijE;dv=ti ds+/fiii dl (9) It reads: The work done by the ecternal tractions and body forces on an ad missible(differentiable and satisfying the homogeneous boundary conditions but otherwise arbitrary) displacement field is equal to the work done by the equilibrated stresses(the actual solution of the problem)on the virtual strains ( the strains produced by the virtual field) Example: Consider the bar under a tensile load shown in the figure
¯ ¯ ¯ � � �� � � � � � � � � � and it must also be zero since the integrand is zero everywhere in the domain. σji,ju¯idV + fiu¯idV = 0 (3) � V V � σjiui ,j ¯ − σjiu¯i,j dV + fiu¯idV = 0 (4) V � V � � σjiu¯injdS − σij �ijdV + fiu¯idV = 0 (5) S V V The integral over the surface can be decomposed into two: an integral over the portion of the boundary where the actual external surface loads (tractions) are specified St and an integral over the portion of the boundary where the displacements are specified (supports) Su. This assumes that these sets are disjoint and complementary, i.e., S = Su ∪ St, Su ∩ St = ∅ (6) tiu¯idS + σjiu¯injdS − σij �ijdV + fiu¯idV = 0 (7) St Su V V We will require that the virtual displacements u¯i vanish on Su, i.e., that the virtual displacement field satisfy the homogeneous essential boundary conditions: u¯i(xj ) = 0, if xj ∈ Su (8) Then, the second integral vanishes. The resulting expression is a statement of the Principle of Virtual Displacements (PVD): σij �ijdV = tiu¯idS + fiu¯idV (9) V St V It reads: The work done by the external tractions and body forces on an admissible (differentiable and satisfying the homogeneous boundary conditions but otherwise arbitrary) displacement field is equal to the work done by the equilibrated stresses (the actual solution of the problem) on the virtual strains (the strains produced by the virtual field). Example: Consider the bar under a tensile load shown in the figure: 2
EA P L The PVD applied to this case is dv= pi a1=L P证 da ld dx EA,-元1 EA dindar,- pa The second term on the left hand side is zero because we have asked that i1=0 at the support. Note we have not asked for any condition on in at 1=L where the load is applied Ea li,dx The only way this expression can be satisfied for any admissible virtual dis- placement field in is if: P=EA EA 0 which represent the equilibrium conditions at the boundary and inside the bar, respectively P=AE E d 3
� � � L P E, A x1 The PVD applied to this case is: du¯1 � σ11 dV = Pu¯1� V dx1 x1=L � L � du1 du¯1 � A E dx1 = Pu¯1� 0 dx1 dx1 x1=L � L� d �du1 � d2u1 � � � EA ¯ ¯ u1� dx1 u1 − dx2 u1 dx1 = P ¯ 0 dx1 1 x1=L � du1 � � du1 � � L d2u1 � � EAdx1 u1 − EAdx1 u1 − EA dx2 ¯ ¯ u¯1dx1 = Pu¯1� x1=L x1=0 0 1 x1=L The second term on the left hand side is zero because we have asked that u¯1 = 0 at the support. Note we have not asked for any condition on u¯1 at x1 = L where the load is applied. � du1 � � � � � � L d2u1 EAdx1 � x1=L − P u1� = EA dx2 ¯ u¯1dx1 x1=L 0 1 The only way this expression can be satisfied for any admissible virtual displacement field u¯1 is if: du1 � P = EA � dx1 x1=L and d2u1 EA = 0 dx2 1 which represent the equilibrium conditions at the boundary and inside the bar, respectively: � du1 �� � � � P = A E � = Aσ11� dx1 x1=L x1=L and d � du1 � d EA = σ11 = 0 dx1 dx1 dx1 3
The solution of this problem is: the boundary conditions are P ea u EA ∈11 EA P Ee A Example: With the exact solution of the problem of the bar under a tensile load, verify the satisfaction of the PVD for the following virtual AE dc= Palo EA Pal= Pal g.e.d AE EA A EPAL Pal g.e.d
The solution of this problem is: u1(x1) = ax1 + b the boundary conditions are: u1(0) = 0 ⇒ b = 0 P = Ea A u1 = P EAx1 �11 = du1 dx1 = P EA σ11 = E�11 = P A Example: With the exact solution of the problem of the bar under a tensile load, verify the satisfaction of the PVD for the following virtual displacement fields: • u¯1 = ax1: � L 0 P adx1 = PaL(?) EA AE PaL = PaL q.e.d. ¯ 1 • u1 = ax : 2 � L 0 P 2ax1dx1 = PaL2 (?) EA AE L2 �A �EP �2a �2 = PaL q.e.d. �E �A Remarks: 4
Principle of Virtual Displacements enforces equilibrium (in weak form) enforces traction(natural) boundary conditions does NoT enforce displacement(essential) boundary conditions will be satisfied for all equilibrated solutions, compatible or in- compatible Unit dummy displacement method Another application of the PVD: provides a way to compute reactions(or dis- placements) in structures directly from PVD. Consider the concentrated reaction force at point 0 of a structure in equilibrium under a set of loads and supports. We can prescribe an arbitrary admissible displacement field ii and the PVD will hold. The unit dummy displacement method consists of choosing the virtual displacement field such that )=l in the direction of the reaction Ro we are interested in. Then the virtual work of the reaction is to Ro=Ro. The PVD then reads(in the absence of body forces Ro·0=/o;;dV (10) ij∈ij 11) where Ei are the virtual strains produced by the virtual displacement field uo Example: 5
¯ ¯ ¯ � � • Principle of Virtual Displacements: – enforces equilibrium (in weak form) – enforces traction (natural) boundary conditions – does NOT enforce displacement (essential) boundary conditions – will be satisfied for all equilibrated solutions, compatible or incompatible Unit dummy displacement method Another application of the PVD: provides a way to compute reactions (or displacements) in structures directly from PVD. Consider the concentrated reaction force at point ��0�� of a structure in equilibrium under a set of loads and supports. We can prescribe an arbitrary admissible displacement field u¯i and the PVD will hold. The unit dummy displacement method consists of choosing the virtual displacement field such that u¯i(x0) = 1 in the direction of the reaction R0 we are interested in. Then the virtual work of the reaction is u¯0 · R0 = |R0. The PVD then reads (in the absence of body forces): R0 · u¯0 = σij �ijdV (10) V R0 = V σij �ijdV (11) where �ij are the virtual strains produced by the virtual displacement field u¯0. Example: 5