Femap 实体模型简化和网格划分实例

Femap 实体模型简化和网格划分实例
Femap 实体模型简化和网格划分实例

Femap实体模型简化和网格划分实例

NST(Beijing)有限公司Frank Fang

关键词:Femap;NX Nastran;CAE;实体模型简化;网格划分;Meshing Toolbox;

1、出现的问题及解决思路

在使用Femap过程中,当已经划分好网格,同时也设置好边界条件后,进行求解计算时,由于各种原因(参数没有正确设置等),很有可能出现一些致命错误。如下图所示:

上图说明出现一个致命错误,点击Show Details,会弹出错误的详细信息:

然后可以通过上面的提示查找出错误的原因,从而进行参数的修正,使分析能够进行下去。

本分析的出错原因是:只对几何模型Geometry的外表面进行了面网的划分,而没有对其进行实体网格划分,如下图所示:

Model的单元类型是Plot Planar,说明网格是面网格;单元形状是Tri 3-noded,说明网格是3个节点的三角形网格。同时把Properties的属性全部隐藏,右边的网格仍然再显示,都说明当前的网格不是实体网格。

很显然,需要的是实体网格,而现在只有三角形的面网格,因此Femap会出现错误提示。

解决思路:对于实体模型,先对其进行面网格划分是没有错误的,但是当面网划分完成,同时

检查单元没有问题后,就必须对其进行实体网格划分,然后才进行分析计算。如果划分的面网或体网有问题,不能正确划分时,就需要判断是否是几何模型正确,或几何模型细节太多,需要对其简化。

本分析可以在现有面网的基础上(参数暂时不用修改),继续划分实体网格,但是需要检查面网的品质,如下图所示:

面网的长宽比、雅克比值都不是特别理想,值偏大。虽然值偏大,但是可以尝试一下看能否正常划分实体网格,执行Mesh-Geometry-Solid命令,划分体网:

但是仍然会出现问题,如下图所示:

说明划分实体网格出现错误。在Message中会有详细提示:

提示最小角度A=0.57433,说明网格质量太差。

到目前为止,实体网格还是不能正常划分,说明之前的面网格质量不太好,即需要重新划分面网格,或者由于本几何体小特征太多,先简化模型后再划分面网。

2、Femap网格划分

重新划分面网,其实就是对那些不重要的小边、小面进行合并等处理,以为得到较好网格做准备。Femap中有个非常有用的工具-Geometry Preparation,可以高效的处理此类问题,如下图所示:

可以按照其默认值1.55,或稍大值1.6来处理,如下图示:

会得到如下图所示效果:

一些面会自动合并,这就给网格划分及提高网格品质带来便利。然后再划分面网,如下图所示:

再查看面网品质,品质有提高,下图为Mesh Size=1.6是的网格品质:

再划分体网:

可以正常划分体网,在Message中有详细提示:

实体网格能够正常划分后,就可以分析求解了,结果如下图所示:

3、Femap几何模型简化工具

虽然上面可以划分出体网,并能够正常分析计算,但是其网格品质也不够理想,其C值最大为32.768,J值最大为0.906。因此如果想更好的提高网格品质,就必须先对几何模型简化,删除或抑制那些不重要的小边、小面、倒角、圆角、孔等。

Femap中对于解决此类问题有一个很强大的工具:Meshing Toolbox,其中有特征抑制、特征编辑、特征删除等命令,可以很方便地处理此类小特征,使模型得到简化。

当然,如果先在CAD软件中简化则更方便,如在NX、Pre/E、CATIA、Solid Edge、Solidworks等中简化模型,然后导入Femap中。Femap可以接收上述主流CAD的数据,如下图示:

大家有兴趣可以尝试一下,从实体模型简化(在CAD软件简化、或在Femap中简化)和导入、网格划分及修改,到最后得到比较理想的分析结果,其中网格划分的过程可参考上面的方法。

4、小结

本文主要介绍了Femap实体模型简化和网格划分中出现的一些常见问题,及解决问题的常规思路和方法,Femap软件简单易用,用许多非常实用的工具,可以方便解决非常复杂的问题。上述的思路和方法供使用Femap的同仁们参考,希望能抛砖引玉。

2012年10月

