文档库 最新最全的文档下载
当前位置:文档库 › 09 ABAQUS-CAE网格划分

09 ABAQUS-CAE网格划分

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

有限元网格划分的基本原则 划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。为建立正确、合理的有限元模型,这里介绍划分网格时应考虑的一些基本原则。 1 网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。 图1 位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。 2 网格疏密 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图b中网格疏密相差更大,它比图a中的网格少48个,但计算出的孔缘最大应力相差1%,而计算时间却减小了36%。由此可见,采用疏密不同的网格划分,既可以保持相当的计算精度,又可使网格数量减

网格划分的几种基本处理方法

网格划分的几种基本处理方法 贴体坐标法: 贴体坐标是利用曲线坐标,并使其坐标线与燃烧室外形或复杂计算区域边界重合,这样所有边界点能够用网格点来表示,不需要任何插值。一旦贴体坐标生成通过变换,偏微分方程求解可以不在任意形状的物理平面上,而在矩形或矩形的组合(空间问题求解域为长方体或它们的组合)转换平面上进行。这样计算与燃烧室外形无关,也与在物理平面上网格间隔无关。 而是把边界条件复杂的问题转换成一个边界条件简单的问题;这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差,而且还可节省计算时间和内存,使流场计算较准确,同时方便求解,较好地解决了复杂形状流动区域的计算,在工程上比较广泛应用。 区域法: 虽然贴体坐标系可以使坐标线与燃烧室外形相重合,从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域,很难用一种网格来模拟,生成单域贴体网格,即使生成了也不能保证网格质量,影响流场数值求解的效果。因此,目前常采用区域法或分区网格,其基本思想是,根据外形特点把复杂的物理域或复杂拓扑结构的网格,分成若干个区域,分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格,或结构块网格。对区域进行分区时,若相邻两个子域分离边界是协调对接,称为对接网格;若相邻两子域有相互重叠部分,则此分区网格称为重叠网格。根据实际数值模拟计算的需要,把整个区域(燃烧室)分成几个不同的子区域,并分别生成网格。这样不仅可提高计算精度,而且还可节省计算机内存,提高收敛精度。但是计算时,必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法,其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程.设压力校正值在最初迭代时为零,为了保证流量连续各个区域应同时求解,然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区,两个区域都对重叠区进行计算,重叠区一边区域内的值,要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系,实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同,要求重叠区二边流量相等。区域法能合理解决网格生成问题,已被大量用来计算复杂形状区域流动。 区域分解法: 对于复杂几何形状的实际燃烧装置,为了保证数值求解流场质量,目前常采用区域分解法。该法基本要点是:根据燃烧室形状特点和流场计算需要,把计算区域分成一个主区域和若干个子区域,对各个区域(块)分别建立网格,并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立,各个子区域大小也尽可能相同,使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同,通常在变量变化梯度大的区域,可以布置较细网格,并采用高阶紊流模型和描述复杂反应的紊流燃烧模型,以便更合理模拟实际流场。对于变量变化不太大区域,可采用较疏的网格和较简单的数学模型,这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑,相邻子区域连接重叠网格或对接网格来实现,在各子区域交界处通过插值法提供各子域求解变量的信息传递,满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和耦合,从而取得全流场解。 非结构网格法: 上述各方法所生成的网格均属于结构化网格,其共同特点是网格中各节点排列有序,每个节点与邻点之间关系是固定的,在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验,计算技术也较成熟,目前被广泛用来构造复杂外形区域

在ANSYS平台上的复杂有限元网格划分技术

在ANSYS平台上的复杂有限元网格划分技术 1. 网格密度 有限元结构网格数量的多少将直接影响计算结果的精度和计算规模的大小。一般来说,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,怎样在这两者之间找到平衡,是每一个CAE工作者都想拥有的技术。网格较少时,增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高很少,而计算时间却大幅度增加。所以应该注意网格数量的经济性。实际应用时,可以比较两种网格划分的计算结果,如果两次计算结果相差较大,应该继续增加网格,重新计算,直到结果误差在允许的范围之内。 在决定网格数量时还应该考虑分析类型。静力分析时,如果仅仅是计算结构的变形,网格数量可以少一点。如果需要计算应力,则在精度要求相同的情况下取相对较多的网格。同样在结构响应计算中,计算应力响应所取的网格数量应该比计算位移响应的多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选取较少的网格,如果计算的阶数较高,则网格数量应该相应的增加。在热分析中,结构内部的温度梯度不大时,不需要大量的内部单元,否则,内部单元应该较多。 有限元分析原则是把结构分解成离散的单元,然后组合这些单元

