有限元复习资料

1、何谓有限元法，其基本原理是什么？有限元法是求解一类偏微分方程和常微分方程的数值分析方法。其基本原理是：1）将连续的求解区 域离散成有限多个单元，在每个单元中设定有限多个节点，将连续体看作是只在节点处相连接的一 组单元的集合体；2）选定场函数的节点值作为基本未知量，单元内的场函数值通过节点处的值用假 定的插值函数表示；3）通过和原问题数学模型等效的变分原理或加权余量法，建立求解基本未知量 的代数方程组或常微分方程组，从而将一个连续域中的无限自由度问题化为离散域中的有限自由度 问题。求解后利用解得的节点值和设定的插值函数确定单元内部以至整个集合体上的场函数。何谓单元、载荷、边界条件？整体结构离散中

2、划分出的每一小块体称为单元。常见的单元类型有线性单元、三角形单元、四边形 单元、四面体单元和六面体单元。单元顶点及插值点称为节点。求解区域所受外界施加的力称为载荷，包括集中载荷和分布载荷等。在不同学科中，载荷的含义有 所不同。边界条件是指求解区域边界上所受到的外加约束。如结构分析中的位移约束，热分析中边界上给定 的温度、热流等条件。写出有限元法列式的一般形式，说明各项（包括元素）的含义。应用有限元法，求解问题的数学模型被转化为一组线性代数方程组，用矩阵表示为K u F，其中 K称为刚度矩阵，u为待求节点场函数向量，F为节点载荷向量。总体刚度矩阵由单元刚度矩阵叠加而成。以结构分析为例，单元刚度矩

3、阵中元素Ke的含义是单元节点第j自由度产生单位位移引起的 j第i自由度方向的力，总体刚度矩阵中元素K.的含义是节点整体编号中各单元在同一节点第j自由 IJ度产生单位位移引起的第i自由度方向力的和，也就是有公共节点的各相关单元对这一节点都有“贡 献”。说明有限元法与有限差分法的异同并列举其它的数值求解方法。相同之处：1）两者都是求解偏微分方程和常微分方程的数值分析法；2）都要对求解区域进行离散；3）最终的形式都是一组以节点场函数值为未知变量的线性代数方程组。不同之处：1）形成以节点场函数为未知量的线性代数方程组的理论不同。有限元法是通过变分原理或加权余量 法等方法实现；而差分法是将控制方程直接应

4、用在节点上，用节点处的差商代替导数，得到以节点 处场函数值为未知量的线性代数方程组，结合给定的边界条件进行求解；2）区域离散方法不同。有限差分法用与坐标线（面）平行的线（面）划分区域，对复杂形状区域较 难处理；有限元法划分区域的单元形状灵活，适于处理复杂形状区域；其它的数值求解方法：如有限容积法，有限分析法。ANSYS软件主要有那几部分组成，主要能解决哪些工程领域的问题？ANSYS软件主要由前处理模块、分析计算模块和后处理模块组成。ANSYS软件是大型多物理场分析求解软件，主要应用领域：机械、化工、铁路、电子、航空航天、 能源动力等。主要功能模块包括：1）结构静力分析；2）结构动力学分析；3）

5、结构非线性分析；4） 运动学分析；5）热分析；6）电磁场分析；7）流体动力学分析；8）声场分析；9）压电分析；10） 偶合场分析；11）高级分析模块，如优化设计等。应用ANSYS软件求解物理问题的一般步骤有哪些？1）建立模型指定文件名和分析问题标题；定义单位制；定义单元类型；定义单元实常数；定义材 料特性；创建几何模型；划分网格生成物理模型。2）加载和求解 定义分析问题类型和分析选项；施加约束和载荷；指定载荷步选项；求解。3）查看结果并进行分析应用通用后处理器POST1和时间-历程后处理器检查和分析结果ANSYS软件结构分析模块主要有那几类单元，各有哪些特点？主要有：集中质量单元，杆单元，梁单

