有限元分析理论基础

有限元分析概念有限元法：把求解区域看作由许多小的在节点处相互连接的单元 （子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单 元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地 适应复杂的几何形状、复杂的材料特性和复杂的边界条件有限元模型：它是真实系统理想化的数学抽象由一些简单形状的 单元组成，单元之间通过节点连接，并承受一定载荷有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷 工况）进行模拟并利用简单而又相互作用的元素，即单元，就可以用 有限数量的未知量去逼近无限未知量的真实系统线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在 小变形假设的基础上在这类问题中，材料的应力与应变呈线性关系， 满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求 解线性方程问题，所以只需要较少的计算时间如果采用高效的代数方 程组求解方法，也有助于降低有限元分析的时间线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两 方面非线性问题与线弹性问题的区别：1） 非线性问题的方程是非线性的，一般需要迭代求解；2） 非线性问题不能采用叠加原理；3） 非线性问题不总有一致解，有时甚至没有解。

有限元求解非线性问题可分为以下三类：1） 材料非线性问题材料的应力和应变是非线性的，但应力与应变却很微小，此时 应变与位移呈线性关系，这类问题属于材料的非线性问题由于从 理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力 与应变之间的非线性关系要基于试验数据，有时非线性材料特性可 用数学模型进行模拟，尽管这些模型总有他们的局限性在工程实 际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、 弹塑性、粘塑性及蠕变等2） 几何非线性问题几何非线性问题是由于位移之间存在非线性关系引起的当物体的位移较大时，应变与位移的关系是非线性关系研究 这类问题一般都是假定材料的应力和应变呈线性关系它包括大位 移大应变及大位移小应变问题如结构的弹性屈曲问题属于大位移 小应变问题，橡胶部件形成过程为大应变问题3） 非线性边界问题在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视， 接触边界属于高度非线性边界平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、 橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界 相接触时通常要考虑非线性边界条件实际的非线性可能同时出现上述两种或三种非线性问题。

有限元理论基础有限元方法的基础是变分原理和加权余量法，其基本求解思想是把 计算域划分为有限个互不重叠的单元，在每个单元，选择一些合适的节 点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导 数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或 加权余量法，将微分方程离散求解采用不同的权函数和插值函数形式， 便构成不同的有限元方法1. 加权余量法:是指采用使余量的加权函数为零求得微分方程近似解的方法称为 加权余量法Weighted residual method WRM)是一种直接从所需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方法加权 余量法是求解微分方程近似解的一种有效的方法设问题的控制微分方程为:在S边界上B(u) - g = 0(5.1.2)式中 :L、B——分别为微分方程和边界条件中的微分算子;f、g ——为与未知函数u无关的已知函数域值;u——为问题待求的未知函数当才U用 才口权余-3弋去习洪薜钮寸._肖先在美桦t或上?虻夷乙 俨试函救 A—般共可如下形式：fl-ozVc^V； = At1 ;-■ I <■7<式中： c 待宝奈敬，也可称为广义坐标；如 币.白 丸■布 普舟:敬％的线|土无关的扯函瓷％由于占 一 般只且卞y求函我E的*tT以雨\ *^3等式（、1.3）代入式国1.1）和式（5.1.2）后书Fi耳不到洒工，者记：| % ~ L（ f^~ f 在\「域内\ Rb ^B（^-g 在险女界上 成E坦然尺 .建戏J奥了・i式函M攵与一点尖宙年兰」御的偏至，它们攻另府象做内WF才口 .'边界余堆°去衣域5 习I入内 都权而1数 仍,我^查界S上引 入边界枳甬数 电 贝9可地立11个消除余好的泰件，一般可表示为；+ ：/上.氏』$ _ 0 打一 L二L』） ⑴今不同的枳函数 井,才口卬命 区映了不同的涓除余苗的准贝儿 从上 式可坎得到求解祜定京数矩阵C的代救方程绍，一经解得祜冗 系数，由式（％ 1.3）即可得■斯需求解显伍向脸的近似解」由于试函数丹 的不后I. 金蚩氏 和 与 可与第下三村情况，依此加枳余主-法可分为：1.内部法试函数满其边界条件. 也即 ^-^-g-o 此时消除余场的条件成为:[膏湿尸-0 （7-1" .0 仁.].曰二边界法试函数满母窒弟K租，也HP 尺=及跖”。