解得到最终的结果。其结果的精度取决于单元的尺寸和分布,粗的网格往往其结果偏小,甚至结果会发生错误。所以必须保证单元相对足够小,考虑到模型的更多的细节,使得到的结果越接近真实结果。由于粗的网格得到的结果是非保守的,因此要认真查看结果,其中有几种方法可以帮助读者分析计算结果与真实结果之间的接近程度。 最常用的方法是用对结果判断的经验来估计网格的质量,以确定网格是否合理,如通过看云图是否与物理现象相一致,如果云图线沿单元的边界或与实际现象不一致,那么很有可能结果是不正确的。 更多的评价网格误差的方法是通过比较平均的节点结果和不平均的单元结果。如在ANSYS中,提供了两条显示结果的命令:PLNS,PLES。前者是显示平均的节点结果,后者是显示不平均的单元结果。PLNS命令是计算节点结果,它是通过对该节点周围单元结果平均后得到的,分析结果是基于单元高斯积分点值,然后外插得到每个节点,因此在给定节点周围的每个单元都由自己的单元计算得到,所以这些节点结果通常是不相同的。PLNS命令是在显示结果之前将每个节点的所有结果进行了平均,所以看到的云图是以连续的方式从一个单元过渡到另外一个单元。而PLES命令不是对节点结果平均,所以在显示云图时单元和单元之间是不连续的。这种不连续程度在网格足够密(即单元足够小)的时候会很小或不存在,而在网格较粗时很大。由于PLNS结果是一个平均值,所以它得到的结果会比PLES的结果小,他

网格划分的方法

网格划分的方法 1.矩形网格差分网格的划分方法 划分网格的原则: 1)水域边界的补偿。舍去面积与扩增面积相互抵消。2)边界上的变步长处理。 3)水、岸边界的处理。 4)根据地形条件的自动划分。 5)根据轮廓自动划分。

2.有限元三角网格的划分方法 1)最近点和稳定结构原则。 2)均布结点的网格自动划分。 3)逐渐加密方法。 35 30 25 20 15 10 5 05101520253035

距离(m)距 离 (m) 3. 有限体积网格的划分方法 1) 突变原则。 2) 主要通道边界。 3) 区域逐步加密。

距离(100m) 离距(100m )距离(100m)离距(100m )

4. 边界拟合网格的划分方法 1) 变换函数:在区域内渐变,满足拉普拉斯方程的边值问题。 ),(ηξξξP yy xx =+ ),(ηξηηQ yy xx =+ 2) 导数变化原则。 ?????? ??????=?????? ??????-ηξ1J y x ,???? ??=ηηξξy x y x J 为雅可比矩阵,??? ? ??--=-ηηξξy x y x J J 11, ξηηξy x y x J -= )22(1 222233ηηξηξηηξηξξηηηηηξξηηξξξηξy y x y y y x y y x x y y x y y x y J xx +-+-+-= 同理可得yy ξ,xx η,yy η。 变换方程为 020222=+++-=+++-)()(ηξηηξηξξηξηηξηξξγβαγβαQy Py J y y y Qx Px J x x x 其中2222,,ξξηξξηηηγβαy x y y x x y x +=+=+=。

ANSYS 13.0 Workbench 网格划分及操作案例

第 3章 ANSYS 13.0 Workbench网格划分及操作案例 网格是计算机辅助工程(CAE)模拟过程中不可分割的一部分。网格直接影响到求解精 度、求解收敛性和求解速度。此外,建立网格模型所花费的时间往往是取得 CAE 解决方案所 耗费时间中的一个重要部分。因此,一个越好的自动化网格工具,越能得到好的解决方案。 3.1 ANSYS 13.0 Workbench 网格划分概述 ANSYS 13.0 提供了强大的自动化能力,通过实用智能的默认设置简化一个新几何体的网 格初始化,从而使得网格在第一次使用时就能生成。此外,变化参数可以得到即时更新的网 格。ANSYS 13.0 的网格技术提供了生成网格的灵活性,可以把正确的网格用于正确的地方, 并确保在物理模型上进行精确有效的数值模拟。 网格的节点和单元参与有限元求解,ANSYS 13.0在求解开始时会自动生成默认的网格。 可以通过预览网格,检查有限元模型是否满足要求,细化网格可以使结果更精确,但是会增 加 CPU 计算时间和需要更大的存储空间,因此需要权衡计算成本和细化网格之间的矛盾。在 理想情况下,我们所需要的网格密度是结果随着网格细化而收敛,但要注意:细化网格不能 弥补不准确的假设和错误的输入条件。 ANSYS 13.0 的网格技术通过 ANSYS Workbench的【Mesh】组件实现。作为下一代网格 划分平台, ANSYS 13.0 的网格技术集成 ANSYS 强大的前处理功能, 集成 ICEM CFD、 TGRID、 CFX-MESH、GAMBIT网格划分功能,并计划在 ANSYS 15.0 中完全整合。【Mesh】中可以根 据不同的物理场和求解器生成网格,物理场有流场、结构场和电磁场,流场求解可采用 【Fluent】、【CFX】、【POLYFLOW】,结构场求解可以采用显式动力算法和隐式算法。不同的 物理场对网格的要求不一样,通常流场的网格比结构场要细密得多,因此选择不同的物理场, 也会有不同的网格划分。【Mesh】组件在项目流程图中直接与其他 Workbench分析系统集成。 3.2 ANSYS 13.0 Workbench 网格划分 ANSYS 网格划分不能单独启动,只能在 Workbench 中调用分析系统或【Mesh】组件启 动,如图 3-1 所示。 图3-1 调入分析系统及网格划分组件

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

ANSYS有限元网格划分的基本原则 引言 ANSYS中有两种建立有限元模型的方法:实体建模和直接生成。使用实体建模,首先生成能描述模型的几何形状的几何模型,然后由ANSYS程序按照指定的单元大小和形状对几何体进行网格划分产生节点和单元。对于直接生成法,需要手工定义每个节点的位置和单元的连接关系。 一般来说对于规模较小的问题才适于采用直接生成法,常见的问题都需要先通过实体建模生成几何模型,然后再对其划分网格生成有限元模型。随着计算机性能的提高,分析模型的复杂性和规模都越来越大,而直接生成法也因其自身的局限性逐渐的被淘汰,所以正确的理解划分网格的目的和掌握划分网格的方法不论是对ANSYS的学习还是对二次开发都有重要的作用,尤其是当模型复杂度大,对模型的某些部分网格需要特殊处理时,这种对划分网格深度的理解作用更加明显。 2 常用高级网格划分方法 随着ANSYS功能的越来越强大和计算机性能的飞速提高,有限元分析向着大型化、复杂化的方向发展,而划分网格的观念也需要逐渐从二维模型向三维模型上上转变。这里主要描述三种常见的高级划分网格的方法,正确的理解和掌握这些划分网格的思想对于二次开发者来说非常的重要。 1)延伸网格划分 延伸网格划分是指将一个二维网格延伸生成一个三维网格;三维网格生成后去掉二维网格,延伸网格划分的步骤大体包括:先生成横截面、指定网格密度并对面进行网格划分、拖拉面网格生成体网格、指定单元属性、拖拉、完成体网格划分、释放已选的平面单元。 这里通过一个延伸网格划分的简单例子来加深对这种网格划分的理解。 图1 延伸网格划分举例 建立如图1所示的三维模型并划分网格,我们可以先建立z方向的端面,然后划分网格,通过拖拉的方法在z方向按照图中所示尺寸要求的三维模型,只需

