耦合及约束方程.ppt

耦合及约束方程耦合典型的耦合自由度的应用定义耦合自由度cpngen,nset,lab,node1,node2,ninc!修改已有的耦合自由度二、耦合重合节点约束方程的应用建立刚性区。在某些特殊情况下，全刚性区给出了约束方程的另一种应用，全刚性区和部分刚性区的约束方程都可由程序自动生成。
过盈装配。与接触耦合相似，但在两个界面之间过盈量或穿透。定义约束方程cerig,maste,slave,ldof,ldof2,ldof3,ldof4,ldof5

将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,toler,dof1,dof2,dof3,dof4,dof5,dof6,moveto这项操作将不相容风格形式的区域“系”在一起。在两区域的交界处，从网格稠密的区域选择节点A，从网格粗糙的区域选择节点B，用区域B单元的形函数，在相关的区域A和B界面的节点处写约束方程。ANSYS允许这些节点位置使用两公差准则。节点在单元之外超过第一公差就认为节点不在界面上。节点贴近单元表面的距离小于第二公差则将节点移动到表面上，如图所示。对ceintf命令有些限制：应力和热通量可能会不连续地穿过界面。界面区域的节点不能指定位移。可用每节点有六个自由度的单元接合六自由度实体。4.从已有约束方程集生成约束方程集注意事项：
所有的约束方程都以小转动理论为基础。因此，它应用在大转动分析中（nlgeom）应当限制在约束方程所包含的自由度方向无重大变化的情况。
约束方程的出现将产生不可预料的反力和节点力结果。
由于相邻区域网格疏密不同，边界上的相容性仍然存在；但是当网格越密，这种不相容的危害就越小，如图所示。应用实例问题描述
wpoffs,,5
wprota,0,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
mshkey,1
vmesh,all
allsel
/solu
asel,s,loc,x,0
da,all,all

allsel
fk,101,fy,-3.0

!定义耦合及约束方程

cp,1,ux,1,21
cp,2,uy,1,21
cp,3,uz,1,21

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))
allsel
solve

/post1
plnsol,s,eqv,0,1
plnsol,u,y,0,1

etable,zl1,smisc,2
etable,zl2,smisc,8

etable,mz1,smisc,6
etable,mz2,smisc,12
plls,zl1,zl2,