动网格模型
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
动网格模型
更新时间:2013-11-28
大纲:
1. 2. 动网格的相关知识介绍; 以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程;
3. 与动网格应用有关的参考文献; 4. 使用动网格进行计算的一些例子。
##1. 动网格的相关知识介绍 有关动网格基础方面的东西,请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节 ,这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是 预先定义的运动,即可以在计算前指定其速度或角速度;也可以是预先未做定义的运动,即 边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界 的变化情况自动完成。在使用动网格模型时,必须首先定义初始网格、边界运动的方式并指 定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运 动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域,则需要将它 们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组 合到各自的初始网格区域中。不同区域之间的网格不必是正则的,可以在模型设臵中用 FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注:一般来讲,在Fluent中使用动网格,基本上都要使用到UDF,所以你最好具备一定的C语言 编程基础。
2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算,即弹簧近似光滑模型(spring-based smoothing)、动态分层模型(dynamic layering)和局部重划模型(local remeshing)。 弹簧近似光滑模型 在弹簧近似光滑模型中,网格的边被理想化为节点间相互连接的弹簧。移动前的网格间距相当 于边界移动前由弹簧组成的系统处于平衡状态。在网格边界节点发生位移后,会产生与位 移成比例的力,力量的大小根据胡克定律计算。边界节点位移形成的力虽然破坏了弹簧系 统原有的平衡,但是在外力作用下,弹簧系统经过调整将达到新的平衡,也就是说由弹簧 连接在一起的节点,将在新的位臵上重新获得力的平衡。从网格划分的角度说,从边界节 点的位移出发,采用虎克定律,经过迭代计算,最终可以得到使各节点上的合力等于零的 、新的网格节点位臵,这就是弹簧光顺法的核心思想。 原则上弹簧光顺模型可以用于任何一种网格体系,但是在非四面体网格区域(二维非三角形) ,最好在满足下列条件时使用弹簧光顺方法: (1)移动为单方向。 (2)移动方向垂直于边界。 如果两个条件不满足,可能使网格畸变率增大。另外,在系统缺省设臵中,只有四面体网格( 三维)和三角形网格(二维)可以使用弹簧光顺法,如果想在其他网格类型中激活该模型 ,需要在dynamic-mesh-menu 下使用文字命令spring-on-all-shapes?,然后激活该选项即可。 动态层模型 对于棱柱型网格区域(六面体和或者楔形),可以应用动态层模型。动态层模型的中心思想是 根据紧邻运动边界网格层高度的变化,添加或者减少动态层,即在边界发生运动时,如果 紧邻边界的网格层高度增大到一定程度,就将其划分为两个网格层;如果网格层高度降低 到一定程度,就将紧邻边界的两个网格层合并为一个层:
如果网格层j扩大,单元高度的变化有一临界值: H_min>(1+alpha_s)*h_0 式中h_min为单元的最小高度,h_0为理想单元高度,alpha_s为层的分割因子。在满足上述条件 的情况下,就可以对网格单元进行分割,分割网格层可以用常值高度法或常值比例法。在 使用常值高度法时,单元分割的结果是产生相同高度的网格。在采用常值比例法时,网格 单元分割的结果是产生是比例为alpha_s的网格。 若对第j层进行压缩,压缩极限为: H_min 式中alpha_c为合并因子。在紧邻动边界的网格层高度满足这个条件时,则将这一层网格与 外面一层网格相合并。 动网格模型的应用有如下限制: (1)与运动边界相邻的网格必须为楔形或者六面体(二维四边形)网格。 (2)在滑动网格交界面以外的区域,网格必须被单面网格区域包围。 (3)如果网格周围区域中有双侧壁面区域,则必须首先将壁面和阴影区分割开,再用 滑动交界面将二者耦合起来。 (4)如果动态网格附近包含周期性区域,则只能用FLUENT 的串行版求解,但是如果周期性 区域被设臵为周期性非正则交界面,则可以用FLUENT 的并行版求解。
如果移动边界为内部边界,则边界两侧的网格都将作为动态层参与计算。如果在壁面上只有一 部分是运动边界,其他部分保持静止,则只需在运动边界上应用动网格技术,但是动网格 区与静止网格区之间应该用滑动网格交界面进行连接。 局部重划模型 在使用非结构网格的区域上一般采用弹簧光顺模型进行动网格划分,但是如果运动边界的位移 远远大于网格尺寸,则采用弹簧光顺模型可能导致网格质量下降,甚至出现体积为负值的 网格,或因网格畸变过大导致计算不收敛。为了解决这个问题,FLUENT 在计算过程中将 畸变率过大,或尺寸变化过于剧烈的网格集中在一起进行局部网格的重新划分,如果重新 划分后的网格可以满足畸变率要求和尺寸要求,则用新的网格代替原来的网格,如果新的 网格仍然无法满足要求,则放弃重新划分的结果。 在重新划分局部网格之前,首先要将需要重新划分的网格识别出来。FLUENT 中识别不合乎要 求网格的判据有二个,一个是网格畸变率,一个是网格尺寸,其中网格尺寸又分最大尺寸 和最小尺寸。在计算过程中,如果一个网格的尺寸大于最大尺寸,或者小于最小尺寸,或 者网格畸变率大于系统畸变率标准,则这个网格就被标志为需要重新划分的网格。在遍历 所有动网格之后,再开始重新划分的过程。局部重划模型不仅可以调整体网格,也可以调 整动边界上的表面网格。 需要注意的是,局部重划模型仅能用于四面体网格和三角形网格。在定义了动边界面以后,如 果在动边界面附近同时定义了局部重划模型,则动边界上的表面网格必须满足下列条件: (1)需要进行局部调整的表面网格是三角形(三维)或直线(二维)。 (2)将被重新划分的面网格单元必须紧邻动网格节点。 (3)表面网格单元必须处于同一个面上并构成一个循环。 (4)被调整单元不能是对称面(线)或正则周期性边界的一部分。
动网格的实现在FLUENT 中是由系统自动完成的。如果在计算中设臵了动边界,则FLUENT 会 根据动边界附近的网格类型,自动选择动网格计算模型。如果动边界附近采用的是四面体 网格(三维)或三角形网格(二维),则FLUENT 会自动选择弹簧光顺模型和局部重划模 型对网格进行调整。如果是棱柱型网格,则会自动选择动态层模型进行网格调整。在静止 网格区域则不进行网格调整。 动网格问题中对于固体运动的描述,是以固体相对于重心的线速度和角速度为基本参数加以定 义的。既可以用型函数定义固体的线速度和角速度,也可以用UDF 来定义这两个参数。同 时需要定义的是固体在初始时刻的位臵。 注:这一小节主要讲述了动网格的更新方法,最好能掌握,尤其是各种方法的适用范围,通常 来讲,在一个case中,我们使用的更新方法都是根据网格类型以及和要实现的运动来选择的 ,很多时候都是几种更新方法搭配起来使用的。 总结一下: 使用弹簧近似光滑法网格拓扑始终不变,无需插值,保证了计算精度。但弹簧近似光滑法不适 用于大变形情况,当计算区域变形较大时,变形后的网格会产生较大的倾斜变形,从而使 网格质量变差,严重影响计算精度。动态分层法在生成网格方面具有快速的优势,同时它 的应用也受到了一些限制。它要求运动边界附近的网格为六面体或楔形,这对于复杂外形 的流场区域是不适合的。使用局部网格重划法要求网格为三角形(二维)或四面体(三维 ),这对于适应复杂外形是有好处的,局部网格重划法只会对运动边界附近区域的网格起 作用。
3、动网格问题的建立 设臵动网格问题的步骤如下: (1)在Solver(求解器)面板中选择非定常流(unsteady)计算。 (2)设定边界条件,即设定壁面运动速度。 (3)激活动网格模型,并设定相应参数,菜单操作如下: Define -> Dynamic Mesh -> Parameters... (4)指定移动网格区域的运动参数,菜单操作如下: Define -> Dynamic Mesh -> Zones... (5)保存算例文件和数据文件。 (6)预览动网格设臵,菜单操作为: Solve -> Mesh Motion... (7)在计算活塞问题时,设定活塞计算中的事件: Define -> Dynamic Mesh -> Events... 并可以通过显示阀与活塞的运动,检查上述设臵是否正确: Display -> IC Zone Motion... (8)应用自动保存功能保存计算结果。 File -> Write -> Autosave... 在动网格计算中,因为每个计算步中网格信息都会改变,而网格信息是储存在算例文件中的,所以必 须同时保存算例文件和数据文件。 (9)如果想建立网格运动的动画过程,可以在Solution Animation(计算结果动画)面板中进行相关设 臵。 注:在这一步中,需要提醒一下,使用动网格进行正式计算之前,最好养成预览动网格更新的习惯; 就是在正式计算前,浏览一下动网格的更新情况,这样可以避免在计算过程中出现动网格更新本 身的问题。在预览更新时,很多人都说会出现负体积的警告,更新不成功,出现这样的问题时, 最好先把时间步长改的更小点儿试试,一般来讲,排除UDF本身的原因,出现更新出错的原因都 与时间步长有关,这需要结合所使用的更新方法多琢磨。
4、设定动网格参数 为了使用动网格模型,需要在dynamic mesh(动网格)面板中激活Dynamic Mesh(动网格)选项 。如果计算的是活塞运动,则同时激活In-Cylinder(活塞)选项。然后选择动网格模型,并 设臵相关参数。 1)选择网格更新模型 在Mesh Methods(网格划分方法)下面选择Smothing(弹簧光顺模型),Layering(动态层模型 )和(或)Remshing(局部重划模型)。 2)设臵弹簧光顺参数 激活弹簧光顺模型,相关参数设臵位于Smoothing(光顺)标签下,可以设臵的参数包括Spring Constant Factor(弹簧弹性系数)、Boundary Node Relaxation(边界点松弛因子)、 Convergence Tolerance(收敛判据)和Number of Iterations(迭代次数)。 弹簧弹性系数应该在0 到1 之间变化,弹性系数等于0 时,弹簧系统没有耗散过程,在图中算例 中,靠近壁面的网格没有被改变,而是保持了原来的网格形状和密度;在弹性系数等于1 时 ,弹簧系统的耗散过程与缺省设臵相同,从图中可以发现壁面发生变形,壁面附近网格因 为过度加密而质量下降。因此在实际计算中应该在0 到1 之间选择一个适当的值。边界点松 弛因子用于控制动边界上网格点的移动。当这个值为零时,边界节点不发生移动;在这个 值为1 时,则边界节点的移动计算中不采用松弛格式。在大多数情况下,这个值应该取为0 到1 之间的一个值,以保证边界节点以合适的移动量发生移动。 收敛判据就是网格节点移动计算中,迭代计算的判据。迭代次数是指网格节点移动计算的最大 迭代次数。 3)动态层 在Layering(动态层)标签下,可以设臵与动态层模型相关的参数。通过设定Constant Height(常值高度)与Constant Ratio(常值比例)可以确定分解网格的两种方法。Split Factor (分割因子)和Collapse Factor(合并因子)则分别为上面介绍的alpha_s和alpha_c。
4)局部重新划分网格 在Remeshing(重划网格)标签下,设臵与局部重划模型相关的参数。可以设臵的参数包括 Maximum Cell Skewness(最大畸变率)、Maximum Cell Volume(最大网格体积)和Minimum Cell Volume(最大网格体积),其含义如前所述,主要用于确定哪些网格需要被 重新划分。在缺省设臵中,如果重新划分的网格优于原网格,则用新网格代替旧网格;否 则,将保持原网格划分不变。如果无论如何都要采用新网格的话,则可以在Options(选项) 下面选择Must Improve Skewness(必须改善畸变率)选项。如果Options(选项)下面的Size Function(尺寸函数)被激活,则还可以用网格尺寸分布函数标志需要重新划分的网格。假 设在某点附近的理想网格尺寸为L ,而某个网格的尺寸为L' ,如果: L’不属于[0.8*gamma*L,1.25*gamma*L] 则网格被标志为需要重新划分的网格,并在随后的计算中被重新划分。式中的gamma用下面的 公式计算: 当alpha>0时,gamma=1+alpha*d_b^(1+2*beta) 当alpha<0时,gamma=1+alpha*d_b^[(1-beta)^-1] 式中d_b为网格到壁面边界的最小距离,alpha和beta就是需要设臵的Size Function Variation(尺寸 函数增量)和Size Function Rate(尺寸函数变化率)。 以上两式中,alpha以边界网格尺寸为基准控制内部网格的大小,alpha等于0.5表示内部网格的尺 寸至少是边界网格的1.5倍,alpha等于-0.5表示内部网格的尺寸等于边界网格的0.5倍,如果 alpha等于0则表示内部网格与边界网格大小相同;bata的取值在-1到1之间,用于控制网格从 边界到内部区域的变化速率。Beta取正值表示网格变化速率较慢取负值则表示变化速率较快 ,取0 则表示网格从边界到内部区域呈线性变化。 另外一个选项Size Function Resolution(尺寸函数分辨率)用于定义尺寸分布函数对网格大小的 分辨率,其参照值是最小网格尺寸。 点击Use defaults(使用缺省值)按钮,可以恢复系统缺省设臵。
更新时间:2013-11-28
大纲:
1. 2. 动网格的相关知识介绍; 以NACA0012翼型俯仰振荡实例进行讲解动网格的应用过程;
3. 与动网格应用有关的参考文献; 4. 使用动网格进行计算的一些例子。
##1. 动网格的相关知识介绍 有关动网格基础方面的东西,请具体参考FLUENT User’s Guide或FLUENT全攻略的相关章节 ,这里只给出一些提要性的知识要点。 1、简介 动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题。边界的运动形式可以是 预先定义的运动,即可以在计算前指定其速度或角速度;也可以是预先未做定义的运动,即 边界的运动要由前一步的计算结果决定。网格的更新过程由FLUENT 根据每个迭代步中边界 的变化情况自动完成。在使用动网格模型时,必须首先定义初始网格、边界运动的方式并指 定参予运动的区域。可以用边界型函数或者UDF 定义边界的运动方式。FLUENT 要求将运 动的描述定义在网格面或网格区域上。如果流场中包含运动与不运动两种区域,则需要将它 们组合在初始网格中以对它们进行识别。那些由于周围区域运动而发生变形的区域必须被组 合到各自的初始网格区域中。不同区域之间的网格不必是正则的,可以在模型设臵中用 FLUENT软件提供的非正则或者滑动界面功能将各区域连接起来。 注:一般来讲,在Fluent中使用动网格,基本上都要使用到UDF,所以你最好具备一定的C语言 编程基础。
2、动网格更新方法 动网格计算中网格的动态变化过程可以用三种模型进行计算,即弹簧近似光滑模型(spring-based smoothing)、动态分层模型(dynamic layering)和局部重划模型(local remeshing)。 弹簧近似光滑模型 在弹簧近似光滑模型中,网格的边被理想化为节点间相互连接的弹簧。移动前的网格间距相当 于边界移动前由弹簧组成的系统处于平衡状态。在网格边界节点发生位移后,会产生与位 移成比例的力,力量的大小根据胡克定律计算。边界节点位移形成的力虽然破坏了弹簧系 统原有的平衡,但是在外力作用下,弹簧系统经过调整将达到新的平衡,也就是说由弹簧 连接在一起的节点,将在新的位臵上重新获得力的平衡。从网格划分的角度说,从边界节 点的位移出发,采用虎克定律,经过迭代计算,最终可以得到使各节点上的合力等于零的 、新的网格节点位臵,这就是弹簧光顺法的核心思想。 原则上弹簧光顺模型可以用于任何一种网格体系,但是在非四面体网格区域(二维非三角形) ,最好在满足下列条件时使用弹簧光顺方法: (1)移动为单方向。 (2)移动方向垂直于边界。 如果两个条件不满足,可能使网格畸变率增大。另外,在系统缺省设臵中,只有四面体网格( 三维)和三角形网格(二维)可以使用弹簧光顺法,如果想在其他网格类型中激活该模型 ,需要在dynamic-mesh-menu 下使用文字命令spring-on-all-shapes?,然后激活该选项即可。 动态层模型 对于棱柱型网格区域(六面体和或者楔形),可以应用动态层模型。动态层模型的中心思想是 根据紧邻运动边界网格层高度的变化,添加或者减少动态层,即在边界发生运动时,如果 紧邻边界的网格层高度增大到一定程度,就将其划分为两个网格层;如果网格层高度降低 到一定程度,就将紧邻边界的两个网格层合并为一个层:
如果网格层j扩大,单元高度的变化有一临界值: H_min>(1+alpha_s)*h_0 式中h_min为单元的最小高度,h_0为理想单元高度,alpha_s为层的分割因子。在满足上述条件 的情况下,就可以对网格单元进行分割,分割网格层可以用常值高度法或常值比例法。在 使用常值高度法时,单元分割的结果是产生相同高度的网格。在采用常值比例法时,网格 单元分割的结果是产生是比例为alpha_s的网格。 若对第j层进行压缩,压缩极限为: H_min 式中alpha_c为合并因子。在紧邻动边界的网格层高度满足这个条件时,则将这一层网格与 外面一层网格相合并。 动网格模型的应用有如下限制: (1)与运动边界相邻的网格必须为楔形或者六面体(二维四边形)网格。 (2)在滑动网格交界面以外的区域,网格必须被单面网格区域包围。 (3)如果网格周围区域中有双侧壁面区域,则必须首先将壁面和阴影区分割开,再用 滑动交界面将二者耦合起来。 (4)如果动态网格附近包含周期性区域,则只能用FLUENT 的串行版求解,但是如果周期性 区域被设臵为周期性非正则交界面,则可以用FLUENT 的并行版求解。
如果移动边界为内部边界,则边界两侧的网格都将作为动态层参与计算。如果在壁面上只有一 部分是运动边界,其他部分保持静止,则只需在运动边界上应用动网格技术,但是动网格 区与静止网格区之间应该用滑动网格交界面进行连接。 局部重划模型 在使用非结构网格的区域上一般采用弹簧光顺模型进行动网格划分,但是如果运动边界的位移 远远大于网格尺寸,则采用弹簧光顺模型可能导致网格质量下降,甚至出现体积为负值的 网格,或因网格畸变过大导致计算不收敛。为了解决这个问题,FLUENT 在计算过程中将 畸变率过大,或尺寸变化过于剧烈的网格集中在一起进行局部网格的重新划分,如果重新 划分后的网格可以满足畸变率要求和尺寸要求,则用新的网格代替原来的网格,如果新的 网格仍然无法满足要求,则放弃重新划分的结果。 在重新划分局部网格之前,首先要将需要重新划分的网格识别出来。FLUENT 中识别不合乎要 求网格的判据有二个,一个是网格畸变率,一个是网格尺寸,其中网格尺寸又分最大尺寸 和最小尺寸。在计算过程中,如果一个网格的尺寸大于最大尺寸,或者小于最小尺寸,或 者网格畸变率大于系统畸变率标准,则这个网格就被标志为需要重新划分的网格。在遍历 所有动网格之后,再开始重新划分的过程。局部重划模型不仅可以调整体网格,也可以调 整动边界上的表面网格。 需要注意的是,局部重划模型仅能用于四面体网格和三角形网格。在定义了动边界面以后,如 果在动边界面附近同时定义了局部重划模型,则动边界上的表面网格必须满足下列条件: (1)需要进行局部调整的表面网格是三角形(三维)或直线(二维)。 (2)将被重新划分的面网格单元必须紧邻动网格节点。 (3)表面网格单元必须处于同一个面上并构成一个循环。 (4)被调整单元不能是对称面(线)或正则周期性边界的一部分。
动网格的实现在FLUENT 中是由系统自动完成的。如果在计算中设臵了动边界,则FLUENT 会 根据动边界附近的网格类型,自动选择动网格计算模型。如果动边界附近采用的是四面体 网格(三维)或三角形网格(二维),则FLUENT 会自动选择弹簧光顺模型和局部重划模 型对网格进行调整。如果是棱柱型网格,则会自动选择动态层模型进行网格调整。在静止 网格区域则不进行网格调整。 动网格问题中对于固体运动的描述,是以固体相对于重心的线速度和角速度为基本参数加以定 义的。既可以用型函数定义固体的线速度和角速度,也可以用UDF 来定义这两个参数。同 时需要定义的是固体在初始时刻的位臵。 注:这一小节主要讲述了动网格的更新方法,最好能掌握,尤其是各种方法的适用范围,通常 来讲,在一个case中,我们使用的更新方法都是根据网格类型以及和要实现的运动来选择的 ,很多时候都是几种更新方法搭配起来使用的。 总结一下: 使用弹簧近似光滑法网格拓扑始终不变,无需插值,保证了计算精度。但弹簧近似光滑法不适 用于大变形情况,当计算区域变形较大时,变形后的网格会产生较大的倾斜变形,从而使 网格质量变差,严重影响计算精度。动态分层法在生成网格方面具有快速的优势,同时它 的应用也受到了一些限制。它要求运动边界附近的网格为六面体或楔形,这对于复杂外形 的流场区域是不适合的。使用局部网格重划法要求网格为三角形(二维)或四面体(三维 ),这对于适应复杂外形是有好处的,局部网格重划法只会对运动边界附近区域的网格起 作用。
3、动网格问题的建立 设臵动网格问题的步骤如下: (1)在Solver(求解器)面板中选择非定常流(unsteady)计算。 (2)设定边界条件,即设定壁面运动速度。 (3)激活动网格模型,并设定相应参数,菜单操作如下: Define -> Dynamic Mesh -> Parameters... (4)指定移动网格区域的运动参数,菜单操作如下: Define -> Dynamic Mesh -> Zones... (5)保存算例文件和数据文件。 (6)预览动网格设臵,菜单操作为: Solve -> Mesh Motion... (7)在计算活塞问题时,设定活塞计算中的事件: Define -> Dynamic Mesh -> Events... 并可以通过显示阀与活塞的运动,检查上述设臵是否正确: Display -> IC Zone Motion... (8)应用自动保存功能保存计算结果。 File -> Write -> Autosave... 在动网格计算中,因为每个计算步中网格信息都会改变,而网格信息是储存在算例文件中的,所以必 须同时保存算例文件和数据文件。 (9)如果想建立网格运动的动画过程,可以在Solution Animation(计算结果动画)面板中进行相关设 臵。 注:在这一步中,需要提醒一下,使用动网格进行正式计算之前,最好养成预览动网格更新的习惯; 就是在正式计算前,浏览一下动网格的更新情况,这样可以避免在计算过程中出现动网格更新本 身的问题。在预览更新时,很多人都说会出现负体积的警告,更新不成功,出现这样的问题时, 最好先把时间步长改的更小点儿试试,一般来讲,排除UDF本身的原因,出现更新出错的原因都 与时间步长有关,这需要结合所使用的更新方法多琢磨。
4、设定动网格参数 为了使用动网格模型,需要在dynamic mesh(动网格)面板中激活Dynamic Mesh(动网格)选项 。如果计算的是活塞运动,则同时激活In-Cylinder(活塞)选项。然后选择动网格模型,并 设臵相关参数。 1)选择网格更新模型 在Mesh Methods(网格划分方法)下面选择Smothing(弹簧光顺模型),Layering(动态层模型 )和(或)Remshing(局部重划模型)。 2)设臵弹簧光顺参数 激活弹簧光顺模型,相关参数设臵位于Smoothing(光顺)标签下,可以设臵的参数包括Spring Constant Factor(弹簧弹性系数)、Boundary Node Relaxation(边界点松弛因子)、 Convergence Tolerance(收敛判据)和Number of Iterations(迭代次数)。 弹簧弹性系数应该在0 到1 之间变化,弹性系数等于0 时,弹簧系统没有耗散过程,在图中算例 中,靠近壁面的网格没有被改变,而是保持了原来的网格形状和密度;在弹性系数等于1 时 ,弹簧系统的耗散过程与缺省设臵相同,从图中可以发现壁面发生变形,壁面附近网格因 为过度加密而质量下降。因此在实际计算中应该在0 到1 之间选择一个适当的值。边界点松 弛因子用于控制动边界上网格点的移动。当这个值为零时,边界节点不发生移动;在这个 值为1 时,则边界节点的移动计算中不采用松弛格式。在大多数情况下,这个值应该取为0 到1 之间的一个值,以保证边界节点以合适的移动量发生移动。 收敛判据就是网格节点移动计算中,迭代计算的判据。迭代次数是指网格节点移动计算的最大 迭代次数。 3)动态层 在Layering(动态层)标签下,可以设臵与动态层模型相关的参数。通过设定Constant Height(常值高度)与Constant Ratio(常值比例)可以确定分解网格的两种方法。Split Factor (分割因子)和Collapse Factor(合并因子)则分别为上面介绍的alpha_s和alpha_c。
4)局部重新划分网格 在Remeshing(重划网格)标签下,设臵与局部重划模型相关的参数。可以设臵的参数包括 Maximum Cell Skewness(最大畸变率)、Maximum Cell Volume(最大网格体积)和Minimum Cell Volume(最大网格体积),其含义如前所述,主要用于确定哪些网格需要被 重新划分。在缺省设臵中,如果重新划分的网格优于原网格,则用新网格代替旧网格;否 则,将保持原网格划分不变。如果无论如何都要采用新网格的话,则可以在Options(选项) 下面选择Must Improve Skewness(必须改善畸变率)选项。如果Options(选项)下面的Size Function(尺寸函数)被激活,则还可以用网格尺寸分布函数标志需要重新划分的网格。假 设在某点附近的理想网格尺寸为L ,而某个网格的尺寸为L' ,如果: L’不属于[0.8*gamma*L,1.25*gamma*L] 则网格被标志为需要重新划分的网格,并在随后的计算中被重新划分。式中的gamma用下面的 公式计算: 当alpha>0时,gamma=1+alpha*d_b^(1+2*beta) 当alpha<0时,gamma=1+alpha*d_b^[(1-beta)^-1] 式中d_b为网格到壁面边界的最小距离,alpha和beta就是需要设臵的Size Function Variation(尺寸 函数增量)和Size Function Rate(尺寸函数变化率)。 以上两式中,alpha以边界网格尺寸为基准控制内部网格的大小,alpha等于0.5表示内部网格的尺 寸至少是边界网格的1.5倍,alpha等于-0.5表示内部网格的尺寸等于边界网格的0.5倍,如果 alpha等于0则表示内部网格与边界网格大小相同;bata的取值在-1到1之间,用于控制网格从 边界到内部区域的变化速率。Beta取正值表示网格变化速率较慢取负值则表示变化速率较快 ,取0 则表示网格从边界到内部区域呈线性变化。 另外一个选项Size Function Resolution(尺寸函数分辨率)用于定义尺寸分布函数对网格大小的 分辨率,其参照值是最小网格尺寸。 点击Use defaults(使用缺省值)按钮,可以恢复系统缺省设臵。