_基于ANSYS的有限元法网格划分浅析

文章编号:1003-0794(2005)01-0038-02 基于ANSYS的有限元法网格划分浅析 杨小兰,刘极峰,陈 旋 (南京工程学院,南京210013) 摘要:为提高有限元数值的计算精度和对复杂结构力学分析的准确性,针对不同分析类型采用了不同的网格划分方法,结合实例阐述了ANSYS有限元网格划分的方法和技巧,指出了采用ANSYS有限元软件在网格划分时应注意的技术问题。 关键词:ANSYS;有限元;网格;计算精度 中图号:O241 82;TP391 7文献标识码:A 1 引言 ANSYS有限元分析程序是著名的C AE供应商美国ANSYS公司的产品,主要用于结构、热、流体和电磁四大物理场独立或耦合分析的CAE应用,功能强大,应用广泛,是一个便于学习和使用的优秀有限元分析程序。在ANSYS得到广泛应用的同时,许多技术人员对ANSYS程序的了解和认识还不够系统全面,在工作和研究中存在许多隐患和障碍,尤为突出的是有限元网格划分技术。本文结合工程实例,就如何合理地进行网格划分作一浅析。 2 网格划分对有限元法求解的影响 有限元法的基本思想是把复杂的形体拆分为若干个形状简单的单元,利用单元节点变量对单元内部变量进行插值来实现对总体结构的分析,将连续体进行离散化即称网格划分,离散而成的有限元集合将替代原来的弹性连续体,所有的计算分析都将在这个模型上进行。因此,网格划分将关系到有限元分析的规模、速度和精度以及计算的成败。实验表明:随着网格数量的增加,计算精确度逐渐提高,计算时间增加不多;但当网格数量增加到一定程度后,再继续增加网格数量,计算精确度提高甚微,而计算时间却大大增加。在进行网格划分时,应注意网格划分的有效性和合理性。 3 网格划分的有效性和合理性 (1)根据分析数据的类型选择合理的网格划分数量 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格。如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,可划分较少的网格。 (2)根据分析数据的分布特点选择合理的网格疏密度 在决定网格疏密度时应考虑计算数据的分布特点,在计算固有特性时,因为固有频率和振型主要取决于结构质量分布和刚度分布,采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差很大,可减小数值计算误差。同样,在结构温度场计算中也趋于采用均匀的网格形式。在计算数据变化梯度较大的部位时,为了更好地反映数据变化规律,需要采用比较密集的网格,而在计算数据变化梯度较小的部位,为了减小模型规模,则应划分相对稀疏的网格,这样整个结构就表现出疏密不同的网格划分形式。 以齿轮轮齿的有限元分析模型为例,由于分析的目的是求出齿轮啮合传动过程中齿根部分的弯曲应力,因此,分析计算时并不需要对整个齿轮进行计算,可根据圣文男原理将整个区域缩小到直接参与啮合的轮齿。虽然实际上参与啮合的齿数总大于1,但考虑到真正起作用的是单齿,通常只取一个轮齿作为分析对象,这样作可以大大节省计算机内存。考虑到轮齿应力在齿根过渡圆角和靠近齿面处变化较大,网格可划分得密一些。在进行疏密不同网格划分操作时可采用ANSYS提供的网格细化工具调整网格的疏密,也可采用分块建模法设置网格疏密度。 图1所示即为采用分块建模法进行网格划分。图1(a)为内燃机中重要运动零件连杆的有限元应力分析图,由于连杆结构对称于其摆动的中间平面,其厚度方向的尺寸远小于长度方向的尺寸,且载荷沿厚度方向近似均匀分布,故可按平面应力分析处 38 煤 矿 机 械 2005年第1期

