ansys使用技巧|隧道及地下工程 - 领先的岩土技术社区，岩土领域的媒体、社区与应用平台！

离线lg3060631135

助理工程师

发帖: 42

土币: 480

威望: 70

原创币: 0

只看楼主倒序阅读使用道具楼主发表于: 2010-10-20

1、如何定制Beam188/189单元的用户化截面
ANSYS提供了几种通用截面供用户选用，但有时不能满足用户的特殊需求。为此，ANSYS提供了用户创建截面（库）的方法。如果你需要创建一个非通用横截面，必须创建一个用户网格文件。具体方法是，首先创建一个2-D实体模型，然后利用SECWRITE命令将其保存(Main Menu>Preprocessor>Sections> -Beam-Write Sec Mesh)。该过程的细节如下：

1.             创建截面的几何模型（二维面模型）。

2.             对所有线设置单元份数或者单元最大尺寸 (Main Menu>Preprocessor> -Meshing-Size Cntrls>-Lines-Picked Lines或使用MeshTool)。记住：保证模型中的区格（cell）数目不能超过250个。

3.             选择菜单Main Menu>Preprocessor>Sections>-Beam-Write Sec Mesh，弹出一个拾取窗口，单击Pick All拾取包含区格的所有面。

4.             ANSYS自动在所有面上创建区格。在划分网格时，ANSYS可能显示单元形状差的消息，也可以被忽略不显示，但是总能看到一条消息“Unable to mesh area....”。如果已经完成上述工作，你清除所有面上的单元(Main Menu>Preprocessor>-Meshing-Clear> Areas)，并重复第2、3、4步，即必须重新控制网格密度。

5.             如果第4步成功则弹出Write Section Library File对话框，File Name域填入一个未用过SECT文件名，Drives域指定一个截面文件存放驱动器，Directories域指定一个截面文件存放目录，然后单击按钮OK，完成用户截面文件建立。如果在第3步中给线指定太多单元份数，区格和节点的数目可能超过限制范围，必须清除所有面上的单元2~4步，直到获得合适数量的区格和节点。

   一旦完成上述工作，用户可以在以后分析中读取用户网格文件 (Main Menu > Preprocessor > Sections > -Beam-Read Sect Mesh)，定义成适当的截面ID号，就与通用梁完全一致。要想检查用户截面，只要绘制截面（网格）图（Main Menu>Preprocessor>Sections>Plot Section）或者列表截面属性（Main Menu>Preprocessor>Sections>List Sections）。

2、如何正确理解ANSYS的节点坐标系
节点坐标系用以确定节点的每个自由度的方向，每个节点都有其自己的坐标系，在缺省状态下，不管用户在什么坐标系下建立的有限元模型，节点坐标系都是与总体笛卡尔坐标系平行。有限元分析中的很多相关量都是在节点坐标系下解释的，这些量包括：
输入数据：
1 自由度常数
2 力
3 主自由度
4 耦合节点
5 约束方程等
输出数据：
1 节点自由度结果
2 节点载荷
3 反作用载荷等



但实际情况是，在很多分析中，自由度的方向并不总是与总体笛卡尔坐标系平行，比如有时需要用柱坐标系、有时需要用球坐标系等等，这些情况下，可以利用ANSYS的“旋转节点坐标系”的功能来实现节点坐标系的变化，使其变换到我们需要的坐标系下。具体操作可参见ANSYS联机帮助手册中的“分析过程指导手册->建模与分网指南->坐标系->节点坐标系”中说明的步骤实现。

