网格划分的几种基本处理方法 贴体坐标法： 贴体坐标是利用曲线坐标，并使其坐标线与燃烧室外形或复杂计算区域边界重合，这样所有边界点能够用网格点来表示，不需要任何插值。一旦贴体坐标生成通过变换，偏微分方程求解可以不在任意形状的物理平面上，而在矩形或矩形的组合（空间问题求解域为长方体或它们的组合）转换平面上进行。这样计算与燃烧室外形无关，也与在物理平面上网格间隔无关。 而是把边界条件复杂的问题转换成一个边界条件简单的问题；这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差，而且还可节省计算时间和内存，使流场计算较准确，同时方便求解，较好地解决了复杂形状流动区域的计算，在工程上比较广泛应用。 区域法： 虽然贴体坐标系可以使坐标线与燃烧室外形相重合，从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域，很难用一种网格来模拟，生成单域贴体网格，即使生成了也不能保证网格质量，影响流场数值求解的效果。因此，目前常采用区域法或分区网格，其基本思想是，根据外形特点把复杂的物理域或复杂拓扑结构的网格，分成若干个区域，分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格，或结构块网格。对区域进行分区时，若相邻两个子域分离边界是协调对接，称为对接网格；若相邻两子域有相互重叠部分，则此分区网格称为重叠网格。根据实际数值模拟计算的需要，把整个区域（燃烧室）分成几个不同的子区域，并分别生成网格。这样不仅可提高计算精度，而且还可节省计算机内存，提高收敛精度。但是计算时，必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法，其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程．设压力校正值在最初迭代时为零，为了保证流量连续各个区域应同时求解，然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区，两个区域都对重叠区进行计算，重叠区一边区域内的值，要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系，实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同，要求重叠区二边流量相等。区域法能合理解决网格生成问题，已被大量用来计算复杂形状区域流动。 区域分解法： 对于复杂几何形状的实际燃烧装置，为了保证数值求解流场质量，目前常采用区域分解法。该法基本要点是：根据燃烧室形状特点和流场计算需要，把计算区域分成一个主区域和若干个子区域，对各个区域（块）分别建立网格，并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立，各个子区域大小也尽可能相同，使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同，通常在变量变化梯度大的区域，可以布置较细网格，并采用高阶紊流模型和描述复杂反应的紊流燃烧模型，以便更合理模拟实际流场。对于变量变化不太大区域，可采用较疏的网格和较简单的数学模型，这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑，相邻子区域连接重叠网格或对接网格来实现，在各子区域交界处通过插值法提供各子域求解变量的信息传递，满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和耦合，从而取得全流场解。 非结构网格法： 上述各方法所生成的网格均属于结构化网格，其共同特点是网格中各节点排列有序，每个节点与邻点之间关系是固定的，在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验，计算技术也较成熟，目前被广泛用来构造复杂外形区域
向量法解立体几何 用传统的方法解立体几何需要烦琐的分析、复杂的计算。而用向量法解题思路清晰、过程简洁。对立体几何的常见问题都可以起到化繁为简，化难为易的效果。 一. 证明两直线平行 已知两直线a 和b , b D C a B A ∈∈,,,,则?b a //存在唯一的实数λ使CD AB λ= 二. 证明直线和平面平行 1.已知直线αα∈∈?E D C a B A a ,,,,,且三点不共线，则a ∥?α存在有序实数 对μλ,使CE CD AB μλ+= 2.已知直线,,,a B A a ∈?α和平面 α的法向量n ,则a ∥n AB ⊥?α 三.证明两个平面平行 已知两个不重合平面βα,,法向量分别为n m ,,则α∥n m //?β 四．证明两直线垂直 已知直线b a ,。b D C a B A ∈∈,,,，则0=??⊥CD AB b a 五.证明直线和平面垂直 已知直线α和平面a ，且A 、B a ∈,面α的法向量为m ,则m AB a //?⊥α 六.证明两个平面垂直 已知两个平面βα,,两个平面的法向量分别为n m ,,则n m ⊥?⊥βα 七．求两异面直线所成的角 已知两异面直线b a ,，b D C a B A ∈∈,,,，则异面直线所成的角θ 为：CD AB ?=θcos 八．求直线和平面所成的角 A B
已知A,B 为直线a 上任意两点，n 为平面α的法向量，则a 和平面α所成的角θ为: 1． 当??? ? ??2, 0π 时?-=2πθ 2． 当??? ??∈?ππ,2 时2πθ-?= 九．求二面角 1．已知二面角βα--l ，且l CD l AB D C B A ⊥⊥∈∈,,,,且βα，则二面角的平面角θ 的大小为：=θ 2.已知二面角,βα--l n m ,分别为面βα,的法向量，则二面角的平面角θ的 大小与两个法向量所成的角相等或互补。即-=πθ 注:如何判断二面角的平面角和法向量所成的角的关系。 (1)通过观察二面角锐角还是钝角,再由法向量的成的角求之。 (2)通过观察法向量的方向,判断法向量所成的角与二面角的平面角相等还是互补。 十.求两条异面直线的距离 已知两条异面直线b a ,,m 是与两直线都垂直的向量,b B a A ∈∈,则两条 异面直线的距离d = 十一.求点到面的距离 已知平面α和点A,B 且αα∈?B A ,,m 为平面α的法向量,则点A 到平面 α 的距离d =
【分享】复杂几何模型的系列网格划分技术 众所周知，对于有限元分析来说，网格划分是其中最关键的一个步骤，网格划分的好坏直接影响到解算的精度和速度。在ANSYS中，大家知道，网格划分有三个步骤：定义单元属性（包括实常数）、在几何模型上定义网格属性、划分网格。在这里，我们仅对网格划分这个步骤所涉及到的一些问题，尤其是与复杂模型相关的一些问题作简要阐述。 一、自由网格划分 自由网格划分是自动化程度最高的网格划分技术之一，它在面上（平面、曲面）可以自动生成三角形或四边形网格，在体上自动生成四面体网格。通常情况下，可利用ANSYS的智能尺寸控制技术（SMARTSIZE命令）来自动控制网格的大小和疏密分布，也可进行人工设置网格的大小（AESIZE、LESIZE、KESIZE、ESIZE等系列命令）并控制疏密分布以及选择分网算法等（MOPT命令）。对于复杂几何模型而言，这种分网方法省时省力，但缺点是单元数量通常会很大，计算效率降低。同时，由于这种方法对于三维复杂模型只能生成四面体单元，为了获得较好的计算精度，建议采用二次四面体单元（92号单元）。如果选用的是六面体单元，则此方法自动将六面体单元退化为阶次一致的四面体单元，因此，最好不要选用线性的六面体单元（没有中间节点，比如45号单元），因为该单元退化后为线性的四面体单元，具有过刚的刚度，计算精度较差；如果选用二
次的六面体单元（比如95号单元），由于其是退化形式，节点数与其六面体原型单元一致，只是有多个节点在同一位置而已，因此，可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元，减少每个单元的节点数量，提高求解效率。在有些情况下，必须要用六面体单元的退化形式来进行自由网格划分，比如，在进行混合网格划分（后面详述）时，只有用六面体单元才能形成金字塔过渡单元。对于计算流体力学和考虑集肤效应的电磁场分析而言，自由网格划分中的层网格功能（由LESIZE命令的LAYER1和LAYER2域控制）是非常有用的。 二、映射网格划分 映射网格划分是对规整模型的一种规整网格划分方法，其原始概念是：对于面，只能是四边形面，网格划分数需在对边上保持一致，形成的单元全部为四边形；对于体，只能是六面体，对应线和面的网格划分数保持一致；形成的单元全部为六面体。在ANSYS中，这些条件有了很大的放宽，包括： 1 面可以是三角形、四边形、或其它任意多边形。对于四边以上的多边形，必须用LCCAT命令将某些边联成一条边，以使得对于网格划分而言，仍然是三角形或四边形；或者用AMAP命令定义3到4个顶点（程序自动将两个顶点之间的所有线段联成一条）来进行映射划分。 2 面上对边的网格划分数可以不同，但有一些限制条件。
常见几何体的面积、体积求法与应用 要计算某材料的密度、重量，研究某物体性能及其物质结构等，特别对于机械专业的学生，必须要求工件的面积、体积等，若按课本上公式来计算，而课本上公式不统一，不好记住，并且很繁杂，应用时要找公式，对号入座很麻烦。笔者在教学与实践中总结出一种计算常见几何体的面积、体积方法。其公式统一，容易记住，且计算简单。对技校学生来说，排除大部分繁琐的概念、定理，以及公式的推导应用等。 由统计学中的用加权平均数对估计未来很准确。比如，估计某商品下个月销售量，若去年平均销售量为y ，设本月权为4，上月权数为1，下月权数为1，各月权数分别乘销售量相加后除以6等于y 。这样能准确地确定下个月销售量。能不能以这种思想方法用到求几何体的面积、体积呢？通过推导与实践，对于常见的几何体确实可用这种方法来求得其面积、体积。下面分别说明求常见几何体的面积、体积统一公式的正确性与可用性。 常见几何体的面积、体积统一公式： ) 4(6 )4(621002100S S S h V C C C h A ++= ++= （其中A 为几何体侧面积，C 0为上底面周长，C 1为中间横截面周长，C 2 为下底面周长，V 为几何体体积，S 0为上底面面积，S 1为中间横截面面积，S 2为下底面面积，h 为高，h 0为斜高或母线长。注：中间横截面为上、下底等距离的截面。） 一、棱柱、棱锥、棱台、圆柱、圆锥、圆台的面积 、体积用统一公式的正确性 1、棱柱： ⑴据棱柱上底周长、下底周长、中间横截面周长相等，即2 1 C C C ==， 可得： 2020210066 )4(6 C h C h C C C h =?= ++，这与课本中的棱柱侧面积公式等同。 以下每个几何体都能推得与课本中相应公式等同，说明这统一公式的正确性。 ⑵据棱柱上底面、下底面、中间横截面相等，可知：2 1 S S S ==，即： h S S S S h S S S h V 2222210)4(6 )4(6 =++= ++= 。 2、棱锥 ⑴设底边长为a 2，边数为n ，斜高为h 0，侧面三角形中位线为a 1，则
第 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 调入分析系统及网格划分组件
立体几何与平面几何计算公式 初中数学几何中，不论是平面几何还是立体几何，他们的计算公式是我们进行数学试题计算的基础，因此，希望中考考生积极的做好几何计算公式的复习。下面是初中数学几何计算公式，一起了解一下： 1 、正方形 C：周长S：面积：a：边长 周长=边长×4 C=4a 正方形面积=边长×边长S= a a 2 、长方形C：周长S：面积a：边长 周长=(长+宽)×2 C = 2(a+b) 长方形面积=长×宽S = a b 3 、三角形s：面积a：底h：高 三角形面积=底×高÷2 s = ah÷2 4 、平行四边形s：面积a：底h：高 平行四边形面积=底×高s = ah 5、梯形s面积a上底b下底h高 梯形面积=(上底+下底)×高÷2 s = (a+b) h÷2 6 、圆形r:半径d：直径c：周长s：面积 半径=直径÷2 r = d/2 半径=周长÷圆周率÷2 r = c/2π 直径=半径×2 d = 2r 直径=周长÷圆周率d = c/π
周长=圆周率×直径 c = πd 周长=圆周率×半径×2 c = 2πr 圆面积=圆周率×半径×半径s = πr r 圆环面积=圆周率×(大圆半径×大圆半径-小圆半径×小圆半径) s=π(R R-r r) 7 、长方体V:体积s:面积a:长b: 宽h:高 体积=长×宽×高V = abh 8、正方体V:体积a:棱长 总棱长=棱长×12 C = 12a 表面积=棱长×棱长×6 S表= a a6 体积=棱长×棱长×棱长V = a a a 9、圆柱体V:体积s:底面积h:高 圆柱体侧面积=底面周长×高s= c h 圆柱体体积=底面积×高V= sh 圆柱体体积=圆周率×半径×半径×高V =πr r h 圆柱体体积=1/2×侧面积×半径V =1/2s侧r 10、圆锥体V:体积s:底面积h:高 圆锥体体积=1/3×底面积×高V = 1/3sh 圆锥体体积=1/3×圆周率×半径×半径×高V = 1/3×πr r h
第3章网格划分技术及技巧 创建几何模型后，必须生成有限元模型才能分析计算，生成有限元模型的方法就是对几何模型进行网格划分，网格划分主要过程包括三个步骤： ⑴定义单元属性 单元属性包括单元类型、实常数、材料特性、单元坐标系和截面号等。 ⑵定义网格控制选项 ★对几何图素边界划分网格的大小和数目进行设置； ★没有固定的网格密度可供参考； ★可通过评估结果来评价网格的密度是否合理。 ⑶生成网格 ★执行网格划分，生成有限元模型； ★可清除已经生成的网格并重新划分； ★局部进行细化。 3.1 定义单元属性 3.1.1 单元类型 1. 定义单元类型 命令：ET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR ITYPE---用户定义的单元类型的参考号。 Ename---ANSYS单元库中给定的单元名或编号，它由一个类别前缀和惟一的编号组成，类别前缀可以省略，而仅使用单元编号。 KOP1~KOP6---单元描述选项，此值在单元库中有明确的定义，可参考单元手册。也可通过命令KEYOPT进行设置。 INOPR---如果此值为1则不输出该类单元的所有结果。 例如： et,1,link8 !定义LINK8单元，其参考号为1；也可用ET,1,8定义 et,3,beam4 !定义BEAM4单元，其参考号为3；也可用ET,3,4定义 2. 单元类型的KEYOPT 命令：KEYOPT,ITYPE,KNUM,V ALUE ITYPE---由ET命令定义的单元类型参考号。 KNUM---要定义的KEYOPT顺序号。 V ALUE---KEYOPT值。 该命令可在定义单元类型后，分别设置各类单元的KEYOPT参数。例如： et,1,beam4 !定义BEAM4单元的参考号为1 et,3,beam189 !定义BEAM189单元的参考号为3 keyopt,1,2,1 !BEAM4单元考虑应力刚度时关闭一致切线刚度矩阵 keyopt,3,1,1 !考虑BEAM189的第7个自由度，即翘曲自由度
一个网格划分实例的详解 该题目条件如下图所示： 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.
空间角和距离的计算(1) 一 线线角 1．直三棱柱A 1B 1C 1-ABC ，∠BCA=900，点D 1，F 1分别是A 1B 1和A 1C 1的中点，若BC=CA=CC 1，求BD 1与AF 1所成角的余弦值． 2．在四棱锥P-ABCD 中，底面ABCD 是直角梯形，∠BAD=900，AD ∥BC ，AB=BC=a ，AD=2a ，且PA ⊥面ABCD ，PD 与底面成300角． （1）若AE ⊥PD ，E 为垂足，求证：BE ⊥PD ； （2）若AE ⊥PD ，求异面直线AE 与CD 所成角的大小． 二．线面角 1．正方体ABCD-A 1B 1C 1D 1中，E ，F 分别为BB 1、CD 的中点，且正方体的棱长为2． （1）求直线D 1F 和AB 和所成的角； （2）求D 1F 与平面AED 所成的角． F 1D 1B 1 C 1A 1 B A C A B C D P E C D E F D 1 C 1 B 1 A 1 A B
2．在三棱柱A 1B 1C 1-ABC 中，四边形AA 1B 1B 是菱形，四边形BCC 1B 1是矩形，C 1B 1⊥AB ，AB=4，C 1B 1=3，∠ABB 1=600，求AC 1与平面BCC 1B 1所成角的大小． 三．二面角 1．已知A 1B 1C 1-ABC 是正三棱柱，D 是AC 中点． （1）证明AB 1∥平面DBC 1； （2）设AB 1⊥BC 1，求以BC 1为棱，DBC 1与CBC 1为面的二面角的大小． 2．ABCD 是直角梯形，∠ABC=900，SA ⊥面ABCD ，SA=AB=BC=1，AD=0.5． （1）求面SCD 与面SBA 所成的二面角的大小； （2）求SC 与面ABCD 所成的角． 3．已知A 1B 1C 1-ABC 是三棱柱，底面是正三角形，∠A 1AC=600，∠A 1AB=450，求二面角B —AA 1—C 的大小． B 1 C 1 A 1 B A C D B 1 C 1 A 1B A C B A D C S B 1 C 1 B C A 1
3.4 网格划分实例 复杂面模型的网格划分 1. 孔板 钢结构螺栓连接中的节点板，其板上都设有一定数量的螺栓孔，这些栓孔可能对称布置也可能不对称布置。要得到四边形映射网格必须满足其要求的条件，可对板进行适当的切分或连接。本例采用切分命令将面切成多个小面，有些可满足4 边的条件，包含曲线的面则不满足4 边的条件，可分别采用AMESH 和AMAP 命令（如用LCCAT 需要不断连接、划分、删除连接线等操作）进行映射网格划分。 示例： ! EX3.18 孔板网格划分 finish $ /clear $ /prep7 a0=300 $ b0=800 $ r0=15 ! 定义参数 blc4,,,a0,b0 $ cyl4,a0/4,b0/8,r0 ! 创建矩形面和一个圆面 agen,2,2,,,a0/2 $ agen,2,2,3,1,,b0/8 ! 复制生成其它圆面 agen,2,2,5,1,,b0*5/8 $ asel,s,,,2,9,1 ! 选择圆面 cm,a2cm,area $ allsel ! 将所选择圆面定义为组件a2cm
asba,1,a2cm ! 用矩形面减圆面，形成孔板 wprota,,-90 ! 将孔板竖向切分 *do,i,1,5 $ wpoff,,,b0/16 $ asbw,all $ *enddo wpoff,,,b0*5/16 $ *do,i,1,5 $ wpoff,,,b0/16 $ asbw,all $ *enddo wprota,,,90 ! 将孔板横向切分 *do,i,1,3 $ wpoff,,,a0/4 $ asbw,all $ *enddo wpcsys,-1 $ numcmp,all lsel,s,radius,,r0 $ lesize,all,,,8 ! 选择圆孔边界线，定义网分数为8 lsel,inve $ lesize,all,,,4 $ lsel,all ! 其余线网分数为4 et,1,82 $ mshape,0,2d $ mshkey,1 ! 定义单元类型、单元形状及网分类型 asel,u,loc,y,b0/16,b0*5/16 ! 不选择带圆孔的面 asel,u,loc,y,b0*11/16,b0*15/16 ! 不选择带圆孔的面 lsla,s$lsel,r,tan1,x ! 选择竖向线 lesize,all,50,,,,1 ! 修改这些线的网分尺寸 amesh,all $ allsel ! 划分这些面的网格
不四 s = —+ 爲Mu = =￡￡sin B 2 2 边形 不四 平边 行形 a. b. c. d —各边长險、爲rsi s -面积右、必一对角线 [H^hY^bh + cH 2 H, 曰-面枳 € _ a￡K abc % 4」戸（尹_&）〔戸 _&）（尹_亡） P-三边和之半 s-三角形囲积 艮-三角形外接圆半径 外 切 角 形 直 角 角 形 尸=匚石一刁 ■S _ 血 P V F P-三边和之半 2 -三角形面积 r -三角形内切圆半径 以=胪亠阱弘b -直角边 c = 3十戸? _斜边 1 , "尹占-面积 c -J/ ■n?十2&曰a'b^ -各边长
隅 角 0 ]073t a s - 面积 d -短轴D - 长轴匸-短 半轴 R -长半轴 扇 形 ISO* -°01745^ 亠二喫 2 360 半径 圆心角= 0.008727r^* 弓-面积
正 六 E 体 正 十 _____ L 面 体 正 多 边 形 (六个正方形 ) 口 -边 数 a - 一边之长 R -外接圆半径 r 内切圆半径 e-巒 财之 1D 心角 顶 用 官-面 积 D -周良 tzFhj u 〔教目） F=6a 2 棱顶点 12 3 丁 = / C 数 目） 稜腆点 30 20 正 立 方 体 截 头 直 锥 （十二个五甬形）爲 柱 卩二 20.6457^ r= 7.663 la 5 F = 6a 2 C L □ -边 长 d-对角线长 = 7^" = 1732^1 。=扌心1 +比） 尸=#餉+宀） + s i 十巧 衍“2 —两端周 围的长 ￡ L-S 2 —两端的 面积 $二gk 十邑+ J 远”叼） C* P -宜截断面周长 F = ^/ + 2s h - 高 V = sh 目-底面积
网格划分的方法 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 +=+=+=。
在中国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）通常情况下首先进行大量的生成，后面的编辑是比较容易的。 某些区域比较重要的网格划分的质量要求高些，如力的作用区域，边界条件所在的区域。一些设计区域和离设计区域比较远的地方可以适当放宽要求，但是最好是一些网格性能指标要满足。
N维空间几何体质心的计算方法 摘要:本文主要是求一个图形或物体的质心坐标的问题,通过微积分方面的知识来求解,从平面推广到空间,问题也由易到难。首先提出质心或形心问题,然后给出重心的定义,再由具体的例子来求解相关问题。 关键字:质心重心坐标平面薄板二重积分三重积分 一.质心或形心问题: 这类问题的核心是静力矩的计算原理。 1.均匀线密度为M的曲线形体的静力矩与质心: 静力矩的微元关系为 , dMx yudl dMy xudl ==. 其中形如曲线L( (, y f x a x b =≤≤的形状体对x轴与y轴的静力矩分别 为( b
a y f x S = ? , ( b y a M u f x =? 设曲线AB L 的质心坐标为( ,x y,则,, y x M M x y
M M == 其 中( b a M u x d x u l == ? 为AB L 的质量,L为曲线弧长。若在式 y M x M
= 与式 x M y M = 两端同乘以2π,则可得 到22( b a y xl f x S ππ == ? ,
22( b a x yl f x S ππ == ? ,其中x S 与y S 分别表示曲线AB L 绕x轴与y轴旋转而成的旋转体的侧面积。 2.均匀密度平面薄板的静力矩与质心: 设f(x为 [],a b 上的连续非负函数,考虑形如区域 {} (,,0(
D x y a x b y f x =≤≤≤≤ 的薄板质心,设M为其密度,利用微元法,小曲边梯形MNPQ的形心坐标为1 (,(, 2 y f y x y x x ≤≤+? ,当分割无限细化时,可当小曲边梯形MNPQ的质量视为集中于点 1 (,( 2 x f x 处的一个质点,将它对x轴与y轴分别取静力矩微元可有 1 (( 2 x dM u f x f x dx
网格划分的几种基本处理方法 学习2010-01-10 17:13:52 阅读48 评论0 字号：大中小 贴体坐标法： 贴体坐标是利用曲线坐标，并使其坐标线与燃烧室外形或复杂计算区域边界重合，这样所有边界点能够用网格点来表示，不需要任何插值。一旦贴体坐标生成通过变换，偏微分方程求解可以不在任意形状的物理平面上，而在矩形或矩形的组合（空间问题求解域为长方体或它们的组合）转换平面上进行。这样计算与燃烧室外形无关，也与在物理平面上网格间隔无关。 而是把边界条件复杂的问题转换成一个边界条件简单的问题；这样不仅可避免因燃烧室外形与坐标网格线不一致带来计算误差，而且还可节省计算时间和内存，使流场计算较准确，同时方便求解，较好地解决了复杂形状流动区域的计算，在工程上比较广泛应 用。 区域法： 虽然贴体坐标系可以使坐标线与燃烧室外形相重合，从而解决复杂流动区域计算问题。但有时实际流场是一个复杂的多通道区域，很难用一种网格来模拟，生成单域贴体网格，即使生成了也不能保证网格质量，影响流场数值求解的效果。因此，目前常采用区域法或分区网格，其基本思想是，根据外形特点把复杂的物理域或复杂拓扑结构的网格，分成若干个区域，分别对每个子区域生成拓扑结构简单的网格。由这些子区域组合而成的网格，或结构块网格。对区域进行分区时，若相邻两个子域分离边界是协调对接，称为对接网格；若相邻两子域有相互重叠部分，则此分区网格称为重叠网格。根据实际数值模拟计算的需要，把整个区域（燃烧室）分成几个不同的子区域，并分别生成网格。这样不仅可提高计算精度，而且还可节省计算机内存，提高收敛精度。但是计算时，必须考虑各区域连接边界处耦合以及变量信息及时、准确地传递问题。处理各个区域连接有多种方法，其中一个办法是在求解各变量时各区域可以单独求解若干次而对压力校正方程．设压力校正值在最初迭代时为零，为了保证流量连续各个区域应同时求解，然后对各个速度和压力进行校正。或者采用在两个区域交界处有一个重叠区，两个区域都对重叠区进行计算，重叠区一边区域内的值，要供重叠区另一边区域求解时用。或通过在重叠内建立两个区域坐标对应关系，实现数据在重叠区内及时传递。如果两个区采用网格疏密分布不相同，要求重叠区二边流量相等。区域法能合理解决网格生成问题，已被大量用来计算复杂形状区域流动。 区域分解法： 对于复杂几何形状的实际燃烧装置，为了保证数值求解流场质量，目前常采用区域分解法。该法基本要点是：根据燃烧室形状特点和流场计算需要，把计算区域分成一个主区域和若干个子区域，对各个区域（块）分别建立网格，并对各个区域分别进行数值求解。区域分解原则是尽量使每个子区域边界简便以便于网格建立，各个子区域大小也尽可能相同，使计算负载平衡有利于平行计算。各区域的网格间距数学模型以及计算方法都可以不同，通常在变量变化梯度大的区域，可以布置较细网格，并采用高阶紊流模型和描述复杂反应的紊流燃烧模型，以便更合理模拟实际流场。对于变量变化不太大区域，可采用较疏的网格和较简单的数学模型，这样可节省计算时间。各子区域的解在相邻子区域边界处通过耦合条件来实现光滑，相邻子区域连接重叠网格或对接网格来实现，在各子区域交界处通过插值法提供各子域求解变量的信息传递，满足各子域流场计算要求通量和动量守恒条件以便实现在交界面处各子域流场解的匹配和 耦合，从而取得全流场解。 非结构网格法： 上述各方法所生成的网格均属于结构化网格，其共同特点是网格中各节点排列有序，每个节点与邻点之间关系是固定的，在计算区域内网格线和平面保持连续。特别是其中分区结构网格生成方法已积累了较多经验，计算技术也较成熟，目前被广泛用来构造复杂外形区域内网格。但是，若复杂外形稍有改变，则将需要重新划分区域和构造网格，耗费较多人力和时间。为此，近年来又发展了另一类网格——非结构网格。此类网格的基本特点是：任何空间区域都被以四面体为单元的网格所划分，网格节点不受结构性质限制，能较好地处理边界，每个节点的邻点个数也可不固定，因此易于控制网格单元的大小、形状及网格的位置。与结构网格相比，此类网格具有更大灵活性和对复杂外形适应性。在20世纪80年代末和90年代初，非结构网格得到了迅速发展。生成非结构网格方法主要有三角化方法和推进阵面法两种。虽然非结构网格容易适合复杂外形，但与结构网格相比还存在一些缺点：（1）需要较大内存记忆单元节点之
·网格划分包含以下3个步骤: –定义单元属性 –指定网格的控制参数 –生成网格 A. 多种单元属性 ·如前所述, 每个单元有以下与之相关的属性: –单元类型(TYPE) –实常数(REAL) –材料特性(MAT) ·许多FEA模型有多种属性. 例如，下图所示的筒仓有两种单元类型, 三种实常数, 以及两种材料. ·只要您的模型中有多种单元类型(TYPEs), 实常数(REALs) 和材料(MATs), 就必须确保给每一种单元指定了合适的属性. 有以下3种途径: –在网格划分前为实体模型指定属性 –在网格划分前对MAT, TYPE,和REAL进行“总体的”设置 –在网格划分后修改单元属性 ·如果没有为单元指定属性, ANSYS将MAT=1, TYPE=1, 和REAL=1作为模型中所有单元的缺省设置. 注意, 采用当前激活的TYPE, REAL, 和MAT 进行网格操作. 为实体模型指定属性 1.定义所有需要的单元类型,材料, 和实常数. 2.然后使用网格工具的“单元属性”菜单条(Preprocessor > MeshTool): –选择实体类型后按SET键. –拾取您想要指定属性的实体. –在后续的对话框设置适当的属性.或选择需要的实体,使用VATT, AATT, LATT, 或KATT命令. 3.当您为实体划分网格时, 它的属性将自动转换到单元上.
使用总体的属性设置 1.定义所有需要的单元类型,材料, 和实常数. 2.然后使用网格工具的“单元属性”菜单条(Preprocessor > MeshTool): –选择Global后按SET 键. –在“网格划分属性”对话框中激活需要的属性组合. 这些被视为激活的TYPE, REAL,和MAT 设置. 或使用TYPE, REAL, 和MAT命令.
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 体上对应线和面的网格划分数可以不同，但有一些限制条件。
一、自由网格划分 自由网格划分是自动化程度最高的网格划分技术之一，它在面上（平面、曲面）可以自动生成三角形或四边形网格，在体上自动生成四面体网格。通常情况下，可利用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 体上对应线和面的网格划分数可以不同，但有一些限制条件。 对于三维复杂几何模型而言，通常的做法是利用ANSYS布尔运算功能，将其切割成一系列四、五或六面体，然后对这些切割好的体进行映射网格划分。当然，这