网格划分

有限元网格划分 摘要:总结近十年有限元网格划分技术发展状况。首先,研究和分析有限元网格划分的基本原则;其次,对当前典型网格划分方法进行科学地分类,结合实例,系统地分析各种网格划分方法的机理、特点及其适用范围,如映射法、基于栅格法、节点连元法、拓扑分解法、几何分解法和扫描法等;再次,阐述当前网格划分的研究热点,综述六面体网格和曲面网格划分技术;最后,展望有限元网格划分的发展趋势。 关键词:有限元网格划分;映射法;节点连元法;拓扑分解法;几何分解法;扫描法;六面体网格 1 引言 有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。 2 有限元网格划分的基本原则 有限元方法的基本思想是将结构离散化,即对连续体进行离散化,利用简化几何单元来近似逼近连续体,然后根据变形协调条件综合求解。所以有限元网格的划分一方面要考虑对各物体几何形状的准确描述,另一方面也要考虑变形梯度的准确描述。为正确、合理地建立有限元模型,这里介绍划分网格时应考虑的一些基本原则。 2.1 网格数量

网格数量直接影响计算精度和计算时耗,网格数量增加会提高计算精度,但同时计算时耗也会增加。当网格数量较少时增加网格,计算精度可明显提高,但计算时耗不会有明显增加;当网格数量增加到一定程度后,再继续增加网格时精度提高就很小,而计算时耗却大幅度增加。所以在确定网格数量时应权衡这两个因素综合考虑。 2.2 网格密度 为了适应应力等计算数据的分布特点,在结构不同部位需要采用大小不同的网格。在孔的附近有集中应力,因此网格需要加密;周边应力梯度相对较小,网格划分较稀。由此反映了疏密不同的网格划分原则:在计算数据变化梯度较大的部位,为了较好地反映数据变化规律,需要采用比较密集的网格;而在计算数据变化梯度较小的部位,为减小模型规模,网格则应相对稀疏。 2.3 单元阶次 单元阶次与有限元的计算精度有着密切的关联,单元一般具有线性、二次和三次等形式,其中二次和三次形式的单元称为高阶单元。高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且高次插值函数可更高精度地逼近复杂场函数,所以增加单元阶次可提高计算精度。但增加单元阶次的同时网格的节点数也会随之增加,在网格数量相同的情况下由高阶单元组成的模型规模相对较大,因此在使用时应权衡考虑计算精度和时耗。 2.4 单元形状 网格单元形状的好坏对计算精度有着很大的影响,单元形状太差的网格甚至会中止计算。单元形状评价一般有以下几个指标: (1)单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准。 (2)扭曲度:单元面内的扭转和面外的翘曲程度。 (3)节点编号:节点编号对于求解过程中总刚矩阵的带宽和波前因数有较大的影响,从而影响计算时耗和存储容量的大小

网格划分实例详细步骤

一个网格划分实例的详解 该题目条件如下图所示: Part 1:本部分将平台考虑成蓝色的虚线 1. 画左边的第一部分,有多种方案。 方法一:最简单的一种就是不用布置任何初始的2dmesh直接用one volume 画,画出来的质量相当不错。 One volume是非常简单而且强大的画法,只要是一个有一个方向可以 mapped的实体都可以用这个方法来画网格,而事实上,很多不能map的单元也都可以用这个命令来画,所以在对三维实体进行网格划分的时候,收件推荐用one volume来试下效果,如果效果不错的话,就没有必要先做二维单元后再来画。 方法二:先在其一个面上生成2D的mesh,在来利用general选项,这样的优点是可以做出很漂亮的网格。

相比之下:方法二所做出来的网格质量要比一要高。 2. 画第二段的网格,同样演示两种方法: 方法一:直接用3D>solid map>one volume 方法二:从该段图形来看,左端面实际上由3个面组成,右端面由一个部分组成,故可以先将左端面的另两个部分的面网格补齐,再用general选项来拉伸,但是,问题是左面砖红色的部分仅为3D单元,而没有可供拉伸的源面网格,故,应该先用face命令生成二维网格后,再来拉伸,其每一步的结果分见下:

在用general选项时,有个问题需要注意:在前面我们说过,source geom和elemes to drag二选一都可以,但是这里就不一样了,因为source geom选面的话,只能选择一个面,而此处是3个面,所以这里只能选elemes to drag而不能选择source geom.

CATIA有限元高级划分网格教程

