ENDRELEASE命令创建铰接的耦合规则及示例

在创建铰接时，常用的方法有几下几种：
1、BEAM44单元，释放自由度： KEYOPT。
2、BEAM188／189，ENDRELEASE命令自动创建耦合。
3、手动创建耦合,最万能，但是也最复杂。

在8.0以后，新的ENDRELEASE命令给188单元带来了非常的活力！但是关于此命令网上的资料却很少。无论BAIDU 还是GOOGLE 的结果都只有一页，而且内容含糊不清。现在我们就来看一下这个命令的用法。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
一、命令格式
在PREP7下：
ENDRELEASE,--,TOLERANCE,Dof1,Dof2,Dof3,Dof4

--：程序预留位，估计以后还要继续开发
TOLERANCE：容许角度，以“度”为单位，超过此数值的单元结合处才进行计算。默认20，输入－1以释放所有单元
Dof1,Dof2,Dof3,Dof4：需要释放的自由度，如果Dof1为空，默认为翘曲（warping）并且Dof2以及以后参数忽略：
＝》WARP：释放翘曲自由度
＝》ROTX：释放X方向转角
＝》ROTY：释放Y方向转角
＝》ROTZ：释放Z方向转角
＝》UX：释放X方向位移
＝》UY：释放Y方向位移
＝》UZ：释放Z方向位移
＝》BALL：创建球铰（等效于释放WARP、ROTX、ROTY、ROTZ）
菜单路径：
Main Menu>Preprocessor>Loads>Load Step Opts>Other>End Releases>On Selected setMain Menu>Solution>Load Step Opts>Other>End Releases>On Selected set

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
二、范例
这里我们将做一个单层厂房的屋架。包括梁和檩条（由于是范例，适当简化）。从正上方看是一个4＊4个口。
!首先准备单元
/prep7
et,1,188 $mp,ex,1,1e7 $sectype,1,beam,rect $secdata,0.1,0.2
!然后创建模型
k,1 $k,2,5 $k,3,,5
l,1,2 $lgen,4,1,1,1,5 $lgen,5,1,4,1,,5
!X方向的这些梁是不释放自由度的，必须首先创建！先不要急着做L，1，3
!梁从L1 到 L20,下一个线从L21开始
l,1,3 $lgen,4,21,21,1,,5 $lgen,5,21,24,1,5
!这些檩条，每根（线）都是铰接于梁上，需要释放两端自由度。
nummrg,kp
!合并KEYPOINT，默认容许范围1e-4，以较小编号合并
lmesh,all
!类型都只有一个，默认就行了。MESH后就不能MERGEL了
!接下来ENDRELEASE就行了。但是为了观察结果，我们给所有关键点都施加铰支约束，通过在梁和檩条上施加力观察弯曲变形的传递，以确定是否达到要求。
dk,ux $dk,uy $dk,uz
endrelease,,,ball
!通过cplist查看耦合表
!下面我们对55，29，69，85四个NODE施加Z方向的力，观察变形
f,55,fz,-100 $f,29,fz,-100 $f,69,fz,-100 $f,85,fz,-100
!下面我们计算并观察变形
/solu $solve
/post1 $pldisp
!适当转动视角，可以注意到：梁（X）方向上弯矩可以传递，檩条（Y）方向弯矩不能传递。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
三、总结
通过实验，我们注意到，对于四条线交于一点的情况，ENDRELEASE命令对于编号最小的线上的单元不改变。对于编号较大的线上的单元，创建两个新的节点并分别用他们取代两侧单元（本例中的檩条）的节点。然后将原节点（依然连接在梁上）和两个新节点的进行UX，UY，UZ 三个方向上的位移耦合（边线上是2个点耦合）。
由此，连线的创建顺序是很有必要的，将同一类结构一起复制，然后再进行其他结构的复制，LINE的编号上形成了分类，清晰明了，易于以后分析，也不容易出错。本文由MaoYJ原创，欢迎讨论。