此时消除余-昼的条件为：1町*妙=0 （7 = LNL.m 4.顷3 一混合沽试函我不满足疫制方程和边界条件，此时用式（％ 1.5｝来沸除氽*，.混合法对于试函数的选取最方便，但在相同精度条件下，工作量最 大对部法和边界法必须使基函数事先满足一定条件，这对复杂结构分 析往往有一定困难，但试函数一经建立，其工作量较小无论采用何种方法，在建立试函数时均应注意以下几点：(1) 试函数应由完备函数集的子集构成已被采用过的试函数有幂 级数、三角级数、样条函数、贝赛尔函数、切比雪夫和勒让德多项式等 等2) 试函数应具有直到比消除余量的加权积分表达式中最高阶导数 低一阶的导数连续性3) 试函数应与问题的解析解或问题的特解相关联若计算问题具 有对称性，应充分利用它显然，任何独立的完全函数集都可以作为权函数按照对权函数的 不同选择得到不同的加权余量计算方法，主要有：配点法、子域法、最 小二乘法、力矩法和伽辽金法其中伽辽金法的精度最高下面以内 部法为例, 介绍按权网 救分类时加权余的 五种基本 方法, 对内 部法来说，涓 除余墨的统一格式是1. 于域活[SiHxhimain Method)此法曾先润■求帛苇域\「划分成11个子域 V., 在每个子域内令权国唾史二字于1, 而在子坦31寸2冈权两孕攵为弃，也即：修)''[0 外)如果在各个子域里分刎危取试丽效，那么它的求解在形式上将 类1以于孑「限力法.2. 配^v^rCollocatioii Method)W, -S（P-P.）子域法是令余星在一个子域上的总和为专°而西己点法是使 余童在才巳定的n个点上寺于多，这些敏才M为配农」此法的枚 两数为：■：：：■ 0' 0戏；口匝威犹拉豆)函敬,它的定义为:1° ga.宅[a 、 b\ \ e [a , b]P, P]——公易u代表求解域内任一点矛口两己点。

