ANSYS中CERIG命令详解

合集下载

ANSYS命令解释

ANSYS命令解释分享⼀下！BEGINNER'S GUIDE TO ANSYS COMMANDSThe symbol '*' corresponds to the following:* --> k, l, a, v, e, n, cm, et, mp, r where ==>k --> Keypointsl --> Linesa --> Areav --> Volumese --> Elementsn --> Nodescm --> componentet --> element typemp --> material propertyr --> real constant$ --> d, f, sf, bf, ic, where ==>d --> DOF constraint (ux... in Structural, Temp in thermal,f --> Force Load ( Heat in thermal)sf --> Surface load on nodesbf --> Body Force on NodesMore Commands can be generated by sensible combinations of " $* " family of commands. See the following list of $* possible options$* --> dk --> DOF constraints on KP (Vx,Vy,Pres... in CFD)dl --> DOF constraints on Linesda --> DOF constraints on Areasfk --> Force on Keypointssfl --> Surface load on Linessfa --> Surface load on Areassfe --> Surface load on element faces中国热模⽹⾸发bfk --> Body Force on Keypointsbfl --> Body Force on Linesbfa --> Body Force on Areabfv --> Body Force on Volumesbfe --> Body Force on Elementsic --> Initial Conditions ",p" --> If ",p" was issued at the end of the Command(in Input Window) the GUI based picking menu will be activated. Useful for listing, plotting, meshing, deleting, etc..**********************************************************1. Listing of picked Entities:COMMAND SYNTAX: *LIS,p & $*LIS,pA few Combinations of this command are:klis,p --> List KPllis,p --> Lists Linesalis,p --> Lists Areasvlis,p --> Lists Volumeselis,p --> Lists Elementsnlis,p --> Lists Nodescmlis,p --> Lists componentscslis,p --> Lists user created local co-ordinate systemsdlis,p --> Lists DOF constraints specfied on nodesdalis,p --> Lists DOF constraints applied on Areasflis,p --> Lists force on Nodes中国热模⽹⾸发sfllis,p --> Lists Surface Load on linesbfalis,p --> Lists body force load applied on Areasiclis,p --> Lists Initial condition on NodesIf ",p" was not issued, all entites currently selected will be listed.For certain commands ",p" cannot be issued. See the below mentioned commandsetlis --> Lists the different element types definedmplis --> Lists whatever Material propertiesrlis --> Lists whatever real constantscslis --> Lists all co-ordinate systemscmlis --> Lists all components*********************************************************2. Plotting of Entities: COMMAND SYNTAX: *plo KPLO / LPLO / APLO / VPLO / EPLO / NPLO / CMPLO / **********************************************************3. deleting of Entities:COMMAND SYNTAX: *DEL,p & $*DEL,pKDEL,p / LDEL,p / ADEL,p / VDEL,p / EDEL,p / NDEL,p / CMDEL,p / DDEL,p /DKDEL,p / DADEL,p / FDEL,p / SFDEL,p / SFEDEL,p / SFADEL,p / SFLDEL,p / 热点模具⽹BFADEL,p / ......The syntax for this command is very similar to LISTING command.**********************************************************4. distance between two entities:COMMAND SYNTAX: *DIS,pndis,p --> Distance between two nodeskdis,p --> Distance between two KPs**********************************************************5. Meshing of geometries:COMMAND SYNTAX: *MES,pKMES,p / LMES,p / AMES,p / VMES,p**********************************************************6. Size settings for Lines and Areas before meshing :COMMAND SYNTAX : *size,,p Lesiz,p / Aesize,p*********************************************************7. Clearing Meshes of already meshed geometries:COMMAND SYNTAX: *CLE,p KCLE,p / LCLE,p / ACLE,p / VCLE,p**********************************************************8. BOOLEAN Operations: IntersectCOMMAND SYNTAX : *IN* AINA,p / VINV,p / LINL,p / AINV,p / LINV,p / LINA,p 热点模具⽹**********************************************************9. BOOLEAN Operations: GLUECOMMAND SYNTAX : *GLUE VGLUE,p / AGLUE,p / LGLUE,p**********************************************************10. Boolean Operations: SUBTRACT/DIVIDE:COMMAND SYNTAX: *sb*,p See the following examples to understand how this works:asba,p --> Subtract Area from Areaasbl,p --> Divide Area by linevsba,p --> Divide volume by Arealsbw,p --> Divide line by Workplanevsbw,p --> Divide volume by Workplaneasbw,p --> Divide area by Workplanevsbv,p --> subtract Volume by another volumeMore combinations exist. The user needs to explore them for themselves --> forms a part of learning**********************************************************11. Boolean Operations: Overlap:COMMAND SYNTAX: *OVLAP,p AOVLAP,p / VOVLAP,p**********************************************************12. Concatenation of Lines / Areas --> for map meshing热点模具⽹COMMAND SYNTAX : *ccat,pLCCAT,p --> Concatenation of Lines for Map meshing AreaACCAT,p --> Concatenation of Areas for Map meshing Volume*********************************************************13. Dragging operationCOMMAND SYNTAX : *drag,pvdrag,p --> Drag areas along a line to create a new volumeadrag,p --> Drag line along a line to create a new arealdrag,p --> Drag KP along a line to create a new line**********************************************************14. Copy Geomtric entitiesCOMMAND SYNTAX : *GEN,,pKGEN,,p / LGEN,,p / AGEN,,p / VGEN,,pPlease note that *GEN commands are also used for MOVE operations. The difference lies in the value specified in the 10th field of these *GEN commands. By default it is 0 --> which does the COPY operation. If specfied as 1 --> it does the MOVE operation**********************************************************15. Bottom -to- Top modeling commands:COMMAND SYNTAX : *,p & **,p 热点模具⽹k,p ---> Allows user to pick KP in the Workplanel,p ---> Create lines from existing KPak,p ---> Create area from KPal,p ---> Create areafrom linesv,p ---> Create Volume from KPva,p ---> Create Volume from Arease,p ---> Create Elem from existing nodesen,p ---> Create Elem from nodes**********************************************************16. To apply common Boundary Conditions such as DOF constraint, Forces, Surface Loads, Bodyforce Loads and Initial conditions* --> is meant for the KLAVE entities only (KLAVEN stands for KP, Lines, Area, Volumes & ELem )16a. DOF constraint :COMMAND SYNTAX : $*,p ( Please Note: NOT all * are valid)See the valid combinations below:D,p --> To apply DOF on nodesDK,p --> To apply DOF on KeypointsDL,p --> Apply DOF on LinesDA,p --> Apply DOF on Areas ( symmetry or Anti-symmetry will be prompted)****************16b. FORCE Loading:COMMAND SYNTAX : $*,p 中国热模⽹⾸发See the valid combinations below:f,p --> Forces on nodesfk,p --> Force on Keypoints(fa,p or FV,p or FL,p ----> Since force cannot be applied on Lines or Area & volumes... this command does not exist.)****************16c. Surface Loads:COMMAND SYNTAX : $*,pSee the valid combinations below:sf,p --> Surface Load on a set of Nodessfl,p --> Surface Load on Linessfa,p --> Surface Load on Areasfe,p --> Surface Load on Element(SFk,p and SFV,p do not exist since pressure cannot be applied on a single Kp and neither can it be applied on a volume) ****************16d. BodyForce Load: COMMAND SYNTAX : bf*,pSee the valid combinations below:bf,p --> Bodyforce Load on a set of Nodesbfk,p --> Bodyforce Load on KPbfl,p --> Bodyforce Load on Linesbfa,p --> Bodyforce Load on Areasbfv,p --> Bodyforce Load on Volumesbfe,p --> Bodyforce Load on Element****************16e. Initial conditions:热点模具⽹ic,p --> Initial Conditions on Nodes(P.S: Initial Conditions can be applied only to nodes. )***********************************************************17. To refine a mesh :COMMAND SYNTAX : *ref,pkREF,p / kREF,p / aREF,p / eREF,p / nREF,p***********************************************************18. To TRANsfer loads from the Solid model to the FE model:COMMAND SYNTAX : $TRANdtran / ftran / sftran / bftran & SBCTRAN(SBCTRAN --> Transfers all solid model loads to FE model)***********************************************************19. Writing / Reading information to a file (ASCII)COMMAND SYNTAX : *read, & *write,NWRITE / MPWRITE / ETWRITE / RWRITE / EWRITE / CDWRITENREAD / MPREAD / ETREAD / RREAD / EREAD / CDREAD / LDREAD(Some of these commands ETWRITE/ETREAD , RWRITE/RREAD are undocumented. But they do work) The Commands CDWRITE and CDREAD are used to write/read all FE model related info (w or w/o geometry to ASCII files) Its recommended the user read the online help on these two commands before using them 热点模具⽹The LDREAD commands are used to read loads (LD) from other analysis types. For example: Tempfrom Thermal results file (*.rth) are applied onto Structural elements.好，我来补充⼀下楼上师兄的命令。