附:如大家有Femap技术问题,欢迎交流:fangfenghao@https://www.360docs.net/doc/281549931.html,。

ABAQUS常用技巧归纳(图文并茂).

ABAQUS学习总结 1.ABAQUS中常用的单位制。-(有用到密度的时候要特别注意) 单位制错误会造成分析结果错误,甚至不收敛。 2.ABAQUS中的时间 对于静力分析,时间没有实际意义(静力分析是长期累积的结果)。对于动力分析,时间是有意义的,跟作用的时间相关。 3.更改工作路径 4.对于ABAQUS/Standard分析,增大内存磁盘空间会大大缩短计算 时间;对于ABAQUS/Explicit分析,生成的临时数据大部分是存储在内存中的关键数据,不写入磁盘,加快分析速度的主要方法是提高CPU的速度。 临时文件一般存储在磁盘比较大的盘符下

提高虚拟内存

5.壳单元被赋予厚度后,如何查看是否正确。 梁单元被赋予截面属性后,如休查看是否正确。 可以在VIEW的DISPLAY OPTION里面查看。 6.参考点 对于离散刚体和解析刚体部件,参考点必须在PART模块里面定义。而对于刚体约束,显示休约束,耦合约束可以在PART ,ASSEMBLY,INTERRACTION,LOAD等定义参考点. PART模块里面只能定义一个参考点,而其它的模块里面可以定义很多个参考点。

7.刚体部件(离散刚体和解析刚体),刚体约束,显示体约束 离散刚体:可以是任意的形状,无需定义材料属性,要定义参考点,要划分网格。 解析刚体:只能是简单形状,无需定义材料属性,要定义参考点,不需要划分网格。 刚体约束的部件:要定义材料属性,要定义参考点,要划分网格。显示体约束的部件:要定义材料属性,要定义参考点,不需要要划分网格(ABAQUS/CAE会自动为其要划分网格)。 刚体与变形体比较:刚体最大的优点是计算效率高,因为它在分析作业过程中不参与所在基于单元的计算,此外,在接触分析,如果主面是刚体的话,分析更容易收敛。 刚体约束和显示体约束与刚体部件的比较:刚体约束和显示体约束的优点是去除约束后,就可以立即变为变形体。 刚体约束与显示体约束的比较:刚体约束的部件会参与计算,而显示约束的部件不会参与计算,只是用于显示作用。 8.一般分析步与线性摄动分析步 一般分析步:每个分析步的开始状态都是前一个分析步结束时刻的模型状态; 如果不做修改的话,前一个分析步所施加的载荷,边界条件,约束都会延续到当前的分析步中;所定义的载荷,边界条件以及得到的分析结果都是总量。

ANSYS网格划分总结大全

有限元分析中的网格划分好坏直接关系到模型计算的准确性。本文简述了网格划分应用的基本理论,并以ANSYS限元分析中的网格划分为实例对象,详细讲述了网格划分基本理论及其在工程中的实际应用,具有一定的指导意义。 1 引言 ANSYS有限元网格划分是进行数值模拟分析至关重要的一步,它直接影响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。从几何表达上讲,梁和杆是相同的,从物理和数值求解上讲则是有区别的。同理,平面应力和平面应变情况设计的单元求解方程也不相同。在有限元数值求解中,单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成,连续体单元以及壳、板、梁单元的面内均采用高斯(Gauss)积分,而壳、板、梁单元的厚度方向采用辛普生(Simpson)积分。辛普生积分点的间隔是一定的,沿厚度分成奇数积分点。由于不同单元的刚度矩阵不同,采用数值积分的求解方式不同,因此实际应用中,一定要采用合理的单元来模拟求解。 2 ANSYS网格划分的指导思想 ANSYS网格划分的指导思想是首先进行总体模型规划,包括物理模型的构造、单元类型的选择、网格密度的确定等多方面的内容。在网格划分和初步求解时,做到先简单后复杂,先粗后精,2D单元和3D单元合理搭配使用。为提高求解的效率要充分利用重复与对称等特征,由于工程结构一般具有重复对称或轴对称、镜象对称等特点,采用子结构或对称模型可以提高求解的效率和精度。利用轴对称或子结构时要注意场合,如在进行模态分析、屈曲分析整体求解时,则应采用整体模型,同时选择合理的起点并设置合理的坐标系,可以提高求解的精度和效率,例如,轴对称场合多采用柱坐标系。有限元分析的精度和效率与单元的密度和几何形状有着密切的关系,按照相应的误差准则和网格疏密程度,避免网格的畸形。在网格重划分过程中常采用曲率控制、单元尺寸与数量控制、穿透控制等控制准则。在选用单元时要注意剪力自锁、沙漏和网格扭曲、不可压缩材

