06有限元分析中的单元性质特征与误差处理_第1页
06有限元分析中的单元性质特征与误差处理_第2页
06有限元分析中的单元性质特征与误差处理_第3页
06有限元分析中的单元性质特征与误差处理_第4页
06有限元分析中的单元性质特征与误差处理_第5页
已阅读5页,还剩46页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、 o6.1单元节点编号与带宽存储 o6.2形状函数矩阵与刚度矩阵的性质 o6.3边界条件的处理与支反力计算 o6.4单元刚度矩阵的缩聚 o6.5为以函数构造与收敛性要求 o6.6C0型单元与C1型单元 o6.7单元的拼片试验 o6.8有限元分析数值解的精度与性质 o6.9单元应力的计算结果的误差与平均处理 o6.10控制误差和提高精度的h方法和p方法 6.1单元节点编号与带宽存储 o计算机进行有限元分析时, 需要存储所有单元和节点信 息,随着所求解问题自由度 的增大,计算规模的增大, 整体刚度矩阵的规模非常巨 大。 o由于整体刚度矩阵中显现出 相邻单元之间的关联性,因 此矩阵中的大部分数据都为

2、 零,反映非零数据的一个指 标就是带宽。 o由于刚度矩阵是对称的,可以看出,若节点的自由度数目为m,则 每一个单元在整体刚度矩阵的半带宽为 di=(第i个单元中节点编号的最大差值+1)m d=max( di ) (i=1,2n) 其中n为整个结构系统的单元数。显然 对于二维问题,m=2 对于三维问题,m=3 88 7877 6766 585655 464544 363433 2322 141211 0 0 00 000 00000 00000 K KK KK KKK KKK KKK KK KKK 88 7877 6766 585655 464544 363433 2322 141211 0 0

3、 0 0 00 0 K KK KK KKK KKK KKK KK KKK 84838281 74737271 64636261 54535251 44434241 34333231 24232221 14131211 KKKK KKKK KKKK KKKK KKKK KKKK KKKK KKKK 6.2形状函数矩阵形状函数矩阵与刚度矩阵的性质的性质 o以一维杆单元为例,杆单元 的位移场为 12 ( )1 ij ee iuijuj xx u xaa xuu ll N uN u iuju NNN 形函数矩阵形函数矩阵 1、左端发生单位位移,右端固定 2、右端发生单位位移,左端固定 3、发生刚体位移

4、 6.2形状函数矩阵与刚度矩阵的性质刚度矩阵的性质 o仍然以一维杆单元为例,它的刚度方程为 1、考虑单元左端发生单位位移,右端固定情况 2、考虑单元右端发生单位位移,左端固定情况 3、考察刚体位移 111211 212222 kkup kkup o性质1:单元刚度矩阵的对角元素kii表示要使单元的第i个节点产生 单位位移,而其它的节点位移为0时,需要在i点施加的节点力。 o性质2:单元刚度矩阵的对角元素kij(ij)表示要使单元的第j个节 点产生单位位移,而其它的节点位移为0时,需要在i点施加的节点 力。 o性质3:单元刚度矩阵是对称的。这可以由功的互等定理得到。对 于线弹性体,力所做的功跟加

5、载次序无关,这可以利用上面的性 质1和2得到。 111211 212221 1 0 kkk kkk 111212 212222 0 1 kkk kkk 第一种加载状态第二种加载状态 第一种加载状态下的外力在第二种加载状态下移动相应位移做的功为 112121 0+1=kkk 第二种加载状态下的外力在第一种加载状态下移动相应位移做的功为 122212 1+0=kkk 根据功的互等定理,可以得到结论:刚度矩阵刚度矩阵是对称对称的。 o性质4:单元刚度矩阵是半正定的。 o性质5:单元刚度矩阵是奇异的。 o性质6:单元刚度矩阵的任意行或列代表一个平衡力系,当节点位 移全部为线位移时,任意行或列的代数和应