3、巧用ANSYS的Toolbar
ANSYS在图形界面方面的最大特点是用户可以根据自己的目的定制自己的图形界面。尽管ANSYS提供了专门的菜单开发工具UIDL语言，但毕竟只有少数熟练的用户能很好地使用。但有一种方法特别适合初学者使用，那就是通过ANSYS提供的Toolbar，把ANSYS常用的命令和宏定制成按钮，这样需要反复地点取菜单来执行的常用命令，通过点击一次Toolbar按钮就可以完成了。
   定制按钮的方法有两种:菜单方式与命令方式。菜单方式就是直接点菜单Utility Menu>Macro>Edit Abbreviations …或Utility Menu>Macro>Edit Toolbar …;而命令方式就是直接执行命令*ABBR(Abbreviation的缩写)。ANSYS的db文件自动保存你所增加的按钮，但是每次重新开始新的分析，都需要重新定义每个按钮。如果想要在每次启动ANSYS或开始新的分析时都能自动加载这些快捷按钮，通常是把这些按钮对应的*ABBR命令放在ANSYS的启动文件star56.ans(对ANSYS5.6)中，ANSYS在每次启动时会自动加载这些命令。ANSYS的启动文件star56.ans在ANSYS安装目录下的DOCU目录中，如c:\ansys56\docu\star56.ans，它是一个文本文件，可以用通用的文本编辑器打开和编辑它。其实该文件提供了一个内容很丰富的模板，它收集了许多ANSYS常用命令的缩略，只不过都把它们注释掉了，你只需去掉前面的注释符(!)，就可以在进入ANSYS后在Toolbar中看到它们。你也可以把所有的*ABBR命令放在一个单独的文件中，然后在启动文件中加入一条命令ABBRES，用ABBRES命令的好处是你在任何时候都可以加载这些快捷按钮，这对于你如果使用由别人建的模型，而他用的是不同的Toolbar，这是因为Toolbar是保存在模型文件中(jobname.db)，你可以发现采用这样的方法给你带来很大的方便。
   我承认我是一个Toolbar迷，把许多命令和宏放在Toolbar中使我的工作更快和更有效率。在我的Toolbar中，大多数按钮是直接执行ANSYS命令和宏，其中也有一些按钮是用来调入新的Toolbar实现一些新的功能。通过将一个Toolbar嵌入到另一个Toolbar中，可以将Toolbar分成前处理和后处理等大类，也可以列出一些特别的信息有助于工作。不过我个人的倾向不赞成用太多的嵌套Toolbar，因为点三到四层按钮与点三到四级菜单差不多。宁可在一个Toolbar中有许多按钮，也不要有许多嵌套Toolbar。

包含43个按钮的toolbar

减少了11个按钮的toolbar

   尽管最好能把所有需要的按钮一次全部调进来，但有时发现这也不是好事。有一天我数了一下Toolbar中的按钮，竟然发现有43个按钮，尽管这不是一个世界记录，但这确实太多了，此时我才意识到我这个toolbar迷困扰于过多的按钮混乱中。得想办法结束这种状态。
我决定摆脱这种按钮迷的状态，开始精简我的toolbar，去掉那些不太常用的按钮。我开始研究toolbar，结果发现没有一个按钮在工作中不太常用。我几乎每天都要用到所有这些命令，少了任何一个，我都要花时间到菜单堆里去找这些命令。我想保持原有toolbar中的所有功能，但我要把按钮数减下来，我能做得到吗？
   再瞧一下toolbar，我发现打开和关闭六种实体（keypoints,lines,areas,volumes,和elements）的编号显示共用了12个按钮，此外，还有许多按钮用来打开和关闭显示节点和单元坐标系，线的方向和局部坐标系等；还有我用单独的按钮控制矢量模式和光栅模式显示。结果发现有22个按钮用来决定实体的编号和符号以及显示模式的开关。如果我能把每组开和关的按钮用一个按钮来代替，那么我可以把按钮数减少11个，这是一个很可观的数字。
我决定写一个宏，它能够确定当前的设置，然后翻转该设置，我给它取名叫numtog.mac。在确定了程序的框架后，我开始查看关于*GET命令的帮助文档，看看是否能得到/PNUM,/PSYMB和显示模式的设置状态。令我失望的是，我发现尽管*GET命令可以设置显示模式，但不能得到关于/PNUM和/PSYMB的设置，没有它我就无法实现我的宏，失望之感油然而生。
非常时候会有非常手段，文档中没有说明由*GET得到的量并不表示一定不可以获得，也许通过我的努力可以得到。我决定从其它地方找，在这以前我经常用此办法找到其它信息，我想到了控制ANSYS菜单系统的文件，因为ANSYS在建立关于/PNUM和/PSYMB的设置时也需要相关的信息，也许它以未公开的方式使用。我开始查找UIFUNC1.GRN和UIFUNC2.GRN，果然在UIFUNC2.GRN文件中找到了我要的信息，在该文件中我找到了控制编号和符号显示的命令，它也是由*GET命令获得，只不过时以一种特殊的方式，我就把相关的命令拷贝和粘帖到我的宏文件中，以节省时间。
   尽管不推荐使用此类*GET命令，但有些高级的APDL编程人员经常使用。问题是未公开的*GET命令在5.6中能使用，但不能保证在5.7中也一定能使用，这只有等你试了以后才知道。
   让我们来看一个例子，它用来确定keyponits编号的设置状态，格式如下：
   *get,ar20,common,,d3com,,int,104
