有限元分析理论基础

1、有限元分析理论基础 有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元 （子域）所构成，其模型给出基本方程的分片（子域）近似解，由于 单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很 好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型： 它是真实系统理想化的数学抽象。 由一些简单形状 的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析： 是利用数学近似的方法对真实物理系统 （几何和载 荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可 以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的

2、， 所考虑的变形建立 在小变形假设的基础上。 在这类问题中， 材料的应力与应变呈线性关 系，满足 xx 胡克定律；应力与应变也是线性关系，线弹性问题可归 结为求解线性方程问题， 所以只需要较少的计算时间。 如果采用高效 的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析 两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类： 1）材料非线性问题 材料的应力和应变是非线性的，但应力与应

3、变却很微小，此 时应变与位移呈线性关系，这类问题属于材料的非线性问题。由 于从理论上还不能提供能普遍接受的本构关系，所以，一般材料 的应力与应变之间的非线性关系要基于试验数据，有时非线性材 料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。 在工程实际中较为重要的材料非线性问题有：非线性弹性（包括 分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研 究这类问题一般都是假定材料的应力和应变呈线性关系。它包括 大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于 大位移小应

4、变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、 xx 、撞击等问题中，接触和摩擦的作用不可忽视， 接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成 型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外 部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。 有限元理论基础 有限元方法的基础是变分原理和加权余量法， 其基本求解思想是 把计算域划分为有限个互不重叠的单元， 在每个单元内， 选择一些合 适的节点作为求解函数的插值点， 将微分方程中的变量改写成由各变 量或其导数的节点值与所选用的插值函数组

5、成的线性表达式， 借助于 变分原理或加权余量法， 将微分方程离散求解。 采用不同的权函数和 插值函数形式，便构成不同的有限元方法。 1. 加权余量法： 是指采用使余量的加权函数为零求得微分方程近似解的方法称 为加权余量法。 (Weighted residual method WRM )是一种直接从所 需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方 法。加权余量法是求解微分方程近似解的一种有效的方法。 设问题的控制微分方程为 在V域内 在S边界上 式中： L、B分别为微分方程和边界条件中的微分算子； f、g 为与未知函数u无关的已知函数域值； u为问题待求的未知函数 出干！J用加枳余

