耦合及约束方程讲座二、约束方程.docxVIP

耦合及约束方程讲座二、约束方程 约束方程提供了比耦合更通用的联系自由度的方法。有如下形式: CoMStaiU 二刀(C阻儁尿加(7) P(7)) 7-1 这里U (I )是自由度，N是方程中项的编号。深一 如何生成约束方程 1.直接生成约束方程 o 直接生成约束方程: 命令：CE GUI: Main Men uPreprocessorCoupli ng / CeqnCon stra int Eqn BEAM3单元与 BEAM3单元与PLANE42单元 (PLANE42单元无平面转动自由度)的连接来完成的： 如果不用约束方程则节点 2处表现为一个铰链。下述方法可在梁和平面应力单元之间传 递力矩，自由度之间满足下面的约束方程： ROTZ2 = (UY3 - UY1)/10 0 = UY3 - UY1 - 10*ROTZ2 相应的ANSYS命令为： CE,1,O,3,UY,1,1,UY,-1,2,ROTZ,-1O o 修改约束方程 在PREP7或 SOLUTION中修改约束方程中的常数项： 命令：CECMOD GUI: Main Me nuPreprocessorCoupli ng / Ceq nModify Con strEq n Mai n Men uPreprocessorLoadsOtherModify Con strE qn Main Men uSolutio n OtherModify Con strEq n 如果要修改约束方程中的其它项，必须在求解前在 PREP7中用使CE命令(或相应 GUI 途径)。 2.自动生成约束方程 o 生成刚性区域 CERIG命令通过写约束方程定义一个刚性区域。 通过连接一主节点到许多从节点来定义 刚性区。（此操作中的主自由度与减缩自由度分析的主自由度是不同的） 命令：CERIG GUI: Main Men uPreprocessorCoupli ng / Ceqn Rigid Regi on 将CERIG命令的Ldof设置为ALL （缺省），此操作将为每对二维空间的约束节点生成 三个方程。这三个方程在总体笛卡尔空间确定三个刚体运动（ UX、UY ROTZ。为在二 维模型上生成一个刚性区域，必须保证 X-Y平面为刚性平面，并且在每个约束节点有 UX、UY和ROTZ三个自由度。类似地，此操作也可在三维空间为每对约束节点生成六个 方程，在每个约束节点上必须有（ UX、UY UZ、ROTX ROY和ROTZ六个自由度。 输入其它标记的 Ldof域将有不同的作用。如果此区域设置为 UXYZ,程序在二维（X，Y） 空间将写两个约束方程，而在三维空间（ X、丫、Z）将写三个约束方程。这些方程将写 成从节点的平移自由度和主节点的平移和转动自由度。类似地， RXYZ标记允许生成忽 略从节点的平移自由度的部分方程。其它标记的 Ldof将生成其它类型的约束方程。 总之，从节点只需要由 Ldof标记的自由度，但主节点必须有所有的平移和转动自由度 （即二维的 UX、UY和ROTZ 三维的 UX、UY、UZ、ROTX ROTY ROTZ。对由没有转动 自由度单元组成的模型，应当考虑增加一个虚拟的梁单元以在主节点上提供旋转自由 度。 o 将疏密不同的已划分网格区域连在一起 可将一个区域（网格较密）的已选节点与另一个区域 （网格较稀）的已选单元用 CEINTF 命令（菜单途径 Main Men uPreprocessorCoupli ng / Ceq nAdjace nt Regio ns ）连起 来生成约束方程。 Fine Mesh Coarse \j1esh (selecl nodes) (select elerre^ts) 域选择节点A，从网格粗糙区域选择单元 B，用区域B单元的形函数，在相关的区域 A 和B界面的节点处写约束方程。 ANSYS允许这些节点位置使用两公差准则。节点在单元 之外超过第一公差就认为节点不在界面上。 节点贴近单元表面的距离小于第二公差则将 节点移到表面上，见下图。 对CEINTF 对CEINTF命令有些限制：应力或热通量可能会不连续地穿过界面。界面区域的节点不 能指定位移。可用每节点有六个自由度的单元接合 6自由度实体。 o 从已有约束方程集生成约束方程集 可用CESGEN命令从已有约束方程集生成约束方程。那么已有约束方程集内的节点编号 将增加以生成另外的约束方程集。另外约束方程集的标记和系数保持与原集的一致。 命令：CESGEN GUI: Main Men uPreprocessorCoupli ng / Ceqn Gen w/same DOF 使用约束方程的注意事项 ? 所有的约束方程都以小转动理论为基础。因此，它应用在大转动分析中〔 NLGEOM应当限制在 约束方程所包含的自由度方