从这个命令格式我们可以推断一些信息：其中的“Common”是指储存在公共数据块中的数据，“D3com”是该公共数据块的名字，“Int”指数值数据以整数格式存储，“104”是我们所要数据存储的位置。所以上面这条命令是告诉ANSYS提取在公共数据块“D3com”104位置的整型数据，并把它赋给变量ar20，通过试验知道，当keypoints编号显示为关闭时ar20为0，而开时为1，即用1和0表示开和关。
   我的新toolbar已用了几星期了，很方便，它比旧的toolbar减少了25%的按钮。现在我的toolbar小多了，我还可以增加其它一些命令，如控制/EDGE的开和关等等。相信你会得到一些收益。

4、如何使用用户定义用户自定义矩阵
Matrix 27用户定义用户自定义矩阵，由单元选项控制定义质量、刚度或阻尼矩阵，你只要在同一组接点，分别定义三次MATRIX27单元（KEYOPT（2）分别为2，4，5）即可，然后在定义实常数时，分别定义三种单元对应的质量、刚度、阻尼矩阵系数。

5 ANSYS 查询函数（Inquiry Function）
在ANSYS操作过程或条件语句中，常常需要知道有关模型的许多参数值，如选择集中的单元数、节点数，最大节点号等。此时，一般可通过*GET命令来获得这些参数。现在，对于此类问题，我们有了一个更为方便的选择，那就是查询函数 — Inquiry Function。
    Inquiry Function类似于ANSYS的 *GET 命令，它访问ANSYS数据库并返回要查询的数值，方便后续使用。ANSYS每执行一次查询函数，便查询一次数据库，并用查询值替代该查询函数。
    假如你想获得当前所选择的单元数，并把它作为*DO循环的上界。传统的方法是使用*GET命令来获得所选择的单元数并把它赋给一个变量，则此变量可以作为*DO循环的上界来确定循环的次数
*get, ELMAX,elem,,count
*do, I, 1, ELMAX
…
…
*enddo
现在你可以使用查询函数来完成这件事，把查询函数直接放在*DO循环内，它就可以提供所选择的单元数
*do, I, ELMIQR(0,13)
…
…
*enddo
这里的ELMIQR并不是一个数组，而是一个查询函数，它返回的是现在所选择的单元数。括弧内的数是用来确定查询函数的返回值的。第一个数是用来标识你所想查询的特定实体（如单元、节点、线、面号等等），括弧内的第二个数是用来确定查询函数返回值的类型的（如选择状态、实体数量等）。
    同本例一样，通常查询函数有两个变量，但也有一些查询函数只有一个变量，而有的却有三个变量。
查询函数的种类和数量很多，下面是一些常用、方便而快速快捷的查询函数
1 AREA—arinqr(areaid,key)
areaid—查询的面，对于key=12,13,14可取为0；
key—标识关于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，节点数，如果已分网
=-7，单元数，如果已分网
kpinqr(kpid,key)的返回值
对于key=1
=-1,未选择
=0，未定义
=1，选择

