2-有限元静力分析基本原理

2001年10月1日Lb-1二、结构静力有限元分析1、材料力学的局限2、弹性力学基础3、平面应力问题的有限元分析原理入门4、ANSYS静力分析实践入门2001年10月1日Lb-21、材料力学的局限带孔板单向拉伸问题的应力确定。图1-3a图1-3bntKt应力集中系数minAPn名义应力2001年10月1日Lb-32、弹性力学基础2.1弹性力学基本方程2.2弹性力学求解2.3主应力2.4强度理论及相当应力2.5弹性力学问题简化2001年10月1日Lb-42.1弹性力学基本方程（1）平衡微分方程—应力体力关系（2）静力边界条件—应力面力关系（3）几何方程—应变位移关系（4）物理方程—应力应变关系2001年10月1日Lb-5Z总和后整理便得到X方向的平衡方程:同理得到x、y方向的平衡方程:XYZ其中X、Y、Z为物体的体力分量。（1）平衡微分方程——应力体力关系xyzoxyxzyxyyzyyyzyzdyyyydyyyxyxdxzyzxdzzyzyzdzzxzxzzzzzzdxxxxdxxxzxzdxxxyxydxzyyyzyzdzzzzzdxxxzxzdyz2001年10月1日Lb-6设斜面ACD为边界面，其外法线n的方向为(l1,l2,l3)，面积为ΔS，边界外力分量为(px，py，,pz)，则三角形ABC、ABD、BCD的面积分别为ΔS在各相应方向上的投影。l1ΔS,l2ΔS,l3ΔS,NxyzoNXNZNY（2）静力边界条件——应力面力关系由x方向的平衡得到：XNΔS=l1ΔSσx+l2ΔSτyx+l3ΔSτzx即XN=l1σx+l2τyx+l3τzx注意，这里边界上的外力是坐标轴方向上的分量。XN=l1σx+l2τyx+l3τzxYN=l1τxy+l2σy+l3τzyZN=l1τxz+l2τyz+l3σz边界条件为：2001年10月1日Lb-7vudxdyABCDdxxuu?dxxvvdyyuu??dyyvv??'A'B'C'DDBbaxy0图1-5A点在X方向的位移分量为u；B点在X方向的位移：ABCD---A’B’C’D’求正应变，用位移分量来表示：xy、xdxxuuuuxudxudxxuux)(线素AB的转角为：A点在Y方向的位移分量为v；B点在Y方向的位移分量：dxxvvBABBtgaaxuxvdxxudxvdxxvv1)(线素AB的正应变为：由于变形是微小的，所以上式可将比单位值小得多的略去，得xuxvayubyuxvxyba因此，剪应变为：同理，Y向线素AD的转角yuxvxwzuzvywzwyvxuxyzxyzzyx,,,,空间问题的位移分量为：u、v、w位移边界条件：wwvvuusss,,（3）几何方程——应变位移关系2001年10月1日Lb-8（4）物理方程——应力应变关系虎克定律ExxxEσμ，εEσμεzxzyxy波松效应EG)1(2yzyzxzxzxyxyyxzzzxyyzyxxGGGEEE1,1,1)]([1)]([1)]([1)]([zyxzyxxzxyxxxE1EσμEσμE图1-42001年10月1日Lb-92.2弹性力学求解yuxvxwzuzvywzwyvxuxyzxyzzyx,,,,TzxyzxyzyxTzxyzxyzyxTzyxwzyxvzyxuf应力应变位移.3.2)],,(),,,(),,,([}{.1XYZyzyzxzxzxyxyyxzzzxyyzyxxGGGEEE1,1,1)]([1)]([1)]([1平衡微分方程（3个）静力边界条件（3个）物理方程（6个）几何方程（6个）一点处自由度（3个）应变（6个）应力（6个）弹性体内部和边界上任意一点处有15各方程可解15个未知数，但很少有解析解。XN=l1σx+l2τyx+l3τzxYN=l1τxy+l2σy+l3τzyZN=l1τxz+l2τyz+l3σz2.3主应力在单元体上的三个相互垂直的平面上即可能是正应力，剪应力或者二者的组合。•主单元体：各个侧面上的剪切应力均为零的单元体；•主平面：单元体上没有切应力的面称为主平面；•主应力：作用在主平面上的正应力。),,,,,(321zxyzxyzyxf2001年10月1日Lb-11相当应力:强度条件中直接与许用应力[σ]比较的量，称为相当应力σr(形状改变比能理论)(最大剪应力理论)(最大拉应力理论)2132322214)()()(21r313r11r3212r（最大伸长线应变理论）2.4强度理论及相当应力(3)强度条件的一般形式r≤[]强度理论是判断材料在复杂应力状态下是否破坏的理论。材料在外力作用下有两种不同的破坏形式：一是在不发生显著塑性变形时的突然断裂,称为脆性破坏;二是因发生显著塑性变形而不能继续承载的破坏，称为塑性破坏。2001年10月1日Lb-12（1）平面应力问题的特点（2）平面应力问题求解（3）二维实体问题总结2.5弹性力学问题简化2001年10月1日Lb-131、弹性力学平面问题模型所谓平面问题指弹性力学的平面应力和平面应变问题。(1)平面应力问题（薄板拉压）A.几何条件：结构是等厚度薄板，即厚度远远小于截面尺寸；B.载荷条件：载荷均平行与板面且沿板厚度方向均匀分布。（1）平面应力问题的特点0zyzxz20tz板面中心板面C.特点：2001年10月1日Lb-14（2）平面应力问题求解yuxvyvxuxyyx,,TxyyxTxyyxTzyxvzyxuf应力应变位移.3.2)],,(),,,([}{.1xyxyzxyyzyxxGEE1)]([1)]([1平衡微分方程（2个）静力边界条件（2个）物理方程（3个）几何方程（3个）一点处自由度（2个）应变（3个）应力（3个）弹性体内部和边界上任意一点处有8各方程可解8个未知数。比三维实体问题减少近1倍方程。XN=l1σx+l2τyxYN=l1τxy+l2σy00YyxXyxyxyxyx2001年10月1日Lb-151、弹性力学平面问题模型所谓平面问题指弹性力学的平面应力和平面应变问题。(1)平面应力问题（薄板拉压）A.几何条件：结构是等厚度薄板，即厚度远远小于截面尺寸；B.载荷条件：载荷均平行与板面且沿板厚度方向均匀分布。C.特点：(2)平面应变问题（长体均剪）A.几何条件：结构是等截面长体，即长度远远大于截面尺寸；B.载荷条件：载荷均垂直于长度方向且沿长度方向均匀分布。C.特点：2、薄板弯曲问题几何条件：结构是等厚度薄板，即厚度远远小于截面尺寸；B.载荷条件：载荷垂直于板面。C.特点：3、轴对称问题A.几何条件：结构是旋转体，即通过子午面绕对称轴旋转而成；B.载荷条件：载荷关于结构对称轴满足轴对称。C.特点：4、杆系结构A.几何条件：长度远大于宽度和厚度的构件。杆、梁、柱、轴等；B.载荷条件：任意；C.特点：符合材料力学，杆单元和梁单元。（3）二维实体问题总结0zyzxz0yzzxz0,0zrzr0||,000zzzwuPPointALH2.6结构有限元法分析的目的及适用条件•进行结构的最优方案设计在进行机械结构设计时，可以通过对可能的结构方案进行有限单元法计算,根据对方案计算结果的分析和比较，按强度、刚度和稳定性要求，对原方案进行修改补充，以便得到较合理的应力、变形分布，并且经济性又较好的设计方案。•分析结构损坏原因，寻找改进途径当结构件在工作中发生故障如裂纹、断裂、磨损过大等时，可应用有限单元法进行计，研究结构损坏的原因，找出危险区域和部位，提出改进设计的方案，并进行相应的计算分析，直至找到合理的结构为止。（2）结构静力分析适用条件结构静力分析用于确定在静载荷引起的结构位移、应力和应变等效应。静载荷是一种假定，即假定载荷随时间的变化非常缓慢。静力分析既可以是线性的也可以是非线性的，非线性静力分析的类型包括：大变形、塑性、蠕变、应力刚化、接触（间隙）单元、超弹性单元等。（1）目的2001年10月1日Lb-173、平面应力问题的有限元分析原理入门问题描述：为说明解平面应力问题的有限元方法，我们给出一个问题详细的解。图5-20所示薄板受表面T=1000MPa的拉力作用，确定节点位移和单元应力。板厚为1mm，E=30×106MPa，μ=0.30。图5-20受拉薄板2001年10月1日Lb-183.1结构的离散化Procedure1......2......3......用有限元法对结构进行应力分析时，首先要将结构进行离散化。即将一个连续体看成由有限个单元组成的体系。弹性力学平面问题中最常见的单元是三角形单元。所有作用在单元上的载荷都按静力等效的原则移置到结点上，并在受几何约束的结点处设置相应的铰支座。这样就得到了用以代替原来弹性体的有限单元计算模型。平面应力问题单元是均质等厚度薄板，常用的单元由三节点三角形单元、六节点三角形单元、四节点矩形单元和八节点矩形单元等2001年10月1日Lb-193.2单元分析3.2.1位移模式(位移插值函数)3.2.2单元中的应变3.2.3单元中的应力3.2.4节点位移与节点力的关系2001年10月1日Lb-203.2.1位移函数（1）位移函数形式（2）待定系数确定（3）形函数2001年10月1日Lb-21（1）位移函数形式在有限单元位移法中，假设结点上的位移是基本未知量。为了能用单元的结点位移表示单元中的应变和应力分量，必须假定一个位移模式，也就是说根据单元的结点位移去构造单元上的位移插值函数。位移插值函数采用多项式形式（泰勒展开式或麦克劳林展开式）待定系数个数=单元自由读数取一个典型的三角形单元进行力学分析。采用线性插值，即假定单元上的位移分量是坐标的线性函数：Tyxvyxuf)],(),,([}{一点处的自由度2121121098726524321yaxyaxayaxaavyaxyaxayaxaau2001年10月1日Lb-22其中：而：是三角形ijm的面积。（2）待定系数确定它们可以由结点位移确定如下：联立求解上述方程，可得：2001年10月1日Lb-23于是可以得到：其中：同理得：（3）形函数mmjjiimjimjivuvuvuNNNNNNvu000000e][Nf2001年10月1日Lb-243.2.2单元中的应变ModuleObjective一点处的应变yuxvyvxuxyyxvuxyyxyuxvyvxuxyxx00eeBNf2001年10月1日Lb-253.2.3单元中的应力ModuleObjective物理方程：弹性矩阵xyxyzxyyzyxxGEE