6、垃-法求用1丛鮮日寸，昔先在求辭域上吐立一于试函皱；H , 一般兵冇如下形式: n = NO(5.1 31 z 式中： g彳寺比系劉j 竝可譯尔为广义坐标； A;职白完侖函故耳的纟钏生无关的基雷数。 由于円一舟丈只绘彳寺衣函劉hi的近1以为军. 因 m旁式Q 1 3) 代入式(A I I冲口式G I 4后将得不到勺垢足. 若iG : | R =(刃$- J在 或 Rb = B(R、-g在百边界上(5 握踽马 、尽反映-丁彳戎朋i爭乂与H实稱之甬1的偏釜，它TH分另M尔 T故内那利边界余量、 若在域：内弓入内部权函数町，在边界S上引入边界权函数 W8 则可址立n个消除余量的条1牛. _般可农示为

7、： L肥R/7F+ L%足於=0(/= L2丄 a0丄为 * v* S 不同的权函数；和打 反映了不同的消除余量的准则:，从上 式可以得到求解待定 系数矩阵C的代数方程组o 经解衢侍定 系数.由式(5,1,3)即可得所需求解边血问趣的近似解 由于试函数的不同，氽鸯鸟和 可有如下三种怖况， 依此加权余於法可分为： 1 +内郡法 试函苗攵満足边界条件，也即 尺=e娜-小 从匕甘寸消除余進的条1牛戚之r WhB.V 0 (i =1,21 )(5 1.6) IB F 2. 边挤法 试两敖满兄控制方程，施即 R=L(-/=a 此时消除余塚的条件为： WSiRS = O (/-L2,L j?) * s 3

8、. 混合法 试函数不满足控制方程在边界条T牛. 此时用 式(5 L5)来消除氽量:J 混合法对于试函数的选取最方便， 但在相同精度条件下，工作量 最大。对内部法和边界法必须使基函数事先满足一定条件，这对复杂 结构分析往往有一定困难，但试函数一经建立，其工作量较小。 无论采用何种方法，在建立试函数时均应注意以下几点： (1)试函数应由完备函数集的子集构成。已被采用过的试函数有 幕级数、三角级数、样条函数、贝赛尔函数、切XX和勒让XX多项式 等等。 试函数应具有直到比消除余量的加权积分表达式中最高阶导 数低一阶的导数连续性。 (3)试函数应与问题的解析解或问题的特解相关联。若计算问题 具有对称性，

9、应充分利用它 显然，任何独立的完全函数集都可以作为权函数。按照对权函数 的不同选择得到不同的加权余量计算方法， 主要有：配点法、子域法、 最小二乘法、力矩法和XXXX法。其中XXXX法的精度最高。 下面以內岂肉去为例, 介绍 抉权函攻攵分类时力廿杈余量的五种 基本 方法:，对内部法来说，消除余壁的统一格式是： 3肿T 二 Q 二 12L a 1, 子域法（Snbdomam Metliotl） 此法首先将求解域划分成n个子域 匚+ 在每个子 域内令权函数 寻于1” 而在子域之夕卜职权函数为马 也即： w =|1厲内） _10 （於M 如巣在各个子域里分别选取试函数.那么 它的求解在形式上将 类似

10、于 冇限元法口 2 配点法（Collocation Method； 子域法兄令余辽在一个子域上的总和为 知 而否己点法兄T丸 余呈在扌旨定的11个点上寻于寺, 这些点滁为配点Q此法的权 函数为： WSiP-P) :瓷皿驭（犹拉克）两数，它的定义为: I 0 X 工 A. 曲 7)=1 5 X = X L* 4(5 - ix = IJ- %叮1号迂c P、P分别代农求解域内任一点霹口西已点。 由于此法只在配点上保证余蚤为冬.因此不需要作积分吃十算. 所以心 3简单的力口杈余垃法 3*姒小二乘法(L亡鬧Square Method) 本法逓过使在整个求解域上余甘的平方和取极小来建立消除余 鈕的条1牛

11、. 若记余蚤平右和为1(C),即J(C)= I* dv= f R：丽卩 则极值条件为：叮八雰“ = 由此可见”本法权函效为、% =学 O=1JX ,M) 亠 伽辽金QgxkiU Method) 本法是使余量与毎一个基函数正交. 也艮P以基函数作为权函数 肌=N (z = 12X “) 斗戎函数 涉包含址个克备函:心厂h亿 用本法必可求彳号粉确 5.矩法(Method of Moment ) 本法与伽辽金法相似，也兢用完备函数弘作权雨数门 但本去的权函数与伽辽金去又冇区别, 它与试函 扳无关。 消除余呈的条T牛是从如开始的各阶矩为冬 因u匕 对一九问趣 阵=0 (J =L2. L 出) 对;细问

12、迥，叫=xy (jj=1,2, l ,塚) 其余类推 这五种基本方2去在待定 系芋丈足鎖多以尔做高阶近T以)日寸. 其榆 度彼此相近。但对低阶近1以(11狡小)t/r况下” 后三种的桥度 要高于前两种_. 基本方法举例 为说明上述基本槪念.以图所示寻截面忿背粱.燮满跨均布荷 :我作用，求悬臂端B的竖向位移亠为例. 说明基本方法的应用。 若取试函数为：C(A5 +7a4 - L4Z2.? + 26/3a3 ) (*) 图示二育的扌空制右程为 咚十。 dx 勝普。 (龙=0) ( = 0 其边弊条件为 不堰验证其濒足卫生界条T牛. 牝即 碍=% 而控制方程的內那余41? R为: Rt = EIc(

13、120 xr241)-(j 因此本问题腿内解法,下面分别用基本方三去进行衣解， 子域法解 由于试函数仅一个待定常数，因此只需职一个子域(昔于全域)即 可.消除氽莹白勺条件为: |： EIc (120a 十 24/ ) - 巒=0 由此可解得：“ q S4/ 代回(#)式可得: 配点法解 同上所连，只需立盘一个配点来建立消除余壁的条件o 若令: Rl x=0.75/ 可得: CAEIl 腹)唧 J 51EI 若令： jo 则得 C- 1445/ (*水) 可贝不冋的配点结果丸刁、一扌斗的Q 此时消 除余量 的 条1牛为: /Jci：120 x+24/ 一？卜国门2加亠 247 寸必二 0.010

14、17 可得： 伽辽金法解 此时. A4 =x5 + /x4-14/2x3+26/3x2 涓除余呈的条件为： JO NRdx = Q 由此可谆2=誉 = 0.1262g/4 El 矩法解 由于只冇一个讨定沿数”因此消除余辽条件只卅零次矩即可, U匕呂寸显然与子域三去完全相同 本例各方法的楙度比校 本问遐的楙研解由粱位移计坪可叫为 ql4 _ 0.1257 SEI EI 由此可得，上述各方法对本例计算的误冬侬次为 -33.3%; 1.75% (22.2%) : 13.9%; 096%; -33.3% 上而22 2%为式(村)结果“ 2、虚功原理 平衡方程和几何方程的等效积分“弱”形式 虚功原理包含

15、虚位移原理和虚应力原理，是虚位移原理和虚应力原 理的总称。他们都可以认为是与某些控制方程相等效的积分“弱”形 式。虚功原理：变形体中任意满足平衡的力系在任意满足协调条件的变 形状态上作的虚功等于零，即体系外力的虚功与内力的虚功之和等于 J | A 零。 虚位移原理是平衡方程和力的边界条件的等效积分的“弱”形式； 虚应力原理是几何方程和位移边界条件的等效积分“弱”形式。 虚位移原理的力学意义：如果力系是平衡的，则它们在虚位移和虚 应变上所作的功的总和为零。反之，如果力系在虚位移（及虚应变）上 所作的功的和等于零，则它们一定满足平衡方程。所以，虚位移原理表 述了力系平衡的必要而充分条件。 一般而言

16、，虚位移原理不仅可以适用 于线弹性问题，而且可以用于非线性弹性及弹塑性等非线性问题。 虚应力原理的力学意义：如果位移是协调的，则虚应力和虚边界约 束反力在他们上面所作的功的总和为零。 反之，如果上述虚力系在他们 上面所作的功的和为零，则它们一定是满足协调的。所以，虚应力原理 表述了位移协调的必要而充分条件。 虚应力原理可以应用于线弹性以及非线性弹性等不同的力学问题。 但是必须指出，无论是虚位移原理还是虚应力原理，他们所依赖的几何 方程和平衡方程都是基于小变形理论的，他们不能直接应用于基于大变 形理论的力学问题。 3、最小总势能法 应变能：作用在物体上的外载荷会引起物体变形， 变形期间外力 所做

17、的功以弹性能的形式储存在物体中，即为应变能。 由n个单元和m个节点组成的物体的总势能为总应变能和外力所 做功的差： nm =上(e) _ 7 FiU i e 4i总 最小势能原理：对于一个稳定的系统，相对于平衡位置发生的位 移总会使系统的总势能最小，即： ，匸1,2,3, ,n 有限元法的收敛性 有限元法是一种数值分析方法，因此应考虑收敛性问题。 有限元法的收敛性是指：当网格逐渐 xx时，有限元解答的序列 收敛到精确解；或者当单元尺寸固定时，每个单元的自由度数越多， 有限元的解答就越趋近于精确解。 有限元的收敛条件包括如下四个方面： 1) 单元内，位移函数必须连续。多项式是单值连续函数，因此

18、选择多项式作为位移函数，在单元内的连续性能够保证 2）在单元内，位移函数必须包括常应变项。每个单元的应变状 态总可以分解为不依赖于单元内各点位置的常应变和由各点位置决 定的变量应变。当单元的尺寸足够小时， 单元中各点的应变趋于相等， 单元的变形比较均匀， 因而常应变就成为应变的主要部分。 为反映单 元的应变状态，单元位移函数必须包括常应变项。 3）在单元内，位移函数必须包括刚体位移项。一般情况下，单 元内任一点的位移包括形变位移和刚体位移两部分。 形变位移与物体 形状及体积的改变相联系， 因而产生应变；刚体位移只改变物体位置， 不改变物体的形状和体积， 即刚体位移是不产生变形的位移。 空间一

19、个物体包括三个平动位移和三个转动位移，共有六个刚体位移分量。 由于一个单元牵连在另一些单元上， 其他单元发生变形时必将带 动单元做刚体位移，由此可见，为模拟一个单元的真实位移，假定的 单元位移函数必须包括刚体位移项。 4）位移函数在相邻单元的公共边界上必须协调。对一般单元而 言，协调性是指相邻单元在公共节点处有相同的位移， 而且沿单元边 界也有相同的位移， 也就是说， 要保证不发生单元的相互脱离开裂和 相互侵入重叠。 要做到这一点， 就要求函数在公共边界上能由公共节 点的函数值唯一确定。 对一般单元， 协调性保证了相邻单元边界位移 的连续性。 但是，在板壳的相邻单元之间，还要求位移的一阶导数连

20、续，只 有这样，才能保证结构的应变能是有界量 总的说来，协调性是指在相邻单元的公共边界上满足连续性条 件。 前三条又叫完备性条件，满足完备条件的单元叫完备单元；第四 条是协调性要求，满足协调性的单元叫协调单元；否则称为非协调单 元。完备性要求是收敛的必要条件，四条全部满足，构成收敛的充分 必要条件。 在实际应用中，要使选择的位移函数全部满足完备性和协调性要 求是比较困难的，在某些情况下可以放松对协调性的要求。 需要指出的是，有时非协调单元比与它对应的协调单元还要好， 其原因在于近似解的性质。假定位移函数就相当于给单元施加了约束 条件，使单元变形服从所加约束，这样的替代结构比真实结构更刚一 些。

21、但是，这种近似结构由于允许单元分离、重叠，使单元的刚度变 软了，或者形成了（例如板单元在单元之间的绕度连续，而转角不连 续时，刚节点变为铰接点）对于非协调单元，上述两种影响有误差相 消的可能，因此利用非协调单元有时也会得到很好的结果。在工程实 践中，非协调元必须通过“小片试验后”才能使用。 应力的单元平均或节点平均处理方法 最简单的处理应力结果的方法是取相邻单元或围绕节点各单元 应力的平均值 ? 1. 取相邻单元应力的平均值 这种方法最常用于 3 节点三角形单元中。这种最简单而又相当 实用的单元得到的应力解在单元内是常数。可以将其看作是单元 内应力的平均值，或是单元形心处的应力。由于应力近似解

22、总是 在精确解上下振荡，可以取相邻单元应力的平均值作为此两个单 元合成的较大四边形单元形心处的应力。 如 2 单元的情况下，取平均应力可以采用算术平均， 即平均应力 =（单元 1的应力+单元 2的应力） /2。 也可以采用精确一些的面积加权平均， 即平均应力=单元1应力X单元1的面积+ 单元2应力x单 元 2面积/ （单元 1 面积+单元 2面积） 当相邻两单元面积相差不大时，两者的结果基本相同。在单元 划分时应避免相邻两单元的面积相差太多，从而使求解的误差相 近。 一般而言， 3 节点三角形单元的最佳应力点是单元的中心点， 此点的应力具有 1 阶的精度。 ? 2. 取围绕节点各单元应力的平均

23、值 首先计算围绕该节点（ i ）周围的相关单元在该节点出的应力 值，然后以他们的平均值作为该节点的最后应力值 ，即 其中，1m围绕在i节点周围的全部单元。取平均值时也可进 行面积加权。 有限元法求解问题的基本步骤 1. 结构离散化 对整个结构进行离散化，将其分割成若干个单元，单元间彼此通 过节点相连； 2. 求出各单元的刚度矩阵 K(e) K(e)是由单元节点位移量(e)求单元节点力向量F(e)的 转移矩阵，其关系式为:F(e)= K(e) (e) 3. 集成总体刚度矩阵 K 并写出总体平衡方程： 总体刚度矩阵K是由整体节点位移向量求整体节点力向量 的转移矩阵，其关系式为F= K ，此即为总体

24、平衡方程。 4. 引入支撑条件，求出各节点的位移 节点的支撑条件有两种：一种是节点 n 沿某个方向的位移为零， 另一种是节点 n 沿某个方向的位移为一给定值。 5. 求出各单元内的应力和应变。 对于有限元方法，其基本思路和解题步骤可归纳为 (1) 建立积分方程，根据变分原理或方程余量与权函数正交化原 理，建立与微分方程初边值问题等价的积分表达式， 这是有限元法的 出发点。 (2) 区域单元剖分， 根据求解区域的形状及实际问题的物理特点， 将区域剖分为若干相互连接、 不重叠的单元。 区域单元划分是采用有 限元方法的前期准备工作， 这部分工作量比较大， 除了给计算单元和 节点进行编号和确定相互之间

25、的关系之外，还要表示节点的位置坐 标，同时还需要列出自然边界和本质边界的节点序号和相应的边界 值。 (3) 确定单元基函数，根据单元中节点数目及对近似解精度的要 求，选择满足一定插值条件的插值函数作为单元基函数。 有限元方法 中的基函数是在单元中选取的，由于各单元 具有规则的几何形状， 在选取基函数时可遵循一定的法则。 (4) 单元分析：将各个单元中的求解函数用单元基函数的线性组 合表达式进行逼近；再将 近似函数代入积分方程，并对单元区域进 行积分，可获得含有待定系数 (即单元中各节点 的参数值 )的代数方 程组，称为单元有限元方程。 (5) 总体合成：在得出单元有限元方程之后，将区域中所有单

26、元 有限元方程按一定法则进行累加，形成总体有限元方程。 (6) 边界条件的处理：一般边界条件有三种形式，分为本质边界 条件(狄里克雷边界条件 ) 、自然边界条件 (xx 边界条件 )、混合边 界条件(XX边界条件)。对于自然边界条件， 一般在积分表达式中可 自动得到满足。 对于本质边界条件和混合边界条件， 需按一定法 则 对总体有限元方程进行 XX 满足。 (7) 解有限元方程： 根据边界条件 XX 的总体有限元方程组， 是含 所有待定未知量的封闭 方程组，采用适当的数值计算方法求解，可 求得各节点的函数值。 单元刚度矩阵的特性 单元刚度矩阵无论在局部坐标系中还是在整体坐标系中都具有相 同的三

27、个特性： 1) 对称性 由材料力学中的位移互等定理可知，对一个构件，作用在点 j 的 力引起点 i 的绕度等于有同样大小而作用于点 i 的力引起的点 j 的绕 度，即 kij(e) = kji(e) ，表明单元刚度矩阵是一个对称矩阵。 2) 奇异性 无逆阵的矩阵就叫做奇异矩阵，其行列式的值为 0，即 |k(e)|=0 ， 这一点可以从例题直接得到验证。其物理意义是引入支撑条件之前， 单元可平移。 3) 分块性 有前面所讲的内容可以看出，矩阵k(e)可以用虚线分成四块, 因此可写成如下的分块形式， “1耐阿广 式中kmn(e)局部坐标系中单元(e)按局部码标记的节点 m n 之间的刚度子矩阵 刚

28、架结构xx节点载荷的处理的方法 在刚架结构以及其他较复杂的结构上，他们所受的载荷可以直接 作用在节点上，又可以不直接作用在节点上而作用于单元节点间的其 他位置上。后一种情况下的载荷称为非节点载荷。有限元分析时，总 体刚度方程中所用到的力向量是节点力向量。因此在进行整体分 析前应当进行载荷的移植，将作用于单元上的力移植到节点上。移植 时按静力等效的原则进行。 处理非节点载荷一般可直接在整体坐标系内进行，其过程为： 1) 将各杆单元看成一根两端固定的 xx，分别求出两个固定端的 约束反力。其结果可直接利用材料力学的公式求得； 2) 将各固定端的约束反力变号，按节点进行集成，获得各节点的 等效载荷

