耦合和约束方程
ANSYS中不同单元之间的连接问题
ANSYS 中不同单元之间的连接问题中不同单元之间的连接问题论坛里常有人问不同单元之间的连接问题，我自己也一直被这个问题所困绕，最近从ANSYS工程分析进阶实例上知道了ANSYS 中不同单元之间的连接原则。
中不同单元之间的连接原则。
感觉收收获不小，现把它感觉收收获不小，现把它上传与大家共享。
上传与大家共享。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，全包含前一种单元的自由度，则只要有公共节点即可，则只要有公共节点即可，则只要有公共节点即可，不需要约束方程，不需要约束方程，否则需要耦合自由度与约事方程。
例如： （1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz 是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
一般来说，beam 与beam 的连接方式默认为刚接，link 与link 的默认链接方式为铰接，link与beam 的默认连接方式为铰接。
那么实际情况中经常遇到beam 与beam 的铰接,如何处理理 举例：举例：有一长为100mm 的矩形截面梁，截面为10X1mm ，与一规格为20mmX7mmX10mm 的实体连接，约束实体的端面，在梁端施加大小为3N 的y 方向的压力，梁与实体都为一材料，弹性模量为30Gpa ，泊松比为0.3。
本例主要讲解梁与实体连接处如何利用耦合及约束方程进行处理。
命令流如下： FINI /CLE /FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名定义工作名/PREP7 ET,1,SOLID95 !定义实体单元类型为SOLID95 ET,2,BEAM4 !定义梁单元类型为BEAM4MP ,EX,1,3E4 !定义材料的弹性模量定义材料的弹性模量 ,PRXY MP ,PRXY,1,0.3 ,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.5 WPROTA,0,90 !将工作平面绕X 轴旋转90度 VSBW,ALL !将实体沿工作平面剖开将实体沿工作平面剖开WPOFFS,0,5 !将工作平面向Y 方向移动5 WPROTA,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 !选择梁线选择梁线 LA TT,1,2,2 !指定梁的单元属性指定梁的单元属性LESIZE,ALL,,,10 !指定梁上的单元份数上的单元份数 LMESH,ALL !划分梁单元划分梁单元 VSEL,ALL !选择所有实体选择所有实体V A TT,1,1,1 !设置实体的单元属性ESIZE,1 !指定实体单元尺寸 MSHAPE,0,2D !设置实体单元为2D MSHKEY ,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 方向自由度方向自由度,2,UY CP ,2,UY,1,21 !,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 PLNSOL, U,Y, 0,1.0 !, 0,1.0 !显示Y 方向位移向位移 PLNSOL, S,EQV PLNSOL, S,EQV, 0,1.0 !, 0,1.0 !显示等效应力显示等效应力ETABLE,ZL1,SMISC,1 !读取梁单元上I 节点X 方向的力 ETABLE,ZL2,SMISC,7 !读取梁单元上J 节点X 方向的力ETABLE,MZ1,SMISC,6 !读取梁单元上I 节点Z 方向的力矩方向的力矩 带格式的: 突出显示批注突出显示批注[微软系统1]: 1是转动的中心点！！！ 带格式的: 突出显示突出显示ETABLE,MZ2,SMISC,12 !读取梁单元上J 节点Z 方向的力矩方向的力矩PLETAB,ZL1 !显示梁单元X 方向的力方向的力 PLETAB,MZ1 !显示梁单元Z 方向力矩方向力矩上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的，要是用耦合自由度和约束方程来实现的，有一定的局限性，有一定的局限性，只适用于小位移，只适用于小位移，下面介绍一种下面介绍一种支持大位移算法的方法，MPC 法。
ANSYS中CERIG命令详解
CERIG命令详解2011-12-21 21:27:42| 分类：ANSYS | 标签：|字号大中小订阅ansys学习－耦合与约束方程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节点单元来连接它们。
ANSYSLSDYNA_经验手册
ANSYS/LSDYNA 经验手册显式与隐式方法对比：隐式时间积分不考虑惯性效应（[C]and[M]）。
在t＋△t时计算位移和平均加速度：{u}={F}/[K]。
线性问题时，无条件稳定，可以用大的时间步。
非线性问题时，通过一系列线性逼近（Newton－Raphson）来求解；要求转置非线性刚度矩阵[k]；收敛时候需要小的时间步；对于高度非线性问题无法保证收敛。
显式时间积分用中心差法在时间t求加速度：{a}=([F(ext)]-[F(int)])/[M]。
速度与位移由：{v}={v0}+{a}t,{u}={u0}+{v}t新的几何构型由初始构型加上{X}={X0}+{U}非线性问题时，块质量矩阵需要简单的转置；方程非耦合，可以直接求解；无须转置刚度矩阵，所有的非线性问题（包括接触）都包含在内力矢量中；内力计算是主要的计算部分；无效收敛检查；保存稳定状态需要小的时间步。
关于文件组织：jobname.lsdyna输入流文件，包括所有的几何，载荷和材料数据jobname.rst后处理文件主要用于图形后处理（post1），它包含在相对少的时间步处的结果。
jobname.his在post26中使用显示时间历程结果，它包含模型中部分与单元集合的结果数据。
时间历程ASCII文件，包含显式分析额外信息，在求解之前需要用户指定要输出的文件，它包括：GLSTAT全局信息，MATSUM材料能量，SPCFORC节点约束反作用力，RCFORC接触面反作用力，RBDOUT刚体数据，NODOUT节点数据，ELOUT单元数据……在显式动力分析中还可以生成下列文件：D3PLOT类似ansys中jobname.rstD3THDT时间历程文件，类似ansys中jobname.his关于单元：ANSYS/LSDYNA有7种单元（所有单元均为三维单元）:LINK160:显式杆单元；BEAM161：显式梁单元；SHELL163：显式薄壳单元；SOLID164：显式块单元；COMBI165：显式弹簧与阻尼单元；MASS166：显式结构质量；LINK167:显式缆单元显式单元与隐式单元不同：每种单元可以用于几乎所有的材料模型。
LS-DYNA常见问题集锦
1 如何处理LS-DYNA中的退化单元在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。
在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。
这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。
其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。
所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。
2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。
该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。
除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。
3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。
特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。
具体操作方法和过程参见下列命令流文件的演示。
4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。
X13ANSYS-耦合和约束方程
• 在一个模型中可以定义任意多个约束方程。 在一个模型中可以定义任意多个约束方程。 • 一个约束方程可以包含任意数量的节点和自由度。约束方程的一般 一个约束方程可以包含任意数量的节点和自由度。 形式是： 形式是：
Coef1 * DOF1 + Coef2 * DOF2 + Coef3 * DOF3 + ... = Constant
– 先选择网格划分较好的部分实体的节点和 另一部分的单元。 另一部分的单元。 – 自动计算所有必要的系数和常数。 自动计算所有必要的系数和常数。 – 适合于实体单元对实体单元 2-D 或 3-D 适合于实体单元对实体单元, 。
November 3, 2003 Inventory #001970 3-15
...耦合
一般应用: 一般应用 • 施加对称条件 • 无摩擦界面 • 铰接
November 3, 2003 Inventory #001970 3-4
耦合和约束方程
...耦合
施加对称条件
• 自由度耦合，常被用来施加平移或者循环对称条件，以保证截面依然是平面。 自由度耦合，常被用来施加平移或者循环对称条件，以保证截面依然是平面。
Objective
一个耦合设置是一组被约束在一起，有着相同大小， 一个耦合设置是一组被约束在一起，有着相同大小，但值未知 的自由度
Definition
耦合设置的特点： 耦合设置的特点： • • • • • 只有一个自由度卷标－如：ux,uy或temp 只有一个自由度卷标－ 或 可含有任意节点数 任意实际的自由度方向－ 在不同的节点上可能是不同的 任意实际的自由度方向－ux在不同的节点上可能是不同的 主、从自由度的概念 加在主自由度上的载荷
第 13 章 耦合和约束方程
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 法。
约束和耦合方程
5. 单击OK
4. 在约束方程中将要使用的自由度
练习－在蜗轮叶片上建立约束方程
在此练习中，将使用约束方程将具有不同单元类型和不同网格的两部分连接起来。这两部分分别是涡轮叶片段及叶片连接的基座
1. 恢复数据库文件（eblade.db1）并在图形窗口中显示单元.2. 选择基座上的单元（mat2）3. 选择叶片底面上的节点 a.首先，unselect附在底座单元上的节点（接第2步）
3. 输入耦合设置参考号，选择自由度卷标.
4. 单击OK.
建立耦合关系（续）
在零偏移量的一组节点之间生成附加耦合关系:Main Menu: Preprocessor > Coupling / Ceqn > Gen w/Same Nodes
3. 单击OK
1. 输入现存耦合设置的参考号.
2. 对每个设置指定新的自由度卷标.
耦合的一般性应用（续）
2. 无摩擦的界面 如果满足下列条件，则可用耦合自由度来模拟接触面:表面保持接触, 此分析是几何线性的（小变形）忽略摩擦在两个界面上，节点是一一对应的.通过仅耦合垂直于接触面的移动来模拟接触. 优点:分析仍然是线性的无间隙收敛性问题
3. 铰接 耦合可用来模拟力耦松弛松，例如铰链、无摩擦滑动器、万向节
约束方程的应用
1. 连接不同的网格:实体与实体的界面2－D或3－D相同或相似的单元类型单元面在同一表面上，但结点 位置不重合
Objective
3-5. 3－5说明约束方程的四种应用.
约束方程的应用（续）
建立转动自由度和移动自由度之间的关系
2. 连接不相似的单元类型:壳与实体垂直于壳或实体的梁.
3. 建立刚性区在某些特殊情况下，全刚性区给出了约束方程的另一种应用全刚性区和部分刚性区的约束方程都可由程序自动生成
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
以现有的约束方程为基础生成约束方程：
1. 生成第一个约束方程: Main Menu: Preprocessor > Coupling / Ceqn > Constraint Eqn. 2. 生成其余的约束方程: Main Menu: Preprocessor > Coupling/Ceqn > Gen w/Same DOF.
中南大学
3
工程软件应用
耦合的一般应用
.
Objective
耦合的一般性应用.
1.
施加对称性条件： • 耦合自由度常被用来实施移动或循环对称条件. • 考虑在均匀轴向压力下的空心长圆柱体，此3－D结构可用下面右图 所示的2－D轴对称模型表示.
y y
11 x 1
12 2
13 3
14 4
15 5 x
由于结构的对称性，上面的一排结点在轴向 上的位移应该相同
1. 输入现存耦合 设置的参考号.
2. 对每个设置指定 新的自由度卷标 . 3. 单击OK
中南大学
8
工程软件应用
建立耦合关系
1. ..... 2. ..... 3. ..... Procedure
在同一位置的节点之间自动生成耦合关系:
Main Menu: Preprocessor > Coupling / Ceqn > Coincident Nodes 1. 指定自由度卷标. 2. 指定节点位置的 容差
人工建立约束方程的菜单路径： Main Menu: Preprocessor > Coupling / Ceqn > Constraint Eqn
1. 输入常数项， 节点号，自由 度卷标和方程 系数.
2. 单击OK.
中南大学
16
工程软件应用
建立约束方程的过程
1. ..... 2. ..... 3. ..... Procedure
Constant = Coef1 * DOF1 + Coef2 * DOF2 + ...
中南大学
11
工程软件应用
约束方程的应用
说明约束方程的应用.
Objective
1.连接不同的网格:
• 实体与实体的界面 • 2－D或3－D • 相同或相似的单元类型 • 单元面在同一表面上，但结点 位置不重合
中南大学
中南大学
4
工程软件应用
耦合的一般性应用
2. 无摩擦的界面 • 如果满足下列条件，则可用耦合自由度来模拟接触面: –表面保持接触, –此分析是几何线性的（小变形） –忽略摩擦 –在两个界面上，节点是一一对应的. •通过仅耦合垂直于接触面的移动来模拟接触. •优点: –分析仍然是线性的 –无间隙收敛性问题
12
工程软件应用
约束方程的应用
2. 连接不相似的单元类型: 壳与实体 垂直于壳或实体的梁.
• 建立转动自由度和移动自由度之间的关系
中南大学
13
工程软件应用
约束方程的应用
3.建立刚性区
• 在某些特殊情况下， 全刚性区给出了约 束方程的另一种应 用
• 全刚性区和部分刚 性区的约束方程都 可由程序自动生成
3. 单击OK
中南大学
9
工程软件应用
约束方程
2015-6-15
中南大学
10
工程软件应用
约束方程
定义“约束方程”
Objective
约束方程定义节点自由度之间的线性关系
Definition
约束方程的特点 • 自由度卷标的任意组合. • 任意节点号. • 任意实际的自由度方向――在不同的节点上ux可能不同. 例
4. 单击OK. 2. 单击OK
3. 输入耦合设置参考 号，选择自由度卷 标.
中南大学
7
工程软件应用
建立耦合关系
1. ..... 2. ..... 3. ..... Procedure
在零偏移量的一组节点之间生成附加耦合关系:
Main Menu: Preprocessor > Coupling / Ceqn > Gen w/Same Nodes
拾取将要连在一起的结点，然后单击OK
1. 选择将要使用的
刚性区的类型（ 自由度设置）
2. 单击OK
中南大学
18
工程软件应用
建立约束方程的过程
1. ..... 2. ..... 3. ..... Procedure
在相邻的区域生成约束方程：
1. 从网格较密的区域中选择节点 2. 从网格较稀的区域中选择单元. Main Menu: Preprocessor > Coupling / Ceqn > Adjacent Regions
中南大学
5
工程软件应用
耦合的一般性应用
考虑一个2D的梁模型，每个节点上有三个 自由度ux、uy和rotz，A点为一铰链连接 。将同一位置节点的自由度ux、uy耦合起 来。
3. 铰接
耦合可用来模拟力 耦松弛，例如铰链、 无摩擦滑动器、万 向节
1
2
A
节点1和节点2 处于同一位置 ，但为于清楚 起见，在图上 分开显示。.
工程软件应用
耦合和约束方程
2015-6-15
中南大学
1工程软件应用源自耦合2015-6-15
中南大学
2
工程软件应用
耦合设置
定义“耦合设置”
Objective
Definition
一个耦合设置是一组被约束在一起，有着相同大小， 但值未知的自由度 耦合设置的特点： • 只有一个自由度卷标－如：ux,uy或temp • 可含有任意节点数 • 任意实际的自由度方向－ux在不同的节点上可能是 不同的
为了模拟铰接，将同一位置两个节点 的移动自由度耦合起来，而不耦合转 动自由度
中南大学
6
工程软件应用
建立耦合关系
通过三种不同的办法建立耦合关系
Objective
1. ..... 2. ..... 3. ..... Procedure
进入人为地创建耦合关系的菜单路径: Main Menu: Preprocessor > Coupling / Ceqn > Couple DOFs 1. 拾取将要耦合的结点
3.
指定容差，此容 差作为单元区域 中最小单元长度 的比率.
4.
在约束方程中 将要使用的自 由度
5. 单击OK
中南大学
19
生成的约束方 程数. 现存约束方 程中的节点 增量 生成的约束方 程的起始序号 ，终止序号和 增量
3. 选择 OK
中南大学
17
工程软件应用
建立约束方程的过程
1. ..... 2. ..... 3. ..... Procedure
通过“刚性区”来建立约束方程
Main Menu: Preprocessor > Coupling / Ceqn > Rigid Region >
中南大学
14
工程软件应用
约束方程的应用
4.
过盈装配
• 与接触耦合相似，但在两个界面之间允许有过盈量 或穿透 • 典型方程: 0.01 = UX (node 51) - UX (node 251)
中南大学
15
工程软件应用
建立约束方程的过程
采用四种不同的方法建立约束方程.
Objective
1. ..... 2. ..... 3. ..... Procedure