有限元二点高斯积分(高斯积分与减缩积分)
有限元二点高斯积分(高斯积分与减缩积分)
2024-11-06 01:43:31  作者:摇头玩  网址:https://m.xinb2b.cn/know/fsi248765.html
0. 有限元积分方法简介

在有限元分析中,为了计算各个单元的刚度矩阵和荷载向量,需要对单元上的积分进行数值近似。两种常见的积分方法是高斯积分和减缩积分。

高斯积分是一种数值积分方法,基于高斯-勒让德公式,用于近似计算函数在给定区间上的积分。在有限元分析中,高斯积分用于计算单元上的刚度矩阵和荷载向量。具体而言,对于每个单元,将其内部积分区域划分为一系列高斯积分点,然后根据高斯积分点的权重和位置,对积分表达式进行数值近似。通过对所有高斯积分点的贡献进行累加,可以得到单元上的刚度矩阵和荷载向量的近似值

减缩积分(Reduced Integration)是一种在有限元分析中用于减少计算量和存储需求的技术,它采用了较少的高斯积分点来近似计算单元的刚度矩阵。减缩积分的核心思想是,通过在单元内部使用较少的高斯积分点来减少计算量,同时尽可能保持良好的数值稳定性。传统的高斯积分在每个单元上使用一组完整的高斯积分点,这可能会导致计算量过大,尤其是对于高阶元素或复杂几何形状的单元。减缩积分通过选择适当的积分点子集,例如在三维问题中选择积分点子集的体心点,来近似计算刚度矩阵。减缩积分的优点是减少了高斯积分点的数量,从而降低了计算量和存储需求。然而,它也可能导致计算结果的不精确性。由于减少了积分点的数量,减缩积分可能无法捕捉到一些细节信息,特别是在处理某些特殊情况(如奇异问题)时。总结来说,减缩积分是一种在有限元分析中用于减少计算量和存储需求的技术,它通过使用较少的高斯积分点来近似计算刚度矩阵。它可以在一定程度上降低计算开销,但也需要注意其对计算结果的影响

1. 高斯积分与减缩积分的精度探究

以平面应力问题为例,使用mathematical计算了四边形单元的刚度矩阵,对比了精确解,高斯积分解,减缩积分解的差异性。

Clear["Global`*"](*Parameters*)t=1;Em=1;v=0.2; (*mathematical计算默认都是精确解,但是我们提供的数据必须是分数;其次,对于大规模矩阵运算,精确解的求解计算非常慢*)Coord={{0,0},{3,2},{4,4},{2,5}};(*Constitutive Law-Plain Stress*)Cc=(Em/(1 v))*{{(1/(1-v)),(v/(1-v)),0},{(v/(1-v)),(1/(1-v)),0},{0,0,1/2}};(*Shape Functions*)N1=(1-eta)*(1-xi)/4;N2=(1-eta)*(1 xi)/4;N3=(1 eta)*(1 xi)/4;N4=(1 eta)*(1-xi)/4;Ni={N1,N2,N3,N4};Head[eta]; (*没有赋值直接拿来用的变量属于symbol数据类型*)Head[N4]; (*symbol数据类型的组合是Times数据类型,实际上还是symbol类型,Times代表symbol类型之间运算后赋值给N1*)(*Coordinate Mapping*)CoordMap={x,y};x=Sum[Ni[[i]]*Coord[[i,1]],{i,1,4}];y=Sum[Ni[[i]]*Coord[[i,2]],{i,1,4}];(*Jacobian Matrix*)NatCoord={xi,eta};J=Table[D[CoordMap[[i]],NatCoord[[j]]],{i,1,2},{j,1,2}];Jinv=Transpose[Inverse[J]];H=Jinv//MatrixForm ; (*这里H仅仅是矩阵格式,不是原来的Jinv,这种格式就是用来显示的,一种代表形式的数据类型*)(*B Matrix*)B1={{1,0,0,0},{0,0,0,1},{0,1,1,0}};B2={{Jinv[[1,1]],Jinv[[1,2]],0,0},{Jinv[[2,1]],Jinv[[2,2]],0,0},{0,0,Jinv[[1,1]],Jinv[[1,2]]},{0,0,Jinv[[2,1]],Jinv[[2,2]]}};B2//MatrixForm;B3=Table[0,{i,1,4},{j,1,8}]; (*如果是四边形八节点的二次单元,这里改为{j,1,8}*)Do[B3[[j,2i-1]]=B3[[j 2,2i]]=D[Ni[[i]],NatCoord[[j]]],{i,1,4},{j,1,2}];B3//MatrixForm; (*不要写成B=B1.B2.B3//MatrixForm*)B=B1.B2.B3;(*Stiffness Matrix*)IntFunc=(Transpose[B].Cc.B)*Det[J]; ke=t*NIntegrate[IntFunc,{xi,-1,1},{eta,-1,1}];ke//MatrixForm; (*如果计算精确解,用Integrate,即,kn=t*Integrate[IntFunc,{xi,-1,1},{eta,-1,1}],但是数据必须都是分数*)(*下面是用高斯积分,又算了一次刚度矩阵中的积分*)xiset={1/Sqrt[3],-1/Sqrt[3]};etaset={1/Sqrt[3],-1/Sqrt[3]};w={1,1}; (*权重系数列表*)kfull=t*Sum[w[[i]]*w[[j]]*(IntFunc/.{xi->xiset[[i]],eta->etaset[[j]]}),{i,1,2},{j,1,2}];(*下面是用减缩积分,又算了一次刚度矩阵中的积分*)kreduced=t*2*2 IntFunc/.{xi->0,eta->0};(*对比三种积分方法的差异*)ke//MatrixForm ;kfull//MatrixForm;kreduced//MatrixForm;



