1、有限元复习一、选择题(每题1分,共10分)二、判断题(每空1分,共10分)三、填空题(每空1分,共10分)三、简答题(共44分)共6题四、综述题(共26分)两题一基本概念1. 平面应力/平面应变问题;空间问题/轴对称问题;杆梁问题;线性与非线性问题平面应力问题(1) 均匀薄板(2)载荷平行于板面且沿厚度方向均匀分布在六个应力分量中,只需要研究剩下的平行于 XOY 平面的三个应力分量,即 ( ) 。xyxyx、 、 000zzxzy, ,一般 , 并不一定等于零,但可由 及 求得,在分析问题时不必0zz 考虑。于是只需要考虑 三个应变分量即可。xyx、 、平面应变问题(1) 纵向很长,且横截面沿
2、纵向不变。 (2)载荷平行于横截面且沿纵向均匀分布只剩下三个应变分量 。也只需要考虑0zyzxxyx、 、三个应力分量即可x、 、轴对称问题物体的几何形状、约束情况及所受外力都对称于空间的某一根轴。轴对称单元的特点(与平面三角形单元的区别):轴对称单元为圆环体,单元与单元间为节圆相连接;节点力与节点载荷是施加于节圆上的均布力;单元边界是一回转面;应变不是常量。 在轴对称问题中,周向应变分量 是与有关。r板壳问题 一个方向的尺寸比另外两个方向尺寸小很多,且能承受弯矩的结构称为板壳结构,并把平分板壳结构上下表面的面称为中面。如果中面是平面或平面组成的折平面,则称为平板;反之,中面为曲面的称为壳。杆
3、梁问题杆梁结构是指长度远大于其横断面尺寸的构件组成的系统。在结构力学中常将承受轴力或扭矩的杆件称为杆,而将承受横向力和弯矩的杆件称为梁。平面(应力应变)问题与板壳问题的区别与联系平面应力问题是指很薄的等厚度薄板,只在板边上受有平行于板面并且不沿厚度变化的面力,同时,体力也平行于板面并且不沿厚度变化。而平面应变问题是指很长的柱形体,在柱面上受有平行于横截面并且不沿长度变化的面力,同时体力也平行于横截面并且不沿长度变化。板壳问题的弹性体受垂直于板面的力的作用,板将变成有弯有扭的曲面。 线性问题/非线性问题线性问题:基于小变形假设,应力与应变方程、应力与位移关系方程、平衡方程都是线性的。非线性问题:
4、材料非线性(非线性弹性、非线性弹塑性),几何非线性(大变形大应变如金属橡胶,小应变大位移如薄壁结构)2. 不同类型单元的节点自由度的理解:3. 有限元法的基本思想与有限元分析的基本步骤(5 步)有限元法的基本思想:离散、分片插值;其中离散的思想吸收了差分法的启示。有限元分析的基本步骤:数学建模(问题分析) ,结构离散(第一次近似) ,单元分析(位移函数,单刚方程) (第二次近似) ,整体分析与求解(总刚度方程,引入约束,解方程组求节点位移,根据节点位移求应力) ,结果分析及后处理。单元类型 节点数 节点自由度杆单元 2 1梁单元 2 3平面单元 3 2平面四边形 4 2轴对称问题 3 2板壳单
5、元 4 3四面体单元 4 34. 里兹法的基本思想及与有限元法区别里兹法的基本思想:先根据描述问题的微分方程和相应定解条件构造等价的泛函变分形式,然后在整个求解区域上假设一个试探函数(或近似函数),通过求解泛函极值来获得原问题的近似解。与有限元法的区别:里兹法是整体场函数用近似函数代替,有限元法是离散求解域,分片连续函数来近似整体未知场函数。5. 有限元法的基本定义(节点、单元、节点力、节点载荷) 单元:即原始结构离散后,满足一定几何特性和物理特性的最小结构域 节点:单元与单元间的连接点。 节点力:单元与单元间通过节点的相互作用力 节点载荷:作用于节点上的外载(等效) 。6. 位移函数的构造方
6、法及满足的基本条件构造方法:(1)广义坐标法,按照帕斯卡三角形选择多项式,项数多少由单元的自由度数决定。 (2)插值函数法,表示为形函数和节点位移的乘积表示。基本条件:(1)位移函数在单元节点的值应等于节点位移(即单元内部是连续的) ;(2)所选位移函数必须保证有限元的解收敛于真实解。7. 位移函数的收敛性条件(协调元、非协调元)及单元协调性的判断位移函数的收敛性条件(1)位移函数应包含刚体位移(2)位移函数应包含常量应变(反映单元的常应变状态)(3)位移函数在单元内连续,在单元之间的边界上要协调满足 1 和 2 称为完备单元,满足 1,2,3 称为协调单元。单元协调性的判断以 3 节点三角形
7、单元为例,位移分量在每个单元中都是坐标的线性函数的话,在公共边界上也会是线性变化的,那么相邻单元在公共边界上的任意一点都具有相同的位移,也就是协调单元。有限元法中,假设一种位移函数近似表达单元内部的真实位移分布,该位移函数可表示为位移函数和节点位移的线性插值。8. 有限元解的性质有限元解具有下限性质,即有限元的解小于实际的精确解。这是因为实际结构本来是具有无限自由度的,当用有限元求解时,结构被离散为有限个单元的集合后,便只有有限个自由度了。由无限自由度变为有限自由度可以认为是对真实位移函数增加了约束,限制了结构的变形能力,从而导致结构的刚度增大、计算的位移减小。9. 虚功原理、最小势能原理及变
8、分法(里兹法)虚功原理:在力的作用下处于平衡状态的体系,当发生与约束条件相符合的任意微小的虚刚体位移时,体系上所有的主动力在虚位移上所作的总功(各力所作的功的代数和)恒等于零。最小势能原理:表明在满足位移边界条件的所有可能位移中,实际发生的位移使弹性体的势能最小。10.形函数特性1)形函数 Ni 为 x、y 坐标的函数,与位移函数有相同的阶次。2)形函数 Ni 在 i 节点处的值等于 1,而在其他节点上的值为 0。3)单元内任一点的形函数之和恒等于 1。4)形函数的值在 0-1 间变化。11.单元刚度矩阵的性质及元素的物理意义单元刚度矩阵的性质特点:(1)对称性(2 )奇异性,|K|=0(3)
9、主对角线元素恒为正值( 4)奇偶行元素之和分别为零(各行或各列元素之和为零)物理意义:单元刚阵K 的物理意义是单元受节点力作用后抗变形的能力。其中分块矩阵K ij的物理意义为:当在 j 节点处产生单位位移而其他节点位移为零时,在 i 节点上需要作用力的大小。其中元素 Kij 表示在第 j 号自由度上产生单位位移时,其他自由度位移为零时,在 i 号自由度上所需要施加的力的大小。单元刚度矩阵的元素表示该单元的各节点沿坐标方向发生单位位移时引起的节点力,它决定于该单元的形状、大小、方位和弹性常数,而与单元的位置无关,即不随单元或坐标轴的平行移动而改变。12.边界条件处理(载荷等效移置 集中力/均布力
10、/ 线性分布力 边界位移约束处理 固定/指定位移等)载荷等效移置连续弹性体离散为单元组合体时,为简化受力情况,需把弹性体承受的任意分布的载荷都向节点移置(分解),而成为节点载荷。载荷移置的原则:能量等效(或静力等效原则) ,即单元的实际载荷与移置后的节点载荷在相应的虚位移上所做的虚功相等。集中力,移置到两端节点,使得 F1 L1 =F2 L 2, F1 +F2=F均布力,移置到两端节点,F 1 =F2=0.5qL线性分布力, F1=1/3 0.5qL ,F 2=2/3 0.5qL边界位移约束一绝对位移约束刚性支座(活动铰支,固定铰支,固接支座) 固定位移弹性支座(线弹性制作,非线性支座) 可变
11、位移强迫约束 指定位移 用载荷等效,装配应力+整体应力二相对位移约束 (如两接触面)1.约束等式2.耦合约束(连接重合节点,模拟滑动边界连接,施加周期对称边界条件)常见的位移约束问题处理约束不足的处理(1)利用对称性引进约束(取 1/n 后,在对称面上施加位移约束)(2)转换载荷为位移约束(受平衡载荷作用,将一部分载荷用位移约束代替)(3)人为增加约束(约束点应尽量远离重要部位,约束点变形要相对小)其他,杆离散为多个杆单元时,须在连接节点增加约束,只允许产生轴向位移。轴对称结构,施加轴向约束。过约束的处理有时需要施加过约束,有时需要释放过约束。引入位移边界条件是为了消除有限元整体刚度矩阵 的奇
12、异性。K13.总体刚度矩阵组装原则及总刚阵特点总体刚度矩阵组装原则:1.在整体离散结构变形后,应保证各单元在节点处仍然协调地相互连接,即在该节点处所有单元在该节点上有相同位移。2.整体离散结构各节点应满足平衡条件。即环绕每个节点的所有单元作用其上的节点力之和应等于作用于该节点上的节点载荷 Ri。总刚度矩阵特点:除了具有单元刚阵的特点外,还有1.稀疏性,是指总刚矩阵的绝大多数元素都是零,非零子块只占一小部分。2.带状性,是指总刚矩阵中非零子块集中在主对角线两侧,呈带状分布。(附,半带宽 B=(相关节点号最大差值+1)*节点自由度数)二 建模与结果分析1. 影响有限元分析精度和成本的因素影响有限元
13、解的误差:1)离散误差 2)位移函数误差分析精度:A、单元阶次 B、单元数量 C、划分形状规则的单元D、建立与实际相符的边界条件 E、减小模型规模 F、避免出现“病态”方程组,当总刚矩阵元素中各行或各列的值相差较大时,则总刚近似奇异。2. 有限元模型的基本构成(节点数据、单元数据、边界条件等)节点数据:节点编号、坐标值、坐标参考系代码、位移参考系代码、节点数量、单元编号单元数据:单元节点、编号单元、材料特性码、单元物理特性值码、单元截面特性、相关几何数据边界条件数据:位移约束数据、载荷条件数据、热边界条件数据、其他边界条件数据归纳起来,网格划分生成节点和单元的过程主要包括定义单元属性、定义网格
14、生成控制和生成网格三个步骤。3. 有限元建模的常用方法理解及应用(如细节处理、分步计算、局部计算、子结构法、对称性简化等)细节处理也称为小特征处理,即删除或抑制对结构力学性能影响不大的细小结构。分步计算,如果结构的局部存在相对尺寸非常小的细节,且又不能进行细节处理,可采用分步计算来控制有限元模型的规模。局部处理就是从所建立的力学模型中抽取一部分出来进行分析,该部分通常是研究者最关心的的危险区域。子结构法是先将大型结构分解为若干个结构区域,每个区域作为一个子结构。子结构被进一步细分为单元,并人为地将子结构上的节点划分为边界节点和内部节点两类.对称性简化,对称性分为反射对称和周期对称(1)反射对称
15、,受对称载荷作用则对称面上的位移条件为垂直于对称面的移动位移分量为零。绕平行于对称面的两相互垂直的轴的转动位移分量均为零。(2)反射对称,受反对称载荷作用则对称面上的位移条件为平行于对称面的移动位移分量为零;绕方向矢量垂直于对称面的轴的转动位移分量为零。(3)对称结构受任意载荷作用(迭加原理)(4)周期对称的位移条件,周期对称边界上的对应点有相同的位移状态4. 边界约束条件的处理(见前) 。5. 单元类型选择的一般原则选择原则:同一问题所选单元应使计算精度高、收敛速度快、计算量小。1、杆系结构:a 、铰接连接时,选杆单元; b、刚性连接时,选刚架单元。2、平面结构:a 、外载平行于平面内,选平
16、面单元 b、外载不在平面内,选弯曲板壳单元3、空间结构:a 、结构和受力具有轴对称性,选轴对称单元;b、一般实体,选三维实体单元。三其他1. 有限单元法首先求出的解是位移,单元应力、应变可由它求得。2.模型简化力学模型简化就是借助对称原理、圣维南原理、叠加原理、等效准则等,在保证考察问题的精度不受影响或具有足够近似精度的要求下,对所研究问题的力学模型进行简化,从而达到减小计算时间和存储空间的目的。常用方法有降维处理、等效变换、对称性利用、组合分析、场耦合分析。3.对称性利用注意事项(1)不仅结构的形状、载荷要对称,还要求位移约束也对称(2)若对称面上作用有载荷,则应取载荷的 1/2 进行分析(
17、3)若对称面上存在板或梁,这时板或梁单元的刚度应取整个单元刚度的 1/2(4)用对称面剖分结构时,应尽量使剖分面不在结构的最大应力位置。其他略。4.整体刚度方程中的约束条件引入引入约束的方法常有:1)降阶法(打乱KR的储存顺序)2)对角元素置 1 法 3)对角元素乘大数法四 ANSYS 操作部分1 ANSYS 有两种模式:一种是 图形交互模式,另一个是 命令流模式。2 ANSYS 软件的基本分析步骤ANSYS 软件基本分析步骤大致可分为三个模块:前处理、加载求解和后处理前处理(1)建立有限元模型所需输入的资料。如节点、坐标资料、元素内节点排列次序。 (2)材料属性 (3)模型的单元划分求解处理
18、:(1)负载条件 (2)边界条件及求解后处理:POST1 用于静态结构分析、屈曲分析及模态分析,将解题部分所得的解答如变形、应力、反力等资料,通过图形接口以及各种不同表示方式把等位移图、等应力图等显示出来。POST26 仅用于动态结构分析,用于与时间相关的时域处理。3 Ansys前处理问题启动ANSYS进入ANSYS交互界面环境,最初的默认激活坐标系(当前坐标系)总是总体直角坐标系。ANSYS 中实体建模的方式有:(1)简单模型,在 ANSYS 中创建几何模型,分自底向上和自顶向下的建模方式;(2)对于复杂模型,导入在其他 CAD 系统创建的模型;(3)直接生成。ANSYS 两种创建有限元模型
19、的方法:(1) 直接法:创建节点,由节点创建单元,适用于简单的结构系统杆梁。(2) 间接法:先几何模型,再网格划分 ,适用于几何外形结构复杂的系统。布尔运算就是对生成的实体模型进行诸如交、并、减等的逻辑运算,给用户快速生成复杂实体模型提供了极大的方便。单元与文件数据类型:在 ANSYS 的单元库中,SOLID45 单元属于结构实体单元,beam188 是 3D 梁单。利用ANSYS进行结构分析时, 载荷步文件为Jobname.sn,结果文件为 jobname.rst 。约束: 对称与反对称约束对称约束限制对称面内所有节点的两个方向旋转自由度,同时限制了垂直于对称面的位移自由度。反对称约束限制了
20、对称面内所有节点在对称面内的两个方向位移自由度,同时限制了垂直于对称面的旋转自由度 。载荷问题:Ansys中的载荷(Loads)包括所有边界条件 以及外部或内部的作用力。具体包括:DOF ,集中载荷 Force ,面载荷 Surface load ,体积载荷 Body load ,惯性载荷 Inertia load,耦合场载荷 Coupled-field loads 。ANSYS划分有限元网格的方法及其区别Ansys 提供两种网格划分方法:自由网格划分和映射网格划分。所谓自由网格划分,体现在没有特定的准则,对于单元形状无限制,生成的单元不规则,基本适用于所有的模型,自由网格生成的内部节点位置比较随意,用户无法控制。映射网格划分则要求满足一定的规则,生成的单元形状比较规则,适用于形状规则的体或面。且映射面网格只包含四边形或三角形单元,而映射体网格只包含六面体单元。Ansys单位制问题Ansys坐标问题Ansys定义材料属性方法4 ANSYS 后处理问题ANSYS中的两个后处理分别是 通用后处理 和 时间历程后处理 。