29、总体刚度矩阵的集成法 使用刚度矩阵获得的方法获得总体刚度矩阵。在此将其扩展到由 整体坐标系中的单元刚度矩阵的子矩阵集成总体刚度矩阵。步骤如 下： 1）对一个有n个节点的结构，将总体刚度矩阵K划分为nxn各 子区间，然后按节点总码的顺序进行编号； 2）将整体坐标系中单元刚度矩阵的各子矩阵根据其下标的两个总 码对号入座，写在总体刚度矩阵相应的子区间； 3）同一子区间内的子矩阵相加， 成为总体刚度矩阵中的相应的子 矩阵。 总体刚度矩阵的特性 1 ）对称性：因为由此特性，在计算机中只需存储其上三角部分； 2）奇异性：物理意义仍为在无约束的情况下，整个结构可做刚体 运动； 3）稀疏性： K 中有许多零子

30、矩阵，而且在非零子矩阵中还有大 量的零元素， 这种矩阵称为稀疏矩阵。 大型结构的总体刚度矩阵一般 都是稀疏矩阵； 4）分块性： 平面问题离散化时的规定 1）单元之间只在节点处相连； 2）所有的节点都为铰接点； 3）单元之间的力通过节点传递； 4）外载荷都要移植到节点上； 5）在节点位移或某一分量可以不计之处， 就必须在该节点安置一 个铰支座或相应的连杆支座。 通过以上的规定来建立平面有限元分析模型。 结构对称性的利用规律 一般来说，作用在对称结构上的载荷系统分为对称的、反对称的 和一般的三种情况。 1. 结构对称，载荷对称或反对称 这种情况下，对称面上的边界条件可按以下规则确定： A. 在不同

