有限单元法基本思想原理数值计算过程

数码资讯 2019-09-10106未知admin

  有限单元法基本思想,原理,数值计算过程_物理_自然科学_专业资料。通过最简单的三结点三角形单元在平面问题中的应用介绍有限单元法基本思想,原理,数值计算过程。

  有限单元法学习报告 在对力学问题分析求解过程中,方法可以概括为两种方法,一种为解析法,对具体问题具 体分析,通过一定的推导用具体的表达式获得解答,由于实际工程中结构物的复杂性,此方 法在处理工程问题是十分困难的;另一种是数值法,有限元法是其中一种方法,其数学逻辑 严谨,物理概念清晰,又采用矩阵形式表达基本公式,便于计算机编程,因此在工程问题中 获得广泛的应用。 有限元法基本原理是, 将复杂的连续体划分为简单的单元体; 将无限自由度问题化为有限 自由度问题, 因为单元体个数是有限的; 将偏微分方程求解问题化为有限个代数方程组的求 解问题。通常以位移为基本未知量,通过虚功原理和最小势能原理来求解。 基本思想是先化整为零, 即离散化整体结构, 把整体结构看作是由若干个通过结点相连的 单元体组成的整体;再积零为整,通过结点的平衡来建立代数方程组,最后计算出结果。 我将采用最简单的三结点三角形为基本单元体, 解决弹性力学中的平面问题为例, 解释有 限单元法的基本原理、演示数值计算过程和一般性应用结论。 一、离散化 解决平面问题时,主要单元类型包括三角形单元(三结点、六结点)和四边形单元(四结 点矩形、四结点四边形、八结点四边形)等。选用不同的单元会有不同的精度,划分的单元 数越多,精度越高,但计算量也会越大。因此在边界曲折,应力集中处单元的尺寸要小些, 但最大与最小单元的尺寸倍数不宜过大。 在集中力作用点及分布力突变的点宜选为结点, 不同厚度,不同材料不能划分在同一单元中。三角形单元以内角接近 60°为最好。充分利 用对称性与反对称性。 二、单元分析 将一个单元上的所有未知量用结点位移表示,并将分布在单元上的外力等效到结点上。 1、位移函数选取: 根据有限元法的基本思路,将连续体离散为有限的单元集合后,此时单元体满足连续性、 均匀性、各向同性、完全线弹性假设。单元与单元之间通过结点连接并传递力,位移法(应 用最广)以结点位移δ i=(ui vi) 为基本未知量,以离散位移场代替连续位移场。单元体内 的位移变化可以用位移函数(位移模式)来表示,因为有限元分析所得结果是近似结果,为 了保证计算精度和收敛性,x 位移函数应尽可能反应物体中的真实位移,即满足完备性和连 续性的要求: T ?位移模式必须能反映单元的刚体位移。 ?位移模式必须能反映单元的常量应变。 ?位移模式应尽可能反应位移的连续性。 设三角形单元三个结点编号为 i、j、m。平面三角形单元位移函数选取为 u=α 1+α 2x+α 3y v=α 4+α 5x+α 6y 可以写成 ? ?u ? u0 ? ? y 的形式,u0、v0 反映了单元的刚体平动,? 反映了单元的刚体转动, v ? v ? ? y 0 ? e 满足完备性和连续性的要求?。 采用插值法由单元结点位移列阵δ = ? ui 0 Nj vi uj vj um vm ? T 计算α 1、α 2、α 3、α 4、 e α 5、α 6.,求出位移 d=[u(x,y), v(x,y)]。6 个未知量,6 个代数方程,得 de=Nδ de= ? ? = ? ? u ? ? Ni ?v? ? 0 0 Ni Nj 0 Nm 0 0 ? ? ?u Nm ? i vi uj vj um vm ? T 式中 Ni=(ai+bix+ciy)/2A, ai = xj xm 1 yi 1 xj yj b i= ci = (i、j、m 轮换)A 为 1 ym 1 xm ym 三角形面积,为避免 A0,i、j、m 按逆时针排列。N 为形函数矩阵,形函数 Ni 的性质有: ?Ni(xi,yi)=1 Ni(xj,yj)=0 Ni(xm,ym)=0 ?Ni(x,y)+Nj(x,y)+Nm(x,y)=1 可推出三个形函数中,两个是独立的,反映了刚体平移。 令 z=Ni,在直接坐标系中画出 Ni、Nj、Nm 的函数图形是以 Ni(xi,yi)=1 为高的四面 体,所以结点位移影响单元的位移场,单元的位移场是线性分布的,相邻单元在公共边上的 位移是连续的, 单元相邻边的位移只取决于单元相邻公共边上的结点而与其他结点无关, 无 论以哪个单元计算相邻边的位移,结果一定相同。 形函数 Nie 决定了单元内的位移模式,反映了 i 结点位移对单元内任意点位移的贡献率。 2、根据几何方程用单元结点位移表示单元应变: ? ?u ? ? ? ? ?x ? ? Bi ? ?v ? 1 ? ? ?? ?? ?0 ? y 2 A ? ? ?c ? i ? ?u ?v ? ? ? ? ? ?y ?x ? 0 bj 0 ci 0 cj bi c j b j 0? ? 0 cm ? ? ui c mb ? ?m bm vi uj vj um ?v T m ? e ? B? e B 为几何矩阵 ? bi ? B 可写为分块矩阵 B=(Bi Bj Bm) ,Bi= ? 0 ?c ? i T 0? ? ci ? ,B 内所有元素与 x,y 无关,所以该单元 bi ? ? 内应变是常量, 反映单元的常量应变, 满足完备性和连续性的要求?, 这是一种常应变单元。 3、根据物理方程用单元结点位移表示单元应力: ? e ? D? e ? ?1 E ? D? ?? 1? ? 2 ? ?0 ? ? 1 0 ? 0 ? ? 0 ? 1? ? ? ? 2 ? D 为弹性矩阵 ? e ? D? e ? DB? e ? S? e e S 为应力矩阵 S=DB 中,每一个元素都是常数,所以 ? 的每一个分量与单元内 x,y 位置无关,这是一种 常应力单元。 应变是常量,其误差为 o ? ?x ? ,比位移精度低。 2 因为在三结点三角形单元中,位移函数中含有坐标的一次项,其误差为 o ?x ,而应力、 ? ? 4、根据虚功原理用单元结点位移表示单元结点力 单元在结点处受力, 单元会发生变形, 因此单元在结点处所受到的力与单元结点位移肯定 有关系。单元间通过结点的相互作用成为整体,因此每一单元的受力——位移关系找出来, 整体的受力——位移关系也就出来了。 e 记单元节点力为 F ? Fi e 单元内应力为 ? ? ? x ? Fj Fm ? T,单元结点虚位移为 ? = ? *i *e ? ? * j ? *m ? T ? ? y ? xy ? T, 单元内虚应变 ? * ? ? ? * x ? * y ? * xy ? T F e ? ?? ?? * ? ? e dxdy ? t ,可得 T A 根据虚功原理, ?? ? *e T F e ? ?? BT DBdxdy ? t ? ? e A 因为 B、D 中元素都是常数, F ? B DBtA? ? K? ,K=BTDBtA 为单元刚度矩阵。 e T e e ? BiT ? ? T? K 为 6 行 6 列矩阵可写为 K ? tA ? B j ? D ? Bi T ? ? Bm ? ? ? k xx kij ? Bi DB j At ? ? ijyx ?k ? ij Bj ? kii ? Bm ? ? ? k ji ?k ? mi kij kj j kmj kim ? ? k jm ? , kmm ? ? kijxy ? xy At , kij 表示 j 结点处发生 y 方向的单位位移时所引起的 i 结 yy ? ? kij ? 点处 x 方向的结点力。不同类型不同形式的单元,只有弹性矩阵 D 和几何矩阵 B 不同,计 算子块矩阵的公式相同,平面问题中,影响刚度矩阵 K 的只有几何矩阵 B。 K 的性质有: ?K 中每个元素表示个单元结点沿坐标方向发生单位位移时所引起的结点力。 ?K 为对称矩阵。 ?单元做刚体位移时,单元内不产生应变应力,结点力为 0,所以 K 中每行每列元素之和为 0,所以 K ? 0 ,所以只根据 F ? K? 无法求得唯一解。 e e 5、根据虚功等效原则计算等效结点力 根据有限元的基本方法,单元内任意点的位移、应变、应力等最终都要用结点位移来表示, 所以作用在物体上的外力也要用结点位移表示。为了计算等效结点力,在任意的虚位移上, 使原载荷与等效载荷虚功相等。 设外力为 f p ,结点虚位移为 ? *e ,则任意点虚位移为 d *e ? N? *e ,等效节点载荷为 FLe ,有 d *eT f pt ? ? *eT FLe 同理得 FL ? e FLe ? N T f pt (集中力) A ? S e N T f ds ? t (面力) , FL ? ?? Nfdxdy ? t (体力) 。 三、整体分析 将结构的所有单元通过结点连接起来, 形成一个整体的离散结构以代替实际的连续体, 以 形成以结点位移为未知量的整体结构的有限元代数方程组,最后求得结点位移。 对结点受力分析: 结点受到与之相关的单元给它的反作用力和外载荷的等效结点力, 这两组 力坐标轴方向相反,所以应该相等,即 ?F ? ?F i e e Li ,设有 n 个结点,每个结点建立两个 方向的方程,不考虑外界约束时,共 2n 个方程,2n 个未知量( ?ix , ?iy , ? jx ... ) ,为了建立这 个代数方程组,建立整个弹性体的结点力和结点位移的关系式 K? ? FL ,K(2n×2n)为整 理刚度矩阵,δ 为整体结点位移列阵,FL 整体结点载荷列阵。 为了求整体刚度矩阵,要找到它与已求得的单元刚度矩阵的关系,在整体中对结点编码, 设整体刚度矩阵中某元素为 Kij,意为 j 个结点在 x 或 y 方向发生位移引起 i 个结点 x 或 y 方向的结点力,找到同时用到 i 与 j 结点的单元,并用与之对应的单元刚度矩阵中的元素 ksm 相加得到 Kij,整体刚度矩阵也是奇异矩阵,必须考虑边界约束条件,消除 K 的奇异性, 才能求解结点位移。再由单位结点等效载荷得到整体结点载荷列阵 FL。这样 K、FL 已知,求 解代数方程,解出整体结点位移列阵δ ,得到相应的单元结点位移δ δ 得到了,相应的 d 、σ 、ε 等就得到了。 e e 。 e e e

Copyright © 2010-2020 南湖社区 版权所有  

联系QQ:1352848661