ANSYS中简化模型和划分网格的方法

广州有道资料网https://www.360docs.net/doc/281549931.html, ANSYS中简化模型和划分网格的方法 本文介绍了ANSYS中简化模型和划分网格的相关方法。 使在建立仿真模型时,经验是非常有助于用户决定哪些部件应该考虑因而必须建立在模型中,哪些部件不应该考虑因而不需建立到模型中,这就是所谓的模型简化。此外,网格划分也是影响分析精度的另外一个因素。本文将集中讨论如何简化模型以获得有效的仿真模型以及网格划分需要注意的一些问题。 理想情况下,用户都希望建立尽可能详细的仿真模型,而让仿真软件自己来决定哪些是主要的物理现象。然而,由于有限的计算机资源或算法限制,用户应该简化电磁仿真的模型。 模型简化 模型简化主要取决于结果参数及结构的电尺寸。例如,如果用户希望分析安装在某电大尺寸载体上的天线的远场方向图,那么模型上距离源区超过一个波长的一些小特征和孔径(最大尺度小于/50)就可以不考虑。另一方面,如果用户希望分析从源到用带有小孔的屏蔽面屏蔽的导线之间的耦合,那么必须对小孔、靠近源的屏蔽面以及导线进行精确建模。另外一个常用的简化是用无限薄的面来模拟有限厚度的导体面。一般而言,厚度小于/100的金属面都可以近似为无限薄的金属面。有限导电性和有限厚度的影响可以在SK卡中设置。对于比较厚的导体面,如果这种影响是次要的,那么用户仍然可以采取这种近似。例如,当建立大反射面天线的馈源喇叭模型时,喇叭壁的有限厚度对于反射面天线主波束的影响就是次要的。然而,如果喇叭天线用于校准标准时,那么喇叭壁的有限厚度就不能忽略。 网格划分 一般而言,网格划分的密度设置为最短波长的十分之一。然而,在电流或电荷梯度变化剧烈的区域,如源所在区域、曲面上的缝隙和曲面的棱边等,必须划分得更密。一个实用的指导原则是网格大小应该与结构间的间隔距离(d)相比拟(%26lt;=2d)。同样地,如果需要计算近场分布,那么网格大小应该同场点到源点间距离(d)相比拟。 总之,用户建立的几何模型应该抓住主要的物理现象,而网格划分则需要权衡输出结果相对于网格大小的收敛性。 广州有道资料网https://www.360docs.net/doc/281549931.html,

基于变分网格的曲面简化高效算法