6、元，连续实体单元，壳单元和组合单元等。这些单元均有二 维和三维之分。另外，质量单元，杆单元，梁单元，壳单元，组合单元除有材料属性外还有几何或 其他属性，如质量，梁的截面积，惯性矩等，这些参数须在材料实常数中定义。在ANSYS软件中划分网格应注意哪些问题？1）对由多种材料，多个零件组成的结构，对不同零件指定不同网格属性：如通过单元类型，材料类 型，材料实常数等区分不同的区域，所划分网格与具体的零件（区域）相对应。2）网格疏密与所分析物理问题，网格所处位置有关，在被求量变化速率大的地方应作加密处理；3） 不同材料区域单元网格尺寸差异不应太大。用ANSYS软件建立几何模型通常有哪两种方法，具体内容是

7、什么。说明图1所示轴承座结构几 何建模要点。图1轴承座答：有两种方法。一是自下向上的建模方法，该方法 利用点、线、面基本图形要素建立模型，二是自上向下的建 模方法，该方法利用ANSYS软件提供的基本几何形体，通过 逻辑关系建立几何模型。图1所示轴承座几何建模要点： 轴承座关于过孔轴线的垂直面对称，建模时先建立一半的模型， 利用镜像生成另一半模型，最后通过合并生成轴承座整体。轴 承座底座为长方体，孔上部为1/4圆柱，下部为长方体，孔通 过1/4圆柱加下部长方体减去中间的圆柱生成，筋通过菱形基 本图形生成。因此，轴承座可利用自上而下的建模方法创建。建立圆柱要根据需要旋转局部坐标系， 使z轴沿着圆柱

8、轴线方向。拉筋厚度方向也要沿局部坐标的z轴方向。10.在ANSYS软件中命令“ Add”、“Glue”、“Overlap”的作用是什么，通常在什么情况下使用？答：“Add”、“Glue”、“Overlap”命令是几何建模中的逻辑命令，“Add”将几个相关联的、独立建立 的模型合并成一个几何体，在生成一个完整实体时使用。“Glue”将几个相关联的、独立建立的模 型粘接起来，彼此还是独立的几何体，只是在界面上连接起来，用于把不同性质的几何体关联起来， 模型划分网格后，界面上的节点具有相同的位移或温度。“Overlap ”命令将重叠的几何体在界面处分割，但界面仍粘在一起。“Overlap”命令具有分

