有限元法基础杆系结构力学问题.pptx

《有限元法基础杆系结构力学问题.pptx》由会员分享，可在线阅读，更多相关《有限元法基础杆系结构力学问题.pptx（54页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、9.杆件结构力学问题 本章要点本章要点l结构单元的力学假设结构单元的力学假设l经典梁单元和经典梁单元和Timoshenko梁单元梁单元l平面杆系和空间杆系的单元矩阵和向总体坐标的转换平面杆系和空间杆系的单元矩阵和向总体坐标的转换有限元法基础1第1页/共54页9.杆件结构力学问题 关键概念关键概念 结构单元结构单元 C C1 1连续连续 经典梁理论经典梁理论 Timoshenko梁理论梁理论 平面杆系平面杆系 空间杆系空间杆系 剪切自锁（剪切自锁（shear locking)有限元法基础2第2页/共54页9.1 结构单元概论l结构单元结构单元 杆梁单元杆梁单元 板壳单元板壳单元特点：特点：1 1

2、）在一个方向或二个方向上的几何尺度远小于其他方）在一个方向或二个方向上的几何尺度远小于其他方向；向；2 2）引入一定的假设，相应简化为一维或二维问题。）引入一定的假设，相应简化为一维或二维问题。有限元法基础结构单元结构单元3第3页/共54页9.1 结构单元概论l使用实体单元计算结构问题，理论上是可行的使用实体单元计算结构问题，理论上是可行的引起的问题引起的问题1 1）如果划分网格适应结构特点，使单元的一个方向或）如果划分网格适应结构特点，使单元的一个方向或两个方向的尺度比其他方向小很多，将使单元刚度系两个方向的尺度比其他方向小很多，将使单元刚度系数差异过大，影响求解精度数差异过大，影响求解精度

3、2 2）保持单元各个方向的尺度接近，将导致单元总数过）保持单元各个方向的尺度接近，将导致单元总数过分庞大分庞大有限元法基础4第4页/共54页9.1 结构单元概论l例例：一一由由三三个个弹弹簧簧单单元元组组成成的的系系统统，单单元元刚刚度度为为K K1 1K K2 210105 5，K K3 31 1，施加于节点的轴力为，施加于节点的轴力为P P2 20 0，P P1 1P P3 31 1。利用。利用5 5为有效数字的算法求解位移。为有效数字的算法求解位移。有限元法基础5第5页/共54页9.1 结构单元概论系统平衡方程系统平衡方程有限元法基础6第6页/共54页9.1 结构单元概论代入数据代入数据

4、高斯消去法求解高斯消去法求解有限元法基础求解失败10000011000007第7页/共54页9.1 结构单元概论解决办法解决办法 1）提高计算有效位数）提高计算有效位数 有限元法基础8第8页/共54页9.1 结构单元概论2）相对自由度方法）相对自由度方法(relative DOF)令令代入有限元方程代入有限元方程 有限元法基础9第9页/共54页9.1 结构单元概论3）主从自由度方法（）主从自由度方法（master-slave DOF）分析系统的特点，由于分析系统的特点，由于 ，近似地认为，近似地认为将将u3作为主自由度，其他为从自由度。只需对作为主自由度，其他为从自由度。只需对u3 进进行建立

5、平衡条件行建立平衡条件由于假设产生了误差。由于假设产生了误差。有限元法基础10第10页/共54页9.1 结构单元概论l基于主从自由度原理的梁单元和板壳单元基于主从自由度原理的梁单元和板壳单元 结结构构力力学学中中的的梁梁理理论论和和板板壳壳理理论论引引入入变变形形假假设设实实际际上上正正是是应应用用了了主主从从自自由由度度的的原原理理，将将问问题题归归结结为为中中面位移。面位移。l基于相对自由度的梁单元和板壳单元基于相对自由度的梁单元和板壳单元 这这类类单单元元 本本质质上上是是二二维维、三三维维实实体体单单元元，为为使使中中面面法法线线变变形形后后保保持持直直线线，沿沿厚厚度度方方向向只只设

6、设置置两两个个节节点点，其其中中一一个个保保持持原原来来意意义义上上的的自自由由度度，另另一一个个节节点的自由度用它与上述节点的相对位移来代替。点的自由度用它与上述节点的相对位移来代替。有限元法基础11第11页/共54页9.2 等截面直杆-梁单元一一.杆单元杆单元 承受轴向拉压载荷，并只经历轴向拉压变形的细承受轴向拉压载荷，并只经历轴向拉压变形的细长构件，称为杆（长构件，称为杆（Rod）。）。假设：应力在截面上均匀分布，原来垂直于轴线假设：应力在截面上均匀分布，原来垂直于轴线的截面变形后仍然垂直于轴线。的截面变形后仍然垂直于轴线。有限元法基础12第12页/共54页9.2 等截面直杆-梁单元l直