刚度矩阵精确结果


刚度矩阵高斯积分结果


刚度矩阵减缩积分结果

2. 结果分析与建议

如上图所示,高斯积分的结果与解析解几乎一致,减缩积分的误差略大,可见减缩积分的确会降低精度。进一步改变单元形状发现,误差与单元不规则程度相关,例如,单元越接近矩形,精度越高。其中当单元形状为矩形时,高斯积分解等于精确解,这就是有限元软件中强调网格质量的原因(当然加密网格密度也是一种解决方法),或许很多人都因为单元扭曲严重而被终止或者警告分析。当然如果因为网格问题影响收敛时,不妨考虑以下方法:

重新网格划分:如果单元扭曲严重,可能需要重新划分网格,以改善单元的形状和质量。使用更合适的网格生成算法或软件工具,可以生成更好的网格布局,减少扭曲现象的发生。使用更适合的单元类型:对于特定类型的问题,某些单元类型可能更适合,并能够减少单元扭曲。局部修复或调整:对于少数扭曲严重的单元,可以尝试通过局部调整来修复它们。这可能包括重新定义节点位置或移动单元的某些节点,以改善其形状。增强积分(Enhanced Integration):增强积分是一种用于提高有限元分析精度的技术。通过在高应变或高梯度区域增加高斯积分点,可以提高对这些区域的数值近似。增强积分可以帮助解决扭曲导致的数值不稳定性问题。错误估计与自适应网格:使用错误估计技术来检测扭曲严重的单元,并基于估计的错误信息进行自适应网格细化或粗化。这可以帮助优化网格,减少扭曲并提高计算结果的准确性。3.结语

