运动微分方程弹性体体积V ,表面积S ,密度ρ,单位质量所受的体力为f,体力场为f(x,t),单位向量为n 的面元dS 的面力场为t(n,x,t),x 为原点到受力点的向量,t 为时间。
弹性体在t 时刻的动量P (t)dV v dt ddV f dS t dtdP F f V f m F dVf dS t F F F dVv m v p Vi Vi si ii Vi si i Vi i ⎰⎰⎰⎰⎰⎰=+=⨯=⨯=+=+===ρρρρρ动量定理合力弹性体动量体体面*******************************************************************************散度定理:散度定理是矢量场中体积分与面积分之间的一个转换。
⎰⎰∙=∙∇sVS d F dV F散度:表征矢量场A 产生的体积(三维)或面积(二维)的相对膨胀率,其表达式为▽·A 。
zRy Q x P R Q P z y x F ∂∂+∂∂+∂∂=∙∂∂∂∂∂∂=∙∇),,(),,( ,P,Q ,R 为F 在x,y,z 上的分量。
散度定理的证明:S d F dV F sV∙=∙∇⎰⎰⎰⎰⎰。
令()R Q P F ,,=,假设F =(0,0,R),则需要证明dS n R dV R sVz⎰⎰⎰⎰⎰∙=),0,0( 如下图,投影区为U。
dxdy y x z y x R y x zy x R dxdy dz R dV R Uy x Z y x Z zDz ))],(,,()),(,,([)(),(),(底顶顶底⎰⎰⎰⎰⎰⎰⎰⎰-==S=S 底+S 顶+S 侧面令S 底=S1,S 顶=S2,S 侧面=S3. 对于顶面,则dxdy yZ x Z dS n )1,,(22∂∂-∂∂-=Rdxdy dxdy y Z x Z R dS n R =∂∂-∂∂-=)1,,)(,0,0(),0,0(22dxdy y x z y x R dxdy R dS n R U⎰⎰⎰⎰⎰⎰==)),(,,(),0,0(顶顶顶对于底面,则dxdy yZ x Z dS n )1,,(11-∂∂∂∂=dxdy y x z y x R dxdy R dS n R U⎰⎰⎰⎰⎰⎰=-=)),(,,(-),0,0(底底底对侧面,S3=0。
F 垂直于侧面。
综上所述,即证。
其他复杂情况分解为这种情况,即证。
******************************************************************************* 利用散度定理dVx dS n dS t Vjji sj ji si ⎰⎰⎰∂∂==ττ---------------------------------------------------------------------------------------------------------------------- 证明j ji i n t τ=设在面X1X2X3上,面积为S ;面OX1X2,OX2X3,OX1X3的面积为S1、S2、S3。
Sh V S n n e S S i i ∆∙=∆∆=∙∆=∆31),cos(受力分析:x 1x 3x 2i i j ji i a f S S t ∙∆=∆∙+∆∙-∆∙V V ρρτi i j ji i a S h S h f S n S t ∙∆∙=∆∙∙+∆-∆∙3131ρρτ0→h ,j ji i n t τ=物理意义:如果已知过点P 与三个坐标轴方向相垂直的三个面元上的九个应力分量ji τ,则过该点任意面元(法线方向n)上的应力向量i t 都可用这九个应力分量按此式表示出来。
柯西应力公式:在一点处三个与坐标轴方向相垂直的面元上的九个应力分量可以确定该点的应力。
---------------------------------------------------------------------------------------------------------------------- 由于F=ma ,即有dV t u dV t v dV dt dv dV dt dv dV v dt dVi V i Vi Vi V i ⎰⎰⎰⎰⎰∂∂=∂∂==22ρρρρρ由i F F F =+体面dV t u dV f dV x ViV i Vj ji⎰⎰⎰∂∂=+∂∂22ρρτ 去掉积分号22tu f x i i j ji∂∂=+∂∂ρρτ若弹性体处于静力平衡状态,a=0。
弹性体的平衡微分方程:0=+∂∂i jji f x ρτ应力张量的对称性该部分弹性体在时刻t 对坐标系原点o 的动量矩⎰⨯=VdV v x t Nρ)(在i e方向上的分量为dV v x e N k j Vijk i ρ⎰=作用在弹性体上的体力与面力的力矩⎰⎰⨯+⨯=VsdV f x dS t x t Mρ)(在i e方向上的分量为⎰⎰+=Vk j ijk k sj ijk i dV f x e dS t x e M ρ由i iM dtdN =,即dV v x e dt ddV f x e dS t x e k j Vijk Vk j ijk k sj ijk ρρ⎰⎰⎰=+ 运用散度定理:dV x x e dV x x e dVx x e dS n x e dS t x ellk j Vjk ijk l lk jVlk jl ijk Vllk j ijkl lk sj ijk k sj ijk][][)(∂∂+=∂∂+=∂∂==⎰⎰⎰⎰⎰ττττδττ对于dV t ux v v e dV dt dv x v v e dV dt v x d e dV v x e dt dk j j k Vijk k jj k Vijk k j V ijk k j Vijk ][][)(22∂∂+=+==⎰⎰⎰⎰ρρρρ因为0=j k ijk v v e ,j k ijk v v e 、与垂直。
dV t u x e dV v x e dt dk j V ijk k j Vijk 22∂∂=⎰⎰ρρ 所以,有dV t u x e dV f x e dV x x e k j Vijk V k j ijk l lk j V jk ijk 22][∂∂=+∂∂+⎰⎰⎰ρρττ 0][][2222=∂∂-+∂∂+=∂∂-+∂∂+⎰⎰⎰⎰⎰Vk k l lk j ijk Vjk ijk V k j V ijk k j ijk l lk j V jk ijk dV tuf x x e dV e dV t u x e dV f x e dV x x e ρρττρρττ因为0][22=∂∂-+∂∂⎰V kk l lk j ijk dV t u f x x e ρρτ,所以0=⎰dV e Vjk ijk τ。
所以0=jk ijk e τ于是, i=1, 03223=-ττ i=2, 01331=-ττ i=3, 01221=-ττ 即jk kj ττ=,这就是剪应力互等定理。
应力边界条件 j ji i n t τ=四面体的的表面元的外法线为n,外来作用面力为t,则弹性体的应力边界条件:j ji i n t τ=。
它表明了应力的边界值与边界面上的表面力的关系。
本构方程(应力--应变关系)ij ij ij e μλθδτ2+=证明:当应力小于比例极限时,应力与应变是成正比的。
将上面理论推广:线性弹性体内一点处的应力张量分量为该点应变张量分量的线性齐次函数,反之亦然。
即kl ijkl ij e C =τ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡333231232221131211333333323331332333223321331333123311323332323231322332223221321332123211313331323131312331223121311331123111233323322331232323222321231323122311223322322231222322222221221322122211213321322131212321222121211321122111133313321331132313221321131313121311123312321231122312221221121312121211113311321131112311221121111311121111333231232221131211e e e e e e e e e C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C τττττττττ其中,ijkl C 称为弹性系数张量,它有81个分量,它们的值由所在点处材料的弹性性质决定。
如果弹性体不是均匀的,其弹性性质随点的不同而不同,即ijkl C 是点坐标i x 的函数。
如果弹性体是均匀的,弹性体内各点的弹性性质相同,即ijkl C 是与点坐标i x 无关的函数。