基于变分网格的曲面简化高效算法? 金勇, 吴庆标+, 刘利刚 (浙江大学数学系,浙江杭州 310027) An Efficient Method for Surface Simplification Based On Variational Shape Approximation* JIN Yong, WU Qing-biao+, LIU Li-gang (Department of Mathematics, Zhejiang University, Hangzhou 310027, China) + Corresponding author: E-mail:qbwu@https://www.360docs.net/doc/281549931.html, Abstract:Providing fast and accurate simplification method for large polygon mesh is one of the most important research focuses in computer graphics. Approximating mesh model with a few polygons can improve the rendering speed, and reduce the storage of the model. The paper presents a local greedy algorithm to minimize the energy defined by variational shape approximation. The algorithm simplifies the mesh by controlling the number of the target polygons, while attempting to get ideal effect by adaptive seed triangles selection. The algorithm has intuitive geometric meaning. The method is efficient enough to be efficiently adopted in the geometric modeling system. Key words: Polygon mesh simplification; variational shape approximation; greedy algorithm; geometric modeling 摘要: 为大型的多边形网格模型提供快速、准确的简化算法是计算机图形学中的一个重要的研究方面.以较少的多边形逼近表示网格模型,能够提高模型的绘制速度,减小模型的存储空间.本文根据变分网格逼近表示所定义的全局误差能量,提出一种局部贪心优化算法,该算法通过控制目标网格分片数来简化网格,通过种子的自适应选取以达到理想的简化效果,具有直观的几何意义.本文方法计算量少,效率较高,能够有效应用于几何造型系统中. 关键词:多边形网格简化;变分网格逼近;贪心算法;几何造型 中图法分类号: TP391文献标识码: A 1 引言 三维多边形网格模型,包括三角形网格、四边形网格等,在计算机辅助几何设计、计算机动画、虚拟现实、计算机游戏和医学影像等领域有着大量的应用.随着三维扫描技术的发展,顶点数为数万的模型已经非常常见, ?Supported by the National Natural Science Foundation of China under Grant No.10871178, 60776799 (国家自然科学基金); Technology Department of Zhejiang Province Grant No. 2008C01048-3(浙江省重大科技创新项目) 作者简介: 金勇(1985-),男,上海人,博士研究生,主要研究领域为数字几何处理和计算机辅助几何设计;吴庆标(1963-),男, 浙江台州人,博士,教授,博士生导师,主要研究领域为图形与图像处理,数值计算方法,高性能并行计算和计算机模拟; 刘利刚(1975-),男,江西吉安人,博士,副教授,博士生导师,主要研究领域为数字几何处理,计算机辅助几何设计,计算机图形学和图像处理.

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 调入分析系统及网格划分组件

本人学习abaqus五年的经验总结-让你比做例子快十倍

第二章 ABAQUS 基本使用方法 [2](pp15)快捷键:Ctrl+Alt+左键来缩放模型;Ctrl+Alt+中键来平移模型;Ctrl+Alt+右键来旋转模型。 ②(pp16)ABAQUS/CAE 不会自动保存模型数据,用户应当每隔一段时间自己保存模型以避免意外丢失。 [3](pp17)平面应力问题的截面属性类型是Solid(实心体)而不是Shell(壳)。 ABAQUS/CAE 推荐的建模方法是把整个数值模型(如材料、边界条件、载荷等)都直接定义在几何模型上。载荷类型Pressure 的含义是单位面积上的力,正值表示压力,负值表示拉力。 [4](pp22)对于应力集中问题,使用二次单元可以提高应力结果的精度。 [5](pp23)Dismiss 和Cancel 按钮的作用都是关闭当前对话框,其区别在于:前者出现在包含只读数 据的对话框中;后者出现在允许作出修改的对话框中,点击Cancel 按钮可关闭对话框,而不保存 所修改的内容。 [6](pp26)每个模型中只能有一个装配件,它是由一个或多个实体组成的,所谓的“实体”(instance) 是部件(part)在装配件中的一种映射,一个部件可以对应多个实体。材料和截面属性定义在部件上,相互作用(interaction)、边界条件、载荷等定义在实体上,网格可以定义在部件上或实体上,对求解过程和输出结果的控制参数定义在整个模型上。 [7](pp26) ABAQUS/CAE 中的部件有两种:几何部件(native part)和网格部件(orphan mesh part)。 创建几何部件有两种方法:(1)使用Part 功能模块中的拉伸、旋转、扫掠、倒角和放样等特征来直 接创建几何部件。(2)导入已有的CAD 模型文件,方法是:点击主菜单File→Import→Part。网格部件不包含特征,只包含节点、单元、面、集合的信息。创建网格部件有三种方法:(1)导入ODB 文件中的网格。(2)导入INP 文件中的网格。(3)把几何部件转化为网格部件,方法是:进入Mesh 功能模块,点击主菜单Mesh→Create Mesh Part。 [8](pp31)初始分析步只有一个,名称是initial,它不能被编辑、重命名、替换、复制或删除。在初始分析步之后,需要创建一个或多个后续分析步,主要有两大类:(1)通用分析步(general analysis step)可以用于线性或非线性分析。常用的通用分析步包含以下类型: —Static, General: ABAQUS/Standard 静力分析 —Dynamics, Implicit: ABAQUS/Standard 隐式动力分析 —Dynamics, Explicit: ABAQUS/ Explicit 显式动态分析 (2)线性摄动分析步(linear perturbation step)只能用来分析线性问题。在ABAQUS/Explicit 中 不能使用线性摄动分析步。在ABAQUS/Standard 中以下分析类型总是采用线性摄动分析步。 —Buckle: 线性特征值屈曲。 —Frequency: 频率提取分析。 —Modal dynamics: 瞬时模态动态分析。 —Random response: 随机响应分析。 —Response spectrum: 反应谱分析。 —Steady-state dynamics: 稳态动态分析。 [9](pp33)在静态分析中,如果模型中不含阻尼或与速率相关的材料性质,“时间”就没有实际的物 理意义。为方便起见,一般都把分析步时间设为默认的 1。每创建一个分析步,ABAQUS/CAE 就会自动生成一个该分析步的输出要求。 [10] (pp34)自适应网格主要用于ABAQUS/Explicit 以及ABAQUS/Standard 中的表面磨损过程 模拟。在一般的ABAQUS/Standard 分析中,尽管也可设定自适应网格,但不会起到明显的作用。 Step 功能模块中,主菜单Other→Adaptive Mesh Domain 和Other→Adaptive Mesh Controls 分别 设置划分区域和参数。 [11](pp37)使用主菜单Field 可以定义场变量(包括初始速度场和温度场变量)。有些场变量与分析步有关,也有些仅仅作用于分析的开始阶段。使用主菜单Load Case 可以定义载荷状况。载荷状况由一系列的载荷和边界条件组成,用于静力摄动分析和稳态动力分析。