CATIA有限元高级网格划分教程 盛选禹李明志 1.1进入高级网格划分工作台 (1)打开例题中的文件Sample01.CATPart。 (2)点击主菜单中的【开始】→【分析与模拟】→【Advanced Meshing Tools】(高级网格划分工具),就进入【Advanced Meshing Tools】(高级网格划分工具)工作台,如图1-1所示。进入工作台后,生成一个新的分析文件,并且显示一个【New Analysis Case】(新分析算题)对话框,如图1-2所示。 图1-1【开始】→【分析与模拟】→【Advanced Meshing Tools】(高级网格划分工具)(3)在【New Analysis Case】(新分析算题)对话框内选择【Static Analysis】(静力分析)选项。如果以后打开该对话框的时候均希望是计算静力分析,可以把对话框内的【Keep as default starting analysis case】(在开始时保持为默认选项)勾选。这样,下次进入本工作台时,将自动选择静力分析。 (4)点击【新分析算题】对话框内的【确定】按钮,关闭对话框。 1.2定义曲面网格划分参数 本节说明如何定义一个曲面零件的网格类型和全局参数。 (1)点击【Meshing Method】(网格划分方法)工具栏内的【高级曲面划分】按钮

,如图1-3所示。需要在【Meshing Method】(网格划分方法)工具栏内点击中间按钮的下拉箭头才能够显示出【高级曲 面划分】按钮。 图1-2【New Analysis Case】(新分析算题)对话框图1-3【高级曲面划分】按钮

网格划分方法

网格划分的几种基本处理方法 学习2010-01-10 17:13:52 阅读48 评论0 字号:大中小 贴体坐标法: 贴体坐标是利用曲线坐标,并使其坐标线与燃烧室外形或复杂计算区域边界重合,这样所有边界点能够用网格点来表示,不需要任何插值。一旦贴体坐标生成通过变换,偏微分方程求解可以不在任意形状的物理平面上,而在矩形或矩形的组合(空间问题求解域为长方体或它们的组合)转换平面上进行。这样计算与燃烧室外形无关,也与在物理平面上网格间隔无关。 而是把边界条件复杂的问题转换成一个边界条件简单的问题;这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差,而且还可节省计算时间和内存,使流场计算较准确,同时方便求解,较好地解决了复杂形状流动区域的计算,在工程上比较广泛应 用。 区域法: 虽然贴体坐标系可以使坐标线与燃烧室外形相重合,从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域,很难用一种网格来模拟,生成单域贴体网格,即使生成了也不能保证网格质量,影响流场数值求解的效果。因此,目前常采用区域法或分区网格,其基本思想是,根据外形特点把复杂的物理域或复杂拓扑结构的网格,分成若干个区域,分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格,或结构块网格。对区域进行分区时,若相邻两个子域分离边界是协调对接,称为对接网格;若相邻两子域有相互重叠部分,则此分区网格称为重叠网格。根据实际数值模拟计算的需要,把整个区域(燃烧室)分成几个不同的子区域,并分别生成网格。这样不仅可提高计算精度,而且还可节省计算机内存,提高收敛精度。但是计算时,必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法,其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程.设压力校正值在最初迭代时为零,为了保证流量连续各个区域应同时求解,然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区,两个区域都对重叠区进行计算,重叠区一边区域内的值,要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系,实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同,要求重叠区二边流量相等。区域法能合理解决网格生成问题,已被大量用来计算复杂形状区域流动。 区域分解法: 对于复杂几何形状的实际燃烧装置,为了保证数值求解流场质量,目前常采用区域分解法。该法基本要点是:根据燃烧室形状特点和流场计算需要,把计算区域分成一个主区域和若干个子区域,对各个区域(块)分别建立网格,并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立,各个子区域大小也尽可能相同,使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同,通常在变量变化梯度大的区域,可以布置较细网格,并采用高阶紊流模型和描述复杂反应的紊流燃烧模型,以便更合理模拟实际流场。对于变量变化不太大区域,可采用较疏的网格和较简单的数学模型,这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑,相邻子区域连接重叠网格或对接网格来实现,在各子区域交界处通过插值法提供各子域求解变量的信息传递,满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和 耦合,从而取得全流场解。 非结构网格法: 上述各方法所生成的网格均属于结构化网格,其共同特点是网格中各节点排列有序,每个节点与邻点之间关系是固定的,在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验,计算技术也较成熟,目前被广泛用来构造复杂外形区域内网格。但是,若复杂外形稍有改变,则将需要重新划分区域和构造网格,耗费较多人力和时间。为此,近年来又发展了另一类网格——非结构网格。此类网格的基本特点是:任何空间区域都被以四面体为单元的网格所划分,网格节点不受结构性质限制,能较好地处理边界,每个节点的邻点个数也可不固定,因此易于控制网格单元的大小、形状及网格的位置。与结构网格相比,此类网格具有更大灵活性和对复杂外形适应性。在20世纪80年代末和90年代初,非结构网格得到了迅速发展。生成非结构网格方法主要有三角化方法和推进阵面法两种。虽然非结构网格容易适合复杂外形,但与结构网格相比还存在一些缺点:(1)需要较大内存记忆单元节点之

网格划分的技巧和策略