由于此己去只在百如炽上保证汆童为 审，因 此小岳牛毛作积分计弃, 斯以尾击二简单的加权金也睡去3. *小二求法(L己溶t Square Method;本话灵1T史在至个呀L辟坦I上余-致的 平方和]艮彳及刈、来 比立参宁食余-疆的-条[牛者记余-进平方和为1(C), 即7(C)- irF^dV- [R^R.dV则极伍条件为： 七1(碧)「底减『=0 由此矽见，本法杈网效为： 叽=空 打=1,辽・5)4. 伽句金法Meth(Kl)I本法是使余是与尊一个基函数正 交，也 即以基 羽数作为 权 两 教 许；,.=dV. (/ = 1,2J.. Ji)当试透*攵条包含摇个完备泗数为时，用本5去必可求得精确解」5. 3日去(hRthcid of Moment)本土去与伽金法相［队，也是用无备函:数WK乍枳甬数心T旦本矛r的 杈函 政与伽辽金乂, 区匆，它与试函 数工关泌除余量的条什是从专开始的各阶矩为 安，因此 对一维问戏L * =产 (:=L£. L 的对二维问题咬=—护』(摄= 12L网其余类推W冬五牙中基本方法布彳才定东*攵足够多声聚做有■惭＜±tf以沪十，兰小衍 反彼此和近：_■. TN刈"TiV介近T以(11较刁、，忻况■下，后三牙中的柄反 要高于前两种、•，密本方法举例为说明上述基本概念，以图祈示寺戏面思汁梁，受满跨均布荷 我昨川，美悬听境R的竖向住移圣为例，说明域本方法的应用° I图示梁的抵:制方程■为：嗜寸。

ML 5更界条［中为不艰验证其满民边界条什，也即 再=臼 而控制方程的内鄢余坚妇为:R- - £7r(120.v+ 24/j-<7囚 此本问 炊屈 内部法:、下 面分列 用 基本方三去进行求解:，子域法解由廿h式网救■［又一个待定帝效. 也p匕只升—个子坷I （寺于金坦I） 叩 可，消除金-矗的条什为：i［如1 汕+W）-“ ^--o由此可繇福■: 口_ 4~ S-1EII代回『）式可得：⑴侦配点法解同上所此，只舁卜造一个配求来起互消除余-*的条1牛’各令RiL=o.^ =可得：°=_弓_ 湛NI 14成V » 3宜若令 W=L °州再：■-=—L_ （W）144 2?J * 7；rt.r时见不同的配瓜结果岂矛、一杵的•、最小法解此时泌除余量纹条件为1：马^^人=I：官如12■皈十汕 一9_|・|_血1叫工十2^」亦=0<7-° 0101^ 萨 归/可■（异： EE 、-—页—伽辽金法解以匕甘寸， N］=盘-氐4 -14广疽-否"部除余造的条件为：口 N^dx - 0由此可得："墀弊理=竺|岌矩法解由 于只有一个冷定沱数，因 此泌 除余量条件只希 答次矩对可.此时显然与子域法完全相同a本例各方法的精度比较本川1 的*1■确 解 由M晋1立才多it乒可彳孚为：I、 _ q/4 _ 0 125“尸'~ 8£I ~ ~EI由D匕矽彳手，上』还各方三去对本他H十算的1奖泾依次为：-33.3%： 1.75% （22.2%） ; 13.9%; 0.96%; -33.3%上而22.2%为式。

琲）结果『2、虚功原理——平衡方程和几何方程的等效积分“弱”形式虚功原理包含虚位移原理和虚应力原理，是虚位移原理和虚应力原理 的总称他们都可以认为是与某些控制方程相等效的积分“弱”形式虚 功原理：变形体中任意满足平衡的力系在任意满足协调条件的变形状态上 作的虚功等于零，即体系外力的虚功与力的虚功之和等于零虚位移原理是平衡方程和力的边界条件的等效积分的“弱”形式；虚应力原理是几何方程和位移边界条件的等效积分“弱”形式虚位移原理的力学意义：如果力系是平衡的，则它们在虚位移和虚应 变上所作的功的总和为零反之，如果力系在虚位移（及虚应变）上所作 的功的和等于零，则它们一定满足平衡方程所以，虚位移原理表述了力 系平衡的必要而充分条件一般而言，虚位移原理不仅可以适用于线弹性 问题，而且可以用于非线性弹性及弹塑性等非线性问题虚应力原理的力学意义：如果位移是协调的，则虚应力和虚边界约束 反力在他们上面所作的功的总和为零反之，如果上述虚力系在他们上面 所作的功的和为零，则它们一定是满足协调的所以，虚应力原理表述了 位移协调的必要而充分条件虚应力原理可以应用于线弹性以及非线性弹性等不同的力学问题但 是必须指出，无论是虚位移原理还是虚应力原理，他们所依赖的几何方程 和平衡方程都是基于小变形理论的，他们不能直接应用于基于大变形理论 的力学问题。

3、最小总势能法应变能：作用在物体上的外载荷会引起物体变形，变形期间外力所 做的功以弹性能的形式储存在物体中，即为应变能由n个单元和m个节点组成的物体的总势能为总应变能和外力所做 功的差：n=以(e)-珏ui i e=1 i=1最小势能原理：对于一个稳定的系统，相对于平衡位置发生的位移总会使系统的总势能最小，即:an 顶人 合》=A (e)-du du dui i e=1 i i =1F = 0，i=1,2,3,……,n i i有限元法的收敛性有限元法是一种数值分析方法，因此应考虑收敛性问题有限元法的收敛性是指：当网格逐渐加密时，有限元解答的序列收 敛到精确解；或者当单元尺寸固定时，每个单元的自由度数越多，有限 元的解答就越趋近于精确解有限元的收敛条件包括如下四个方面：1） 单元，位移函数必须连续多项式是单值连续函数，因此选择 多项式作为位移函数，在单元的连续性能够保证2） 在单元，位移函数必须包括常应变项每个单元的应变状态总 可以分解为不依赖于单元各点位置的常应变和由各点位置决定的变量 应变当单元的尺寸足够小时，单元中各点的应变趋于相等，单元的变 形比较均匀，因而常应变就成为应变的主要部分。

为反映单元的应变状 态，单元位移函数必须包括常应变项3） 在单元，位移函数必须包括刚体位移项一般情况下，单元任 一点的位移包括形变位移和刚体位移两部分形变位移与物体形状及体 积的改变相联系，因而产生应变；刚体位移只改变物体位置，不改变物 体的形状和体积，即刚体位移是不产生变形的位移空间一个物体包括 三个平动位移和三个转动位移，共有六个刚体位移分量由于一个单元牵连在另一些单元上，其他单元发生变形时必将带动 单元做刚体位移，由此可见，为模拟一个单元的真实位移，假定的单元 位移函数必须包括刚体位移项4） 位移函数在相邻单元的公共边界上必须协调对一般单元而言， 协调性是指相邻单元在公共节点处有相同的位移，而且沿单元边界也有 相同的位移，也就是说，要保证不发生单元的相互脱离开裂和相互侵入 重叠要做到这一点，就要求函数在公共边界上能由公共节点的函数值 唯一确定对一般单元，协调性保证了相邻单元边界位移的连续性但是，在板壳的相邻单元之间，还要求位移的一阶导数连续，只有 这样，才能保证结构的应变能是有界量总的说来，协调性是指在相邻单元的公共边界上满足连续性条件前三条又叫完备性条件，满足完备条件的单元叫完备单元；第四条 是协调性要求，满足协调性的单元叫协调单元；否则称为非协调单元。

完备性要收敛的必要条件，四条全部满足，构成收敛的充分必要条件在实际应用中，要使选择的位移函数全部满足完备性和协调性要比 较困难的，在某些情况下可以放松对协调性的要求需要指出的是，有时非协调单元比与它对应的协调单元还要好，其 原因在于近似解的性质假定位移函数就相当于给单元施加了约束条 件，使单元变形服从所加约束，这样的替代结构比真实结构更刚一些 但是，这种近似结构由于允许单元分离、重叠，使单元的刚度变软了， 或者形成了（例如板单元在单元之间的绕度连续，而转角不连续时，刚 节点变为铰接点）对于非协调单元，上述两种影响有误差相消的可能， 因此利用非协调单元有时也会得到很好的结果在工程实践中，非协调 元必须通过“小片试验后”才能使用应力的单元平均或节点平均处理方法最简单的处理应力结果的方法是取相邻单元或围绕节点各单元应 力的平均值• 1.取相邻单元应力的平均值这种方法最常用于3节点三角形单元中这种最简单而又相当实 用的单元得到的应力解在单元是常数可以将其看作是单元应力的 平均值，或是单元形心处的应力由于应力近似解总是在精确解上 下振荡，可以取相邻单元应力的平均值作为此两个单元合成的较大 四边形单元形心处的应力。

如2单元的情况下，取平均应力可以采用算术平均，即平均应力=（单元1的应力+单元2的应力）/2也可以采用精确一些的面积加权平均，即平均应力=［单元1应力乂单元1的面积+单元2应力乂单元2面积］/ （单元1面积+单元2面积）当相邻两单元面积相差不大时，两者的结果基本相同在单元划o分时应避免相邻两单元的面积相差太多，从而使求解的误差相近一般而言，3节点三角形单元的最佳应力点是单元的中心点，此 点的应力具有1阶的精度• 2.取围绕节点各单元应力的平均值首先计算围绕该节点（i）周围的相关单元在该节点出的应力值 ，然后以他们的平均值作为该节点的最后应力值=上理b ei i m ie=1其中，1~m是围绕在i节点周围的全部单元取平均值时也可进行面 积加权有限元法求解问题的基本步骤1. 结构离散化对整个结构进行离散化，将其分割成若干个单元，单元间彼此通过 节点相连；2. 求出各单元的刚度矩阵[K](e)[K3)是由单元节点位移量{①} (e)求单元节点力向量{F}(e)的转移矩 阵，其关系式为：{F}(e)= [K] (e) {①} (e)3. 集成总体刚度矩阵[K]并写出总体平衡方程：总体刚度矩阵[K]是由整体节点位移向量{①}求整体节点力向量 的转移矩阵，其关系式为{F}= [K] {①},此即为总体平衡方程。

4. 引入支撑条件，求出各节点的位移节点的支撑条件有两种：一种是节点n沿某个方向的位移为零，另 一种是节点n沿某个方向的位移为一给定值5. 求出各单元的应力和应变对于有限元方法，其基本思路和解题步骤可归纳为：(1) 建立积分方程，根据变分原理或方程余量与权函数正交化原理， 建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发 点2) 区域单元剖分，根据求解区域的形状及实际问题的物理特点， 将区域剖分为若干相互连接、不重叠的单元区域单元划分是采用有限 元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点 进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时 还需要列出自然边界和本质边界的节点序号和相应的边界值3) 确定单元基函数，根据单元中节点数目及对近似解精度的要求， 选择满足一定插值条件的插值函数作为单元基函数有限元方法中的基 函数是在单元中选取的，由于各单元 具有规则的几何形状，在选取基 函数时可遵循一定的法则4) 单元分析：将各个单元中的求解函数用单元基函数的线性组合 表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积 分，可获得含有待定系数(即单元中各节点 的参数值)的代数方程组， 称为单元有限元方程。

5) 总体合成：在得出单元有限元方程之后，将区域中所有单元有 限元方程按一定法则进行累加，形成总体有限元方程6 )边界条件的处理：一般边界条件有三种形式，分为本质边界条 件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条 件(柯西边界条件)对于自然边界条件，一般在积分表达式中可自动 得到满足对于本质边界条件和混合边界条件，需按一定法则对总体 有限元方程进行修正满足7)解有限元方程：根据边界条件修正的总体有限元方程组，是含 所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求 得各节点的函数值单元刚度矩阵的特性单元刚度矩阵无论在局部坐标系中还是在整体坐标系中都具有相同 的三个特性:1) 对称性由材料力学中的位移互等定理可知，对一个构件，作用在点［的力引起点i的绕度等于有同样大小而作用于点i的力引起的点j的绕度，即 k..(e)=k.i(e)，表明单元刚度矩阵是一个对称矩阵2) 奇异性无逆阵的矩阵就叫做奇异矩阵，其行列式的值为0，即|k(e)|=0，这 一点可以从例题直接得到验证其物理意义是引入支撑条件之前，单元 可平移3) 分块性有前面所讲的容可以看出，矩阵［k(e)］可以用虚线分成四块，因此可写成如下的分块形式，1221 22{①}12式中k (e)mn局部坐标系中单元(e)按局部码标记的节点m、n之间的刚度子矩阵刚架结构中非节点载荷的处理的方法在刚架结构以及其他较复杂的结构上，他们所受的载荷可以直接作 用在节点上，又可以不直接作用在节点上而作用于单元节点间的其他位 置上。