3 LINE—lsinqr(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，选择的数目
=14，定义的最大数
=-2，超单元标记
=-3，主自由度
=-4，激活的自由度
=-5，附着的实体模型
ndinqr(node,key)的返回值
对于key=1
=-1,未选择
=0，未定义
=1，选择

5 VOLUMES—vlinqr(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,未选择
=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将重力以惯性力的方式施加，所以在输入加速度时，其方向应与实际的方向相反。
对于结构上的设备重量可以用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、如何实现壳单元的偏置
在构建诸如汽车车身、飞机引擎罩、发动机气道以及其他薄壁结构的有限元模型时通常选用壳单元，ANSYS提供的SHELL43、SHELL63、SHELL93以及5.5版开始新增加的SHELL143和SHELL181单元都能很方便地用于薄壁结构的建模，并能得到精确的结果，同时避免了采用实体单元构建此类结构所带来的复杂性和规模的庞大。像SOLID45和SOLID95之类的实体单元,其厚度由节点的位置所决定，而壳单元的厚度是由实常数决定的。由于壳单元的厚度独立于几何模型，因此可以很方便地进行参数化和优化设计。此外，建立一个用于生成壳单元的面几何模型要比建一个三维实体几何模型简单得多。

   尽管使用壳单元有许多优点，但也有一些缺点在使用时需要加以注意。首先，壳单元只适合某几种结构类型，因为它假设弯曲应力沿壳单元厚度方向线性变化；其次，为了得到理想的结果，壳单元必须建在实体结构的中心层或中面上。图1显示了如何用壳单元构建一个工字截面梁，左边是实体模型，右边是对应的壳单元模型。如果壳单元用实体模型的顶面或底面生成，那么有限元模型的刚度与实际结构的刚度将相差很大。

图1. 工字型梁的三维实体和中面壳单元模型

   与模型的复杂程度有关，有的时候很难或不可能得到中面，许多时候分析人员只能用模型的顶面或底面，希望壳单元的偏置影响较小可以忽略，这样做有的场合可以，有时则不行。而且，分析人员从来不知道这样做带来的误差有多大。
   关于这类问题的一种解决办法是用SHELL91层状壳单元替代标准壳单元。SHELL91主要用于复合材料分析，当然也可以用于传统材料的分析。采用SHELL91的好处是可以通过keyopt(11)来控制厚度朝单元的顶面或底面偏置。比如：当SHELL91的keyopt(11)=1时表示朝单元底面偏置，而当keyopt(11)=2时朝单元顶面偏置，而当keyopt(11)=0时则表示不发生偏置。

   由于SHELL91也带中间节点，因此它可以直接替换SHELL93单元。为了评价按这种方式使用SHELL91的效果，我分别用SOLID95、SHELL93和SHELL91来分析不同壁厚的弯管如图2所示，注意只改变管的厚度，而管的直径和弯肘半径保持不变。比值t/r，即壁厚/弯管半径用来表示中面偏置的程度。由SOLID95所得到的结果作为基准来归一化SHELL93和SHELL91模型得到的结果。所有情况都受两个载荷步：第一步在管的顶端受一个分布拉力，方向偏离固定端，使得在肘的内侧产生很大的拉应力，这个载荷条件用于确定ANSYS计算的应力是否正确；第二个载荷步是在同方向施加一个加速度，用它可以确定用带偏置的SHELL91单元时ANSYS能否正确考虑质量。

图2. 弯管实体模型及对应的SOLID95划分的有限元模型

由图3的结果显示可知:当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用在管的内侧和外侧时得到结果的精度显然没有用SHELL91采用合适的偏置所得到的结果好。
在第二载荷步施加加速度所的结果与载荷步一时的情况相同，由此可以得出结论:在中面得不到的情况下，用SHELL91替代SHELL93处理比较好。
由于我经常用CAD软件来得到我使用的几何模型，而很少有CAD软件能直接提取中面，因此使用带偏置的SHELL91单元是我经常采用的办法。

图3. 不同单元产生的变形误差

图4. 不同单元产生的应力误差

9、耦合及约束方程讲座一、耦合
当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。耦合自由度集包含一个主自由度和一个或多个其它自由度。
典型的耦合自由度应用包括：

模型部分包含对称；

在两重复节点间形成销钉、铰链、万向节和滑动连接；

迫使模型的一部分表现为刚体。

如何生成耦合自由度集

1. 在给定节点处生成并修改耦合自由度集
命令：CP
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs
在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。

2. 耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTF
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes

3. 除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：

o        如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>Numbering Ctrls>Merge Items）合并节点。

o        可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements >At Coincid Nd）通在重复节点对之间生成2节点单元来连接它们。

o        用CEINTF命令（GUI：Main Menu>Preprocessor> Coupling/Ceqn >Adjacent Regions）将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。

生成更多的耦合集

一旦有了一个或多个耦合集，可用这些方法生成另外的耦合集：

1. 用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。
命令：CPLGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same Nodes

2. 用下列方法生成与已有耦合集不同（均匀增加的）节点编号但有相同的自由度标记的新的耦合集：
命令：CPSGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same DOF

使用耦合注意事项

1. 每个耦合的节点都在节点坐标系下进行耦合操作。通常应当保持节点坐标系的一致性。

2. 自由度是在一个集内耦合而不是集之间的耦合。不允许一个自由度出现在多于一个耦合集中。

3. 由D或共它约束命令指定的自由度值不能包括在耦合集中。

4. 在减缩自由度分析中，如果主自由度要从耦合自由度集中选取，只有主节点的自由度才能被指定为主自由度。

5. 在结构分析中，耦合自由度以生成一刚体区域有时会引起明显的平衡破坏。不重复的或不与耦合位移方向一致的一个耦合节点集会产生外加力矩但不出现在反力中。

10、耦合及约束方程讲座二、约束方程
约束方程提供了比耦合更通用的联系自由度的方法。有如下形式：

这里U（I）是自由度，N是方程中项的编号。

如何生成约束方程

1. 直接生成约束方程

o        直接生成约束方程：
命令：CE
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Constraint Eqn
下面为一个典型的约束方程应用的例子，力矩的传递是由BEAM3单元与PLANE42单元（PLANE42单元无平面转动自由度）的连接来完成的：

图12-1建立旋转和平移自由度的关系

如果不用约束方程则节点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 Menu>Preprocessor>Coupling / Ceqn>Modify ConstrEqn
Main Menu>Preprocessor>Loads>Other>Modify ConstrEqn
Main Menu>Solution>Other>Modify ConstrEqn
如果要修改约束方程中的其它项，必须在求解前在PREP7中用使CE命令（或相应GUI途径）。

自动生成约束方程

o        生成刚性区域
CERIG命令通过写约束方程定义一个刚性区域。通过连接一主节点到许多从节点来定义刚性区。（此操作中的主自由度与减缩自由度分析的主自由度是不同的）
命令：CERIG
GUI: Main Menu>Preprocessor>Coupling / Ceqn>Rigid Region
将CERIG命令的Ldof设置为ALL（缺省），此操作将为每对二维空间的约束节点生成三个方程。这三个方程在总体笛卡尔空间确定三个刚体运动（UX、UY、ROTZ）。为在二维模型上生成一个刚性区域，必须保证X─Y平面为刚性平面，并且在每个约束节点有UX、UY和ROTZ三个自由度。类似地，此操作也可在三维空间为每对约束节点生成六个方程，在每个约束节点上必须有（UX、UY、UZ、ROTX、ROY和ROTZ）六个自由度。
输入其它标记的Ldof域将有不同的作用。如果此区域设置为UXYZ，程序在二维（X，Y）空间将写两个约束方程，而在三维空间（X、Y、Z）将写三个约束方程。这些方程将写成从节点的平移自由度和主节点的平移和转动自由度。类似地，RXYZ标记允许生成忽略从节点的平移自由度的部分方程。其它标记的Ldof将生成其它类型的约束方程。
总之，从节点只需要由Ldof标记的自由度，但主节点必须有所有的平移和转动自由度（即二维的UX、UY和ROTZ；三维的UX、UY、UZ、ROTX、ROTY、ROTZ）。对由没有转动自由度单元组成的模型，应当考虑增加一个虚拟的梁单元以在主节点上提供旋转自由度。

将疏密不同的已划分网格区域连在一起
可将一个区域（网格较密）的已选节点与另一个区域（网格较稀）的已选单元用CEINTF命令（菜单途径Main Menu>Preprocessor>Coupling / Ceqn>Adjacent Regions）连起来生成约束方程。

这项操作将不相容网格形式的区域“系”在一起。在两区域的交界处，从网格稠密的区域选择节点A，从网格粗糙区域选择单元B，用区域B单元的形函数，在相关的区域A和B界面的节点处写约束方程。ANSYS允许这些节点位置使用两公差准则。节点在单元之外超过第一公差就认为节点不在界面上。节点贴近单元表面的距离小于第二公差则将节点移到表面上，见下图。

o
对CEINTF命令有些限制：应力或热通量可能会不连续地穿过界面。界面区域的节点不能指定位移。可用每节点有六个自由度的单元接合6自由度实体。

o        从已有约束方程集生成约束方程集
可用CESGEN命令从已有约束方程集生成约束方程。那么已有约束方程集内的节点编号将增加以生成另外的约束方程集。另外约束方程集的标记和系数保持与原集的一致。
命令：CESGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn >Gen w/same DOF

使用约束方程的注意事项

·         所有的约束方程都以小转动理论为基础。因此，它应用在大转动分析中〔NLGEOM〕应当限制在约束方程所包含的自由度方向无重大变化的情况。

·         约束方程的出现将产生不可预料的反力和节点力结果。

由于相邻区域网格疏密不同，边界上的相容性仍然存在。但是当网格越密，这种不相容的危害就越小。

11、膜元Shell41是否能作大变形分析？
膜元Shell41可以用于大变形分析。但是膜元shell41在处理大变形问题时必须采用三角形单元，因为四边形单元会产生跷曲（warping），所以在划分网格时请选Tri。

12、如何得到径向和周向的计算结果?
在圆周对称结构中，如圆环结构承受圆周均布压力。要得到周向及径向位移，可在后处理/POST1中，通过菜单General Postproc>Options for Outp>Rsys>Global cylindric 或命令Rsys,1 将结果坐标系转为极坐标，则X方向位移即为径向位移，Y向位移即为周向位移。

13、如何加快计算速度
在大规模结构计算中，计算速度是一个非常重要的问题。下面就如何提高计算速度作一些建议：

1. 充分利用ANSYS MAP分网和SWEEP分网技术，尽可能获得六面体网格，这一方面减小解题规模，另一方面提高计算精度。

2. 在生成四面体网格时，用四面体单元而不要用退化的四面体单元。比如95号单元有20节点，可以退化为10节点四面体单元，而92号单元为10节点单元，在此情况下用92号单元将优于95号单元。

选择正确的求解器。对大规模问题，建议采用PCG法。此法比波前法计算速度要快10倍以上（前提是您的计算机内存较大）。对于工程问题，可将ANSYS缺省的求解精度从1E-8改为1E-4或1E-5即可。

14、ANSYS的几种动画模式
Q: ANSYS中有好几种动画模式，但并非同时可用。那么有哪几种类型，何时可用呢？ sDW"j\

A: ANSYS中有四种动画模式：Bitmap，AVI，Display List，Pixmap。该种模式是否可用取决于运行环境是UINX还是PC,以及是使用二维还是三维显示设备。可以通过选择菜单Utility Menu>PlotCtrls>Device Options(或键入适当的命令)来制作动画。下面是该命令选项的小结：

二维显示设备的PC： Bitmap（命令－/device,anim,bmp）

AVI（命令－/device,anim,avi）

三维显示设备的PC： Bitmap（命令－/dv3d,anim,1）

AVI（命令－/dv3d,anim,2）

Display List（命令－/dv3d,anim,0）

二维显示设备的UNIX：Pixmap（缺省，无相应命令）

三维显示设备的UNIX：Pixmap（命令－/dv3d,anim,1）

Display List（命令－/dv3d,anim,0）

Q: 这些动画模式的区别是什么？ m9\@kA

A: Display List选项（三维设备可用）与其他方式的不同在于：

Display List允许在播放动画的过程中对模型进行动态操作（如放大缩小等），而其它模式不能；

播放器不同：Bitmap，Display List，Pixmap动画模式在ANSYS图形窗口进行播放，可以通过ANSYS动画控制器面板控制动画；AVI（只在PC上）用WINDOWS媒体播放器播放。

需注意在ANSYS5.5中缺省的动画模式为AVI，播放器为媒体播放器。而在ANSYS5.6中二维设备的缺省动画模式为Bitmap，三维设备的缺省动画模式为Display List，动画在ANSYS图形窗口播放。



Q: 如何存储在ANSYS中创建的动画？

A: 对于PC环境，Bitmap或AVI模式的动画会自动存储为jobname.avi（jobname是你指定的分析名称）。Display List动画在PC上不会自动存储，必须通过菜单Utility Menu> PlotCtrls >Animate>Save Animation指定存储，典型格式为jobname.anim。

对于UNIX环境，Pixmap及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上播放.anim或.avi文件。还可以将.anim文件转变为.avi文件。所有的ANSYS产品都提供 ANIMATE程序，该程序位于ANSYS安装盘的根目录中。

Q: ANIMATE是否能读入所有的动画文件？

A: 不能。ANIMATE程序接受UNIX上的Pixmap模式的.anim文件或PC上的Bitmap或AVI模式的.avi文件。