最近刚复习了一遍有限元分析基础,旨在希望加深理解原理,进一步更好的应用有限元软件,有不当之处,欢迎指正,后续会抽时间更新其他笔记,感兴趣的同学可以看看我的上一篇有限元编程之杆单元,希望大家一起学习交流。

  • 如何用烤箱烤面包
  • 2024-11-06如何用烤箱烤面包1、准备好所需原材料:面粉、鸡蛋、黄油、白糖、酵母、温水2、把除黄油之外的其他原材料混合,揉成光滑面团3、将隔水融化的黄油与面团揉在一起,形成一个光滑的面团4、面团盖上保鲜膜或者塑料袋,醒发2小时5、。
  • 七情六欲代表什么意思(七情六欲分别指什么)
  • 2024-11-06七情六欲分别指什么七情六欲是一个心理反应一般来说,七情是指喜、怒、忧、思、悲、恐、惊,感情的表现或心理活动;六欲是指人的眼、耳、鼻、舌、身、意的生理需求或愿望不同学派对七情六欲内容的描述有所不同七情指“喜、怒、忧、思、。
  • 光遇新手任务怎么做(是怎么玩的)
  • 2024-11-06是怎么玩的寻找子民,这里的子民指的是游戏中的人影,找到蓝色人影并触碰它,跟着光点,咱们能看到一个深灰色的石像,点燃石像上的火苗,它就会教咱们一个新的互动姿势,找到黄色人影,与它们互动收集“光之翼”,咱们的飞行技。
  • 爱情三部曲第二部叫什么(爱在三部曲之后)
  • 2024-11-06爱在三部曲之后《爱在》三部曲是指由美国独立电影人理查德·林克莱编剧及执导的三部爱情电影片,它们分别是于1995年上映的《爱在黎明破晓前》、2004年上映的《爱在日落黄昏时》以及2013年上映的《爱在午夜降临前》有人。
  • 子沫女孩名字寓意(子沫女孩名字寓意有什么)
  • 2024-11-06子沫女孩名字寓意有什么女孩叫子沫有光彩照人、细腻轻柔的含义“瀑下一曰沫也”《说文》中有关于沫字的详细解释,本意就是液体形成的许多细泡,也可称作沫子、泡沫等,子有指孩子的代称,所以女孩叫子沫就是父母希望孩子以后如阳光下的泡沫。
  • 春分的古诗四句(春分的古诗四句有哪些)
  • 2024-11-06春分的古诗四句有哪些《春日田家》宋琬野田黄雀自为群,山叟相过话旧闻夜半饭牛呼妇起,明朝种树是春分《七绝》徐铉春分雨脚落声微,柳岸斜风带客归时令北方偏向晚,可知早有绿腰肥《隋堤柳》翁承赞春半烟深汴水东,黄金丝软不胜风轻笼行。
  • 如何避免热水器漏电(有效避免热水器漏电)
  • 2024-11-06有效避免热水器漏电在家享受热水澡的时候,很多小伙伴可能有过这样的担忧——这电热水器......不会漏电吧!?不怪大家想太多,电热水器漏电的事件每隔一段时间就会出现一例,实在是人心惶惶图片来源:网络而且,这种用电安全事故。
  • 今日全国最新猪价走势图官网(全国猪价再涨一天)
  • 2024-11-06全国猪价再涨一天今日猪评-8.13今日猪价二分天下,一半地区以稳为主,一半地区以上涨为主其中上海、浙江、福建和广东四个地区已经突破7块5,上海和广东两地更是突破7块7,离破8也是越来越近了另外还有个不错的消息是辽宁猪。
  • 本科毕业后的档案怎么处理(档案应该怎么处理)
  • 2024-11-06档案应该怎么处理这篇文章本来应该是每年的6月份给大家呈现的,虽然现在的学生已经都毕业了,但是准毕业生们还是可以用到的,咱们来了解一下吧其中,我们应该了解一下毕业生的身份档案的处理要分为两种情况,一种就是毕业时没有找到。
  • 吃石榴需要吐籽吗(吃石榴到底要不要吐籽)
  • 2024-11-06吃石榴到底要不要吐籽吃石榴要不要吐籽是个国际性话题在互联网上你可以看到不同文明在不同时期对此展开的大量争议其中还有人提出应该连皮一起啃当然如果你从来只是拿它榨汁,那么可以无视这个问题对很多石榴食客来说,要不要吐籽是句废话。
  • 安吉丽娜朱莉谍战电影(阿富汗版花木兰)
  • 2024-11-06阿富汗版花木兰有这样一部动画片,它讲述的是一个残酷的故事,但是人们看完后却心生温暖这部影片深深吸引了好莱坞巨星安吉丽娜·朱莉,她不仅为影片担任出品人,帮着找钱找人,还亲自去电影节走红毯为其站台影片上映后一片叫好之声。
  • 教你如何辨别真假燕窝 教你如何1分钟鉴别燕窝优劣
  • 2024-11-06教你如何辨别真假燕窝 教你如何1分钟鉴别燕窝优劣评价燕窝好不好的关键在发头吗?影响燕窝质量的指标很多,比如说品相,纤维粗细,洁净程度,燕窝酸和蛋白质含量,炖煮之后的味道等都是鉴别燕窝品质的关键,只能说发头只是众多鉴别指标中的一种买回来的燕窝有一种氨。