后一种情况下的载荷称为非节点载荷有限元分析时，总体刚度 方程中所用到的力向量 是节点力向量因此在进行整体分析前应当 进行载荷的移植，将作用于单元上的力移植到节点上移植时按静力等 效的原则进行处理非节点载荷一般可直接在整体坐标系进行，其过程为：1） 将各杆单元看成一根两端固定的梁，分别求出两个固定端的约束 反力其结果可直接利用材料力学的公式求得；2） 将各固定端的约束反力变号，按节点进行集成，获得各节点的等 效载荷总体刚度矩阵的集成法使用刚度矩阵获得的方法获得总体刚度矩阵在此将其扩展到由整体坐标系中的单元刚度矩阵的子矩阵集成总体刚度矩阵步骤如下：1） 对一个有n个节点的结构，将总体刚度矩阵［K］划分为nXn各子 区间，然后按节点总码的顺序进行编号；2） 将整体坐标系中单元刚度矩阵的各子矩阵根据其下标的两个总码 对号入座，写在总体刚度矩阵相应的子区间；3） 同一子区间的子矩阵相加，成为总体刚度矩阵中的相应的子矩阵总体刚度矩阵的特性1） 对称性：因为由此特性，在计算机中只需存储其上三角部分；2） 奇异性：物理意义仍为在无约束的情况下，整个结构可做刚体运动；3） 稀疏性：［K］中有许多零子矩阵，而且在非零子矩阵中还有大量的零元素，这种矩阵称为稀疏矩阵。