ANSYS各类型单元连接专题讲解(一)之连接总则

ANSYS各类型单元连接专题讲解（一）之连接总则一直以来，有不少同学咨询水哥关于ANSYS中杆单元、梁单元、壳单元、实体单元的连接问题。

之所以要用到各单元的连接，主要是由于我们在实际项目中，常常需要各种单元组合模拟，例如框架结构计算中的框架柱、框架梁采用梁单元模拟，楼板采用壳单元模拟，如此便会产生各类型单元之间的连接问题。

为解决部分朋友们的疑问，水哥依自己的理解将从以下几个方面系统讲解下ANSYS中杆单元、梁单元、壳单元、实体单元的连接，其中若有不合理之处，还望各位朋友批评指正。

本系列讲解目录如下：1、单元连接总原则。

2、杆与梁、壳、体单元的连接。

3、梁单元与实体单元铰接。

4、2D梁单元与2D实体单元刚接。

5、3D梁单元与3D实体单元刚接。

6、壳单元与实体单元连接。

7、单元连接综合实例。

本篇推文为该系列文章的首篇，主要说下ANSYS中单元连接总的原则以及简单介绍两个概念。

一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约束方程。

例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。

（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。

（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。

（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。

从上述也可见，ANSYS无非是通过三种方法来实现单元之间的连接：共用节点、耦合、约束方程。

