《ANSYS使用技巧24则-.pdf》由会员分享,可在线阅读,更多相关《ANSYS使用技巧24则-.pdf(35页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1、如何定制 Beam188/189单元的用户化截面2、如何正确理解 ANSYS 的节点坐标系3、巧用 ANSYS 的 Toolbar 4、如何使用用户定义用户自定义矩阵5、ANSYS 查询函数(Inquiry Function)6、ANSYS 是否具有混合分网的功能?7、如何实现壳单元的偏置8、如何考虑结构分析中的重力9、耦合及约束方程讲座一、耦合10、耦合及约束方程讲座二、约束方程11、膜元 Shell41 是否能作大变形分析?12、如何得到径向和周向的计算结果?13、如何加快计算速度14、ANSYS 的几种动画模式15、如何提取模态质量16、利用 ANSYS 随机振动分析功能实现随机疲劳
2、分析17、为什么在用 BEAM188 和 189 单元划分单元时会有许多额外的节点?可不可以将它们删除?18、用 ANSYS 分析过整个桥梁施工过程19、用单元死活模拟浇铸过程中的温度分布20、在 ANSYS5.6 中如何施加函数变化的表面载荷21、在 ANSYS 中怎样给面施加一个非零的法向位移约束?22、在任意面施加任意方向任意变化的压力1、如何定制 Beam188/189单元的用户化截面ANSYS 提供了几种通用截面供用户选用,但有时不能满足用户的特殊需求。为此,ANSYS 提供了用户创建截面(库)的方法。如果你需要创建一个非通用横截面,必须创建一个用户网格文件。具体方法是,首先创建一个
3、 2-D 实体模型,然后利用SECWRITE 命令将其保存(Main MenuPreprocessorSections-Beam-Write Sec Mesh)。该过程的细节如下:1.创建截面的几何模型(二维面模型)。2.对所有线设置单元份数或者单元最大尺寸(Main MenuPreprocessor-Meshing-Size Cntrls-Lines-Picked Lines或使用 MeshTool)。记住:保证模型中的区格(cell)数目不能超过 250 个。3.选择菜单 Main MenuPreprocessorSections-Beam-Write Sec Mesh,弹出一个拾取窗口,
4、单击Pick All拾取包含区格的所有面。4.ANSYS 自动在所有面上创建区格。在划分网格时,ANSYS 可能显示单元形状差的消息,也可以被忽略不显示,但是总能看到一条消息“Unable to mesh area.”。如果已经完成上述工作,你清除所有面上的单元(Main MenuPreprocessor-Meshing-Clear Areas),并重复第 2、3、4 步,即必须重新控制网格密度。5.如果第 4步成功则弹出 Write Section Library File对话框,File Name 域填入一个未用过SECT 文件名,Drives 域指定一个截面文件存放驱动器,Directo
5、ries域指定一个截面文件存放目录,然后单击按钮OK,完成用户截面文件建立。如果在第3 步中给线指定太多单元份数,区格和节点的数目可能超过限制范围,必须清除所有面上的单元24步,直到获得合适数量的区格和节点。一旦完成上述工作,用户可以在以后分析中读取用户网格文件(Main Menu Preprocessor Sections -Beam-Read Sect Mesh),定义成适当的截面ID 号,就与通用梁完全一致。要想检查用户截面,只要绘制截面(网格)图(Main MenuPreprocessorSectionsPlot Section)或者列表截面属性(Main MenuPreprocess
6、orSectionsList Sections)。2、如何正确理解ANSYS 的节点坐标系节点坐标系用以确定节点的每个自由度的方向,每个节点都有其自己的坐标系,在缺省状态下,不管用户在什么坐标系下建立的有限元模型,节点坐标系都是与总体笛卡尔坐标系平行。有限元分析中的很多相关量都是在节点坐标系下解释的,这些量包括:输入数据:1 自由度常数2 力3 主自由度4 耦合节点5 约束方程等输出数据:1 节点自由度结果2 节点载荷3 反作用载荷等但实际情况是,在很多分析中,自由度的方向并不总是与总体笛卡尔坐标系平行,比如有时需要用柱坐标系、有时需要用球坐标系等等,这些情况下,可以利用ANSYS 的“旋转节
7、点坐标系”的功能来实现节点坐标系的变化,使其变换到我们需要的坐标系下。具体操作可参见 ANSYS 联机帮助手册中的“分析过程指导手册-建模与分网指南-坐标系-节点坐标系”中说明的步骤实现。3、巧用 ANSYS 的 Toolbar ANSYS 在图形界面方面的最大特点是用户可以根据自己的目的定制自己的图形界面。尽管ANSYS 提供了专门的菜单开发工具UIDL语言,但毕竟只有少数熟练的用户能很好地使用。但有一种方法特别适合初学者使用,那就是通过ANSYS 提供的 Toolbar,把 ANSYS 常用的命令和宏定制成按钮,这样需要反复地点取菜单来执行的常用命令,通过点击一次 Toolbar 按钮就可
8、以完成了。定制按钮的方法有两种:菜单方式与命令方式。菜单方式就是直接点菜单 Utility MenuMacroEdit Abbreviations,或Utility MenuMacroEdit Toolbar,;而命令方式就是直接执行命令*ABBR(Abbreviation的缩写)。ANSYS 的 db 文件自动保存你所增加的按钮,但是每次重新开始新的分析,都需要重新定义每个按钮。如果想要在每次启动 ANSYS 或开始新的分析时都能自动加载这些快捷按钮,通常是把这些按钮对应的*ABBR命令放在 ANSYS 的启动文件 star56.ans(对ANSYS5.6)中,ANSYS 在每次启动时会自动
9、加载这些命令。ANSYS 的启动文件 star56.ans在 ANSYS 安装目录下的 DOCU 目录中,如c:ansys56docustar56.ans,它是一个文本文件,可以用通用的文本编辑器打开和编辑它。其实该文件提供了一个内容很丰富的模板,它收集了许多 ANSYS 常用命令的缩略,只不过都把它们注释掉了,你只需去掉前面的注释符(!),就可以在进入 ANSYS 后在 Toolbar 中看到它们。你也可以把所有的*ABBR 命令放在一个单独的文件中,然后在启动文件中加入一条命令ABBRES,用 ABBRES 命令的好处是你在任何时候都可以加载这些快捷按钮,这对于你如果使用由别人建的模型,而
10、他用的是不同的 Toolbar,这是因为 Toolbar 是保存在模型文件中(jobname.db),你可以发现采用这样的方法给你带来很大的方便。我承认我是一个 Toolbar 迷,把许多命令和宏放在Toolbar 中使我的工作更快和更有效率。在我的Toolbar 中,大多数按钮是直接执行ANSYS 命令和宏,其中也有一些按钮是用来调入新的Toolbar 实现一些新的功能。通过将一个Toolbar 嵌入到另一个 Toolbar 中,可以将Toolbar 分成前处理和后处理等大类,也可以列出一些特别的信息有助于工作。不过我个人的倾向不赞成用太多的嵌套Toolbar,因为点三到四层按钮与点三到四级
11、菜单差不多。宁可在一个 Toolbar 中有许多按钮,也不要有许多嵌套Toolbar。包含 43个按钮的toolbar 减少了11 个按钮的toolbar 尽管最好能把所有需要的按钮一次全部调进来,但有时发现这也不是好事。有一天我数了一下Toolbar 中的按钮,竟然发现有43 个按钮,尽管这不是一个世界记录,但这确实太多了,此时我才意识到我这个 toolbar迷困扰于过多的按钮混乱中。得想办法结束这种状态。我决定摆脱这种按钮迷的状态,开始精简我的toolbar,去掉那些不太常用的按钮。我开始研究toolbar,结果发现没有一个按钮在工作中不太常用。我几乎每天都要用到所有这些命令,少了任何一个
12、,我都要花时间到菜单堆里去找这些命令。我想保持原有 toolbar中的所有功能,但我要把按钮数减下来,我能做得到吗?再瞧一下 toolbar,我发现打开和关闭六种实体(keypoints,lines,areas,volumes,和 elements)的编号显示共用了 12 个按钮,此外,还有许多按钮用来打开和关闭显示节点和单元坐标系,线的方向和局部坐标系等;还有我用单独的按钮控制矢量模式和光栅模式显示。结果发现有22 个按钮用来决定实体的编号和符号以及显示模式的开关。如果我能把每组开和关的按钮用一个按钮来代替,那么我可以把按钮数减少11 个,这是一个很可观的数字。我决定写一个宏,它能够确定当前
13、的设置,然后翻转该设置,我给它取名叫 numtog.mac。在确定了程序的框架后,我开始查看关于*GET命令的帮助文档,看看是否能得到/PNUM,/PSYMB 和显示模式的设置状态。令我失望的是,我发现尽管*GET命令可以设置显示模式,但不能得到关于/PNUM 和/PSYMB 的设置,没有它我就无法实现我的宏,失望之感油然而生。非常时候会有非常手段,文档中没有说明由*GET得到的量并不表示一定不可以获得,也许通过我的努力可以得到。我决定从其它地方找,在这以前我经常用此办法找到其它信息,我想到了控制 ANSYS 菜单系统的文件,因为 ANSYS 在建立关于/PNUM 和/PSYMB 的设置时也需
14、要相关的信息,也许它以未公开的方式使用。我开始查找 UIFUNC1.GRN 和UIFUNC2.GRN,果然在 UIFUNC2.GRN 文件中找到了我要的信息,在该文件中我找到了控制编号和符号显示的命令,它也是由*GET命令获得,只不过时以一种特殊的方式,我就把相关的命令拷贝和粘帖到我的宏文件中,以节省时间。尽管不推荐使用此类*GET命令,但有些高级的APDL 编程人员经常使用。问题是未公开的*GET命令在 5.6 中能使用,但不能保证在5.7 中也一定能使用,这只有等你试了以后才知道。让我们来看一个例子,它用来确定keyponits编号的设置状态,格式如下:*get,ar20,common,d
15、3com,int,104 从这个命令格式我们可以推断一些信息:其中的“Common”是指储存在公共数据块中的数据,“D3com”是该公共数据块的名字,“Int”指数值数据以整数格式存储,“104”是我们所要数据存储的位置。所以上面这条命令是告诉ANSYS 提取在公共数据块“D3com”104 位置的整型数据,并把它赋给变量ar20,通过试验知道,当keypoints编号显示为关闭时ar20 为 0,而开时为 1,即用 1 和 0 表示开和关。我的新 toolbar已用了几星期了,很方便,它比旧的toolbar减少了 25%的按钮。现在我的toolbar小多了,我还可以增加其它一些命令,如控制/
16、EDGE 的开和关等等。相信你会得到一些收益。4、如何使用用户定义用户自定义矩阵Matrix 27 用户定义用户自定义矩阵,由单元选项控制定义质量、刚度或阻尼矩阵,你只要在同一组接点,分别定义三次MATRIX27 单元(KEYOPT(2)分别为2,4,5)即可,然后在定义实常数时,分别定义三种单元对应的质量、刚度、阻尼矩阵系数。5 ANSYS 查询函数(Inquiry Function)在 ANSYS 操作过程或条件语句中,常常需要知道有关模型的许多参数值,如选择集中的单元数、节点数,最大节点号等。此时,一般可通过*GET命令来获得这些参数。现在,对于此类问题,我们有了一个更为方便的选择,那就
17、是查询函数 Inquiry Function。Inquiry Function类似于 ANSYS 的*GET 命令,它访问ANSYS数据库并返回要查询的数值,方便后续使用。ANSYS 每执行一次查询函数,便查询一次数据库,并用查询值替代该查询函数。假如你想获得当前所选择的单元数,并把它作为*DO循环的上界。传统的方法是使用*GET命令来获得所选择的单元数并把它赋给一个变量,则此变量可以作为*DO循环的上界来确定循环的次数*get,ELMAX,elem,count*do,I,1,ELMAX*enddo 现在你可以使用查询函数来完成这件事,把查询函数直接放在*DO循环内,它就可以提供所选择的单元数
18、*do,I,ELMIQR(0,13)*enddo 这里的 ELMIQR 并不是一个数组,而是一个查询函数,它返回的是现在所选择的单元数。括弧内的数是用来确定查询函数的返回值的。第一个数是用来标识你所想查询的特定实体(如单元、节点、线、面号等等),括弧内的第二个数是用来确定查询函数返回值的类型的(如选择状态、实体数量等)。同本例一样,通常查询函数有两个变量,但也有一些查询函数只有一个变量,而有的却有三个变量。查询函数的种类和数量很多,下面是一些常用、方便而快速快捷的查询函数1 AREA arinqr(areaid,key)areaid 查询的面,对于key=12,13,14 可取为 0;key标
19、识关于 areaidr 的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数=-1,材料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数arinqr(areaid,key)的返回值对于 key=1=0,areaid未定义=-1,areaid 未被选择=1,areaid被选择2 KEYPOINTS kpinqr(kpid,key)kpid 查询的关键点,对于key=12,13,14 为 0 key 标识关于 kpid 的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数目=-1,数料号=-2,单元类型=-3,实常数=-4,节点数
20、,如果已分网=-7,单元数,如果已分网kpinqr(kpid,key)的返回值对于 key=1=-1,未选择=0,未定义=1,选择3 LINElsinqr(lsid,key)lsid 查询的线段,对于key=12,13,14 为 0 key标识关于 lsid的返回信息=1,选择状态=2,长度=12,定义的数目=13,选择的数目=14,定义的最大数=-1,材料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数4 NODE ndinqr(node,key)node节点号,对于key=12,13,14 为 0 key标识关于 node的返回信息=1,选择状态=12,定义的数目=13,选择的
21、数目=14,定义的最大数=-2,超单元标记=-3,主自由度=-4,激活的自由度=-5,附着的实体模型ndinqr(node,key)的返回值对于 key=1=-1,未选择=0,未定义=1,选择5 VOLUMESvlinqr(vnmi,key)vnmi查询的体,对于key=12,13,14 为 0 key标识关于 vnmi 的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数目=-1,数料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数=-8,单元形状=-9,中节点单元=-10,单元坐标系vlinqr(vnmi,key)的返回值对于 key=1=-1,未选
22、择=0,未定义=1,选择6 ANSYS是否具有混合分网的功能?ANSYS 具有混合网格剖分的功能。例如两个粘在一起的面,可以对一个面进行三角形划分,再对另一个面进行四边形划分。过程见下列命令:/prep7 et,1,42 rect,1,1 rect,1,2,1 aglue,all mshape,0,2d amesh,1 mshape,1,2d amesh,3 7 如何考虑结构分析中的重力在结构分析中,如何模拟结构自重和设备重量是一个经常遇到的问题,对于结构自重有两点要注意:1 在材料性质中输入密度,如果不输入密度,则将不会产生重力效果。2 因为 ANSYS 将重力以惯性力的方式施加,所以在输入
23、加速度时,其方向应与实际的方向相反。对于结构上的设备重量可以用MASS21 单元来模拟,该单元为一个空间“点”单元。设备重量可通过单元实常数来输入。下面附上一个小例子(设重力方向向下)。/prep7 et,1,42 et,2,21 r,2,10,10,10 mp,ex,1,2e5 mp,nuxy,1,0.3 mp,dens,1,1 rect,10,1 esize,.5 amesh,all type,2 real,2 e,node(5,1,0)fini/solu dk,1,all dk,2,uy,acel,10 solve fini/post1 plnsol,u,sum,2 8、如何实现壳单元的
24、偏置在构建诸如汽车车身、飞机引擎罩、发动机气道以及其他薄壁结构的有限元模型时通常选用壳单元,ANSYS 提供的 SHELL43、SHELL63、SHELL93 以及 5.5 版开始新增加的SHELL143 和 SHELL181 单元都能很方便地用于薄壁结构的建模,并能得到精确的结果,同时避免了采用实体单元构建此类结构所带来的复杂性和规模的庞大。像 SOLID45 和SOLID95之类的实体单元,其厚度由节点的位置所决定,而壳单元的厚度是由实常数决定的。由于壳单元的厚度独立于几何模型,因此可以很方便地进行参数化和优化设计。此外,建立一个用于生成壳单元的面几何模型要比建一个三维实体几何模型简单得多
25、。尽管使用壳单元有许多优点,但也有一些缺点在使用时需要加以注意。首先,壳单元只适合某几种结构类型,因为它假设弯曲应力沿壳单元厚度方向线性变化;其次,为了得到理想的结果,壳单元必须建在实体结构的中心层或中面上。图 1 显示了如何用壳单元构建一个工字截面梁,左边是实体模型,右边是对应的壳单元模型。如果壳单元用实体模型的顶面或底面生成,那么有限元模型的刚度与实际结构的刚度将相差很大。图 1.工字型梁的三维实体和中面壳单元模型与模型的复杂程度有关,有的时候很难或不可能得到中面,许多时候分析人员只能用模型的顶面或底面,希望壳单元的偏置影响较小可以忽略,这样做有的场合可以,有时则不行。而且,分析人员从来不
26、知道这样做带来的误差有多大。关于这类问题的一种解决办法是用SHELL91 层状壳单元替代标准壳单元。SHELL91 主要用于复合材料分析,当然也可以用于传统材料的分析。采用SHELL91 的好处是可以通过keyopt(11)来控制厚度朝单元的顶面或底面偏置。比如:当SHELL91 的 keyopt(11)=1时表示朝单元底面偏置,而当keyopt(11)=2时朝单元顶面偏置,而当keyopt(11)=0时则表示不发生偏置。由于 SHELL91 也带中间节点,因此它可以直接替换SHELL93 单元。为了评价按这种方式使用SHELL91 的效果,我分别用SOLID95、SHELL93 和 SHEL
27、L91 来分析不同壁厚的弯管如图2 所示,注意只改变管的厚度,而管的直径和弯肘半径保持不变。比值 t/r,即壁厚/弯管半径用来表示中面偏置的程度。由SOLID95所得到的结果作为基准来归一化SHELL93 和SHELL91 模型得到的结果。所有情况都受两个载荷步:第一步在管的顶端受一个分布拉力,方向偏离固定端,使得在肘的内侧产生很大的拉应力,这个载荷条件用于确定 ANSYS 计算的应力是否正确;第二个载荷步是在同方向施加一个加速度,用它可以确定用带偏置的SHELL91 单元时 ANSYS 能否正确考虑质量。图 2.弯管实体模型及对应的SOLID95划分的有限元模型由图 3 的结果显示可知:当
28、SHELL93 和 SHELL91 放在中面上时,对不同的 t/r值所得结果的误差与用SOLID95单元相同。当壳单元用管的外侧生成时,SHELL93 预计的误差随t/r的增加而减少,而由SHELL91 在激活 keyopt(11)使向内偏置时得到与SOLID95相同的结果,同样,当壳单元用管的内侧生成时,SHELL93 预计的误差随 t/r的增加而增加,而由SHELL91 在激活 keyopt(11)使向外偏置时得到与 SOLID95相同的结果。最大的平均Von Mises 应力(PLNS,S,EQV)发生在肘管的内侧,如图4 所示,随着 t/r的增大,SHELL93 用在管的内侧和外侧时得
29、到结果的精度显然没有用SHELL91 采用合适的偏置所得到的结果好。在第二载荷步施加加速度所的结果与载荷步一时的情况相同,由此可以得出结论:在中面得不到的情况下,用SHELL91 替代 SHELL93处理比较好。由于我经常用CAD 软件来得到我使用的几何模型,而很少有 CAD 软件能直接提取中面,因此使用带偏置的SHELL91 单元是我经常采用的办法。图3.不同单元产生的变形误差图4.不同单元产生的应力误差9、耦合及约束方程讲座一、耦合当需要迫使两个或多个自由度取得相同(但未知)值,可以将这些自由度耦合在一起。耦合自由度集包含一个主自由度和一个或多个其它自由度。典型的耦合自由度应用包括:模型部
30、分包含对称;在两重复节点间形成销钉、铰链、万向节和滑动连接;迫使模型的一部分表现为刚体。如何生成耦合自由度集1.在给定节点处生成并修改耦合自由度集命令:CP GUI:Main MenuPreprocessorCoupling/CeqnCouple DOFs 在生成一个耦合节点集之后,通过执行一个另外的耦合操作(保证用相同的参考编号集)将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用 CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集(即增、删节点或改变自由度标记)可用CPNGEN命令。(不能由 GUI直接得到 CPNBGEN命令)。2.耦合重合节点。CPIN
31、TF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点(诸如一条缝处)尤为有用。命令:CPINTF GUI:Main MenuPreprocessorCoupling/CeqnCoincident Nodes 3.除耦合重复节点外,还可用下列替换方法迫使节点有相同的表现方式:o如果对重复节点所有自由度都要进行耦合,常用NUMMRG命 令(GUI:Main MenuPreprocessorNumbering CtrlsMerge Items)合并节点。o可用 EINTF命令(GUI:Main Menu PreprocessorCreate E
32、lements At Coincid Nd)通在重复节点对之间生成2节点单元来连接它们。o用CEINTF 命 令(GUI:Main MenuPreprocessor Coupling/Ceqn Adjacent Regions)将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。生成更多的耦合集一旦有了一个或多个耦合集,可用这些方法生成另外的耦合集:1.用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。命令:CPLGEN GUI:Main MenuPreprocessorCoupling/CeqnGen w/Same
33、 Nodes 2.用下列方法生成与已有耦合集不同(均匀增加的)节点编号但有相同的自由度标记的新的耦合集:命令:CPSGEN GUI:Main MenuPreprocessorCoupling/CeqnGen w/Same DOF 使用耦合注意事项1.每个耦合的节点都在节点坐标系下进行耦合操作。通常应当保持节点坐标系的一致性。2.自由度是在一个集内耦合而不是集之间的耦合。不允许一个自由度出现在多于一个耦合集中。3.由 D或共它约束命令指定的自由度值不能包括在耦合集中。4.在减缩自由度分析中,如果主自由度要从耦合自由度集中选取,只有主节点的自由度才能被指定为主自由度。5.在结构分析中,耦合自由度以
34、生成一刚体区域有时会引起明显的平衡破坏。不重复的或不与耦合位移方向一致的一个耦合节点集会产生外加力矩但不出现在反力中。10、耦合及约束方程讲座二、约束方程约束方程提供了比耦合更通用的联系自由度的方法。有如下形式:这里 U(I)是自由度,N是方程中项的编号。如何生成约束方程1.直接生成约束方程o直接生成约束方程:命令:CE GUI:Main MenuPreprocessorCoupling/CeqnConstraint Eqn 下面为一个典型的约束方程应用的例子,力矩的传递是由 BEAM3 单元与 PLANE42 单元(PLANE42 单元无平面转动自由度)的连接来完成的:图 12-1 建立旋转
35、和平移自由度的关系如果不用约束方程则节点2 处表现为一个铰链。下述方法可在梁和平面应力单元之间传递力矩,自由度之间满足下面的约束方程:ROTZ2=(UY3-UY1)/10 0=UY3-UY1-10*ROTZ2相应的 ANSYS 命令为:CE,1,0,3,UY,1,1,UY,-1,2,ROTZ,-10 o修改约束方程在 PREP7 或 SOLUTION 中修改约束方程中的常数项:命令:CECMOD GUI:Main MenuPreprocessorCoupling/CeqnModify ConstrEqn Main MenuPreprocessorLoadsOtherModify ConstrE
36、qn Main MenuSolutionOtherModify ConstrEqn 如果要修改约束方程中的其它项,必须在求解前在PREP7 中用使 CE命令(或相应 GUI途径)。自动生成约束方程o生成刚性区域CERIG 命令通过写约束方程定义一个刚性区域。通过连接一主节点到许多从节点来定义刚性区。(此操作中的主自由度与减缩自由度分析的主自由度是不同的)命令:CERIG GUI:Main MenuPreprocessorCoupling/CeqnRigid Region 将 CERIG 命令的 Ldof 设置为 ALL(缺省),此操作将为每对二维空间的约束节点生成三个方程。这三个方程在总体笛卡
37、尔空间确定三个刚体运动(UX、UY、ROTZ)。为在二维模型上生成一个刚性区域,必须保证XY平面为刚性平面,并且在每个约束节点有UX、UY和 ROTZ 三个自由度。类似地,此操作也可在三维空间为每对约束节点生成六个方程,在每个约束节点上必须有(UX、UY、UZ、ROTX、ROY 和 ROTZ)六个自由度。输入其它标记的Ldof 域将有不同的作用。如果此区域设置为 UXYZ,程序在二维(X,Y)空间将写两个约束方程,而在三维空间(X、Y、Z)将写三个约束方程。这些方程将写成从节点的平移自由度和主节点的平移和转动自由度。类似地,RXYZ 标记允许生成忽略从节点的平移自由度的部分方程。其它标记的Ld
38、of 将生成其它类型的约束方程。总之,从节点只需要由Ldof 标记的自由度,但主节点必须有所有的平移和转动自由度(即二维的 UX、UY和 ROTZ;三维的 UX、UY、UZ、ROTX、ROTY、ROTZ)。对由没有转动自由度单元组成的模型,应当考虑增加一个虚拟的梁单元以在主节点上提供旋转自由度。将疏密不同的已划分网格区域连在一起可将一个区域(网格较密)的已选节点与另一个区域(网格较稀)的已选单元用 CEINTF命令(菜单途径Main MenuPreprocessorCoupling/CeqnAdjacent Regions)连起来生成约束方程。这项操作将不相容网格形式的区域“系”在一起。在两区
39、域的交界处,从网格稠密的区域选择节点A,从网格粗糙区域选择单元B,用区域B单元的形函数,在相关的区域 A和B界面的节点处写约束方程。ANSYS允许这些节点位置使用两公差准则。节点在单元之外超过第一公差就认为节点不在界面上。节点贴近单元表面的距离小于第二公差则将节点移到表面上,见下图。o对 CEINTF命令有些限制:应力或热通量可能会不连续地穿过界面。界面区域的节点不能指定位移。可用每节点有六个自由度的单元接合6 自由度实体。o从已有约束方程集生成约束方程集可用 CESGEN 命令从已有约束方程集生成约束方程。那么已有约束方程集内的节点编号将增加以生成另外的约束方程集。另外约束方程集的标记和系数
40、保持与原集的一致。命令:CESGEN GUI:Main MenuPreprocessorCoupling/Ceqn Gen w/same DOF 使用约束方程的注意事项所有的约束方程都以小转动理论为基础。因此,它应用在大转动分析中,NLGEOM?应当限制在约束方程所包含的自由度方向无重大变化的情况。约束方程的出现将产生不可预料的反力和节点力结果。由于相邻区域网格疏密不同,边界上的相容性仍然存在。但是当网格越密,这种不相容的危害就越小。11、膜元 Shell41 是否能作大变形分析?膜元 Shell41 可以用于大变形分析。但是膜元 shell41 在处理大变形问题时必须采用三角形单元,因为四边
41、形单元会产生跷曲(warping),所以在划分网格时请选Tri。12、如何得到径向和周向的计算结果?在圆周对称结构中,如圆环结构承受圆周均布压力。要得到周向及径向位移,可在后处理/POST1中,通过菜单 General PostprocOptions for OutpRsysGlobal cylindric 或命令 Rsys,1 将结果坐标系转为极坐标,则X方向位移即为径向位移,Y向位移即为周向位移。13、如何加快计算速度在大规模结构计算中,计算速度是一个非常重要的问题。下面就如何提高计算速度作一些建议:1.充分利用 ANSYS MAP 分网和 SWEEP 分网技术,尽可能获得六面体网格,这一
42、方面减小解题规模,另一方面提高计算精度。2.在生成四面体网格时,用四面体单元而不要用退化的四面体单元。比如 95 号单元有 20节点,可以退化为 10 节点四面体单元,而 92 号单元为 10 节点单元,在此情况下用92 号单元将优于95 号单元。选择正确的求解器。对大规模问题,建议采用 PCG 法。此法比波前法计算速度要快10 倍以上(前提是您的计算机内存较大)。对于工程问题,可将 ANSYS 缺省的求解精度从1E-8 改为 1E-4 或 1E-5 即可。14、ANSYS 的几种动画模式Q:ANSYS中有好几种动画模式,但并非同时可用。那么有哪几种类型,何时可用呢?A:ANSYS 中有四种动
43、画模式:Bitmap,AVI,Display List,Pixmap。该种模式是否可用取决于运行环境是UINX还是 PC,以及是使用二维还是三维 显 示 设 备。可 以 通 过 选择 菜 单 Utility MenuPlotCtrlsDevice Options(或键入适当的命令)来制作动画。下面是该命令选项的小结:二维显示设备的PC:Bitmap(命令/device,anim,bmp)AVI(命令/device,anim,avi)三维显示设备的PC:Bitmap(命令/dv3d,anim,1)AVI(命令/dv3d,anim,2)Display List(命令/dv3d,anim,0)二维显
44、示设备的UNIX:Pixmap(缺省,无相应命令)三维显示设备的UNIX:Pixmap(命令/dv3d,anim,1)Display List(命令/dv3d,anim,0)Q:这些动画模式的区别是什么?A:Display List选项(三维设备可用)与其他方式的不同在于:Display List允许在播放动画的过程中对模型进行动态操作(如放大缩小等),而其它模式不能;播放器不同:Bitmap,Display List,Pixmap 动画模式在 ANSYS 图形窗口进行播放,可以通过ANSYS动画控制器面板控制动画;AVI(只在PC上)用 WINDOWS媒体播放器播放。需注意在 ANSYS5.
45、5中缺省的动画模式为AVI,播放器为媒体播放器。而在 ANSYS5.6 中二维设备的缺省动画模式为Bitmap,三维设备的缺省动画模式为 Display List,动画在 ANSYS 图形窗口播放。Q:如何存储在 ANSYS 中创建的动画?A:对于 PC环境,Bitmap 或 AVI 模式的动画会自动存储为jobname.avi(jobname 是你指定的分析名称)。Display List动画在 PC上不会自动存储,必须通过菜单 Utility Menu PlotCtrls AnimateSave Animation指定存储,典型格式为jobname.anim。对于 UNIX环境,Pixma
46、p及 Display List都只能按 ANSYS 指定的格式进行存储,方法同上。Q:A 用 Bitmap 模式存储的.avi文件与用 AVI模式存储的.avi文件有何不同?A:Bitmap模式的文件比 AVI 模式的文件要小,因为采用Bitmap 模式存储的.avi文件不包括回放的画面,它在 ANSYS 动画控制器中播放,可以向前或向后显示。而WINDOWS媒体播放器没有forward/backward选项,回放的画面必须包括在AVI 动画文件中,因此文件较大。Q:UNIX 上制作的动画能否在PC上播放?A:可以。可以使用 ANSYS ANIMATE 程序(ANIMATE.exe)在 PC上
47、播放.anim或.avi文件。还可以将.anim 文件转变为.avi文件。所有的ANSYS 产品都提供 ANIMATE 程序,该程序位于ANSYS 安装盘的根目录中。Q:ANIMATE是否能读入所有的动画文件?A:不能。ANIMATE 程序接受 UNIX上的 Pixmap模式的.anim 文件或 PC上的 Bitmap 或 AVI模式的.avi 文件。15、如何提取模态质量模态分析过程中打开振型型则化开关(MODOPT 命令的 Nrmkey设置为ON),ANSYS 程序将自动将每阶模态的最大位移单位化,就可以提取模态质量。计算方法如下:1、利用 SSUM 对 ETABLE 动能数据求和获得结构
48、总动能();2、将结构总动能除以得到,其中是系统的角频率。下面是 ANSYS Verification Manual中 VM89.DAT 稍加修改后提取模态质量的例子:/PREP7/TITLE,VM89,NATURAL FREQUENCIES OF A TWO-MASS-SPRING SYSTEM C*VIBRATION THEORY AND APPLICATIONS,THOMSON,2ND PRINTING,PAGE 163,EX 6.2-2 ET,1,COMBIN14,2 ET,2,MASS21,4 R,1,200!SPRING CONSTANT=200 R,2,800!SPRING CO
49、NSTANT=800 R,3,.5!MASS=.5 R,4,1!MASS=1 N,1 N,4,1 FILL E,1,2!SPRING ELEMENT(TYPE,1)AND K=200(REAL,1)TYPE,2 REAL,3 E,2!MASS ELEMENT(TYPE,2)AND MASS=.5(REAL,3)TYPE,1 REAL,2 E,2,3!SPRING ELEMENT(TYPE,1)AND K=800(REAL,2)TYPE,2 REAL,4 E,3!MASS ELEMENT(TYPE,2)AND MASS=1(REAL,4)TYPE,1 REAL,1 E,3,4!SPRING EL
50、EMENT(TYPE,1)AND K=200(REAL,1)M,2,UX,3 OUTPR,BASIC,1 D,1,UY,4 D,1,UX,4,3 FINISH/SOLU ANTYPE,MODALMODOPT,subspa,2,2,ON MXPAND,2,YES SOLVE FINISH/post1 set,1,1 etabl,kene,kene ssum*get,keneval1,ssum,item,kene*get,freqval1,mode,1,freq eigen1=(2*3.14159*freqval1)*2 pmass1=2*keneval1/eigen1 set,1,2 etabl