《弹性力学平面问题有限元法》由会员分享,可在线阅读,更多相关《弹性力学平面问题有限元法(73页珍藏版)》请在金锄头文库上搜索。
1、3弹性力学平面问题有限元法材料力学主要研究杆、梁、柱结构力学主要研究杆系(或梁系) 弹性力学主要研究实体和板得受力和变形 弹性力学假设所研究的物体:连续的、完全弹性的 均匀的、各向同性的、微小变形的和无初应力的在这假设基础上研究受力物体一点上的应力、应变、 变形和平衡关系。线性:(非线性)家构的应力占应变的关系(本构关系)呈线性变化。弹性:(塑性)结构在外力拆除后能够完全恢复原有形状的特性。静力分析:(动态分析) 结构所受外力是不随时间变化的恒力。一、弹性力学中的物理量载荷、应力、应变、位移1 载荷载荷是外界作用在弹性体上的力,又称为外力。它包括体力、面力和 中力二无中无呂壬弋体力鬼分話于整个
2、弹性体体积内的外力,如重力和惯性力。在弹性体 内任一点,单位体积的体力用人表示,它可分解为给定坐标系x、y和z 三个坐标轴上的投影称为体力分量。TP = PvX Pry PVZ面力是作用于弹性体表面上的外力,如流体压力和接触压力。TPs = PSX Psy Psz如果外力作用面很小,或者说外力作用在某一点上,则这种外力称为集中力。TPc = PcX Pcy Pcz无论那个位置的体力、那一边界面上的面力,均以正 标向为正,且斜面上的面力是以单位斜面面积上的作用 力数值来表示。内力A定义:物体本身不同部分之间相互作用的力。求解方法:截面法2.应力:内力集度。反映内力分布情况(应力场)AFAtO W
3、矢量P方向沿沿截面切向和法向分解为厂和b量纲:匸MT ?AF的极限方向应力的两种不同分解方法a)沿坐标轴分解b)沿截面法向和切向分解厂除了在推导公式过程中沿坐标轴分解外,通常 釆用沿截面法向和切向分解的方式,即分解为正应力b和切应力因为与物体形变和材料强度之间相关的是应力在其作用截面的法线方向和切线方向的分量。dcr、T-dy正六面单元体的取法*经过物体内任一点如P点取出一个微小的正六面 体,它的棱边分别平行于三个坐标轴而长度分别 为:PA = Ax,PB = Ay,PC = Az 将每个面上的应力分 解为一个正应力和两个切应力。正应力用CF表 示,切应力用厂表示。*应力下标的含意:A.作用面
4、的外法线方向B.力的指向切应力互等定理在受力物体相互垂直的希个平面上,切应力必然成对存在,且数值相等;两者都垂直于两平面的交线,方向共同指向或背离这一交线。弹力规定/oyT一-TTT乂,材力规定TTTTw3 形变A定义:形状的改变(长度的改变和角度的改变)厂线应变(正应变):线段单位长度的伸缩。记号: 正负:伸长为正,压缩为负厂切应变(剪应变):两方向线段夹角的改变。记号:7(以弧而非角度表示)正负:X矗一点的形变状态的概念几何规律:过空间一点有无数根直线。力学特点:即使过同一点,不同方向线段的伸长也同;任两根直线之间夹角的改变也不相同。同一点的应力状态情况一样,可证明,在物体内任意一点,若已
5、知&八8八&乙、丫小丫0即可求得经过该点的任意截面上(方向余弦已知)的正应变和切应变。故这六个应变分量完全确定了该点的应变状态。4 位移定义:位置的改变。记号:况、V、W正负:沿坐标轴正向为正,负向为负。分类:与形变有关的位移和与形变无关位移(刚体位移)八弹性力学基本方程弹性力学基本方程描述弹性体内任一点应力、应变、位移以及外 力之间的关系,它包括平衡方程、儿何方程和物理方程三类。1 平衡方程(应力和体力之间关系)(do)cr +x dxYdydz - o dxdy +(drt +w dyV Ydxdy)z、dz dxdy - t . dxdy + Xdxdydz -1z在X方向有工耳=do+
6、dx dy应力和体力在三个坐标方向上 满足一下平衡方程C + 6 + dx dy+ X = 0 dzdr+Y = 0 zdo+Z = 0 dz平衡方程是弹性体内部必须满足的条件,它说明六个应力分量不是独立的, 它们通过三个平衡方程相互联系。2儿何方程(儿何量位移和应变)X7 du t、 u Hdx+ dx - u-dx(丿dxdudx dxdu dvdududv8 -7 二二+dx9 xydxdvdwdv8 -人-=+yzdydwdudwE. dz/ 二/ xz+dx人= + 0 q tan + tan 0 =dv 竺. dydvdu1 + S 1 + edxdvxydxdx 1dx丿dydy
7、 1 +dv3物理方程(应力分量与应变分量;与材料的物理特性有关)从静力学角度导出了平衡微分方程和静力边界条件; 从几何学角度导出了几何方程和应变协调条件;在推导 过程中并没有涉及到弹性体本身材料的固有特性,故这 些方程适用于一切连续介质。从也 的角度分析可知,不同材料的弹性体其应 力应变关系即本构关系是不同的,对于对于理想弹性体 ,在小变形情况下,应力应变关系服从广义胡克定律物理方程的表达形式以应力表示应变以应变表示应力1- 二一(7 V(CT + CT )E-r-v7=Ga = X6 + 2Gs t = Gyxxxy/ xy1iqa = X9 + 2Gs t = Gvyy克/ yz = o
8、 一 v(cr + (7 y厂y*zEF Ga = A0 + 2Gs t = GyZZxz M1 - = 7 J70. +J)| /.=:0 = S + 8 + 8E 1-1GXyZb二咏为泊松比E为材料的切变弹性模量G =2(1 + )由上可见,三类基本方程中包括15个方程,含6个应力分量、6个应变分量和3个位移分量共15个未知量。实际求解时并不是同时求出全部未知量,而是先求出一部 分(称为基本未知量),再通过基本方程求出其他未知量。 位移法、应力法、混合法选取基本未知量不同四、平面问题1平面应力问题工程中链传动中的链片、发动机中的连杆.内燃机的飞轮、轧机的机架和齿宽较小的 直齿圆柱齿轮等a
9、)弹性体是等厚的薄板(潘 向等厚廃),厚度尺寸远远小于截面尺寸,tL/15;b)体力、面力和约束都只有平面内的量即(人,人;人,人;*),且都不沿z向变化;厂应力边界(面力和约束只作用于板边,在板面上 ()没任何面力和约束的作用。应力边界为:CT 1= 0z5z 二土一2T 1= 0ZXOZ 二土一 22板很薄,外力不沿厚度方向变化,因应力沿厚度 方向连续分布,故可认为所有各点:CT 二二 0Z由切应力互等定律得:T XZ O P W _ 0计只有平行于xy面的平面应力分量平面应力xbyxy yx严由于物体形状、外力和约束沿Z向均不变化,应 力分量和应变分量均只是兀、y的函数;从几何 方程积分
10、求位移可知位移与z有关。矗平面应力问题只有平面应力分量和Txy = Tyx,且仅为兀、歹的函数的弹性力学问题CT = 0 T T 0 T T = 0 入=J 二 cry(x,y) 乙厂為(兀,” Tyx = Tyx(x9y)物理方程几何方程cr二Ds物理方程式中1010001 - “EID=1 “称为平面应力问题的弹性矩阵2)平面应变问题工程中滚针轴承的滚针、轧钢机的轧辐、水坝.受内压管道、齿宽较大的直齿轮等2条件、a)弹性体为常截面的很长柱体。b)体力、面力和约束都只有q平面内的量,且都不沿Z向变化;假想柱体无限长,则任一Z截面均为对称面,即w = 0,只有平面位移况和u存在,即平面位移。A
11、由于截面形状.外力和约束沿z向均不变化,位移分量只是X、y的函数u - y) v = v(x, y) w = 0A假想柱体无限长,则任一 Z截面均为对称面,即w = o,只有平面位移和V存在,即平面位移。由于截面形状、外力和约束沿Z向均不变化,位移分量只是X、y的函数.u - y) v v(x, y) w = 0从数学和几何学角度推导dudvu = y)dw dv+ =0dy dzdw=0Szrxzdu dw+ =0 dz dx只有平行于xy面的平面应变分量平面应变从力学角度推导由对称性(对称结构承受对称荷载,反对称力为 零)可知: =0、T =0zxzy歼由剪切互等定律可知: T x_ =
12、0、 Ty_=0厂由胡克定律可知:冷二丫口 = 、乙厂乙尸*平面应变问题只有平面应变分量、)和/rv = Zvr,且仅为兀、歹的函数的弹性力学问题rzx = rxz心=N = Sx=sx(x,y) y = Wy (兀,歹) rxy =乙),(兀,刃 ryx = 物理方程6=云-“汁+%=“(6 + by)6的存在说明了沿Z向无限长的柱体的假设限 制了每一个横截面的纵向位移。当柱体受到垂直于 轴的外力作用时,这些衡截面之间必然产生挤压 应力。冬 可直接由、6计算得到,故不作为独立的 未知量。物理方程物理方程= &式中D =E(l-“)(! + /)(!- 2“)1 01 2“2(1 “)称为平面应变问题的弹性矩阵名称平面应力问题平面应变问题未知量已知量未知量已知量位移况、Vw工0%、 Vw = 0应变、 F 、 兀yrXyY = y 0As =-(cr + a )乙厂兀E、 F 、兀yV xyv = y =w=0应力CT 、 CT 、 兀yTT = T = CT = 0CT、 CT、xyTxyT - T