这里简单介绍下耦合与约束方程的基本概念。

一、耦合所谓耦合，其实是一种比较特殊的约束方程，只不过为了区别于普通一般的约束方程，方便用户操作，特定提出来的一个概念。

他具体指当我们需要迫使两个或多个自由度取得相同值（值未知）时，可以将这类自由耦合在一起。

ANSYS中加载力矩

现在比较常用的加矩方法一般有以下3种：1．将矩转换成一对的力偶,直接施加在对应的节点上面.2．在构件中心部位建立一个节点,定义为MASS21单元,然后跟其他受力节点藕荷,形成刚性区域,就是用CERIG命令.然后直接加转矩到主节点,即中心节点上面3 使用MPC184单元.是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,,从而形成刚性面.最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。

第一种方法一般可以加扭矩，但是在节点上很容易出现应力集中现象。

后两种情况都是形成刚性区域,但是CERIG命令是要在小变形或者小旋转才能用,只支持静力,线形分析；而第三种方法适用多种情况,不仅支持大应变,还支持非线形情况。

（以上内容来自网络）第一种方法就不介绍了，网上介绍很多，看看就会了。

后2种方法说的较多，但没有具体的例子和验证。

现在做的一个项目上要加弯距和扭矩，我是采用的第二种方法，同时自己对改方法进行过验证，把过程写下来，大家共同探讨。

模型简单取1个空心梁，采用实体建模，单元用95，内半径10mm，外半径15mm，长度200mm，矩10E6 Nmm，mass21单元实常数取零。

(1)梁建模划网(2)在梁轴线上生成1节点，偏离端面10mm。