6、该为0。 同样,由单元刚度矩阵所组装的整体刚度矩阵也有以下性质: 1)对称性 2)奇异性 3)半正定性 4)稀疏性 5)非零元素呈现带状分布 11 22 11 112121121 1221221 2 1 12121 11 0 22 11 . 22 1 . 2 nn t ijij ij nnnn nnnnn Uq Kqk u u k uk u uk u uk uk u uk u u k uk u uk u u 6.3边界条件的处理与支反力的计算 o位移边界条件在大多数情况下有两种类型。 1、零位移边界条件 2、给定具体数值的位移边界条件 根据上述两类边界条件,刚度方程的求解有以下几种方法: 1、

7、直接法 2、置“1”法 3、乘大数法 4、罚函数法 o直接法直接法 abba bbbb K qP K qP aaabaa babbbb KKqP KKqP 0 a q aaabba babbbb K uK qP K uK qP aaabaa babbbb KKqP KKqP a qu -1 bbbbba qKPK u 1、既可以处理零约束,又可以处理非零约束的情况。 2、处理过程直观。 3、待求矩阵的规模变小(维数变小),适合于手工处理。 4、矩阵的节点编号及排序改变,不利于计算机的规范化处理。 o置置“1”法法 bbbb K qP 100 0 a bbbb q KqP 0 a q 1、只能处

8、理零约束情况。 2、待求矩阵的规模不变,不需重新排列,适合于计算机处理。 3、保持整体刚度矩阵的对称性,利于计算机的规范化处理。 abba bbbb K qP K qP aaabaa babbbb KKqP KKqP 0 a q 直接法 o乘大数法乘大数法 aaabaaa babbbb MKKqMKu KKqP 1、 既可以处理零约束,又可以处理非零约束的情况。 2、待求矩阵的规模不变,不需重新排列。 3、保持整体刚度矩阵的对称性,利于计算机的规范化处理。 abba bbbb K qP K qP aaabaa babbbb KKqP KKqP 0 a q 直接法 aaaabbaa babbbb

9、 MKqK qMKu K uK qP o罚函数法罚函数法 罚函数法的最大好处是可以直接求出位移边界上的支反力。 支反力的计算:支反力的计算: 除了罚函数法能够求出支反力以外,其它的方法都需要求解一定的 方程得到。 6.4单元刚度矩阵的缩聚 o采用高次位移函数的单元也常被称为高阶单元。对于高次单元来 说,除了几何端点以外,其余的那些节点可能与其它的单元不发 生关系,当中间的节点与其它单元无关时,我们称作是内部节点。 而其余的节点是外部节点。既然内部节点与其他单元无关,那么 在组成整体刚度之前,就可以把他们消去,也就是把内部节点的 位移用外部节点的位移来表示。 o以一维三节点杆单元为例 1 122

10、33 ( )u xN uN uN u 11121311 21222322 31323333 kkkup kkkup kkkup o以一维三节点杆单元为例 11121311 21222322 31323333 kkkup kkkup kkkup aaabaa babbbb KKqP KKqP aaaabba baabbbb K qK qP K qK qP 1 1 aaaabbbbbaaa bbbbbaa K qK KPK qP qKPK q 其中a代表的是外部节点外部节点,b代表的是内部节点内部节点。 6.5位移函数构造与收敛性要求 o单元中的位移模式一般采用设有待定系数的有限多项式作为近似 函

11、数,优先多项式的选取原则应该考虑以下几个方面: 1、待定系数是由节点位移条件确定的,因此它的个数应该与节点 位移DOF个数相等。 2、在选取多项式时,必须选择常数项和完备的一次项。单元位移 模式中的常数项和一次项可以反映单元的刚体位移合唱应变的特 性。这是因为当划分的单元数趋于无穷时,即单元缩小趋于一点, 此时单元应变趋于常数。 3、选择多项式应该由低到高,尽量选取完全多项式以提高单元的 精度。 o因此,在构造一个单元的位移函数时,应该参考由多项式函数构 成的Pascal三角形和上述原则进行函数项次的选取与构造。 o收敛性问题收敛性问题 在有限元分析中,当节点数目或单元插值函数的项数趋于无穷大