7、杆的基本方程直杆的基本方程 几何关系几何关系 本构方程本构方程 平衡方程平衡方程 端部条件端部条件 有限元法基础13第13页/共54页9.2 等截面直杆-梁单元l最小势能原理最小势能原理 有限元法基础14第14页/共54页9.2 等截面直杆-梁单元l单元列式单元列式位移插值函数位移插值函数 其中其中NN1,N2,Nn,qeu1,u2,unT坐标变换坐标变换2 2节点单元节点单元3 3节点单元节点单元 有限元法基础15第15页/共54页9.2 等截面直杆-梁单元单元矩阵单元矩阵 由变分原理变分取驻值，得有限元方程由变分原理变分取驻值，得有限元方程其中其中 有限元法基础16第16页/共54页9.2

8、 等截面直杆-梁单元二二.扭转单元扭转单元 几何关系几何关系 本构关系本构关系 平衡方程平衡方程 端部条件端部条件 一般一般 m(x)0，在端部给定扭矩，在端部给定扭矩M。有限元法基础17第17页/共54页9.2 等截面直杆-梁单元最小势能原理最小势能原理设转角的插值函数为设转角的插值函数为单元刚度矩阵单元刚度矩阵 有限元法基础18第18页/共54页9.2 等截面直杆-梁单元l扭转单元的限制扭转单元的限制 1 1）采用自由扭转理论，截面在变形后依然保持平）采用自由扭转理论，截面在变形后依然保持平面，这只适用于圆截面情形；面，这只适用于圆截面情形；2 2）在约束扭转时，截面将发生翘曲，需引入翘曲

9、）在约束扭转时，截面将发生翘曲，需引入翘曲修正因子；修正因子；3 3）只有在截面有两个对称轴时，如圆、椭圆矩形，）只有在截面有两个对称轴时，如圆、椭圆矩形，截面才是绕形心转动。截面才是绕形心转动。有限元法基础19第19页/共54页9.2 等截面直杆-梁单元三三.梁单元梁单元 有限元法基础经典梁理论假设经典梁理论假设经典梁理论假设经典梁理论假设：变形前垂直于中心线的截面，变：变形前垂直于中心线的截面，变形后仍保持为平面，且垂直于中心线。形后仍保持为平面，且垂直于中心线。20第20页/共54页9.2 等截面直杆-梁单元l基本方程基本方程 几何方程几何方程 本构关系本构关系 平衡方程平衡方程 端部条

10、件端部条件 有限元法基础21第21页/共54页9.2 等截面直杆-梁单元l最小势能原理最小势能原理本质边界条件：位移端部条件本质边界条件：位移端部条件自然边界条件：力的端部条件自然边界条件：力的端部条件 有限元法基础端部条件本质边界条件自然边界条件简支固定自由22第22页/共54页9.2 等截面直杆-梁单元l单元列式单元列式 梁弯曲是梁弯曲是C C1 1类问题，采用类问题，采用Hermite插值插值 有限元法基础23第23页/共54页9.2 等截面直杆-梁单元l将插值函数代入泛函，由将插值函数代入泛函，由 得到有限元方程得到有限元方程l单元刚度矩阵为单元刚度矩阵为 有限元法基础24第24页/共

11、54页9.2 等截面直杆-梁单元4.4.考虑剪切变形的梁单元考虑剪切变形的梁单元 经典梁理论只适用于细长梁，经典梁理论只适用于细长梁，即梁高远小于跨度。即梁高远小于跨度。有限元法基础l 变形假设变形假设 变形前垂直于中线的截面，变形前垂直于中线的截面，变形后仍然保持平面，但不再变形后仍然保持平面，但不再垂直于中线。垂直于中线。25第25页/共54页9.2 等截面直杆-梁单元l剪切变形剪切变形l曲率与转角的关系曲率与转角的关系l最小势能原理最小势能原理 有限元法基础折剪系数折剪系数26第26页/共54页9.2 等截面直杆-梁单元l关于折剪系数关于折剪系数 由于变形假设，剪应力在截面上为常数，由于