(3)选定mass21单元属性，通过第二步的节点直接生成单元(E或EN)。

(4)指定刚性化区域（CERIG），先选取mass21对应的节点，再选取梁端面节点。

(5)加约束和载荷，梁另一端全约束，在mass21对应节点上加扭矩。

(6)求解后验证结果。

(7)在mass21上加弯距，求解验证结果。

加扭矩，按公式计算得到最大剪应力235.179Mpa。

应力分布连续，无应力集中现象。

加弯距，按公式计算轴向应力为470.357Mpa，但在应力云图上梁端部明显出现端部效应，取中间截面轴向应力对比，结果基本吻合。

所以加弯距时考虑结构加长以减小端部效应的影响。

第三种方法有空在试试。

ANSYS命令流总结(全)

ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。

s就是select，选择。

DIM是定义数组的意思。

array 数组。

MP命令用来定义材料参数。

K是建立关键点命令。

K,关键点编号,x坐标,y坐标，z坐标。

K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。

NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。

这个命令对于重复的线面都可以用。

这个很简单，压缩关键。

Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。

LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。

SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108，超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81，SOLID98，FLUID129，INFIN110，111，FLUID116，130界面单元:INTER192，193，194，195显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性E-弹性(Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent)等。

ansys加转矩

ANSYS中如何加弯矩或扭矩A转矩一般有三种施加的方法:第一种,将矩转换成一对一对的力偶,直接施加在对应的节点上面.第二种,在构件中心部位建立一个节点,定义为MASS21单元,然后跟其他受力节点藕荷,形成刚性区域,就是用CERIG命令.然后直接加转矩到主节点,即中心节点上面第三种,使用MPC184单元.是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,,从而形成刚性面.最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷.上面三种方法计算的结果基本一致,我做过实验的.只不过是后两种情况都是形成刚性区域,但是CERIG命令是要在小变形或者小旋转才能用,只支持静力,线形分析.而第三种方法适用多种情况,不仅支持大应变,还支持非线形情况.如果你需要例子,我下次在发给你看A思路1：矩或扭矩说白了就是矩，所谓矩就是力和力臂的乘积。

施加矩可以等效为施加力；思路2：直接施加弯矩或扭矩，此时需要引入一个具有旋转自由度的节点；可以选择单元21，或者184Q：模态分析得到的结果是不是某个方向上的各阶频率啊，我要得到各个方向的一阶频率能做到吗？A：模态分析得到的结果是你所选择的自由度内的振型。

如果需要得到某一个方向内的振型，用RUDECED法选择你所希望的主自由度如UX，将所求的频率值设置为一即可。

Q：不过大型复杂结构的振型一般都不是一个方向的应该是一个方向为主其它方向为次，如果把所要求的方向设为主自由度会不会跟真实的实际结果有出入？Q：对于复杂的大型实体其主自由度无法给出，此时就无法应用reduce法，可不可以改约束还是就认为一次得到前3阶频率就是x,y,z方向的一阶频率？A：正是因为复杂结构的主自由度不好确定，所以倾向于用SUBSPACE和BLOCK LANCZOS 法，约束不要随便改。

它应该描述结构真实的状态，不能说前3阶频率即为XYZ方向的第一阶频率。

你可以好好想想所用方程中[K]和[M]的意义。

ansys常用命令-自编