大型结构的总体刚度矩阵一般都是稀疏矩阵；4） 分块性：平面问题离散化时的规定1） 单元之间只在节点处相连；2） 所有的节点都为铰接点；3） 单元之间的力通过节点传递；4） 外载荷都要移植到节点上；5） 在节点位移或某一分量可以不计之处，就必须在该节点安置一个铰支座或相应的连杆支座通过以上的规定来建立平面有限元分析模型结构对称性的利用规律一般来说，作用在对称结构上的载荷系统分为对称的、反对称的和一般的三种情况1. 结构对称，载荷对称或反对称这种情况下，对称面上的边界条件可按以下规则确定：A. 在不同的对称面上，将位移分量区分为对称分量和反对称分量；B. 将载荷也按不同的对称面分别区分为对称分量和反对称分量；C. 对于同一个对称面，如载荷是对称的，则对称面上位移的反对称 分量为零，如载荷是反对称的，则对称面上位移的对称分量为零如果所分析的结构对称，但载荷是不对称的，也不是反对称的，这 时可以将这种结构系统简化成载荷为对称和/或反对称情况的组合，仍 可以简化分析过程，提高分析的综合效率如图a所示，结构对称，载荷一般，可将其载荷分解为图b和图c 的组合图b为对称结构，载荷对x、y轴均为对称，图c为结构对称， 载荷对x轴反对称、对y轴对称，此时可取相同的四分之一进行研究， 分别施加对称面上节点的边界条件，进行两次分析计算，并将计算结果 迭加起来，即可得到原结构四分之一的解答，进而得出整个结构的解答。