在中国CAE论坛上看到这个,挺不错的 壳体单元网格划分时,如果能了解一些网格划分的技巧和策略,将会事半功倍。壳体网格划分可以从3个方面入手:几何模型、划分方法和解决策略。 1 几何模型 可以从以下几个方面了解和处理几何模型问题 (1)了解部件的形状,主要集中在尺寸小的部分。 (2)什么样的特征可以被忽略,例如小的倒角和圆孔。 (3)何种特征对分析是关键的特征,这些特征对确保好的单元质量是需要的。 2 划分方法(自动+手工) 可以采用如下方法 (1)将部件分割为不同的区域。 (2)每个区域必须有可能只使用一种三维网格模式。 (3)寻找下述特点区域:大量生成区域、对称性区域、产生困难的区域。 (4)寻找大量不同区域和方法。 (5)注意什么样的二维网格模式被要求。 (6)观察周围区域:什么功能可以在那里使用。 (7)二维网格模式是否可以延伸到相邻区域中。 (8)寻找对网格模式不能处理位置进行网格划分的方法:如果这样做了,寻找网格可以触及的曲面;注意周围网格将与此模式相融合。 (9)小特征融入大特征中;大特征划分网格时必须考虑到小特征。 (10)注意网格模式。 3 解决策略 壳体网格划分的主要策略如下 (1)内部特征衔接外部特征: l 不能变成被限制的。 l 网格模式需要一个面流入以便它们可以停止 l 从内到外划分网格可以避免此问题。 (2)小特征融入到大特征中:注意模式、大特征划分网格时必须考虑到小特征。 (3)硬特征应当先处理,否则它们会变得难于处理。 (4)通常情况下首先进行大量的生成,后面的编辑是比较容易的。 某些区域比较重要的网格划分的质量要求高些,如力的作用区域,边界条件所在的区域。一些设计区域和离设计区域比较远的地方可以适当放宽要求,但是最好是一些网格性能指标要满足。

ANSYS结构有限元分析中的网格划分技术及其应用实例

一、前言 有限元网格划分是进行有限元数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。同理,平面应力和平面应变情况设计的单元求解方程也不相同。在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。 CAD软件中流行的实体建模包括基于特征的参数化建模和空间自由曲面混合造型两种 方法。Pro/E和SoildWorks是特征参数化造型的代表,而CATIA与Unigraphics等则将特征参数化和空间自由曲面混合造型有机的结合起来。现有CAD软件对表面形态的表示法已经大大超过了CAE软件,因此,在将CAD实体模型导入CAE软件的过程中,必须将CAD 模型中其他表示法的表面形态转换到CAE软件的表示法上,转换精度的高低取决于接口程序的好坏。在转换过程中,程序需要解决好几何图形(曲线与曲面的空间位置)和拓扑关系(各图形数据的逻辑关系)两个关键问题。其中几何图形的传递相对容易实现,而图形间的拓扑关系容易出现传递失败的情况。数据传递面临的一个重大挑战是,将导入CAE程序的CAD模型改造成适合有限元分析的网格模型。在很多情况下,导入CAE程序的模型可能包含许多设计细节,如细小的孔、狭窄的槽,甚至是建模过程中形成的小曲面等。这些细节往往不是基于结构的考虑,保留这些细节,单元数量势必增加,甚至会掩盖问题的主要矛盾,对分析结果造成负面影响。 CAD模型的“完整性”问题是困扰网格剖分的障碍之一。对于同一接口程序,数据传递的品质取决于CAD模型的精度。部分CAD模型对制造检测来说具备足够的精度,但对有限元网格剖分来说却不能满足要求。值得庆幸的是,这种问题通常可通过CAD软件的“完整性检查”来修正。改造模型可取的办法是回到CAD系统中按照分析的要求修改模型。一方面检查模型的完整性,另一方面剔除对分析无用的细节特征。但在很多情况下,这种“回归”很难实现,模型的改造只有依靠CAE软件自身。CAE中最直接的办法是依靠软件具有的“重构”功能,即剔除细部特征、缝补面和将小面“融入”大曲面等。有些专用接口在模型传递过程中甚至允许自动完成这种工作,并且通过网格剖分器检验模型的“完整性”,如发现“完整性”不能满足要求,接口程序可自动进行“完整性”修复。当几何模型距CAE分析的要求相差太大时,还可利用CAE程序的造型功能修正几何模型。“布尔运算”是切除细节和修理非完整特征的有效工具之一。 目前数据传递一般可通过专用数据接口,CAE程序可与CAD程序“交流”后生成与CAE 程序兼容的数据格式。另一种方式是通过标准图形格式如IGES、SAT和ParaSolid传递。现有的CAD平台与通用有限元平台一般通过IGES、STL、Step、Parasolid等格式来数据

器件网格划分方法的教学总结与归纳