12、时, 即单元尺寸趋于零时,最后的解答如果能够无线的逼近准确解, 那么这样的位移函数或形函数是逼近于真实的,这就称为收敛。 为使有限元分析的解答收敛,位移函数必须满足一些收敛准则,这 些准则都经过过严密的理论验证。主要包括以下三个方面。 o收敛性准则收敛性准则 定义:当单元尺寸趋于零时,有限元的解趋于真实解。 准则准则1:完备性准则(针对单元内部)。如果在势能泛函中所出现 的位移函数的最高阶导数是m阶,则有限元解答收敛性的条件之一 是选取单元内的位移场函数至少是m阶完全多项式。 准则准则2:协调性准则(针对单元之间)。如果在势能泛函中位移函 数所出现的最高阶导数是m阶,那么位移函数在单元交界面上

13、必须 具有直至m-1阶的连续导数,即Cm-1连续性。 o以一般的梁问题为例 21 dd 2 ll UW EIz vxp x v xx 从上式可以看出,所出现的物理量是关于位移的最高阶导数为2,因 此假定形状函数的时候,形函数至少应该包含完整的二次多项式。 由准则2可知,位移函数为C1连续,即在单元之间的位移函数至少要 求一阶导数连续。 o以一般的平面问题为例 1 d 2 dd p xxyyzzxyxyyzyzzxzx xyzxyz s UW b ub vb wp up vp wA 从上式可以看出,所出现的物理量是关于位移的最高阶导数为1,因 此假定形状函数的时候,形函数至少应该包含完整的一次多

14、项式。 由准则2可知,位移函数为C0连续,即在单元之间的位移函数要求零 阶导数连续。亦即函数的本身连续,而其一阶导数可以不连续。 o协调元与非协调元协调元与非协调元 当单元的位移函数满足完备性要求时,称单元是完备的(一般都比 较容易满足),当单元的位移函数满足协调性条件时,称单元是 协调的(在单元与单元之间的公共边界上对于高阶连续性要求较 难满足)。 当单元的位移函数即完备又协调时,则有限元分析的解答是收敛的, 即当单元尺寸趋于零时,有限元分析的解答趋于真实解。我们称 这种单元为协调单元。 o一般情况下,当泛函中的导数高于一阶时,则要求许可函数在单 元交界面上具有C1或更高的连续性,这时构造单

15、元的插值函数往 往比较困难。如果在单元之间的交界面上位移或导数不连续,将 在交界面上引起无限大的变形,这时必须产生附加应变能,而我 们建立泛函时,并没有考虑这种情况。因此,基于最小势能原理 得到的有限元分析解答就不可能收敛于正确解。 在某些情况下,可以放松对协调性的要求,只要这种单元能够通过 拼片试验,有限元分析的解答仍然可以收敛于正确的解答。这样 的单元称为非协调性单元。 6.6C0和C1型单元 oC0型单元 在泛函中位移函数的最 高阶导数为1,在交界面 上具有0阶的连续导数, 即节点上仅仅要求位移 连续。 杆单元、平面问题单元、 空间问题单元等 6.6C0和C1型单元 oC1型单元 在泛函

16、中位移函数的最高 阶导数为2,在交界面上具 有1阶的连续导数,即节点 上除要求位移连续外,还 要求1阶导数连续。 梁单元、板单元、壳单元 等 6.7单元的拼片试验 o由于非协调单元之间的位移不能保 证位移协调,可以通过拼片试验来 考证是否能描述常应变和刚体位移, 若能通过拼片试验,则解得收敛性 就能得到保证。 o如图所示的单元状况,其中至少一 个节点被单元所完全包围,若节点i 完全被单元所包围,节点i的平衡方 程为 1 0 m ee ijji e K uP o对于非协调单元,需要考察它的收敛性,即考察它是否具有常应 变的能力,因此,我们设计这样一个试验(拼片试验拼片试验): o当对单元片中的各

