预览加载中,请您耐心等待几秒...
1/5
2/5
3/5
4/5
5/5

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

有限元分析中的网格划分好坏直接关系到模型计算的准确性。本文简述了网格划 分应用的基本理论,并以ANSYS限元分析中的网格划分为实例对象,详细讲述了 网格划分基本理论及其在工程中的实际应用,具有一定的指导意义。 1引言 ANSYS有限元网格划分是进行数值模拟分析至关重要的一步,它直接影 响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、 单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素.从几何 表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。同理,平面 应力和平面应变情况设计的单元求解方程也不相同。在有限元数值求解中,单元 的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、 板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采 用辛普生(Simpson)积分。辛普生积分点的间隔是一定的,沿厚度分成奇数积分 点。由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应 用中,一定要采用合理的单元来模拟求解。 2ANSYS网格划分的指导思想 ANSYS网格划分的指导思想是首先进行总体模型规划,包括物理模型的 构造、单元类型的选择、网格密度的确定等多方面的内容。在网格划分和初步求 解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。为提高 求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴 对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度.利 用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则 应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精 度和效率,例如,轴对称场合多采用柱坐标系。有限元分析的精度和效率与单元 的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免 网格的畸形。在网格重划分过程中常采用曲率控制、单元尺寸与数量控制、穿透 控制等控制准则。在选用单元时要注意剪力自锁、沙漏和网格扭曲、不可压缩材 料的体积自锁等问题 ANSYS软件平台提供了网格映射划分和自由适应划分的策略。映射划分 用于曲线、曲面、实体的网格划分方法,可使用三角形、四边形、四面体、五面 体和六面体,通过指定单元边长、网格数量等参数对网格进行严格控制,映射划 分只用于规则的几何图素,对于裁剪曲面或者空间自由曲面等复杂几何体则难以 控制。自由网格划分用于空间自由曲面和复杂实体,采用三角形、四边形、四面 体进行划分,采用网格数量、边长及曲率来控制网格的质量。 3ANSYS网格划分基本原则 3。1网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小.一般来讲,网 格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数 量时应权衡两个因数综合考虑。 图1位移精度和计算时间随网格数量的变化 图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2 代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计 算精度明显提高,而计算时间不会有大的增加.当网格数量增加到一定程度后,再 继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网 格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果 相差较大,可以继续增加网格,相反则停止计算。 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是 计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同 的情况下应取相对较多的网格.同样在响应计算中,计算应力响应所取的网格数 应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态, 可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分 析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网 格。 3。2网格疏密 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算 数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好 地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的 部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏 密不同的网格划分形式。下面通过实例给出网格疏密对计算精度的影响. 图2较粗网格的有限元模型图3图2网格对应得环向应力云图 图4缺口处较细网格图5较密网格所得的环向应力云图 图2是中心带圆孔方板的对称模型,其网格划分反映了疏密不同的划分 原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小, 网格分得较稀。其中图3中在缺口处网格划分较疏;而图4种在缺口处的网格划 分较密。其应力计算结果:图4