网格划分实例详细步骤

一个网格划分实例的详解 该题目条件如下图所示: 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.

一种新的边折叠网格模型简化算法

Computer Engineering and Applications 计算机工程与应用 2013,49(1)1引言降水粒子的智能识别,一直是大气探测领域的热门研究方向之一,降水粒子的有效识别和信息的提取,对于降水天气现象的自动化识别和数据分析至关重要。获取降水粒子的尺度谱和速度谱,并计算和统计降水强度、降水粒子总数、累积降水量以及雷达反射率因子等物理量,将为分析降水过程,了解地面雨滴谱微物理量变化提供基础资料[1]。对雨滴粒径的传统观测方法有面粉球法和吸水纸斑迹法等。面粉球法将雨滴收集在盛有面粉的容器中,当下降雨滴与面粉接触后,每个雨滴就产生一个小小的湿面球。每次测量都必须让面粉球在取样器中自然风干一天,然后放在烘箱内105℃条件下烘48h ,再用高精度电子天 平逐级称量,过程繁琐实时性差,且不适用于测量直径太小的雨滴。斑迹法是历史悠久、应用最广泛的一种雨滴粒径测量方法。该方法基于水滴在同一材料上形成的斑迹大小与水滴的粒径大小成正比的假定,通过实验测量预先设定好水滴粒径与斑迹粒径之间的比例关系,然后通过测量雨滴在相同材料上形成的斑迹大小推知相应的雨滴粒径;该方法雨滴粒径与斑迹粒径之间比例关系的确定以及滤纸和涂料的选取,对测量精度的影响都很大,并且无法避免雨滴溅射对观测结果的影响。虽然图像自动识别方法在利用斑迹法进行雨滴粒径观测中的应用能够减少资Snake 模型在雨滴边缘检测中的应用 卞真稳1,2,吕伟涛2,杨俊2,马颖2,马明1 BIAN Zhenwen 1,2,LV Weitao 2,YANG Jun 2,MA Ying 2,MA Ming 1 1.中国科学技术大学地球和空间科学学院,合肥230026 2.中国气象科学研究院大气探测研究所,北京100081 1.School of Earth and Space Sciences,University of Science and Technology of China,Hefei 230026,China 2.Institute of Atmospheric Sounding,Chinese Academy of Meteorological Sciences,Beijing 100081,China BIAN Zhenwen,LV Weitao,YANG Jun,et al.Application of snake model on raindrop edge https://www.360docs.net/doc/281549931.html,puter Engineering and Applications,2013,49(1):186-190. Abstract :Digital photography has very good application prospect in automatic observation of precipitation phenomenon,in which how to accurately detect the edge of precipitation particle from the digital image is a key technology.Snake model has the ability to merge prior knowledge and image processing algorithms,which can be used to accurately identify the outline of target.Considering the characteristics of raindrop image,an improved method to automatically select the initial contour of Snake model based on target shape heart automatic calibration method is presented,and the iteration processing of the greedy algorithm is used in the detection process.The raindrop detection method based on Snake model can be used to accurately detect the raindrop edge profile and has good https://www.360docs.net/doc/281549931.html,pared with the traditional edge detection operators,the proposed method has better edge detection effect. Key words :raindrop image;edge detection;Snake model;greedy algorithm 摘要:数字摄像技术在降水粒子的自动观测中具有非常好的应用前景,如何在数字图像中准确地进行降水粒子的边缘检测是其中的一项关键技术。Snake 模型具有很好的融合图像上层知识和底层特征的能力,能够实现目标轮廓的准确定位。结合雨滴图像自身的特点,提出了目标形心的自动标定方法,在此基础上改进了Snake 模型初始轮廓点的选取方法,并通过贪婪算法进行迭代处理,实现了基于Snake 模型的雨滴边缘检测算法。算法能够准确地对数字图像中的雨滴边缘轮廓进行检测,且具有较好的稳定性。与传统的边缘检测算子相比,该方法对雨滴图像获得了更好的边缘检测效果。关键词:雨滴图像;边缘检测;Snake 模型;贪婪算法 文献标志码:A 中图分类号:TP391doi :10.3778/j.issn.1002-8331.1112-0245 基金项目:国家科技部科研院所技术开发研究专项(No.NCSTE-2006-JKZX-303)。 作者简介:卞真稳(1982—),男,硕士研究生,主要研究领域:大气探测;吕伟涛,男,博士,研究员;杨俊,男,博士,副研究员;马颖,女,高级 工程师;马明,男,博士,副教授。E-mail :zhwbian@https://www.360docs.net/doc/281549931.html, 收稿日期:2011-12-16修回日期:2012-02-10文章编号:1002-8331(2013)01-0186-05 CNKI 出版日期:2012-05-21https://www.360docs.net/doc/281549931.html,/kcms/detail/11.2127.TP.20120521.1142.058.html 186