17、个节点赋予对应于常应变状态的位移和载荷值 时,核对对i点平衡方程的正确性,如果能够满足,也就是单元满 足常应变要求,因此当单元尺寸不断减小时,有限元解能够收敛 于真正解。 o以平面问题为例 012 012 jjj jjj uaa xa y vbb xb y 由片面问题的平衡方程可知,当单元内的应变或应力都为常数时, 则对应的体积力为零。对应于图中的i点,它的边界力也为零, 因此 。所以此时,通过拼片试验的前提是,当赋予各节 点以上位移模式的位移时,i点的平衡方程变为 0 e i P 1 0 m e ijj e K u 即必须在节点i施加附加约束,该约束力所作的功等于单元交界 面上位移不协调引起

18、的附加应变能。 o仍以平面问题为例 012 012 jjj jjj uaa xa y vbb xb y 由片面问题的平衡方程可知,当单元内的应变或应力都为常数时, 则对应的体积力为零。对应于图中的i点,它的边界力也为零, 因此 。所以i节点以外节点有以上位移模式的位移时,对 于i点的平衡方程 0 e i P 如果求解上式得到的位移值和常应变状态下的位移相一致,则认 为通过拼片试验。否则认为不能通过拼片试验。 1 0 m e ijj e K u 6.8有限元数值解的精度与性质 o求解精度估计 以平面问题为例,单元的位移场可以展开成以下形式 . i i i uv uuxy xy 如果单元尺寸为h,

19、则上式中的x和y都是h量级,若单元的位 移函数采用p阶完全多项式,即它能逼近上述泰勒级数的前p阶多 项式,那么位移解u的误差将是O(hp+1)量级。 3节点3角形单元(p次多项式): 量级 位移应变应变能 误差收敛速度误差误差 h/1量级: O(h2)O(h2)O(h)O(h2) h/2量级: O(h2/4) 。 h/3量级: O(h2/9) 。 h/4量级: O(h2/16) 。 o这里讨论的都是仅仅局限于网格的离散误差,即当一个连续的求 解域被离散成有限个子域,由单元的试函数来逼近整体的域的场 函数所引起的误差。另外,实际误差还应该包括计算机的数值运 算误差。 2 1 2 2 4 /4 i

20、i ii O h uu uuO h 21 1 4 3 iii uuu精确解与不同网格计算结果之间的关系 o有限元分析的下限性质 有限元是把结构无限多的自由度简化为有限多的自由度,结构的 刚度被夸大了,即使是用无限多个自由度来描述,也必然使得原 系统刚度增加,变得更加刚硬,即刚度矩阵的总体数值变大,由 刚度方程知,计算出的位移结果偏小。 由于位移函数的收敛性准则包含完备性和协调性两方面的要求, 而完备性要求比较容易满足,而协调性则较难满足,因此这往往 是研究的重点。 o位移解的下限性质是基于协调单元单调收敛的前提得到的,在有 些情况下,使用非协调单元也可以得到工程上的满意解答,有时 甚至更好,这

21、是由于位移不协调引所造成的误差与其它误差相抵 消的缘故。 6.9单元应力计算结果的误差和平均 o应力结果的误差性质 对于弹性问题,其三大变量 对于一个具体问题,成了求 2关于的极值问题。它是一 个误差泛函。 true true true uuu 1 ddd 2 p T TT trueijitrueitrue s trueUWb up uA true uu o可见,对于求近似解极值近似解极值的问题 从力学上看,是求位移变分引起的总势能为极小值的问题。 从数学上看,是求应变差和应力差在弹性矩阵加权意义下的最小 二乘问题。 因此,应变和应力的近似解的性质,是在加权残值最小二乘意义 上对真实应变和真实应力的逼

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论