利用结构的对称性取某一部分建立有限元模型时，往往会产生约束不足现象例如，若取上例中图c的四分之立有限元时，根据上述分析，在两 对称面上应加水平放置的滚动铰支座，因此模型在垂直方向存在刚体位 移对这种约束不足问题，利用有限元分析时，必须增加附加约束，以 消除模型的刚体位移在本例中，垂直方向可以用刚度很小的杆单元或 边界弹簧单元连接到模型某节点上，使得既消除了模型的刚体位移，又 不致于因附加的杆单元或边界弹簧单元刚度太大而影响结构原有的变 形状态单元形态的选择原则单元形态包括单元形状、边中节点的位置、细长比等，在结构离散 化过程中必须合理选择一般来说，为了保证有限元分析的精度，必须 是单元的形态尽可能的规则对于三角形单元，三条边长尽量接近，不应出现大的钝角、大的边 长这是因为根据误差分析，应力和位移的误差都和单元的最小角的正 弦成反比因而，等边三角形单元的形态最好，它与等腰直角三角形单 元的误差之比为sin45°:sin60° =1:1.23但是为了适应弹性体边界， 以及单元由小到大逐渐过渡，不可能是所有的三角形单元都接近等边三 角形实际上，常常使用等腰直角三角形对于矩形单元来说，细长比不宜过大。