ANSYS网格划分的一些例子

虽然做出来了.但是我还是有一个问题想请教大家: vsweep和mapp分网后形成网格各有什么规律?如何结合两种方法划分出整齐规则的网格呢. 比如:为什么图中的(1)部分用MAPP划分,(2)部分用SWEEP划分呢就可以出现上图中的那种整齐规则的网格?反过来(1)部分用SWEEP,(2)部分用MAPP划分就不会出现整齐规则的网格呢? 部分(1)和部分(2)不可看成一个整体划分吗? 我试了一个,如果把两个部分看成整体,可以分网但是不会出现那种整齐的网格. 只有掌握了生成网格规律才容易得到合理,整齐,规则的网格,总不能分网时把各种方法都试一遍吧. 恳请各位谈点自己的在分网方面的经验.谢谢 1的三个边如果都设了分段数则sweep和map是一样的 et,1,42 et,2,45 cyl4,,,20 lsel,all lesize,all,,,10

esize,,10 vext,1,,,,,20 aclear,all amesh,1 不过好象中间不大好的!!!还望高手指点! 命令流; et,1,42 et,2,45 blc4,,,10,5 lesize,1,,,5 lesize,2,,,10 mshape,0,2d mshkey,1 amesh,1 esize,,5 vrotat,1,,,,,,1,4 aclear,all

用map也可以,,取其四分之一,单元大小可控制!做了一个!

/PREP7 CYL4, , ,5 RECTNG,-1,1,-1,1, FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,-2 AOVLAP,P51X wpro,,90.000000, wpro,,,45.000000 ASBW, 3 wpro,,,-45.000000 wpro,,,-45.000000 FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,4 ASBW,P51X wpro,,,45.000000 wpro,,90.000000, ET,1,PLANE42 ESIZE,1,0, AMAP,6,12,9,7,8 WPSTYLE,,,,,,,,0 AMAP,7,12,10,5,8 AMAP,3,10,11,6,5 AMAP,5,6,7,9,11 AMAP,2,8,5,6,7 TYPE, 1 EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,_Z2 REAL,_Z4 ESYS,0 ET,2,SOLID45 TYPE, 2 EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,_Z2 REAL,_Z4 ESYS,0 VOFFST,2,5, , VOFFST,6,5, , VOFFST,7,5, ,

