有限元分析的基本步骤.docx

一个典型的ANSYS分析过程可分为以下6个步骤：1 定义参数2 创建几何模型3 划分网格4 加载数据5 求解6 结果分析1定义参数1.1指定工程名和分析标题启动ANSYS软件，选择Jobname命令选择Title菜单命令1.2定义单位(2)设置计算类型ANSYS Main Menu: Preferencef Material Props f Material Models f Structural f OK(3)定义分析类型ANSYS Main Menu: Preprocessor f Loads f Analysis Type f New Analysisf STATIC f OK1.3定义单元类型选择 Main Menuf Preprocessorf Element Typef Add/Edit/Delete 命令单击Options按钮，在Element behavior下拉列表中选择Plane strs w/thk选项，单击确 定1.4定义单元常数在 ANSYS 程序主界面 中选择 Main Menu f Preprocessor f Real Constants f Add/Edit/Delete 命令单击Add按钮，进行下一个Choose Element Type对话框1.5定义材料参数在 ANSYS 程序主界面，选择 Main Menuf Preprocessorf Material Propsf Material Models命令(1)选择对话框右侧 Structuralf Linearf Elasticf Isotropic 命令，并单击Isotropic选项， 接着弹出如下所示Linear Isotropic Properties for Material Number 1对话框。在EX文本框中输入弹性模量“200000”，在PRXY文本框中输入泊松比“0.3”，单击 OK2创建几何模型在 ANSYS 程序主界面，选择 Main Menuf Preprocessorf ModelingfCreatfAreasf Rectangle f By 2Corners 命令选择 Main Menuf Preprocessorf Modelingf Creatf Areasf Circlef Solid Circle 命令3网格划分(之前一定要进行材料的定义和分配)选择 Main Menuf Preprocessorf Modelingf Operatef Booleansf Subtractf Arears Circle 命 令选择 Main Menuf Preprocessorf Meshingf Meshf Areasf Free 命令，弹出实体选择对话框，单击Pick All按钮，得到如下所示网格4加载数据(1) 选择 Main Menu f Preprocessor f Loads f Define Loads f Apply f Structural f Displacementf On Lines 命令，出现如下所示对话框，选择约束ALL DOF选项，并设置Displacement value伪0，单击OK。(2)选择 Main MenuPreprocessorLoadsDefine LoadsApplyStructuralPressure On Lines 命令5求解选择 Main Menu SolutionSolveCurrent LS 命令，弹出如下所示窗口6结果分析6.1显示变形图(1)选择 Main MenuGeneral PostProcRead ResultsFirst Set 命令，读入最初结果文 件选择 Main MenuGeneral PostProcPlot ResultsDeformed Shape 命令，弹出Plot Deformed Shape对话框6.2显示Von Mises等效应力1) 选择 Main MenuGeneral PostProcPlot ResultsContour PlotNodal Solu 命令2) 选择 FavoritesNodal SolutionVon Mises stress 命令，并展开Additional Options 折叠菜单，设置Interpolation Nodes为All applicable，单击OK,出现如下所示等效应 力图。单元类型的选择：单元类型的选择，跟你要解决的问题本身密切相关。在选择单元类型前，首 先你要对问题本身有非常明确的认识，然后，对于每一种单元类型，每个节点有 多少个自由度，它包含哪些特性，能够在哪些条件下使用，在ANSYS的帮助文 档中都有非常详细的描述，要结合自己的问题，对照帮助文档里面的单元描述来 选择恰当的单元类型。1. 该选杆单元(Link)还是梁单元(Beam)?这个比较容易理解。杆单元只能承受沿着杆件方向的拉力或者压力，杆单元 不能承受弯矩，这是杆单元的基本特点。梁单元则既可以承受拉，压，还可以承受弯矩。如果你的结构中要承受弯矩， 肯定不能选杆单元。对于梁单元，常用的有beam3,beam4,beam188这三种，他们的区别在于：1) beam3是2D的梁单元，只能解决2维的问题。2) beam4是3D的梁单元，可以解决3维的空间梁问题。3) beam188是3D梁单元，可以根据需要自定义梁的截面形状。2. 壁结构，是选实体单元还是壳单元？对于薄壁结构，最好是选用shell单元，shell单元可以减少计算量，如 果你非要用实体单元，也是可以的，但是这样计算量就大大增加了。实际工程中常用的shell单元有shell63,shell93。shell63是四节点的shell单元(可以退化为三角形)，shell93是带中间节点的四边形shell 单元(可以退化为三角形),shell93单元由于带有中间节点，计算精度比 shell63更高，但是由于节点数目比shell63多，计算量会增大。对于一般 的问题，选用shell63就足够了。3. 实体单元的选择。体单元类型也比较多，实体单元也是实际工程中使用最多的单元类型。常用的实体单元类型有solid45, solid92,solid185,solid187这几种。其中把solid45,solid185可以归为第一类，他们都是六面体单元，都可以退 化为四面体和棱柱体，单元的主要功能基本相同,(SOLID185还可以用于不 可压缩超弹性材料）。Solid92, solid187可以归为第二类，他们都是带中 间节点的四面体单元，单元的主要功能基本相同。如果所分析的结构比较简单，可以很方便的全部划分为六面体单元，或 者绝大部分是六面体，只含有少量四面体和棱柱体，此时，应该选用第一类 单元，也就是选用六面体单元；如果所分析的结构比较复杂，难以划分出六 面体，应该选用第二类单元，也就是带中间节点的四面体单元。前面把常用的实体单元类型归为2类了，对于同一类型中的单元，应该 选哪一种呢？通常情况下，同一个类型中，各种不同的单元，计算精度几乎 没有什么明显的差别。选取的基本原则是优先选用编号高的单元。比如第一 类中，应该优先选用solid185。第二类里面应该优先选用solid187。对于实体单元，总结起来就一句话：复杂的结构用带中间节点的四面体， 优选solid187，简单的结构用六面体单元，优选solid185。结构静力学中常用的单元类型类别形状和特性单元类型杆普通双线性LINK1 LINK8 LINK10梁普通截面渐变塑性考虑剪切变形BEAM3,BEAM4BEAM54,BEAM44BEAM23,BEAM24BEAM188,BEAM189管普通浸入塑性PIPE16,PIPE17,PIPE18PIPE59PIPE20,PIPE602-D实体四边形三角形 超弹性单元 粘弹性 大应变 谐单元 P单元PLANE42,PLANE82,PLANE182 PLANE2HYPER84,HYPER56,HYPER74VISCO88VISO106,VISO108PLANE83,PPNAE25PLANE145,PLANE1463-D实体块四面体层各向异性 超弹性单元 粘弹性 大应变 P单元SOLID45,SOLID95,SOLID73,SOLID185SOLID92,SOLID72SOLID46SOLID64,SOLID65HYPER86,HYPER58,HYPER158VISO89VISO107SOLID147,SOLID148壳四边形轴对称层剪切板P单元SHELL93,SHELL63,SHELL41,SHELL43,SHELL181SHELL51,SHELL61SHELL91,SHELL99SHELL28 ,HELL150一、典型分析过程前处理一一创建有限元模型1）2）3）4）1.单元属性定义（单元类型、实常数、材料属性） 创建或读入几何实体模型有限元网格划分2.施加载荷进行求解1）定义分析选项和求解控制2）定义载荷及载荷步选项2）求解 solve后处理1）查看分析结果2）检验结果|分析的三个主要步骤可在主菜单中得到明确体现.1.建立有限元模型2.施加载荷求解3.查看结果3.找 AN SYS Main MenuPrefevencesPreprocessor1 > Solut ion > General Postproc > T±meH ist Postpro > Design Opt > Radiation Matvix > Bun-Time Stats > Finish施加约束条件、载荷条件1、ANSYS GUI中的功能排列按照一种动宾结构，以动词开始（如Create）.随后是一个名 词（如3）.菜单的排列，按照由前到后、由简单到复杂的顺序，与典型分析的顺序相同2、ANSYS的数据库，是指在前处理、求解及后处理过程中，ANSYS保存在内存中的数据。 数据库既存储输入的数据，也存储结果数据。3、（左侧）主菜单包含ANSYS的主要功能，分为前处理、求解、后处理等。4、（左上角）应用菜单包含例如文件管理、选择、显示控制、参数设置等功能.三、前处理实体建模参数化建模体素库及布尔运算拖拉、旋转、拷贝、蒙皮、倒角等多种自动网格划分工具，自动进行单元形态、求解精度检杳及修正自由/映射网格划分、智能网格划分、自适应网格划分复杂几何体Sweep映射网格生成六面体向四面体自动过渡网格：金字塔形边界层网格划分在几何模型或FE模型上加载：点载荷、分布载荷、体载荷、函数载荷可扩展的标准梁截面形状库四、几何实体模型并不参与有限元分析.所有施加在几何实体边界上的载荷或约束必须最终传递到有限元模型上（节点或单元上）进行求解.五、由几何模型创建有限元模型的过程叫作网格划分几何模型有限元模型六、-体（3D模型）由面围成，代表三维实体.-面（表面）由线围成.代表实体表面、平面形状或壳（可以是三维曲面）.-线（可以是空间曲线）以关键点为端点，代表物体的边.-关键点（位于3D空间）代表物体的角点.七、从最低阶到最高阶，模型图元的层次关系为：关键点（Keypoints） 线（Lines）面（Areas）体（V olumes）提示：如果低阶的图元连在高阶图元上，则低阶图元不能删除.八、加载、求解ANSYS中的载荷可分为：自由度DOF -定义节点的自由度（DOF）值（结构分析_位移、热分析_温度、 电磁分析_磁势等）集中载荷-点载荷（结构分析_力、热分析_热导率、电磁分析_ magnetic current segments）面载荷-作用在表面的分布载荷（结构分析_压力、热分析_热对流、电磁分析_magnetic Maxwell surfaces 等）体积载荷-作用在体积或场域内（热分析_体积膨胀、内生成热、电磁分析_ magnetic current density 等）惯性载荷-结构质量或惯性引起的载荷（重力、角速度等）直接在实体模型加载的优点：+几何模型加载独立于有限元网格.重新划分网格或局部网格修改不影响载荷.+加载的操作更加容易，尤其是在图形中直接拾取时.Main Menu: Solution > -Loads- Apply > Pressure > On Lines说明：压力数值为正表示其方向指向表面在求解初始化前，应进行分析数据检查，包括下面内容:-统一的单位-单元类型和选项-材料性质参数-考虑惯性时应输入材料密度-热应力分析时应输入材料的热膨胀系数-实常数（单元特性）-单元实常数和材料类型的设置, 实体模型的质量特性（Preprocessor > Operate > Calc Geom Items）-模型中不应存在的缝隙,壳单元的法向-节点坐标系集中、体积载荷, 面力方向,温度场的分布和范围,热膨胀分析的参考温度（与ALPX材料特性协调?）ANSYS有两个后处理器：,通用后处理器（即“POST1”）只能观看整个模型在某一时刻的结果（如：结果的照 相 "snapshot” ）.,时间历程后处理器（即“POST26”）可观看模型在不同时间的结果。但此后处理器 只能用于处理瞬态和/或动力分析结果静力分析结果后处理的步骤主要包括：1. 绘变形图2. 变形动画3. 支反力列表4. 应力等值线图5. 网格密度检查ANSYS补充说明一、坐标系1. 工作平面坐标系叩喝：类似于绘图图板，缺省时总与总体坐标系重合，能以网格捕捉形 式显示，并可相对当前激活总体坐标系移动或旋转，其编号永远为“4”2. 总体坐标系global cs：包括三种形式总体直角坐标系(x , y , z)编号为“0”总体柱坐标系(r , 0, z)编号为“1”总体球坐标系(r , 0, 8)编号为“2”3. 局部坐标系local cs：局部坐标系是在任意位置的用户定义坐标系，即不一定与总体坐标 系平行或重合，可以是任意方向，编号为大于等于“11”4. 节点坐标系node cs：所有的力及其他方向的与节点相关的载荷都是在节点坐标系下进行 的，例如力的方向等只与节点坐标系相关节点坐标系上可以输入力和力矩；位移约束；耦合及约束过程5. 单元坐标系element cs：即材料坐标系，例如弹性模量在材料为各向异性时每一方向将 不同，此时则根据单元坐标系输入不同方向的E6. 结果坐标系：结果的输出形式位移，支反力，力矩等都是与结果坐标系相关的，结果坐 标系即当前激活坐标系，同节点坐标系一样，二者可以是任何一种当前激活坐标系二、CAD模型建模原则应考虑多少细节：如倒角和孔处，对分析无用时可忽略，但对分析目标有用，而且此处将 会出现最大应力则不能忽略是否具有对称性：包括轴，旋转，平面或镜面，重复或平移对称等。但下列因素必须对称- 几何形状；材料属性；载荷工况。此时可取一部分分析，而后叠加即可应力奇异：指在有限元模型中那些应力值无限大的点处，如点载荷的集中力和力矩作用处;孤立的约束点；尖角处等。建模时最好避免之三、网格划分器自由式free：对复杂的拓扑结构无限制，形状不定映射式mapped:拓扑结构有限制，只适用规则的体形状，如四，六面体等，可通过global set进行密度设置扫略sweep:适用于柱体形状，同mapped 一样可控制密度Smart size :智能尺寸是根据几何模型的形状，确定网格密度，适于free划分，可通过滑杆 确定网格密度四、划分原则网格划分的单元形状四方和六方的没有可比性Sweep扫略网格须上下面即对应面完全一致能用mapped,sweep划分网格最好先用之，不行再用自由式free网格划分最好按线，面，体的顺序分配单元属性千万不能分配错误面尽量用四边形的网格，体尽量用六面体的网格关心应力结果的区域须进行详细网格划分仅关心位移结果的地方网格可以粗糙些高效率建模技术在建立分析模型之前必须制订好建模方案：-必须考虑那些细节问题？-对称/反对称/轴对称？-模型中存在应力奇异？-选用那种类型的单元？. 线单元. 壳单元.X-Y平面单元-平面应力或应变单元-轴对称单元-谐单元.实体单元. 专用单元. 线性单元/高阶单元/P单元.四边形单元/三角形单元，块单元/四面体单元线弹性是基于材料的应力和应变关系是常数关系的假设一“弹性模量EX”或“杨氏模量 PERO”为常数。物理系统对称分析要求具有以下对称性条件：-几何结构对称-材料特性对称-具有零位移约束-存在非零位移约束主要单元类型举例-平面应力假定在Z方向上的应力为零，主要有以下特点：4当Z方向上的几何尺寸远远小于X和Y方向上的尺寸才有效。4所有的载荷均作用在XY平面内。4在Z方向上存在应变。4 运动只在XY平面内发生。4 允许具有任意厚度（Z方向上）。平面应力分析是用来分析诸如承受面内载荷的平板、承受压力或远离中心载荷的薄圆盘等 结构。-平面应变假定在Z方向的应变为零，主要具有以下特点：4当Z方向上的几何尺寸远远大于X和Y方向上的尺寸才有效。4所有的载荷均作用在XY平面内。4在Z方向上存在应力。4 运动只在XY平面内发生。平面应变分析是用于分析那种一个方向的尺寸（指定为总体Z方向）远远大于其它两个方 向的尺寸，并且垂直于Z轴的横截面是不变的。步骤一、建立实体模型1. 建立和修改工作文件名2. 创建标题名3. 定义单元类型BEAM （杆单元-一维）/PLANE （平面单元-二维）/SOLID （体单元-三维）4. 定义材料属性MURX相对磁导率/DENS材料密度/ EX杨氏弹性模量/ KXX热导率/PRXY泊松比步骤二、建立有限元模型步骤三、加载和求解1、定义分析类型和设置分析选项ANSYS 一共可以求解7种不同类型的分析，它们是：静态分析、瞬态分析、谐振态分 析、模态分析、频谱分析、屈曲分析和子结构分析。选择合适的求解方法：GUI: Main MenuOSolutionO Analysis TypeOSolh Controls2、求解GUI: Main MenuOSolutionOSolveOCurrent LS实例一：问题描述：悬臂梁杆一端固定，另一端为自由端。从零时刻起，给自由度施加随时间变化的应变，确 定不同时刻的应力分布。模型描述要求：该悬臂梁长5米，厚度为0.1米，弹性模量为3.08e6，泊松比为0.3。操作步骤：1、Utility MenuO New注：清除当前数据库并开始新的分析2、Utility MenuO Jobname注：输入新的工作文件名为Cantilever，并选 中 New log and error files 复选框。3、Utility MenuO FileO Change Title 注：输入 Changing pressure on one point of a pole4、Main MenuOPreprocessorOElement TypeOAdd/Edit/Delete注：选择单元类型 Structural-Solid-Quad 4node 425、Main MenuOPreprocessorOMaterial propsOMaterial Models注：在右侧材料库列表中依次双击Structural（机构） OLinear（线性）OElastic（弹性）OIsotropic （各向同性）6、Main MenuOPreprocessorOModelingOCreateO AreasORectangleOBy Dimensions注：输入悬臂梁的尺寸7、Main MenuOPreprocessorOMeshing OMeshAttributesOAll Areas注：为实体模型分配单元属性8、Main MenuOPreprocessorOMeshing OSizeCntrlsOSmartSizeOBasic注：设置智能网格划分水平，网格尺寸选为19、Main MenuOPreprocessorOMeshing OMeshOAreaOFree注：ANSYS便开始进行网格划分，划分后生成 有限元模型10、Main MenuOSolutionOAnalysis TypeONewAnalysis注：定义分析类型，选择Transient （瞬态）11、Main MenuOSolutionODefine Load （定义载荷）O Apply （加载）OStructural （结构） ODisplacement （位移）OOn Lines （在线上） 注：选择 ALL DOF,在 Displacement Value一栏中输入012、Main MenuOSolutionODefine Load （定义载荷）O Apply （加载）OStructural （结构） OForce/Moment （力 /力矩）OOn Keypoint （在关键点上）注：施加集中力载荷，拾取悬臂梁右上端的节点。在 Direction of force/mom 一栏中选择 FY， VALUE中输入-5e313、Main MenuOSolutionOLoad Step Opts （载荷步设置）OTime/Frequence （时间/频率）O Time and Substps （时间-时间步b长）注 1 ：先单击 Main MenuOSolutionOUnabridgedMenu （完整菜单），才能调出Time/Frequnc 菜单项注2：在Time at end of load step （载荷步终止时间） 中输入5；在Stepped or ramped b.c （阶跃或 斜坡）选择 Ramped;在 Maximum no. of substeps（最大子步数）一栏中输入5。14. Main MenuOSolutionOLoad Step Opts （载荷步设置）OWrite LS File （写载荷步文件） 注：创建第一个载荷步文件，在Load step filenumber n 中输入 1。15. 同理，第二次载荷步的创建，同12，13，14步。注：载荷施加方式应改为stepped （阶跃）16. 同理，第三次载荷步的创建。注：载荷大小变为-5e4。17. Main MenuOSolutionOLoad Step Opts （载荷步设置）OOutput Ctrls （输出控制）OSoluPrintout （求解应答）注：弹出Solution Printout Control （求解打印对话 框），在 Item for printout control 中选 All items （输出所有计算结果）； 在FREQ中选Every substep （每一个子步）18. 读取载荷步文件，开始求解Main MenuOSolutionOSolve （求解）OFromLS File （根据载荷步文件求解）注：LSMIN （起始载荷步）：1LSMAX （结束载荷步）：3LSINC （载荷步间隔）：119. 求解结束后，进入后处理器，查看求解结果。Main MenuOGeneral Postproc （通用后处理器）ORead Result （读入结果文件）OBy Time/Freq（根据时间/频率）注：TIME中输入5，则第5秒的求解结果读入结 果后处理器中。20. Main MenuOGeneral PostprocOPlot Results（绘制结果）OContour Plot （绘制等值图）ONodal Solu （节点解）注：选择 Nodal Solution （节点解）OElastic Strain（弹性应变）OY-Component of elastic strain （弹性应变Y分量）此求解值的大小为第5秒Y方向上的应变，应 记住，以便和之后第15秒时刻的求解结果对 比。