Ansys自由度耦合_不同单元之间的连接
ANSYS程序变截面梁与耦合自由度的应用
一、前言变截面梁，即两端截面不同的梁。
在工程中通常将梁受力吃紧的一端采用较大的截面，另一端采用较小的截面，以实现等强度设计，节省原材料。
这种变截面梁生活中普遍存在，远至古代娶亲用的花轿轿杆、农村大车车辕，直至当今建筑钢结构中各种采光大蓬横梁（包括弯曲横梁）、各种连杆摇臂结构等。
特别是悬臂梁几乎处处可见，就连输电用的水泥电杆也做成根粗尖细的。
由于这种变截面梁在工程中普遍存在，在ANSYS程序结构分析中，专门设置了变截面梁单元的功能，使用起来特别方便。
具体操作过程在下面的实例分析中详细介绍。
自由度耦合即构件连接处两个节点的自由度（包括移动自由度和转动自由度）变化是一致的，主节点如何变化，从节点随着同样变化。
自由度耦合在静力分析时常用在连接件上，特别用在具有转动的连接件上。
例如：汽车挂挡手柄连接端的球铰，各种销钉与耳环的连接，各种转盘与转轴的连接等。
自由度耦合的概念与自由度释放的概念正好相反。
在ANSYS程序中没有自由度释放功能，只有自由度耦合功能，但用自由度耦合功能完全可以达到自由度释放的目的。
这种功能具体操作过程，也在下面的实例分析中详细介绍。
二、雨蓬计算分析该雨蓬结构是由11根变截面工字钢梁、1根等截面工字钢梁、3根圆管钢梁、1根槽钢梁和6根拉杆钢组成。
其结构如下：图1雨蓬结构1．原始数据（1）坐标数据关键点号XYZ100020063003045000 （第一拉杆上端点）40-10000 （变截面梁Z轴方向）5010000 （第一排工字钢Z轴方向）60-10006300（第六排槽钢Z轴方向）（2）材料数据主钢梁变截面工字钢：大端300×150×8×6小端150×150×8×6第一横梁工字钢：150×100×8×6第二、第三、第四横梁圆管钢：Φ102×4 Ri=47 R0=51第五横梁槽钢：160×80×8×6拉杆：Φ102×4A=1231.5 =0.0001(3) 载荷由玻璃均布载荷计算而得:中间节点: F=5560N边节点: F=2780N角节点: F=1690N(4) 单元梁单元: BEAM188杆单元: LINK8(5) 边界条件主梁固定于墙上:位移和转角全约束；拉杆与主梁连接: 位移耦合,转角自由；拉杆与墙连接: 位移约束, 转角自由。
ansys中耦合的作用以及详细解释(经典)
节点耦合可以模拟螺栓连接，他表示两者并不是一体，但某一方向的运动是一致的当需要迫使两个或多个自由度（DOFs）取行相同（介未知）值，可以将这些自由度耦合在一起。
耦合自由度集包含一个主自由度和一个或多个其他自由度。
耦合只将主自由度保存在分析的矩阵方程里，而将耦合集内的其他自由度删除。
计算的主自由度值将分配到耦合集内的所有其他自由度中支。
耦合的用途主要有以下几种：●在两重复节点间形成万向节、铰链、销钉以及滑动连接。
其原理是仅仅耦合三个平动自由度（ux,uy,uz）为铰接，耦合三个平动（ux,uy,uz）和两个转动（如RotX,RotY）则等于释放一个转动自由度为销接，其他情形如此类推。
●耦合自由度用于施加循环对称约束条件或重复循环对称约束条件，保证截面始终保持原始形状。
例如在循环对称模型中，将圆盘扇区模型的两个对称边界上的对应节点，在各个自由度上耦合；在锯齿形模型的半齿形模型（重复循环对称）中，需要将一侧边上所有节点的每个自由度进行耦合处理。
●实现小位移条件下的无摩擦接触面模型，仅仅耦合接触面在垂直于接触变面方向上的节点自由度，切线方向自由度不耦合●如果将模型中局部区域内的一部分节点都耦合起来，等于在该局部区域形成一个局部刚体（类似于约束方程中的刚性区）。
Couple DOFs：耦合节点间的单个自由度，用该命令生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。
也可用选择逻辑来耦合所选节点的全部耦合。
可用CP命令输入负的节点号来删除耦合集中的节点。
要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令（不能由GUI直接得到CPNBGEN命令）。
操作方法是：选择该子菜单，弹出拾取节点对话框，用鼠标选取参与耦合的节点，至少两个节点以上，单击OK按钮弹出图所示定义耦合对话框，在Set reference number 项输入唯一的没有占用的新耦合序列号，在Degree-of-freedom label下拉列表中选中某个自由度或者ALL（表示所有自由度），然后单击OK按钮执行耦合操作。
ANSYS中耦合自由度的方法
标题: ANSYS中耦合自由度的方法- dongyijun123 2010-03-02 19:36 阅读:3- 评论:0当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。
耦合自由度集包含一个主自由度和一个或多个其它自由度。
典型的耦合自由度应用包括：模型部分包含对称；在两重复节点间形成销钉、铰链、万向节和滑动连接；迫使模型的一部分表现为刚体。
生成耦合自由度集步骤：1.在给定节点处生成并修改耦合自由度集命令：CPGUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。
也可用选择逻辑来耦合所选节点的相应自由度。
用CP命令输入负的节点号来删除耦合集中的节点。
要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。
（不能由GUI直接得到CPNBGEN命令）。
2.耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。
此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTFGUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>Numbering Ctrls>Merge Items）合并节点。
可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements >At Coincid Nd）通在重复节点对之间生成2节点单元来连接它们。
耦合和约束方程ansys教程课件
16.1.8 练习：耦合循环对称边界
❖ 在此练习中，由生成耦合DOF 设置来模拟有循环对称性的模 型的接触问题
1.建模并在图形窗口中画单元 2.在总体柱坐标系下，生成具有
Y的增量为30的节点复制件 a.将当前坐标系变为总体柱坐
标系 b. 在当前坐标系中，以Y=30
的增量拷贝所有的结点
例如, cp,,ux,all 是把所有选择节点在UX方 向上耦合
输入耦合设置参考号 ，选择自由度卷标
耦合和约束方程ansys教程课件
16.1.6 创建耦合设置(续)
❖②在零偏移量的一组节点之间生成附加耦合关系: ❖Main Menu: Preprocessor > Coupling / Ceqn > Gen w/Same Nodes
耦合和约束方程ansys教程课件
16.2.1 约束方程的特点
约束方程的特点 ❖ 自由度卷标的任意组合 ❖ 任意节点号 ❖ 任意实际的自由度方向――在不同的节点上
ux可能不同
耦合和约束方程ansys教程课件
16.2.2 一般应用
❖ 连接不同的网格 • 实体与实体的界面 • 2－D或3－D • 相同或相似的单元类型 • 单元面在同一表面上，但结点位置不重合
如：用耦合施加循环对称 性，在循环对称切面上的 对应位置实施自由度耦合
耦合和约束方程ansys教程课件
16.1.3 施加对称条件
❖ 耦合自由度常被用来施加移动或循环对称性条件。 这可以保证平面截面依然是平面。例如：
❖ -对圆盘扇区模型 (循环对称)，应使两个对称边界上 的对应节点在各个自由度上耦合。
性区的其它节点上 ❖ 使用CERIG 命令(或 Preprocessor > Coupling/Ceqn
ANSYS各类型单元连接专题讲解(一)之连接总则
ANSYS各类型单元连接专题讲解（一）之连接总则一直以来，有不少同学咨询水哥关于ANSYS中杆单元、梁单元、壳单元、实体单元的连接问题。
之所以要用到各单元的连接，主要是由于我们在实际项目中，常常需要各种单元组合模拟，例如框架结构计算中的框架柱、框架梁采用梁单元模拟，楼板采用壳单元模拟，如此便会产生各类型单元之间的连接问题。
为解决部分朋友们的疑问，水哥依自己的理解将从以下几个方面系统讲解下ANSYS中杆单元、梁单元、壳单元、实体单元的连接，其中若有不合理之处，还望各位朋友批评指正。
本系列讲解目录如下：1、单元连接总原则。
2、杆与梁、壳、体单元的连接。
3、梁单元与实体单元铰接。
4、2D梁单元与2D实体单元刚接。
5、3D梁单元与3D实体单元刚接。
6、壳单元与实体单元连接。
7、单元连接综合实例。
本篇推文为该系列文章的首篇，主要说下ANSYS中单元连接总的原则以及简单介绍两个概念。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约束方程。
例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
从上述也可见，ANSYS无非是通过三种方法来实现单元之间的连接：共用节点、耦合、约束方程。
这里简单介绍下耦合与约束方程的基本概念。
一、耦合所谓耦合，其实是一种比较特殊的约束方程，只不过为了区别于普通一般的约束方程，方便用户操作，特定提出来的一个概念。
他具体指当我们需要迫使两个或多个自由度取得相同值（值未知）时，可以将这类自由耦合在一起。
总结一下ANSYS中不同单元之间的连接问题
论坛里常有人问不同单元之间的连接问题，我自己也一直被这个问题所困绕,最近从ANＳＹＳ工程分析进阶实例上知道了AＮSYS中不同单元之间的连接原则。
感觉收收获不小，现把它上传与大家共享。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可,不需要约束方程,否则需要耦合自由度与约事方程。
例如：(1)杆与梁、壳、体单元有公共节点即可，不需要约束方程。
ﻫ（2)梁与壳有公共节点怒可,也不需要约束写约束方程;壳梁自由度数目相同,自由度也相同，尽管壳的ｒoｔz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
(3)梁与体则要在相同位置建立不同的节点,然后在节点处耦合自由度与施加约束方程。
(４)壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的,有一定的局限性,只适用于小位移，下面介绍一种支持大位移算法的方法，ＭＰC法。
MPC即Mｕｌtiｐoint Conｓtraｉｎt,多点约束方程,其原理与前面所说的方程的技术几乎一致,将不连续、自由度不协调的单元网格连接起来,不需要连接边界上的节点完全一一对应。
MPC能够连接的模型一般有以下几种。
ﻫsolid 模型-ｓoｌid模型ﻫsｈelｌ模型-ｓheｌｌ模型ｓolｉｄ模型-shell 模型ｓolｉd 模型-beaｍ模型shelｌ模型-bｅam模型ﻫ在ＡNSＹS中,实现上述MPC技术有三种途径。
ﻫ(1）通过MPC18４单元定义模型的刚性或者二力杆连接关系。
定义MＰC1８4单元模型与定义杆的操作完全一致，而MPC单(2）利用约束元的作用可以是刚性杆（三个自由度的连接关系)或者刚性梁(六个自由度的连接关系)。
ﻫ方程菜单路径Ｍain Meｎｕ>pｒeproceｓsor>Coupliｎｇ/Cｅqn>shelｌ／ｓｏｌid Ｉｎｔerface创建壳与实体模型之间的装配关系。
ansys之节点耦合专题
/viewthread.php?tid=916386&highlight=%F1%EE%BA%CF /viewthread.php?tid=914150&highlight=%F1%EE%BA%CF/viewthread.php?tid=786833&highlight=%CD%E4%BE%D8（这个很有用！！）前言做一个总结意义！！！原因：最近网上有较多的朋友在咨询关于实体加载的方法目的：希望这个问题不再成为大家的疑惑的一部分一、说说施加方法思路1：矩或扭矩说白了就是矩，所谓矩就是力和力臂的乘积。
施加矩可以等效为施加力；思路2：直接施加弯矩或扭矩，此时需要引入一个具有旋转自由度的节点；二、在ANSYS中实现的方法这里说说3个基本方法，当然可以使用这3个方法的组合方法，组合方法就是对3个基本方法的延伸，但原理仍不变。
方法1：引入mass21，利用cerig命令Ex1:/prep7block,0,1,0,1,0,2k,9,0.5,0.5,2.5mp,ex,1,2e10mp,prxy,1,0.2mp,prxy,1,0.3r,2,1e-6et,1,45et,2,21keyopt,2,3,0lesize,all,0.2vmesh,allksel,s,,,9type,2real,2kmesh,allallselnsel,s,loc,z,2,3NPLOTCERIG,node(0.5,0.5,2.5),ALL,ALL, , , ,allsel/SOLUf,node(0.5,0.5,2.5),my,100e3FINISH/SOLnsel,s,loc,z,0d,all,allallselsolve方法2：利用mpc184单元/prep7block,0,1,0,1,0,2mp,ex,1,2e10mp,prxy,1,0.2mp,prxy,1,0.3et,1,45et,2,184keyopt,2,1,1lesize,all,0.2vmesh,alln,1000,0.5,0.5,2.5type,2mat,2*do,i,1,36e,1000,36+i*enddoallselallsel/SOLUf,node(0.5,0.5,2.5),my,100e3 FINISH/SOLnsel,s,loc,z,0d,all,allallselsolve方法3：使用rbe3命令/prep7block,0,1,0,1,0,2k,9,0.5,0.5,2.5mp,ex,1,2e10mp,prxy,1,0.2mp,prxy,1,0.3r,2,1e-6et,1,45et,2,21keyopt,2,3,0lesize,all,0.2vmesh,allksel,s,,,9type,2real,2kmesh,allallsel*dim,sla,array,36*do,i,1,36sla(i)=i+36*enddo*dim,sla2,array,36*do,i,1,36sla2(i)=i+36*enddoallselrbe3,node(0.5,0.5,2.5),all,sla,sla2allsel/SOLUf,node(0.5,0.5,2.5),my,100e3FINISH/SOLnsel,s,loc,z,0d,all,allallselsolve三、使用结论方法1和方法2的结果一致，方法3偏大。
ANSYS-coupple
1 耦合当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。
耦合自由度集包含一个主自由度和一个或多个其它自由度。
典型的耦合自由度应用包括：•模型部分包含对称；•在两重复节点间形成销钉、铰链、万向节和滑动连接；•迫使模型的一部分表现为刚体。
如何生成耦合自由度集1.在给定节点处生成并修改耦合自由度集命令：CPGUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。
也可用选择逻辑来耦合所选节点的相应自由度。
用CP命令输入负的节点号来删除耦合集中的节点。
要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。
（不能由GUI直接得到CPNBGEN命令）。
2.耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。
此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTFGUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：o如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：MainMenu>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）将两个有不相似网格模式的区域连接起来。
Ansys经典界面耦合场分析单元转换
维度：3D
耦合场：热-应力；热-电-磁；热-应力-电-磁
节点数：8
自由度：温度；结构位移；电势；磁矢势
单元几何形状
PLANE13
维：2D
耦合场：热-应力；热-电
节点数：4
自由度：温度；结构位移；磁矢势
单元几何形状
PLANE67
维度：2D
耦合场：热-电
节点数：4
自由度：温度；电势
单元几何形状
单元几何形状
PLANE223
维度：2D
耦合场：热-应力；热-电；热-应力-电；热-压电
节点数：8
自由度：结构位移；温度；电势
单元几何形状
SOLID226
维度：3D
耦合场：热-应力；热-电；热-应力-电；热-压电
节点数：20
自由度：结构位移；温度；电势
单元几何形状
SOLID227
维度：3D
耦合场：热-应力；热-电；热-应力-电；热-压电
④4节点四边形
⑤8节点四边形
⑥圆锥体
⑦圆柱体
⑧球体
CONTA171
维度：2D
耦合场：热-应力；热-电；热-应力-电
节点数：2
自由度：结构位移；温度；电势
单元几何形状
CONTA172
维度：2D
耦合场：热-应力；热-电；热-应力-电
节点数：3
自由度：结构位移；温度；电势
单元几何形状
CONTA173
维度：3D
LINK68
维度：3D
耦合场：热-电
节点数：2
自由度：温度；电势
单元几何形状
SOLID69
维度：3D
耦合场：热-电
节点数：8
ANSYS 杆单元、壳单元的单元耦合问题学习
ANSYS中杆单元、壳单元的单元耦合问题关键词：单元耦合ANSYS在比较复杂的结构的有限元分析中，不同的结构部件通常使用不同类型的单元来模拟。
通常情况下，不同类型的单元的各个节点的自由度数目是不同的，不同类型单元的连接节点处的自由度的耦合问题，是一个比较令人头疼的问题。
在ANSYS中通常可以用耦合命令CP来耦合不同类型单元在连接节点处的自由度(DOF)。
也可以用CE命令来认为添加自由度之间的约束方程来达到耦合的目的。
下面是一个简单的算例，使用了CE命令来耦合连接节点处的自由度。
模型是航天器的机翼的一个Section的某一个隔框。
上下表皮是薄壳结构，用Shell63单元来模拟，在上下表皮之间有起支撑作用的杆件，用link8单元来模拟。
建模的时候，link8单元和shell63单元在连接有各自独立的节点。
即：link8单元和shell63单元的节点在连接处是重合的，但是，节点编号是各自独立的。
link8单元在每个节点有ux,uy,uz3个平动自由度；shell63在每个节点有ux,uy,uz这3个平动自由度和rotx,roty,rotz这3个转个自由，共6个自由度。
在耦合节点处，两个耦合节点的ux,uy,uz自由度应该是相等的。
这个等式可以用CE命令来描述。
完整的命令流如下：finish/clear,start/prep7!定义第一种材料属性；mp,ex,1,30e6mp,prxy,1,0.3!定义shell63单元和实常数；et,1,shell63r,1,1e-3!建立几何模型；rectng,31.8,33.2,0,0.3556agen,2,1,1,1,0,0,1a,1,4,8,5a,6,7,3,2KL,7,0.5, ,KL,3,0.5, ,在关键点处生成节点；nkpt,100,4 !与编号为117的节点耦合nkpt,101,9 !与编号为169的节点耦合nkpt,102,10 !与编号为120的节点耦合nkpt,103,7 !与编号为160的节点耦合mat,1type,1real,1lesize,1,,,6lesize,3,,,6lesize,5,,,6lesize,7,,,6lesize,9,,,6lesize,10,,,6lesize,11,,,6lesize,12,,,6lesize,2,,,6lesize,4,,,6lesize,6,,,6lesize,8,,,6MSHAPE, 0, 2DMSHKEY, 1allselamesh,all!定义第二种材料属性；mp,ex,2,30e4mp,prxy,2,0.3!定义link8单元和实常数；et,2,link8r,2,28.26e-6mat,2type,2real,2e,101,102e,100,101e,102,103!CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, !Lab3, C3ce,1,0,100,ux,1,117,ux,-1 !节点100的ux＝节点117的ux;ce,2,0,100,uy,1,117,uy,-1 !节点100的uy＝节点117的uy;ce,3,0,100,uz,1,117,uz,-1 !节点100的uz＝节点117的uz;ce,4,0,101,ux,1,169,ux,-1 !同上；ce,5,0,101,uy,1,169,uy,-1ce,6,0,101,uz,1,169,uz,-1ce,7,0,102,ux,1,120,ux,-1ce,8,0,102,uy,1,120,uy,-1ce,9,0,102,uz,1,120,uz,-1ce,10,0,103,ux,1,160,ux,-1ce,11,0,103,uy,1,160,uy,-1ce,12,0,103,uz,1,160,uz,-1!施加约束，底面约束所有的自由度；DA,1,all!施加分布载荷；SFA,2,1,PRES,-1e4finish/solusolve !求解；finish/post1PLNSOL, S,X, 0,1.0finishPLNSOL, S,X, 0,1.0 对应的结果云图PLNSOL, S,EQV, 0,1.0对应的结果云图。
ansys在处理体壳连接问题中的应用
3.1 模型一（小规模连接模型）
模型一：如图 2 所示的悬臂梁，由一个短厚的块体（20m×10m×7m ）和一块非常长 的薄板（100m×10m×1m ）在块体的中部连接而成，块体和板的弹性模量都是 3×1011pa, 泊松比为 0.3，在板的自由端施加垂直于板且大小为 600N/m 的均布荷载，求最大的位移与 应力。
全部壳单元
8.0171
0.00
352800
2.00
解析解
8.0170
360000
从表 1 可以看出，对于小规模的，简单的体壳组合结构模型，三种方法都能够处理实体
单元和壳单元的连接问题，且处理结果与解析解的结果偏差较小, 一般处于误差范围之内。
3.2 模型二（大规模连接模型）
当对大型的体壳组合结构模型进行有限元分析时，如苏通大桥索塔锚固区的钢锚箱与混 凝土塔壁组合结构，由于其结构的特殊性，在进行网格划分时，需要采用自由划分，此时在 接触面上实体单元和壳单元的节点就可能不重合，且在接触面上生成的节点数量很大，属于 大规模的接触问题，如图 3 所示。
其次定义接触单元和目标单元以组成“接触对”。“接触对”的接触单元必须建立在壳单 元侧，目标单元建立在实体单元侧，接触单元和目标单元利用 ESURF 命令创建，它们通过 共享实常数设置组成“接触对”。
-2-
2.3 用复合单元处理接触面
在实际生活中, 我们经常遇到图 1 所示结构的连接情况，即实体与板采用搭接的方式连 接。
ansys中如何进行刚性连接
ansys中如何进行刚性连接一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。
例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
举例：有一长为100mm的矩形截面梁，截面为10X1mm，与一规格为20mmX7mmX10mm的实体连接，约束实体的端面，在梁端施加大小为3N的y方向的压力，梁与实体都为一材料，弹性模量为30Gpa，泊松比为0.3。
本例主要讲解梁与实体连接处如何利用耦合及约束方程进行处理。
命令流如下：FINI/CLE/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名/PREP7ET,1,SOLID95 !定义实体单元类型为SOLID95ET,2,BEAM4 !定义梁单元类型为BEAM4MP,EX,1,3E4 !定义材料的弹性模量MP,PRXY,1,0.3 !定义泊松比R,1 !定义实体单元实常数R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数BLC4,,,20,7,10 !创建矩形块为实体模型WPOFFS,0,3.5 !将工作平面向Y方向移动3.5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPOFFS,0,5 !将工作平面向Y方向移动5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致K,100,20,3.5,5 !创建关键点K,101,120,3.5,5 !创建关键点L,100,101 !连接关键点生成梁的线实体LSEL,S,LOC,X,21,130 !选择梁线LATT,1,2,2 !指定梁的单元属性LESIZE,ALL,,,10 !指定梁上的单元份数LMESH,ALL !划分梁单元VSEL,ALL !选择所有实体VATT,1,1,1 !设置实体的单元属性ESIZE,1 !指定实体单元尺寸MSHAPE,0,2D !设置实体单元为2DMSHKEY,1 !设置为映射网格划分方法VMESH,ALL !划分实体单元ALLS !全选FINI !退出前处理!------------------------/SOLU !进入求解器ASEL,S,LOC,X,0 !选择实体的端面DA,ALL,ALL !约束实体端面ALLS !全选FK,101,FY,-3.0 !在两端施加Y向压力CP,1,UX,1,21 !耦合节点1和节点21X方向自由度CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) !设置约束方程CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283)) !设置约束方程CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283)) !设置约束方程ALLS !全选SOLVE !保存FINI !退出求解器!------------------------/POST1 !进入通用后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, S,EQV, 0,1.0 !显示等效应力ETABLE,ZL1,SMISC,1 !读取梁单元上I节点X方向的力ETABLE,ZL2,SMISC,7 !读取梁单元上J节点X方向的力ETABLE,MZ1,SMISC,6 !读取梁单元上I节点Z方向的力矩ETABLE,MZ2,SMISC,12 !读取梁单元上J节点Z方向的力矩PLETAB,ZL1 !显示梁单元X方向的力PLETAB,MZ1 !显示梁单元Z方向力矩!**********************************************上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的，有一定的局限性，只适用于小位移，下面介绍一种支持大位移算法的方法，MPC 法。
ansys中nummrg与glue命令及自由度耦合的比较ZZ
ansys中nummrg与glue命令及自由度耦合的比较ZZCP, nset, lab,node1,node2,……node17定义或改变耦合节点自由度PREP7: Coupled DOFnset：耦合组编号设置如下：n：随机设置数量HIGH：使用最高定义的耦合数量（如果Lab=all，此为默认值）。
该选项用于在已有组中增加节点。
NEXT：将定义的最高耦合数量增加1。
该项用于在现有组未改变时自动定义耦合组。
Lab: 耦合节点的自由度。
定义类型随NSET所选类型改变：结构类：UX, UY, or UZ (位移); ROTX, ROTY, or ROTZ (角度)；热分析类： TEMP, TBOT, TE2, TE3, . . ., TTOP (温度)；流体分析类: PRES (压力); VX, VY, or VZ (速率)；电子类: VOLT (电压); EMF (电场耦合值); CURR (电流).磁分析类: MAG (标量磁位差); AX, AY, or AZ (矢量磁位差); CURR (电流).Explicit analysis labels: UX, UY, or UZ (位移)。
node1~node17: 待耦合的节点号。
输入相同的节点号会被忽略。
如果某一节点号为负，则此节点从该耦合组中删去。
如果node1=all,则所有选中节点加入该耦合组。
注意：1，不同自由度类型将生成不同编号2，不可将同一自由度用于多套耦合组耦合自由度的结果是耦合组中的一个元素与另一个元素有相同的属性。
耦合可以用于模型不同的结点和联结效果。
一般定义耦合可以使用约束公式（CE）。
对结构分析而言，耦合节点由节点方向定义。
耦合的结果是，这些节点在指定的结点坐标方向上有相同的位移。
对于一组没有定义位移的耦合节点，可能会产生应力弯矩，这些弯矩不是由作用力产生的。
对特定节点的实际自由度是由元素类型（ET）所指定的。
例如，BEAM3的自由度是UX，UY和ROTZ。
(整理)ANSYS单元节点自由度耦合.
ANSYS自由度耦合当生成模型时，典型地是用单元去连接节点以建立不同自由度间的关系，但是，有时需要能够刻划特殊细节（刚性区域结构的铰链连接，对称滑动边界，周期条件，和其他特殊内节点连接等），这些用单元不足以来表达，可用耦合和约束方程来建立节点自由度间的特殊联系，利用这些技术能进行单元做不到的自由度连接。
1、什么是耦合当需要迫使两个或多个自由度（DOFs）取得相同（但未知）值，可以将这些自由度耦合在一起，耦合自由度集包含一个主自由度和一个或多个其他自由度。
耦合只能将主自由度保存在分析的矩阵方程里，而将耦合集内的其他自由度删除。
计算的主自由度值将分配到耦合集内的所有其他自由度中去。
典型的耦合自由度应包括：部分模型包含对称；在两个重复节点间形成销钉，铰链，万向节和滑动连接；迫使模型的一部分表现为刚体。
2、如何生成耦合命令：CPGUI：Preprocessor——Coupl/Ceqn——Couple DOF在生成一个耦合节点之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中。
也可用选择逻辑来耦合所选节点的全部耦合。
可用CP命令输入负的节点号来删除耦合集合中的节点。
要修改一耦合自由度（即增、删节点或改变自由度标记）用CPNGEN命令（无GUI）。
CPINTF命令通过在对每对重合节点上定义自由度标记生成一耦合集而实现对模型重合节点的耦合。
此操作对“扣紧”几对节点（诸如一条缝）尤为有用。
命令:CPINTFGUI：Preprocessor——Coupl/Ceqn——Coincident nodes除耦合重复节点外，还可用下列替换方法迫使节点有相同的变现方式：（1）如果对重复节点所有自由度都要耦合，通常用NRMMRG （numbering——mergeit）将这些节点合并起来更方便；（2）可用EINTF命令（create——element——at coincident）在重复节点生成2节点单元连接；（3）用EINTF(preprocessor——couple/ceqn——adjacent rejoins)将两个不相似网格模式的区域连接起来，这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程；（4）用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。
[转载]ansys中实体单元与壳单元的连接处理方法
![[转载]ansys中实体单元与壳单元的连接处理方法](https://img.taocdn.com/s3/m/61e6e4024531b90d6c85ec3a87c24028905f8541.png)
[转载]ansys中实体单元与壳单元的连接处理⽅法原⽂地址：ansys中实体单元与壳单元的连接处理⽅法作者：埃及⽂字2010为简化模型，在有些模型中采⽤壳和实体混合的单元类型，由于壳单元节点表⽰的是⼀个截⾯，因此除平动⾃由度外，⼀般具有转动⾃由度；⽽实体单元节点表⽰的是⼀个点，因此只有平动⾃由度。
因此在这些模型中，壳单元与实体单元过渡处的节点由于⾃由度的不同，不能单纯的进⾏节点耦合。
⽬前实体单元与壳单元连接常⽤的处理使⽤MPC⽅法(SHSD命令)定义两者之间的装配关系，这种⽅法⼀般可理解为绑定接触形式。
尽管采⽤了接触模块，采⽤MPC⽅法在⼩变形时不需要平衡迭代，在⼤变形中在每个平衡迭代中不断进⾏更新，⼜克服了传统约束⽅程只适⽤于⼩应变的限制。
该⽅法是处理节点⾃由度耦合较为理想的形式。
本⽂分别采⽤混合单元与纯实体单元两种⽅式进⾏建模，分析，结果及命令流如下：采⽤混合单元进⾏分析时的应⼒分布云图（第三应⼒强度理论）：单纯采⽤实体单元进⾏分析时应⼒分布云图（第三应⼒强度理论）：采⽤混合单元分析的命令流如下：/CLEAR/FILNAME, EXAMPLE26/PREP7ET, 1, SOLID95ET, 2, SHELL63ET, 3, TARGE170KEYOPT, 3, 5, 1ET, 4, CONTA175KEYOPT, 4, 2, 2KEYOPT, 4, 12, 5R, 1, 0.02R, 2R, 3R, 4R, 5MP, EX, 1, 2E11MP, PRXY, 1, 0.3/VIEW, 1, 1, 1, 1BLOCK, -0.14, 0.14, -0.14, 0.14, 0,0.98VDELE, 1,,,0ADELE, 1, 2, 1, 1BLOCK, -0.15, 0.15, -0.15, 0.15, 0.98, 1K, 20, 0, 0, 0.98K, 21, 0, 0.1, 0.98K, 22, 0, 0.1, 1K, 23, 0, 0.075, 1K, 24, 0, 0.075, 1.3 K, 25, 0, 0, 1.3 LSEL, NONEL, 20, 21L, 21, 22L, 22, 23L, 23, 24L, 24, 25L, 25, 20LFILLT, 27, 28, 0.025 AL, ALL VROTAT, 11,,,,,,20, 25 ALLS VOVLAP, ALL AATT, 1, 1, 2 ESIZE, 0.02 MSHAPE, 0 MSHKEY, 1 AMESH, 3, 6, 1 VATT, 1, 1, 1 ESIZE, 0.0175 SMRTSIZE, 5 MSHAPE, 1 MSHKEY, 0 VMESH, ALLALLSASEL, S,,,47 NSLA, S, 1TYPE, 3REAL, 2ESURFALLSLSEL, S,,,5ESURF ALLS ASEL, S,,,47 NSLA, S, 1 TYPE, 3 REAL, 3 ESURF ALLS LSEL, S,,,6 NSLL, S, 1 TYPE, 4 REAL,3 ESURF ALLS ASEL, S,,,47 NSLA, S, 1 TYPE, 3 REAL, 4 ESURF ALLS LSEL, S,,,7 NSLL, S, 1 TYPE, 4 REAL, 4 ESURF ALLS ASEL, S,,,47 NSLA, S, 1 TYPE, 3 REAL, 5 ESURF ALLSREAL, 5ESURFALLSSHSD, 2, CREATE SHSD, 3, CREATESHSD, 4, CREATESHSD, 5, CREATEFINISH/SOLUASEL, S,,,21, 27, 6ASEL,A,,,15,33,18NSLA, S, 1D, ALL, ALLALLSNSEL, S, LOC, Z, 0D, ALL, UZALLSSFA, 4,2, PRES, -1E6*3/2.804 SFA, 8, 1, PRES, -1E6*3/2.804 SOLVEFINISH/POST1PLNSOL, S, EQV, 0, 1 FINISH。
梁，壳单元与实体单元耦合自由度
梁，壳单元与实体单元耦合自由度将梁或者壳向体单元深入一层.如图所示. 但是不要指望连接处的应力会精确图2. 稳定连接一般用MPC是最好的。
深入两层，只相差1%。
From ANSYS Knowledge BaseQ: Do you have any examples showing how to connect a 6 DOF shell to a 3 DOF solid by constraint equations and by embedding the shell 2 elements deep into the solid?A: The following input files show that there is only a 1% difference in V on Mises stress between using constraint equations and embedding the shell 2 elements deep into the solid.!Example of attaching a 6 DOF shell to a 3 DOF solid!Done by embedding the shell 2 layers deep into the solid FINI/CLEAR/PREP7MP,EX,1,30E6MP,NUXY,1,0.3ET,1,SOLID185BLC4,0,0,10,10,-10BLC4,0,0,10,10,10ESIZE,,10VMESH,ALL!NUMMRG,NODEASEL,U,,,ALLNSEL,U,,,ALLET,2,SHELL181R,2,0.5BLC4,0,8,10,2BLC4,0,10,10,10/VIEW,1,1,1,1/PNUM,TYPE,1TYPE,2REAL,2LESIZE,ALL,1AMESH,ALLALLSNUMMRG,NODEEPLONSEL,S,LOC,Y,0D,ALL,ALLALLSNSEL,S,LOC,Y,20F,ALL,FZ,100ALLS/SOLUNLGEOM,OFFEQSLV,SPARSESOLVE!Unselect embedded shell elements - with a nodal solution!these elements will reduce true stress level in shells!due to averagingASEL,S,,,13ESLA,UPLNSOL,S,EQV***************************************************Method #2!Example of attaching a 6 DOF shell to a 3 DOF solid!Done by using Constraint EquationsFINI/CLEAR/PREP7MP,EX,1,30E6MP,NUXY,1,0.3ET,1,SOLID185BLC4,0,0,10,10,-10BLC4,0,0,10,10,10ESIZE,,10VMESH,ALLNUMMRG,NODEASEL,U,,,ALLNSEL,U,,,ALLET,2,SHELL181R,2,0.5BLC4,0,10,10,10/VIEW,1,1,1,1TYPE,2REAL,2AMESH,ALLALLSEPLONSEL,S,LOC,Y,0D,ALL,ALLALLSNSEL,S,LOC,Y,20F,ALL,FZ,100ALLS!**********Apply Constraint Equations**********! Get Master nodes*DO,I,0,10,1ESEL,S,TYPE,,2 !Select Shell ElementsNSLE !Select Nodes attached to elementsNSEL,R,LOC,Y,10 !Select Shell Element nodes at interfaceNSEL,R,LOC,X,I !Select Master Node*GET,MN%I%,NODE,0,NUM,MAX !Create parameter for node number ALLS*ENDDO! Get Slave nodes*DO,J,0,10,1ESEL,S,TYPE,,1 !Select Solid ElementsNSLE !Select Nodes attached to elementsNSEL,R,LOC,Y,10 !Select Solid Element nodes at top surfaceNSEL,R,LOC,Z,-1,1 !Select Solid Element nodes at interfaceNSEL,R,LOC,X,J !Select Slave NodeCM,TEMPNODE,NODENSEL,A,NODE,,MN%J%NLIST,ALLCERIG,MN%J%,ALL,UXYZALLS*ENDDO/SOLUNLGEOM,OFF !Large Deformations not allowed with CEEQSLV,SPARSE !shown to ensure PCG isn't chosen by userSOLVE/POST1PLNSOL,S,EQV。
第16章 耦合和约束方程(ansys教程)
16.2.5 建立刚性区
约束方程通常被用来模拟刚性区 作用在节点(主节点 主节点)上的载荷将被恰当地分配到刚 作用在节点 主节点 上的载荷将被恰当地分配到刚 性区的其它节点上 使用CERIG 命令 或 Preprocessor > Coupling/Ceqn 命令(或 使用 > Rigid Region) 在某些特殊情况下，全 在某些特殊情况下， 刚性区给出了约束方程 的另一种应用 全刚性区和部分刚性区 的约束方程都可由程序 自动生成
16.1.1 耦合设置的特点
只有一个自由度卷标－ 只有一个自由度卷标－如：ux,uy或temp 或 可含有任意节点数 任意实际的自由度方向－ 在不同的节点上 任意实际的自由度方向－ux在不同的节点上 可能是不同的 主、从自由度的概念 加在主自由度上的载荷
16.1.2 一般应用
施加对称条件 无摩擦界面 铰接 如 ： 用耦合施加循环对称 性 ， 在循环对称切面上的 对应位置实施自由度耦合
16.1.5 铰接
耦合可用来模拟铰接， 耦合可用来模拟铰接，如：万向节、铰链 万向节、 借助力矩释放可模拟铰接： 借助力矩释放可模拟铰接： 只耦合连接节点间的位 移自由度， 移自由度，不耦合旋转自由度 例如： 下图中， 若 A处重合两节点在 例如 ： 下图中 ， 处重合两节点在UX、 UY方向 、 方向 处重合两节点在 上耦合，旋转不耦合， 上耦合，旋转不耦合，则A连接可模拟成铰接 连接可模拟成铰接 节点1和节点 节点 和节点2 和节点 重合， 重合 ， 为了看 清分开显示
16.2.4 连接不同类型的单元
如果需要连接自由度集不同的单元类型， 如果需要连接自由度集不同的单元类型，则要求写出 约束方程以便于从一类单元向另一类单元传递载荷: 约束方程以便于从一类单元向另一类单元传递载荷 梁与实体或垂直于壳的梁 壳与实体 命令： 命令： CE 命令 (Preprocessor > Coupling/Ceqn > Constraint Eqn) 建立转动自由 度和移动自由 度之间的关系
ansys ls-dyna分析全过程(跌落)
R14w笔记本跌落仿真分析一笔记本仿真模型的建立1.1模型的建立1.1.1 模型简化对分析结果无关紧要的一些细节部分常常使模型相当复杂，在实体建模时往往可将这一步略去。
在某些情况下，由于一些很小的局部而破坏了整个结构的对称性.有时可略去这些局部(或将它们作对称处理)以保持对称结构，缩小分析的规模。
必须权衡简化模型的利弊(损失精度以减小花费)审慎地略去不对称部分。
在这里还想进一步说明模型简化中一些方法和技巧。
“子结构“是将一组有限元压缩成为一个用用一个矩阵表示的超单元.采用子结构的原因有：1.减少计算时间.在非线性分析中，可用子结构计算结构的线性部分，以便那部分的单元短阵不必重复计算每一个平衡迭代；对于有重复部分的结构分忻,可以生成一个超单元来表示这部分图形、然后拷贝到不同的位置：2.利用有限的计算机资料解决非常大型的问题。
当一个分析相对于计算机波前空间或磁盘空间来说太大了、用子结构可使每一部分都满足计算机的要求：“子模型”是为了获得模型中某一区域的更精确的解而产生的一种有限元技术。
当整个模型的网格划分相对于某一区域太组时，可不必重新对整个模型进行更纫的划分，只需对这一区域重新划分。
这就大大节约了时间和费用。
“等效结构”的概念为：将原来的复杂结构用一简单结构模拟，新结构的材料和几何特性与原结构有所不同但刚度等效。
其等效结构是指那些具有重复性的均匀结构，如蜂窝结构、晶体结构等。
1.1.2 单元类型选择1.2单元类型的选择Ansys 隐式单元ANSYS的单元库提供了100 多种的单元类型,单元类型选择的工作就是将单元的选择范围缩小到少数几个单元上,通常采用以下方法。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。
例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
ANSYS网格划分的一些例子
虽然做出来了.但是我还是有一个问题想请教大家:vsweep和mapp分网后形成网格各有什么规律?如何结合两种方法划分出整齐规则的网格呢.比如:为什么图中的(1)部分用MAPP划分,(2)部分用SWEEP划分呢就可以出现上图中的那种整齐规则的网格?反过来(1)部分用SWEEP,(2)部分用MAPP划分就不会出现整齐规则的网格呢?部分(1)和部分(2)不可看成一个整体划分吗?我试了一个,如果把两个部分看成整体,可以分网但是不会出现那种整齐的网格.只有掌握了生成网格规律才容易得到合理,整齐,规则的网格,总不能分网时把各种方法都试一遍吧. 恳请各位谈点自己的在分网方面的经验.谢谢1的三个边如果都设了分段数则sweep和map是一样的et,1,42et,2,45cyl4,,,20lsel,alllesize,all,,,10esize,,10vext,1,,,,,20aclear,allamesh,1不过好象中间不大好的!!!还望高手指点! 命令流;et,1,42et,2,45blc4,,,10,5lesize,1,,,5lesize,2,,,10mshape,0,2dmshkey,1amesh,1esize,,5vrotat,1,,,,,,1,4aclear,all用map也可以,,取其四分之一,单元大小可控制!做了一个!/PREP7CYL4, , ,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,_Z2REAL,_Z4ESYS,0ET,2,SOLID45 TYPE, 2 EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,0 EXTOPT,ATTR,0,0,0 MAT,_Z2REAL,_Z4ESYS,0 VOFFST,2,5, , VOFFST,6,5, , VOFFST,7,5, ,VOFFST,5,5, ,VOFFST,3,5, ,总算搞出来了:如何把网格画成这样？请教高手，如何把界面处的网格画成这样？ansys可以画成那样的，要多切割点volume，那个图看上去就是分成几个体单独划分网格，然后再merge起来的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
问题：如下图所示block单元和beam单元如何连接在一起？先看例子：FINI/CLE/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名/PREP7ET,1,SOLID95 !定义实体单元类型为SOLID95ET,2,BEAM4 !定义梁单元类型为BEAM4MP,EX,1,3E4 !定义材料的弹性模量MP,PRXY,1,0.3 !定义泊松比R,1 !定义实体单元实常数R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数BLC4,,,20,7,10 !创建矩形块为实体模型WPOFFS,0,3.5 !将工作平面向Y方向移动3.5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPOFFS,0,5 !将工作平面向Y方向移动5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致K,100,20,3.5,5 !创建关键点K,101,120,3.5,5 !创建关键点L,100,101 !连接关键点生成梁的线实体LSEL,S,LOC,X,21,130 !选择梁线LATT,1,2,2 !指定梁的单元属性LESIZE,ALL,,,10 !指定梁上的单元份数LMESH,ALL !划分梁单元VSEL,ALL !选择所有实体VATT,1,1,1 !设置实体的单元属性ESIZE,1 !指定实体单元尺寸MSHAPE,0,2D !设置实体单元为2DMSHKEY,1 !设置为映射网格划分方法VMESH,ALL !划分实体单元ALLS !全选FINI !退出前处理/SOLU !进入求解器ASEL,S,LOC,X,0 !选择实体的端面DA,ALL,ALL !约束实体端面ALLS !全选FK,101,FY,-3.0 !在两端施加Y向压力CP,1,UX,1,21 !耦合节点1和节点21X方向自由度CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) !设置约束方程CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283)) !设置约束方程CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283)) !设置约束方程ALLS !全选SOLVE !保存FINI !退出求解器/POST1 !进入通用后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, S,EQV, 0,1.0 !显示等效应力ETABLE,ZL1,SMISC,1 !读取梁单元上I节点X方向的力ETABLE,ZL2,SMISC,7 !读取梁单元上J节点X方向的力ETABLE,MZ1,SMISC,6 !读取梁单元上I节点Z方向的力矩ETABLE,MZ2,SMISC,12 !读取梁单元上J节点Z方向的力矩PLETAB,ZL1 !显示梁单元X方向的力PLETAB,MZ1 !显示梁单元Z方向力矩ansys中联系自由度的方法(CE)——耦合和约束方程CE命令1 耦合当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。
耦合自由度集包含一个主自由度和一个或多个其它自由度。
典型的耦合自由度应用包括："模型部分包含对称；"在两重复节点间形成销钉、铰链、万向节和滑动连接；"迫使模型的一部分表现为刚体。
如何生成耦合自由度集1.在给定节点处生成并修改耦合自由度集命令：CPGUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。
也可用选择逻辑来耦合所选节点的相应自由度。
用CP命令输入负的节点号来删除耦合集中的节点。
要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。
（不能由GUI直接得到CPNBGEN命令）。
2.耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。
此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTFGUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：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.用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。
命令：CPLGENGUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same Nodes2.用下列方法生成与已有耦合集不同（均匀增加的）节点编号但有相同的自由度标记的新的耦合集：命令：CPSGENGUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same DOF使用耦合注意事项1.每个耦合的节点都在节点坐标系下进行耦合操作。
通常应当保持节点坐标系的一致性。
2.自由度是在一个集内耦合而不是集之间的耦合。
不允许一个自由度出现在多于一个耦合集中。
3.由D或共它约束命令指定的自由度值不能包括在耦合集中。
4.在减缩自由度分析中，如果主自由度要从耦合自由度集中选取，只有主节点的自由度才能被指定为主自由度。
5.在结构分析中，耦合自由度以生成一刚体区域有时会引起明显的平衡破坏。
不重复的或不与耦合位移方向一致的一个耦合节点集会产生外加力矩但不出现在反力中。
2约束方程约束方程提供了比耦合更通用的联系自由度的方法。
有如下形式：这里U（I）是自由度，N是方程中项的编号。
如何生成约束方程1.直接生成约束方程直接生成约束方程：命令：CEGUI: Main Menu>Preprocessor>Coupling / Ceqn>Constraint Eqn下面为一个典型的约束方程应用的例子，力矩的传递是由BEAM3单元与PLANE42单元（PLANE42单元无平面转动自由度）的连接来完成的：图12-1建立旋转和平移自由度的关系如果不用约束方程则节点2处表现为一个铰链。
下述方法可在梁和平面应力单元之间传递力矩，自由度之间满足下面的约束方程：ROTZ2 = (UY3 - UY1)/100 = UY3 - UY1 - 10*ROTZ2相应的ANSYS命令为：CE,1,0,3,UY,1,1,UY,-1,2,ROTZ,-10修改约束方程在PREP7或SOLUTION中修改约束方程中的常数项：命令：CECMODGUI: Main Menu>Preprocessor>Coupling / Ceqn>Modify ConstrEqnMain Menu>Preprocessor>Loads>Other>Modify ConstrEqnMain Menu>Solution>Other>Modify ConstrEqn如果要修改约束方程中的其它项，必须在求解前在PREP7中用使CE命令（或相应GUI途径）。
2.自动生成约束方程生成刚性区域CERIG命令通过写约束方程定义一个刚性区域。
通过连接一主节点到许多从节点来定义刚性区。
（此操作中的主自由度与减缩自由度分析的主自由度是不同的）命令：CERIGGUI: 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界面的节点处写约束方程。