31、的对称面上， 将位移分量区分为对称分量和反对称分量； B. 将载荷也按不同的对称面分别区分为对称分量和反对称分量； C. 对于同一个对称面，如载荷是对称的，则对称面上位移的反对 称分量为零，如载荷是反对称的，则对称面上位移的对称分量为零。 如果所分析的结构对称，但载荷是不对称的，也不是反对称的， 这时可以将这种结构系统简化成载荷为对称和 /或反对称情况的组 合，仍可以简化分析过程，提高分析的综合效率。 如图a所示，结构对称，载荷一般，可将其载荷分解为图 b和图 c的组合。图b为对称结构，载荷对x、y轴均为对称，图c为结构 对称，载荷对x轴反对称、对y轴对称，此时可取相同的四分之一进 行研究，分

32、别施加对称面上节点的边界条件，进行两次分析计算，并 将计算结果迭加起来，即可得到原结构四分之一的解答，进而得出整 个结构的解答。 利用结构的对称性取某一部分建立有限元模型时，往往会产生约 束不足现象。 例如，若取上例中图 c 的四分之一建立有限元时， 根据上述分析， 在两对称面上应加水平放置的滚动铰支座， 因此模型在垂直方向存在 刚体位移。对这种约束不足问题，利用有限元分析时，必须增加附加 约束，以消除模型的刚体位移。在本例中，垂直方向可以用刚度很小 的杆单元或边界弹簧单元连接到模型某节点上， 使得既消除了模型的 刚体位移，又不致于因附加的杆单元或边界弹簧单元刚度太大而影响 结构原有的变形状态

