求助：方圆k型相贯节点极限承载力分析

Crystalfly

本人用ansys分析方圆k型相贯节点的极限承载力和参数分析时遇到很多问题，现请教各位大师：
模型尺寸：方主管长1200mm，厚10mm，边长200mm
          圆支管长500mm，厚7mm，直径120mm，两支管间隙40mm
1．        在划分网格时，因为主管和支管的单元类型相同但实常数不同，如何用命令流方式将定义好的两种实常数分配到主、支管上。
2．        对主支管端部加约束时，主管一端固定一端辊动支撑，支管两端都是简支。现已将主管约束加上，但支管端部节点的约束却加不成功。
我在支管端部加约束时，先定义了两个局部柱坐标系；然后将两支管端部的节点坐标系分别转动到这两个局部坐标系，之后在支管端部节点加径向和切向约束。
local,12,1,g/2+(d2-t2)/2/cos(jiao),(b1-t1)/2,0,jiao,0,90
local,13,1,-(g/2+(d2-t2)/2/cos(jiao)),(b1-t1)/2,0,jiao,-90,0
！将右端圆支管端部节点的节点坐标系转到局部坐标系12
csys,12            
nsel,s,loc,z,l2
nrotat,all
allsel
！将右端圆支管端部节点的节点坐标系转到局部坐标系13
csys,13
nsel,s,loc,z,l2
nrotat,all
allsel
！求解
/solu
antype,static
！主管两端加约束
csys,0
nsel,s,loc,x,-l1/2
d,all,,0,,,,all,,,,,
allsel
nsel,s,loc,x,l1/2
d,all,,0,,,,uy,uz,rotx,roty,rotz,
allsel
！支管两端加约束
csys,12
nsel,s,loc,z,l2
d,all,,0,,,,ux,,,,,
csys,13
nsel,s,loc,z,l2
d,all,,0,,,,ux,,,,,
我个人感觉步骤上没有什么问题，但总是加不成功，ansys不执行上面的某些命令流，我就用GUI方式加上去，但总也加不成功（用肉眼看上去，总感觉方向不对劲儿）。不知哪里出了问题，请大师们多多指教。附图如下：
我的QQ：386511902。