几何建模命令：/VIEW, WN, XV, YV, ZV选择一组线：命令：LSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP倒圆角命令：LFILLT, NL1, NL2, RAD, PCENTVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP创建圆弧或圆命令：CIRCLE, PCENT, RAD, PAXIS, PZERO, ARC, NSEG创建面：命令：A, P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11, P12, P13,P14, P15, P16, P17, P18 命令：AL, L1, L2, L3, L4, L5, L6, L7, L8, L9, L10（可用all）命令：ADRAG, NL1, NL2, NL3, NL4, NL5, NL6, NLP1, NLP2, NLP3, NLP4, NLP5, NLP6命令：AROTAT, NL1, NL2, NL3, NL4, NL5, NL6, PAX1, PAX2, ARC, NSEG命令：AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE创建体命令：V, P1, P2, P3, P4, P5, P6, P7, P8四边形或四方体命令：BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH四边形或四方体命令：BLC5, XCENTER, YCENTER, WIDTH, HEIGHT, DEPTH工作平面原点创建圆面或环面：PCIRC, RAD1, RAD2, THETA1, THETA2圆或圆柱命令：CYL4, XCENTER, YCENTER, RAD1, THETA1, RAD2, THETA2, DEPTH圆或圆柱命令：CYL5, XEDGE1, YEDGE1, XEDGE2, YEDGE2, DEPTH正多边形命令：RPOLY, NSIDES, LSIDE, MAJRAD, MINRAD正多边形命令：RPR4, NSIDES, XCENTER, YCENTER, RADIUS, THETA, DEPTH通过3个坐标点定义工作平面命令：WPLANE,WN,XORIG,YORIG,ZORIG,XXAX,YXAX,ZXAX,XPLAN,YPLAN,ZPLAN通过3个关键点定义工作平面命令：KWPLAN, WN, KORIG, KXAX, KPLAN将工作平面移动到一组关键点的中间位置命令：KWPAVE, P1, P2, P3, P4, P5, P6, P7, P8, P9将工作平面移动到一组指定坐标的中间位置命令：WPAVE, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3工作平面恢复到ANSYS默认状态的命令：WPSTYL,DEFA将工作平面沿其自身坐标轴移动命令：WPOFFS, XOFF, YOFF, ZOFF工作平面的旋转命令：WPROTA, THXY, THYZ, THZX工作平面设置：WPSTYL,SNAP,GRSPAC,GRMIN,GRMAX,WPTOL,WPCTYP,GRTYPE,WPVIS, SNAPANG 将既有坐标系的XY平面定义为工作平面命令：WPCSYS,WN,KCN（直接用WPCSYS即可恢复到默认位置）线线减运算：LSBL, NL1, NL2, SEPO, KEEP1, KEEP2面面减运算：ASBA, NA1, NA2, SEPO, KEEP1, KEEP2体体减运算：VSBV, NV1, NV2, SEPO, KEEP1, KEEP2创建体：V, P1, P2, P3, P4, P5, P6, P7, P8命令：VA, A1, A2, A3, A4, A5, A6, A7, A8, A9, A10拖拽面命令：VDRAG, NA1, NA2, NA3, NA4, NA5, NA6, NLP1, NLP2, NLP3,NLP4, NLP5, NLP6面旋转命令：VROTAT, NA1, NA2, NA3, NA4, NA5, NA6, PAX1,PAX2, ARC, NSEG面偏移创建体命令：VOFFST, NAREA, DIST, KINC镜像面命令：ARSYM, Ncomp, NA1, NA2, NINC, KINC, NOELEM, IMOVE通过面延伸命令：VEXT, NA1, NA2, NINC, DX, DY, DZ, RX, RY, RZ复制体命令：VGEN, ITIME, NV1, NV2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE镜像体命令：VSYMM, Ncomp, NV1, NV2, NINC, KINC, NOELEM, IMOVE正棱柱体命令：RPR4, NSIDES, XCENTER, YCENTER, RADIUS, THETA, DEPTH球体命令：SPHERE, RAD1, RAD2, THETA1, THETA2任意位置球命令：SPH4, XCENTER, YCENTER, RAD1, RAD2通过直径端点生成球体命令：SPH5, XEDGE1, YEDGE1, XEDGE2, YEDGE2以工作平面原点为圆心创建圆锥体命令：CONE, RBOT, RTOP, Z1, Z2, THETA1, THETA2任意位置圆锥体命令：CON4, XCENTER, YCENTER, RAD1, RAD2, DEPTH/trlcy,lab,value设置透明显示取中间值表示半透明。

ANSYS项目中常用的命令流