12、变形假设，剪应力在截面上为常数，但是，实际情况剪应力不是均匀分布的，但是，实际情况剪应力不是均匀分布的，而是按抛物线分布。因此需引进校正因子，而是按抛物线分布。因此需引进校正因子，剪力和剪应变的关系为剪力和剪应变的关系为k的的值值因因截截面面形形状状的的不不同同而而异异，一一般般按按能能量量等等效效的的方方法法计计算算，即即简简化化后后的的剪剪切切应应变变能能与与实实际际剪剪应应力力和和剪剪应应变分布计算的应变能相等。变分布计算的应变能相等。有限元法基础27第27页/共54页9.2 等截面直杆-梁单元l单元列式单元列式位移插值位移插值单元刚度矩阵单元刚度矩阵 有限元法基础28第28页/共54页

13、9.2 等截面直杆-梁单元载荷项载荷项其他矩阵的定义其他矩阵的定义 有限元法基础29第29页/共54页9.2 等截面直杆-梁单元l剪切自锁（剪切自锁（Shear locking）以以2 2节点单元为例，有节点单元为例，有 有限元法基础精确积分：精确积分：2个积分点个积分点矩阵秩为矩阵秩为230第30页/共54页9.2 等截面直杆-梁单元注意注意:是梁高是梁高h的的3 3次函数次函数 是梁高是梁高h的的1 1次函数次函数当当 时，并且时，并且 是非奇异的是非奇异的 有限元方程有限元方程 零解零解解决办法解决办法 有限元法基础“剪切自锁剪切自锁”现象现象（1 1）减缩积分）减缩积分（2 2）假设应

14、变）假设应变31第31页/共54页9.2 等截面直杆-梁单元（1 1）减缩积分）减缩积分 2 2节点单元对节点单元对 使用使用1 1点积分点积分 使用使用1 1点积分点积分 有限元法基础矩阵秩为矩阵秩为132第32页/共54页9.2 等截面直杆-梁单元产生原因产生原因 2 2节点单元的剪应变节点单元的剪应变当当 ,要求要求 有限元法基础第一式表示转角与挠度的关系第一式表示转角与挠度的关系第二式导致转角相等，不能表示弯曲变形第二式导致转角相等，不能表示弯曲变形 33第33页/共54页9.2 等截面直杆-梁单元（1 1）假设应变）假设应变 由于直接导出的剪切应变可能不正确，因此，重由于直接导出的剪

15、切应变可能不正确，因此，重新假设剪应变新假设剪应变并令并令 有限元法基础 插值点取在应变计算比较准确的点上，一般插值点取在应变计算比较准确的点上，一般在应力最优点上，如在应力最优点上，如2 2节点单元去在形心点处，节点单元去在形心点处，3 3节点单元取节点单元取2 2点点GaussGauss积分点上。积分点上。34第34页/共54页9.2 等截面直杆-梁单元例：计算悬臂梁问题例：计算悬臂梁问题 使用一个单元计算使用一个单元计算（1 1）经典梁单元）经典梁单元 有限元法基础35第35页/共54页9.2 等截面直杆-梁单元（2 2）Timoshenko梁单元梁单元纯弯曲纯弯曲 有限元法基础减缩积分

16、减缩积分精确积分精确积分36第36页/共54页9.2 等截面直杆-梁单元（2 2）Timoshenko梁单元梁单元剪切力剪切力 有限元法基础减缩积分减缩积分精确积分精确积分罚因子罚因子误差误差2537第37页/共54页9.3 平面杆系l特点特点（1）平面杆系受力状态）平面杆系受力状态 a）横向弯曲）横向弯曲 b）沿轴向的拉压）沿轴向的拉压（2）在局部坐标下）在局部坐标下 平面杆单元梁单元受轴向力杆单元平面杆单元梁单元受轴向力杆单元（3）平面杆系单元通过坐标转换）平面杆系单元通过坐标转换 单元坐标单元坐标 总体坐标总体坐标 坐标转换坐标转换 有限元法基础38第38页/共54页9.3 平面杆系l单

17、元刚度矩阵单元刚度矩阵l例：例：2 2节点单元刚度矩阵节点单元刚度矩阵 有限元法基础39第39页/共54页9.3 平面杆系l坐标变换坐标变换 总体坐标系总体坐标系 单元坐标系单元坐标系 单元坐标与总体坐标的单元坐标与总体坐标的夹角为夹角为 ，方向余弦为，方向余弦为 有限元法基础40第40页/共54页9.3 平面杆系节点节点DOF的转换关系的转换关系 有限元法基础41第41页/共54页9.3 平面杆系总体坐标下的刚度矩阵总体坐标下的刚度矩阵 有限元法基础42第42页/共54页9.3 平面杆系l铰节点的处理铰节点的处理 铰接的特点铰接的特点：在铰接处的所有节点位移相同，转：在铰接处的所有节点位移相

