有限元网格划分的基本原则

更新时间:2023-07-26 12:23:32 阅读: 评论:0

有限元网格划分的基本原则
有限元网格划分的基本原则
划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。为建立正确、合理的有限元模型,这里介绍划分网格时应考虑的一些基本原则。
1 网格数量
网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。
图1 位移精度和计算时间随网格数量的变化
在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。
2 网格疏密
网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图b中网格疏密相差更大,它比
图a中的网格少48个,但计算出的孔缘最大应力相差1%,而计算时间却减小了36%。由此可见,采用疏密不同的网格划分,既可以保持相当的计算精度,又可使网格数量减
蒸肉丸子
小。因此,网格数量应增加到结构的关键部位,在次要部位增加网格是不必要的,也是不经济的。1
图2 带孔方板的四分之一模型
徐东明划分疏密不同的网格主要用于应力分析(包括静应力和动应力),而计算固有特性时则趋于采用较均匀的钢格形式。这是因为固有频率和振型主要取决于结构质量分布和刚度分布,不存在类似应力集中的现象,采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大,可减小数值计算误差。同样,在结构温度场计算中也趋于采用均匀网格。
3 单元阶次
许多单元都具有线性、二次和三次等形式,其中二次和三次形式的单元称为高阶单元。选用高阶单元可提高计算精度,因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且高次插值函数可更高精度地逼近复杂场函数,所以当结构形状不规则、应
力分布或变形很复杂时可以选用高阶单元。但高阶单元的节点数较多,在网格数量相同的情况下由高阶单元组成的模型规模要大得多,因此在使用时应权衡考虑计算精度和时间。图3是一悬臂梁分别用线性和二次三角形单元离散时,其顶端位移随网格数量的收敛情况。可以看出,但网格数量较少时,两种单元的计算精度相差很大,这时采用低阶单元是不合适的。当网格数量较多时,两种单元的精度相差并不很大,这时采用高阶单元并不经济。例如在离散细节时,由于细节尺寸限制,要求细节附近的网格划分很密,这时采用线性单元更合适。
图3 不同阶次单元的收敛情况
增加网格数量和单元阶次都可以提高计算精度。因此在精度一定的情况下,用高阶单元离散结构时应选择适当的网格数量,太多的网格并不能明显提高计算精度,反而会使计算时间大大增加。为了兼顾计算精度和计算量,同一结构可以采用不同阶次的单元,即精度要求高的重要部位用高阶单元,精度要求低的次要部位用低阶单元。不同阶次单元之间或采用特殊的过渡单元连接,或采用多点约束等式连接。
美国大雪4 网格质量
网格质量是指网格几何形状的合理性。质量好坏将影响计算精度。质量太差的网格甚至会中止计算。直观上看,网格各边或各个内角相差不大、网格面不过分扭
曲、边节点位于边界等份点附近的网格质量较好。网格质量可用细长比、锥度比、内角、翘曲量、拉伸值、边节点位置偏差等指标度量。划分网格时一般要求网格质量能达到某些指标要求。在重点研究的结构关键部位,应保证划分高质量网格,即
2
粤语教程使是个别质量很差的网格也会引起很大的局部误差。而在结构次要部位,网格质量可适当降低。当模型中存在质量很差的网格(称为畸形网格)时,计算过程将无法进行。图4是三种常见的畸形网格,其中a单元的节点交叉编号,b单元的内角大于180°,c单元的两对节点重合,网格面积为零。
图4 几种常见的畸形网格
5 网格分界面和分界点
结构中的一些特殊界面和特殊点应分为网格边界或节点以便定义材料特性、物理特性、载荷和位移约束条件。即应使网格形式满足边界条件特点,而不应让边界条件来适应网格。常见的特殊界面和特殊点有材料分界面、几何尺寸突变面、分布载荷分界线(点)、集中载荷作用点和位移约束作用点等。图5是具有上述几种界面的结构及其网格划分形式。
图5 特殊界面和特殊点网格划分
6 位移协调性
位移协调是指单元上的力和力矩能够通过节点传递相邻单元。为保证位移协调,一个单元的节点必须同时也是相邻单元的节点,而不应是内点或边界点。相邻单元的共有节点具有相同的自由度性质。否则,单元之间须用多点约束等式或约束单元进行约束处理。图6是两种位移不协调的网格划分,图a中的节点1仅属于一个单元,变形后会产生材料裂缝或重叠。图b中的平面单元和梁单元节点的自由度性质不同,粱单元的力矩无法传递到平面单元。
图6 位移不协调的网格划分3
河对岸7 网格布局
当结构形状对称时,其网格也应划分对称网格,以使模型表现出相应的对称特性(如集中质矩阵对称)。不对称布局会引起一定误差,如在图7中,悬臂粱截面相对y轴对称,在对称载荷作用下,自由端两对称节点1、2的挠度值本应相等。
但若分图b所示的不对称网格,计算出的y1=0.0346,y2=0.0350。若改用图c所示的网格,则y1和y2完全相同。
图7 网格布局对计算结果的影响
阿切尔8 节点和单元编号
节点和单元的编号影响结构总刚矩阵的带宽和波前数,因而影响计算时间和存储容量的大小,因此合理的编号有利于提高计算速度。但对复杂模型和自动分网而言,人为确定合理的编号很困难,目前许多有限元分析软件自带有优化器,网格划分后可进行带宽和波前优化,从而减轻人的劳动强度。
法国餐桌礼仪
45、分析完之后想查看之前的加载情况:PlotCtrl>Symbol,在其中选择要选择选择要显示的项目即可
46、catia装配图导入ansys之前,先通过Tools>Generate CATPart from Product 将Product转化为Part,如果将Part导入ansys后丢失元素则需在Part中进行一下布尔相加运算再往ansys里导入,此外在装配图转零件图之前应对实体进
行一下修改如去掉一些小孔、倒角或圆角,以便于ansys中进合理的网格划分
孙悟空原型提高分析精度;
装配体导入ansys后(多个体)划分网格有三种方式:a、先GLUE,之后对每个体划分网格,粘接导致不能划分网格的利用连接(Concatenate)命令后再划分,此法粘接后对体映射或扫略划分网格有时不是很理想,对于各个体自由划分的装配体此法还是比较方便的;b、先add将所有体合为一体,再切割或不切割以实现映射划分网格;c、先对每个体单独进行网格划分(映射、扫略或自由都可以),之后利用约束耦合将各个体约束耦合在一起,常用到CEINTF命令,
此约束命令比较方便,比较常用,应熟练掌握。
47、对复杂形状的体由面分割为几个立方体进行映射(mapped)网格划分时,注意可能会出现立方体不能映射划分的情况,原因应该是四方体的某些面上不具有相同的性质,可能四方体一个面由于和多个分体接触其实已被分割成个多个面,虽是四方体但并不是面对面而是面对多面,所以映射不了;其中一个比较麻烦的解决方法是把立方体切成很多小块,每块都满足映射网格划分的条件;如果不宜进行Bool“加”运算而有元素丢失的话,可以自己补上丢失的元素,另外转化为了零件但各个体是分离的,可以进行“粘接(Glue)“运算再划分网格;装配图转零件之后(不进行Bool“加”运算)有时会多出很多线(多余的线),如一个立方体

本文发布于:2023-07-26 12:23:32,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/fanwen/fan/82/1118165.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:网格   单元   划分   计算   结构
相关文章
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图