器件网格划分方法的教学总结与归纳 一、前言微电子产业规模和技术水平是衡量国家综合实力的重要指标,在促进国民经济可持续性发展的同时,对国家安全战略的保护也有着重要的贡献。 [1] 积极培养掌握先进半导体知识与集成电路设计技术并符合企业需求的高端人才,是高等学校肩负的不可推卸的重要职责。在微电子相关课程体系教学过程中,引入半导体器件计算机模拟仿真技术,可以帮助学生理解抽象、复杂的基础理论,加强学生半导体技术实际应用能力的培养,实现理论教学与实践教学的紧密结合 [2] ,在一定程度上可以缓解教学投入与学校有限办学经费之间的矛盾。 要顺利开展半导体器件模拟仿真工作,首先面临所谓的网格划分问题。 [3] 网格划分指的是将非线性偏微分方程所描述的几何区域分割成有限个子区域的方法,把非线性偏微分方程的求解,简化为在更小单个子区域内线性方程组的求解。网格划分的优劣决定了方程求解速度的快慢,关系到数值求解是否能收敛及误差大小。在正确划分网格的基础上,越细致的网格,得到的数据与真实值的误差就越小,但仿真任务所需计算时间增加的就越快,所需计算硬件资源就越多越昂贵,甚至超出高等学校实际的办学条件。 半导体工艺及器件仿真工具 Sentaurus TCAD 是由 Synopsys 公司开发的最新软件,可以用来模拟集成器件的工艺制程、器件 物理特性和互连特性等,支持的仿真器件类型包括CMO、S 功率

器件、存储器、太阳能电池和光电探测器等,在高校微电子与半导体相关专业教学中逐渐得到了推广。 [4][5] 本文将以 Sentaurus软件对半导体PN结仿真模拟的任务为例,针对软件中SDE 模块中涉及的网格划分的主要内容与方法进行归纳整理,为相关课程的教学提供参考借鉴。 二、步骤与策略 网格的划分大致分为三个步骤:定义网格划分的策略,定义划分网格的区域,将网格划分的策略施加到相应区域上。这是 SDE中网格划分的基本的方法,当有部分区域没有被定义为网格划分区域时,将自动为该部分区域进行网格划分,但相对划分的部分会粗略许多。 对于网格区域的定义需要根据器件的结构和网格的划分策 略,SDE中提供了三种网格区域定义的方法:自定义窗口区域(Window)、通过选定器件的结构区域(Region)、通过选定器件的材料( Material )。三种网格区域定义的方法各有侧重,需要根据情况得当使用。 在网格区域划分的基础上便需要进行网格划分策略的选择,定义网格划分的策略也是网格划分过程中最核心的部分。软件中网格划分的基本思路是,在三个坐标轴方向上设定最小( Min)和最大划分因子(Max)的值,通过调节比例参数(Ratio ),改变最小因子到最大因子的变化速率(当比例参数为 1 时表示选用最小划分因子进行相应坐标轴上的划分)。按照各坐标轴的正方向由最小因子至最大因子的步长,并由最大因子的步长完成整个网格划分剩余的过程。在这里

有限元网格划分和收敛性