网格划分

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

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

ABAQUS实体单元类型总结

在ABAQUS中,基于应力/位移的实体单元类型最为丰富: (1)在ABAQUS/Sandard中,实体单元包括二维和三维的线性单元和二次单元,均可以采用完全积分或缩减积分,另外还有修正的二次Tri单元(三角形单元)和Tet单元(四面体单元),以及非协调模式单元和杂交单元。 (2)ABAQUS/Explicit中,实体单元包括二维和三维的线性缩减积分单元,以及修正的二次二次Tri单元(三角形单元)和Tet单元(四面体单元),没有二次完全积分实体单元。 ------------------------------------------------------------------------------------------------------------ 按照节点位移插值的阶数,ABAQUS里的实体单元可以分为以下三类: 线性单元(即一阶单元):仅在单元的角点处布置节点,在各个方向都采用线性插值。 二次单元(即二阶单元):在每条边上有中间节点,采用二次插值。 修正的二次单元(只有Tri 或Tet 才有此类型):在每条边上有中间节点,并采用修正的二次插值。 ******************************************************************************* *************** 1、线性完全积分单元:当单元具有规则形状时,所用的高斯积分点的数目足以对单元刚度矩阵中的多项式进行精确积分。 缺点:承受弯曲载荷时,会出现剪切自锁,造成单元过于刚硬,即使划分很细的网格,计算精度仍然很差。 2、二次完全积分单元: 优点: (1)应力计算结果很精确,适合模拟应力集中问题; (2)一般情况下,没有剪切自锁问题(shear locking)。 但使用这种单元时要注意: (1)不能用于接触分析; (2)对于弹塑性分析,如果材料不可压缩(例如金属材料),则容易产生体积自锁(volumetric locking);

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等格式来数据

abaqus网格划分

如何使用3D实体单元? 1 如果不需要模拟非常大的应变或进行一个复杂的、改变接触条件的问题,则应采用二次减缩 积分单元(CAX8R,CRE8R,CPS8R.C3D20R等)。 2 如果存在应力集中,则应在局部采用二次完全积分单元(CAX8,CPE8,CPS8,C3D20等)。它 们可在较低费用下对应力梯度提供最好的解决。尽量不要使用线性减缩积分单元。用细化的二次减缩积分单元与二次完全积分单元求解结果相差不大,且前者时间短。 3 对含有非常大的网格扭曲模拟(大应变分析),采用细网格划分的线性减缩积分单元 (CAX4R,CPE4R.CPS4R,C3D8R等)。 4 对接触问题采用线性减缩积分单元或非协调单元(CAX4I,CPE4I,CPS4II,C3D8I等)的细网 格划分。 5 对以弯曲为主的问题,如能保证所关心部位单元扭曲较小,使用非协调单元(如C3D8I), 求解很精确。 6 对于弹塑性分析,不可压缩材料(如金属),不能使用二次完全积分单元,否则易体积自锁, 应使用修正的二次三角形或四面体单元、非协调单元,以及线性减缩积分单元。若使用二次减缩积分单元,当应变超过20%-40%要划分足够密的网格。 7 除平面应力问题之外,如材料完全不可压缩(如橡胶),应使用杂交单元; 在某些情况下,近似不可压缩材料也应使用杂交单元。 8 当几何形状复杂时,万不得已采用楔形和四面体单元。这些单元的线性形式,如C3D6和 C3D4,是较差的单元(若需要时,划分较细的网格以使结果达到合理的精度),这些单元也应远离需要精确求解的区域。 9 如使用了自由网格划分技术,四面体单元应选二次的,其结果对小位移问题应该是合理的, 但花时间多。在ABAQUS/Standard中选C3D10,ABAQUS/Explicit中选修正的(C3D10M)。 如有大的塑性变形,或模型中存在接触,且使用默认的“硬”接触关系,也应选C3D10M。 10 ABAQUS/Explicit模拟冲击或爆炸,应选线性单元。

各种网格划分方法