18、同，转 动不同，不承受弯矩动不同，不承受弯矩 处理办法处理办法：在集成总体刚度矩阵前将相关转动：在集成总体刚度矩阵前将相关转动 DOF DOF凝聚，不与其他单元的刚度矩阵集成。凝聚，不与其他单元的刚度矩阵集成。有限元法基础43第43页/共54页9.3 平面杆系 有限元法基础44第44页/共54页9.3 平面杆系如图如图2 2的一个端点与刚架铰接的一个端点与刚架铰接 单元刚度矩阵为单元刚度矩阵为 为与为与4 4节点连接的铰接节点连接的铰接DOFDOF 由第二组方程得由第二组方程得 有限元法基础45第45页/共54页9.3 平面杆系凝聚后的方程组凝聚后的方程组有限元法基础为了编程方为了编程方便，凝

19、聚后便，凝聚后的刚度矩阵，的刚度矩阵，仍以零元素仍以零元素补充为补充为6X6的矩阵的矩阵46第46页/共54页9.4 空间杆系空间杆系的受力特点空间杆系的受力特点1 1）两个方向的弯曲）两个方向的弯曲2 2）沿轴向的拉压）沿轴向的拉压3 3）绕轴线的扭转）绕轴线的扭转有限元法基础力和平移力和平移力矩和转动力矩和转动单元坐标和总体坐标单元坐标和总体坐标坐标转换关系坐标转换关系1 1）3 3个角度定位个角度定位2 2）9 9个方向余弦个方向余弦47第47页/共54页9.4 空间杆系l空间杆单元空间杆单元 每节点有每节点有6 6个个DOFDOF，分别对应，分别对应3 3个平移和个平移和3 3个转动个

20、转动单元刚度矩阵中的截面几何常数：单元刚度矩阵中的截面几何常数：截面面积截面面积 A xz面内的惯性矩面内的惯性矩 Iyxy面内的惯性矩面内的惯性矩Iz 绕绕x轴的扭转惯性矩轴的扭转惯性矩 Ix翘曲常数翘曲常数 扭转常数扭转常数有限元法基础48第48页/共54页9.4 空间杆系l空间杆单元空间杆单元 每节点有每节点有6 6个个DOFDOF，分别对应，分别对应3 3个平移和个平移和3 3个转动个转动单元刚度矩阵中的截面几何常数：单元刚度矩阵中的截面几何常数：截面面积截面面积 A xz面内的惯性矩面内的惯性矩 Iyxy面内的惯性矩面内的惯性矩Iz 绕绕x轴的扭转惯性矩轴的扭转惯性矩 Ix有限元法基

21、础49第49页/共54页9.4 空间杆系ANSYSANSYS中中BEAM188BEAM188单元单元截面性质截面性质 有限元法基础yz50第50页/共54页9.4 空间杆系在空间中的位置在空间中的位置 由单元端部节点坐标可确认由单元端部节点坐标可确认单元在总体坐标的位置，但截单元在总体坐标的位置，但截面在总体坐标中的方位需引进面在总体坐标中的方位需引进辅助节点来确定。辅助节点来确定。有限元法基础*辅助节点辅助节点xz51第51页/共54页9.4 空间杆系坐标变换坐标变换 单元坐标单元坐标 x-y-z 与总体坐标系与总体坐标系 的方向余的方向余弦为弦为变换关系为变换关系为 有限元法基础52第52

22、页/共54页9.5 小结 1 1）基于经典梁理论（）基于经典梁理论（Bernoulli梁理论梁理论）和基于）和基于Timoshenko梁理论分别要求梁理论分别要求C1 1类和类和C0 0类连续性；类连续性；2 2）Bernoulli梁单元梁单元采用采用Hermite插值容易插值容易C1 1满足连满足连续性要求；续性要求；3）Timoshenko梁单元降低了连续性要求，出现罚梁单元降低了连续性要求，出现罚因子现象，可能出现因子现象，可能出现“剪切自锁剪切自锁”；4 4）本章讨论只限于等截面直杆，也没有考虑约束）本章讨论只限于等截面直杆，也没有考虑约束扭转的影响；扭转的影响；5 5）对变截面杆，非规则截面等还要对单元刚度矩）对变截面杆，非规则截面等还要对单元刚度矩阵进行修正和补充。阵进行修正和补充。有限元法基础53第53页/共54页54感谢您的观看！第54页/共54页