一、基本有限元网格概念 1.单元概述 几何体划分网格之前需要确定单元类型。 单元类型的选择应该根据分析类型、 形状特征、 计算数据特点、精度要求和计算的硬件条件等因素综合考虑。 为适应特殊的分析对象和边界 条件,一些问题需要采用多种单元进行组合建模。 2?单元分类 选择单元首先需要明确单元的类型,在结构中主要有以下一些单元类型: 平面应力单元、 平面应变单元、轴对称实体单元、空间实体单元、板 单元、壳单元、轴对称壳单元、杆单 元、梁单元、弹簧单元、间隙单元、质量单元、摩擦单元、刚体单元和约束单元等。根据不 同的分类方法,上述单元可以分成以 下不同的形式。 3. 按照维度进行单元分类 根据单元的维数特征,单元可以分为一维单元、二维单元和三维单元。 一维单元的网格为一条直线或者曲线。 直线表示由两个节点确定的线性单元。 曲线代表 由两个以上的节点确定的高次单元, 或者由具有确定形状的线性单元。 杆单元、梁单元和轴 对称壳单元属于一维单元,如图 1?图 3所示。 二维单元的网格是一个平面或者曲面,它没有厚度方向的尺寸。这类单元包括平面单元、 轴对称实体单元、板单元、壳单元和复合材料壳单元等,如图 4所示。二 维单元的形状通 常具有三角形和四边形两种, 在使用自动网格剖分时, 这类单元要求的几何形状是表面模型 图1捋果詰柯与一维杆单无犠型(直豉) &2桁舉第构石一隼杆早死撲型(曲线) B3毀姑构与一纯梁单元除世(直疑和呦疚〕

或者实体模型的边界面。采用薄壳单元通常具有相当好的计算效率。

洞丨伍金哉钩和潯壳社电 三维单元的网格具有空间三个方向的尺寸,其形状具有四面体、五面体和六面体,这类单元 包括空间实体单元和厚壳单元,如图5所示。在自动网格划分时,它要求的是几何模型是实 体模型(厚壳单元是曲面也可以)。 图5三址乙勺久和父侬草无 4. 按照插值函数进行单元分类 根据单元插值函数多项式的最高阶数多少,单元可以分为线性单元、二次单元、三次 单元和更高次的单元。 线性单元具有线性形式的插值函数,其网格通常只具有角节点而无边节点,网格边界为直线或者平面。这类单元的优点是节点数量少,在精度要求不高或者结果数据梯度不太大 的情况下,采用线性单元可以得到较小的模型规模。但是由于单元位移函数是线性的,单元 着应力突变,如图6所示。 S6錢41吕节点点单无fu节庖实体羊元

ANSYS 网格划分方法总结

(1) 网格划分定义:实体模型是无法直接用来进行有限元计算得,故需对它进行网格划分以生成有限元模型。有限元模型是实际结构和物质的数学表示方法。 在ANSYS中,可以用单元来对实体模型进行划分,以产生有限元模型,这个过程称作实体模型的网格化。本质上对实体模型进行网格划分也就是用一个个单元将实体模型划分成众多子区域。这些子区域(单元),是有属性的,也就是前面设置的单元属性。 另外也可以直接利用单元和节点生成有限元模型。 实体模型进行网格划分就是用一个个单元将实体模型划分成众多子区域(单元)。 (2)为什么我选用plane55这个四边形单元后,仍可以把实体模型划分成三角 形区域集合??? 答案:ansys为面模型的划分只提供三角形单元和四边形单元,为体单元只提供四面体单元和六面体单元。不管你选择的单元是多少个节点,只要是2D单元,肯定构成一个四边形或者是三角形,绝对没有五、六边形等特殊形状。网格划分也就是用所选单元将实体模型划分成众多三角形单元和四边形子区域。 见下面的plane77/78/55都是节点数目大于4的,但都是通过各种插值或者是合并的方式形成一个四边形或者三角形。 所以不管你选择什么单元,只要是对面的划分,meshtool上的划分类型设置就只有tri和quad两种选择。 如果这个单元只构成三角形,例如plane35,则无论你在meshtool上划分设置时tri还是quad,划分出的结果都是三角形。

所以在选用plane55单元,而划分的是采用tri划分时,就会把两个点合并为一个点。如上图的plane55,下面是plane单元的节点组成,可见每一个单元上都有两个节点标号相同,表明两个节点是重合的。 。 同样在采用plane77 单元,进行tri划分时,会有三个节点重合。这里不再一一列出。(3)如何使用在线帮助: 点击对话框中的help,例如你想了解plane35的相关属性,你可以

ANSYS 网格划分详细介绍

ANSYS 网格划分详细介绍 2008-09-27 18:01 众所周知,对于有限元分析来说,网格划分是其中最关键的一个步骤,网格划分的好坏直接影响到解算的精度和速度。在ANSYS中,大家知道,网格划分有三个步骤:定义单元属性(包括实常数)、在几何模型上定义网格属性、划分网格。在这里,我们仅对网格划分这个步骤所涉及到的一些问题,尤其是与复杂模型相关的一些问题作简要阐述。 一、自由网格划分 自由网格划分是自动化程度最高的网格划分技术之一,它在面上(平面、曲面)可以自动生成三角形或四边形网格,在体上自动生成四面体网格。通常情况下,可利用ANSYS的智能尺寸控制技术(SMARTSIZE命令)来自动控制网格的大小和疏密分布,也可进行人工设置网格的大小(AESIZE、LESIZE、KESIZE、ESIZE等系列命令)并控制疏密分布以及选择分网算法等(MOPT命令)。对于复杂几何模型而言,这种分网方法省时省力,但缺点是单元数量通常会很大,计算效率降低。同时,由于这种方法对于三维复杂模型只能生成四面体单元,为了获得较好的计算精度,建议采用二次四面体单元(92号单元)。如果选用的是六面体单元,则此方法自动将六面体单元退化为阶次一致的四面体单元,因此,最好不要选用线性的六面体单元(没有中间节点,比如45号单元),因为该单元退化后为线性的四面体单元,具有过刚的刚度,计算精度较差;如果选用二次的六面体单元(比如95号单元),由于其是退化形式,节点数与其六面体原型单元一致,只是有多个节点在同一位置而已,因此,可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元,减少每个单元的节点数量,提高求解效率。在有些情况下,必须要用六面体单元的退化形式来进行自由网格划分,比如,在进行混合网格划分(后面详述)时,只有用六面体单元才能形成金字塔过渡单元。对于计算流体力学和考虑集肤效应的电磁场分析而言,自由网格划分中的层网格功能(由LESIZE命令的LAYER1和LAYER2域控制)是非常有用的。 二、映射网格划分 映射网格划分是对规整模型的一种规整网格划分方法,其原始概念是:对于面,只能是四边形面,网格划分数需在对边上保持一致,形成的单元全部为四边形;对于体,只能是六面体,对应线和面的网格划分数保持一致;形成的单元全部为六面体。在ANSYS中,这些条件有了很大的放宽,包括: 1 面可以是三角形、四边形、或其它任意多边形。对于四边以上的多边形,必须用LCCAT命令将某些边联成一条边,以使得对于网格划分而言,仍然是三角形或四边形;或者用AMAP命令定义3到4个顶点(程序自动将两个顶点之间的所有线段联成一条)来进行映射划分。 2 面上对边的网格划分数可以不同,但有一些限制条件。 3 面上可以形成全三角形的映射网格。 4 体可以是四面体、五面体、六面体或其它任意多面体。对于六面以上的多面体,必须用ACCAT命令将某些面联成一个面,以使得对于网格划分而言,仍然是四、五或六面体。 5 体上对应线和面的网格划分数可以不同,但有一些限制条件。

相关文档
相关文档 最新文档