各种网格划分方法 1.输入实体模型尝试用映射、自由网格划分,并综合利用多种网格划分控制方法 本题提供IGES 文件 1. 以轴承座为例,尝试对其进行映射,自由网格划分,并练习一般后处理的多种技术,包 括等值图、云图等图片的获取方法,动画等。 2. 一个瞬态分析的例子 练习目的:熟悉瞬态分析过程 瞬态(FULL)完全法分析板-梁结构实例 如图所示板-梁结构,板件上表面施加随时间变化的均布压力,计算在下列已知条件下结构的瞬态响应情况。 全部采用A3钢材料,特性: 杨氏模量=2e112/m N 泊松比=0.3 密度=7.8e33 /m Kg 板壳: 厚度=0.02m 四条腿(梁)的几何特性: 截面面积=2e-42m 惯性矩=2e-84m 宽度=0.01m 高度=0.02m 压力载荷与时间的 关系曲线见下图所示。 图 质量梁-板结构及载荷示意图 0 1 2 4 6 时间(s ) 图 板上压力-时间关系 分析过程 第1步:设置分析标题 1. 选取菜单途径Utility Menu>File>Change Title 。 2. 输入“ The Transient Analysis of the structure ”,然后单击OK 。 第2步:定义单元类型 单元类型1为SHELL63,单元类型2为BEAM4 第3步:定义单元实常数 实常数1为壳单元的实常数1,输入厚度为0.02(只需输入第一个值,即等厚度壳)

实常数2为梁单元的实常数,输入AREA 为2e-4惯性矩IZZ=2e-8,IYY =2e-8宽度TKZ=0.01,高度TKY=0.02。 第5步:杨氏模量EX=2e112/m N 泊松比NUXY=0.3 密度DENS=7.8e33 /m Kg 第6步:建立有限元分析模型 1. 创建矩形,x1=0,x2=2,y1=0,y2=1 2. 将所有关键点沿Z 方向拷贝,输入DZ =-1 3. 连线。将关键点1,5;2,6;3,7;4,8分别连成直线。 4. 设置线的分割尺寸为0.1,首先给面划分网格;然后设置单元类型为2,实常数为2, 对线5到8划分网格。 第7步:瞬态动力分析 1. 选取菜单途径Main Menu>Solution>-Analysis Type-New Analysis ,弹出New Analysis 对话框。 2. 选择Transient ,然后单击OK ,在接下来的界面仍然单击OK 。 3. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Damping , 弹出Damping Specifications 窗口。 4. 在Mass matrix multiplier 处输入5。单击OK 。 5. 选取菜单途径Main Menu > Solution > -Loads-Apply > -Structural- Displacement>On Nodes 。弹出拾取(Pick )窗口,在有限元模型上点取节点232、242、252和262,单击OK ,弹出Apply U,ROT on Nodes 对话框。 6. 在DOFS to be constrained 滚动框中,选种“All DOF ”(单击一次使其高亮度显示, 确保其它选项未被高亮度显示)。单击OK 。 7. 选取菜单途径Utility Menu>Select>Everything 。 8. 选取菜单途径Main Menu>Solution>-Load Step Opts-Output Ctrls>DB/Results File ,弹出Controls for Database and Results File Writing 窗口。 9. 在Item to be controlled 滚动窗中选择All items ,下面的File write frequency 中选择Every substep 。单击OK 。 10. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Time – Time Step ,弹出Time – Time Step Options 窗口。 11. 在Time at end of load step 处输入1;在Time step size 处输入0.2;在Stepped or ramped b.c 处单击ramped ;单击Automatic time stepping 为on ;在Minimum time step size 处输入0.05;在Maximum time step size 处输入0.5。单击OK 。 12. 选取菜单途径Main Menu>Solution>-Loads-Apply>-Structure-Pressure>On Areas 。弹出Apply PRES on Areas 拾取窗口。 13. 单击Pick All ,弹出Apply PRES on Areas 对话框。 14. 在pressure value 处输入10000。单击OK 15. 选取菜单途径Main menu>Solution>Write LS File ,弹出Write Load Step File 对 话框。 16. 在Load step file number n 处输入1,单击OK 。 17. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Time – Time Step ,弹出Time – Time Step Options 窗口。

相关文档
最新文档