9、割和粘接的功能。二、图2所示为一受轴向拉力作用的阶梯轴，试用最小势能原理建立该阶梯轴的有限元列式。解：小势能原理可叙述为：在外力和约束作 弹性体，其总势能为最小，或者说，在所有 实位移场的总势能取最用下处于平衡状态的可能的位移场中，真小值。数学描述即总势能的一阶变分为零，而且二阶变分是正定的（大于零）。对阶梯轴：兀=J b dv - F u2外力作用下的真实位移使:兀是该问题的泛函。位移与应变的关系为:dudx应力与应变的关系为:（本构关系）设位移模式：u =a 1 +a2xx = 0, u = u 对第单元：x = l , u = u得：ax1-厂1 e J(11=Bu ；,jJ,jJ对阶梯

10、轴划分单元:V 1 V _兀=乙 J ns dv 一乙 F u ；2对典型单元:兀 e = J !b dv -E Fu 2le JAdx - (u , u )BuAdx 一 uF=j 2(u ,u .)|f 1E-一1e! Adx-(u - u)J真实位移使总势能取极小值（8兀=0 ），对单元也有（8兀e = 0 ）。得，8兀e=0,=8 uj11212ee111212ee1_11212ee110。d u、00, j EA0, j EAf 0uj应使王ui dx -(1,0 ) =dxEA整体刚度阵的装配:扩展单元刚度阵，k(1)k(1)01112k(1)k(1)02122000k(总体刚度阵

11、，k = k+ k（2）整体有限元列式：kx5= F ；=(u,u , u )T , F = (F ,0, F )T。12313三、图3所示为一双层平壁导热问题，试用工程法建立求解该双层平壁温度场的有限元列式。若T = 40 C， T = 15 C， X = 0.35 W/ (m.K) , L = 25mm ,人=0.2 W/ （m.K）, L2 = 5mm ,用建立的有限元列式求平壁1 一半厚度处的温度和 通过平壁传递的热量。解：节点热流与单元节点温度的关系T - T 人A 人A 节点处：。=人A tj = T - fT ; TOC o 1-5 h z ， e e l l i l ， eee

12、T -T 入 A _ 入 A _节点处：Q =人 A T-=-T +T ;图 3j e e ll i l ，入Ale入Ae4IQ 1。用矩阵表示为：，I =Q.JJ单元有限元列式的集结：整体热导矩阵可通过单元热导矩阵扩展或节点热流平衡方程求出。 对节点1列热流平衡方程：le入Ale ；或者：Q e = keTe，k为热导矩阵图。3Qw +Qi +Qi = 0 ；112T（】）=T , T（2）= Ti1 i2 ； T（】）=T , T（2）= T节点 1：k（】）T + k（i）T. =4 ；节点2：k（1）T（1）+ k（1）T（1）+ k（2）T 俱）+ k（2）T（2）= Q ；21 i

13、22 J11 i12 J2k（2）T（2）+ k（2）T（2）=21 i 用矩阵表示为：节点3：22 JI| = 或，-Q1wQ w |-Q w 3节点上热流:Q w未知，1Q w =0，2边界处：x = 0,T = T，x = l + l ,T = Tn12w用置1法求未知节点的温度:k(1)210Vk + k(2)22110k(2)121(T T2TV 3(T(X X )产n + : Tw )V 11TV w12X )+ 了 A,(A = A = A )；热流：如顼12=e w3X A, (T1 n1人A1l1(X/TV1 nX)X X2 T12 A1wJ1 1X :=X 2X21 +21

14、112 JJ12+、(T - T )；n w2 J代入给定数据:(T t2TV 311TV w+ 了Tw七 J2+?2 J(40(0.35)40 +15 1V 0.0250.005J(空+QV 0.025 0.005 J15 u0.2(4021.48e=-(Ti n1,板1(301- 4) = 160W , ik 0.030.3 X 1半厚度处的温度为(Tn+T2) /2=22C。四、用有限元法计算二维平面应力问题。计算区域用三角形网格划分，典型网格如图5所示， 元内位移插值函数采用线性插值函数u (X y) =。+。X + a y。1231)2)3)通过节点位移值求差值函数系数的表达式；导出

15、用单元节点位移表示单元内部位移的计算式;应用虚位移原理导出单元刚度矩阵表达式。x ii(xi, yi)解：1)位移插值函数:;代入节点处坐标和位移值，得u =a*a x + a；u =以+以 x +以 y ； u =a*a x +a；i 12 i 3 i j 12 j 3 j k 12 k 3 k用矩阵表示为：(1xy )a ui11xya =u ,jj21xkyk JaI 3u侦k Ja =一 XTjXk, y求解得a 1如2/akbkuu . 7=12 Aabajbjacccul 3)V ijk JI k)bi1，y.1, y1, y=- j = y y , b，k= j1, y k=yk

16、一 y jb = 1 = y yj1, x1, x1, xci=j = x - x ,，c .=，k一i1, xk=x -x，c = i = x - xj其中：a = x y 一 x y ,i i k k j=x y - x y k i i k,a = x y 一 x yk i j j i2 A= a + a + ai j k1 u (a+ b x + c y ) u +1(a + b x + c y )u +1(a+ b x + c y ) u2 Aiiii2 Ajjjj2 Akkkk2)=N (x, y ) u + N (x, y ) u + N (x, y ) uii jj kk之8 x08 Nr8 x08 =088yu, Lv=08 Ni8 y_8_8 N8 N8 y8 x8 y8 x3)根据虚位移原理，有5 *TF = jj&*T Q tdxdyj08 x8 N0t8 y8 N8 NTj8 y8 xd N8 N1r08 x8 N01r8 y8 N8 N1t1t8 y8 xuiviujvjumv= b b a x 2 A . jbxbyT xyyYxy=DB5 e将以上格式代入虚功方程式，得K = j BTDBtdAA五、耦合场分析常用的方法有哪两种 分析的过程。以图5所示的冷却栅管热应力分析为例说明其中一种耦合