目录提取质量与体积 (2)加弯矩命令流 (2)显示未画网格的体 (2)面单元与体单元耦合 (3)ELEMENT TABLE (4)提取质量与体积TVOLU=0 ！起初单元体积设定值*GET,MaxEleNum,ELEM,,NUM,MAX ！提取对所要求的体对应的单元的最大号*do,i,1,MaxEleNum ！进行do循环*get,volu,elem,i,volu ！提取每个单元的体积TVOLU=TVOLU+volu ！对各个单元体积进行叠加得到最后总的体积*enddo ！循环结束，在parameter中得到TVOLU的值，即为总体积TolWEIGHT=TVOLU*dens ！得到这个体的质量，在parameter中得到TolWEIGHT的值，即为总质量TVOLU=0*GET,MaxEleNum,ELEM,,NUM,MAX*do,i,1,MaxEleNum*get,volu,elem,i,voluTVOLU=TVOLU+volu*enddoTolWEIGHT=TVOLU*1注释：1.在有些结构需要求解质量和体积时要用到这些，如塔器的附加质量，来求密度用到的体积。

2.对于多个体的体积分别提取时，只要把参数变化一下就可以，如把MaxEleNum变为MaxEleNum2等。

需要变化的已经标记。

加弯矩命令流F,671341,FX, ——F,671341,FY, ——F,671341,FZ, ——F,671341,MX, ——F,671341,MY, ——F,671341,MZ, ——显示未画网格的体allsel,all ！首先全选vsel,s,type,,1 ！TYPE为单元类型，这里选取单元类型为1的体组件，vsel,inve ！把刚才选取的体进行反选，反选为整体的体减去上一步的体。

Vplot ！显示体注：（1）上述APDL是在只有一种单元的类型情况，如有多种单元类型情况如下：allsel,all ！首先全选vsel,s,type,,1 ！TYPE为单元类型，这里选取单元类型为1的体组件，vsel,a,type,,2vsel,2,type,,3 ！Also select 同时选择单元类型为2的和3的。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

CERIG命令详解
2011-12-21 21:27:42| 分类：ANSYS | 标签：|字号大中小订阅
ansys学习－耦合与约束方程
1 耦合
当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。

耦合自由度集包含一个主自由度和一个或多个其它自由度。

典型的耦合自由度应用包括：
"模型部分包含对称；
"在两重复节点间形成销钉、铰链、万向节和滑动连接；
"迫使模型的一部分表现为刚体。

如何生成耦合自由度集
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.在结构分析中，耦合自由度以生成一刚体区域有时会引起明显的平衡破坏。

不重复的或不与耦合位移方向一致的一个耦合节点集会产生外加力矩但不出现在反力中。

2约束方程
约束方程提供了比耦合更通用的联系自由度的方法。

有如下形式：
这里U（I）是自由度，N是方程中项的编号。

如何生成约束方程
1.直接生成约束方程
直接生成约束方程：
命令：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
修改约束方程
在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途径）。

2.自动生成约束方程
生成刚性区域
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允许这些节点位置使用两公差准则。

节点在单元之外超过第一公差就认为节点不在界面上。

节点贴近单元表面的距离小于第二公差则将节点移到表面上，见下图。

对CEINTF命令有些限制：应力或热通量可能会不连续地穿过界面。

界面区域的节点不能指定位移。

可用每节点有六个自由度的单元接合6自由度实体。

从已有约束方程集生成约束方程集
可用CESGEN命令从已有约束方程集生成约束方程。

那么已有约束方程集内的节点编号将增加以生成另外的约束方程集。

另外约束方程集的标记和系数保持与原集的一致。

命令：CESGEN
GUI: Main Menu>Preprocessor>Coupling / Ceqn >Gen w/same DOF
使用约束方程的注意事项
"所有的约束方程都以小转动理论为基础。

因此，它应用在大转动分析中〔NLGEOM〕应当限制在约束方程所包含的自由度方向无重大变化的情况。

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

"由于相邻区域网格疏密不同，边界上的相容性仍然存在。

但是当网格越密，这种不相容的危害就越小。