33、。 单元形态的选择原则 单元形态包括单元形状、边中节点的位置、细长比等，在结构离 散化过程中必须合理选择。一般来说，为了保证有限元分析的精度， 必须是单元的形态尽可能的规则。 对于三角形单元，三条边长尽量接近，不应出现大的钝角、大的 边长。这是因为根据误差分析， 应力和位移的误差都和单元的最小内 角的正弦成反比。因而，等边三角形单元的形态最好，它与等腰直角 三角形单元的误差之比为 sin45 :sin60 =1:1.23 。但是为了适应 弹性体边界， 以及单元由小到大逐渐过渡， 不可能是所有的三角形单 元都接近等边三角形。实际上，常常使用等腰直角三角形。 对于矩形单元来说，细长比不宜过大。细长

34、比是指单元最大尺寸 和最小尺寸之比。 最优细长比在很大程度上取决于不同方向上位移梯 度的差别。梯度较大的方向，单元尺寸要小些，梯度小的方向，单元 尺寸可以大一些； 如果各方向上位移梯度大致相同， 则细长比越接近 1，精度越高。有文献推荐，一般情况下，为了得到较好的位移结果， 细长比不应超过 7；为了获得较好的应力结果，细长比不应超过 3。 一般情况下，正方形单元的形态最好。 对于一般的四边形单元应避免过大的边长比，过大的边长比会导 致病态的方程组。 边界条件的确定 确定边界条件是建立有限元模型的重要一环，合理确定有限元模 型的边界条件是成功地进行结构有限元分析的基本要求。 一般情况下，建模对象的边界条件是明确的。根据力学模型的边 界条件可以很容易确定其有限元模型的边界条件。 例如电线杆插入地 基的一端为固定端，桥梁一端为固定铰支座，另一