细长比是指单元最大尺寸和 最小尺寸之比最优细长比在很大程度上取决于不同方向上位移梯度的 差别梯度较大的方向，单元尺寸要小些，梯度小的方向，单元尺寸可 以大一些；如果各方向上位移梯度大致相同，则细长比越接近1，精度 越高有文献推荐，一般情况下，为了得到较好的位移结果，细长比不 应超过7；为了获得较好的应力结果，细长比不应超过3一般情况下，正方形单元的形态最好对于一般的四边形单元应避免过大的边长比，过大的边长比会导致 病态的方程组边界条件的确定确定边界条件是建立有限元模型的重要一环，合理确定有限元模型 的边界条件是成功地进行结构有限元分析的基本要求一般情况下，建模对象的边界条件是明确的根据力学模型的边界 条件可以很容易确定其有限元模型的边界条件例如电线杆插入地基的 一端为固定端，桥梁一端为固定铰支座，另一端为滚动较支座但是，在机械工程中，建模对象往往是整个结构中的一部分，在建 立有限元模型，确定其边界条件时，必须考虑其余部分的影响这方面 主要考虑如下两类问题1.边界位置的确定在建立连续弹性体局部区域的有限元模型时，往往取该局部区域为 隔离体，取其隔离边界条件为零位移约束，并通过试探校正确定零位移 边界条件的位置。

例如，进行齿轮齿有限元分析时，取一个轮齿的局部 区域为隔离体，如图所示，设定PQRS的边界条件为零位移约束，通过 改变边界深度PQ和边界宽度PS研究边界位置对齿根最大拉应力的影 响，最后确定合理的边界条件2. 边界条件的确定有些分析对象的边界位置是零部件的连接部位在建立有限元模型 时，必须研究如何给定边界位置上的边界条件，以反映相连接结构的影 响确定这种问题的边界条件是用简单支撑连杆替代相连接结构的作 用，使替代后结构的系统刚度等价于原结构的系统刚度如分析机床主 轴和传动轴时，可以利用等刚度的杆单元替代轴承和支座的作用，使轴 的分析中包含有轴承和支座的影响当利用整体刚度矩阵的带状特征进行存贮和求解方程组时，单元节 点编号直接影响系统整体刚度矩阵的半带宽，也就是影响在计算机中存 贮信息的多少、计算时间和计算费用因而，要求合理的节点编号使带 宽极小化半带宽的计算公式：半带宽d二（单元节点号的最大差值+1）X节点自由度由此，进行网格节点编号时应使网格中单元节点号的最大差值最小，这样才能保证半带宽最小试比较下图图所示网格的四种编号方案中，单元节点标号的最大差值分别为5,3, 5, 9显然，图2方案要合理。

由此得出结论：沿着短边方向按列- 列-列-列地顺序编号比沿着长度方向按行-行-行-行地顺序要合理（半 带宽小）平面问题中非节点载荷转换为等效节点载荷由于三角形单元复杂的力学性质，不能像分析刚架时那样简单地利 用材料力学公式来求解，而要用虚功方程将加在结构上的非节点载荷转 换为等效节点载荷掌握以下两种常见的非节点载荷的移植结果1）作用在单元一条侧边上的集中力设Q平行于x方向，如图4-14所示，则等效节点载荷为F = ", F = :Q, F = 0, F = 0, F = 0, F = 0 ix l jx l kx iy jy ky若Q平行于y方向，结果与此相仿x2）作用在单元一条侧边上呈三角形分布的载荷设载荷平行于x方向，如图4-15所示，则等效节点载荷为F =四,F =亚,F = 0, F = 0, F = 0, F = 0 ix 6 jx 3 kx iy jy ky若分布载荷为集度是4的均布载荷，则F= j气其余分量为零求解时模型是否准备就绪?在求解初始化前，应进行分析数据检查，包括下面容：1.统一的单位;2.单元类型和选项;3.材料性质参数:考虑惯性时应 输入材料密度；热应力分析时应输入材料的热膨胀系数；4.实常数（单 元特性）；5.单元实常数和材料类型的设置；6.实体模型的质量特性 (Preprocessor > Operate > Calc Geom Items)；7 .模型中不应存在的 缝隙;8.壳单元的法向;9.节点坐标系；10.集中、体积载荷面力方向；11. 温度场的分布和围；12.热膨胀分析的参考温度。