问题描述：
ansys壳单元与梁单元的铰接问题
一个四边形壳单元与四根3D梁单元的铰接问题
先在同一位置建立2个不同节点,然后耦合他们的ux,uy,uz
在板的中间施加一x方向的力,计算不了啊.显示壳单元发生了
不受约束的刚体位移!
好像这个约束只能传递竖向的力（z方向）,x,y方向的力不能通过耦合节点传到下面的梁上哦!
finish
/clear
/prep7
k,1,1,1,1
k,2,1,0,1
k,3,0,1,1
k,4,1
k,5,1,1
k,6,1,0
k,7,0,1
k,8,
k,9,1,1,1
k,10,1,0,1
k,11,0,1,1
k,12,1
et,1,shell63
et,2,beam44
mp,ex,1,2.1E11
mp,prxy,1,0.25
sectype,1,beam,csolid
secdata,0.2
l,9,5
l,10,6
l,11,7
l,12,8
latt,1,2,1
lsel,s,length,1
esize,0.2
mshape,0
mshkey,1
lmesh,all
r,1,0.05
a,4,3,1,2
aatt,1,1,1
asel,all
esize,0.1
mshape,0
mshkey,1
amesh,all
cp,1,ux,13,26
cp,2,uy,13,26
cp,3,uz,13,26
cp,4,ux,1,36
cp,5,uy,1,36
cp,6,uz,1,36
cp,7,ux,19,25
cp,8,uy,19,25
cp,9,uz,19,25
cp,10,ux,7,46
cp,11,uy,7,46
cp,12,uz,7,46
/sol
dk,5,all
dk,6,all
dk,7,all
dk,8,all
f,105,fx,-2
solve
/post1
pldisp,1
prrsol
壳单元与梁单元的铰接问题啊!
这个是我做的一个实现壳与梁单元连接的例子哦.
是可以运行的,可是变形图不对啊!支座反力没有,铰接没有传递荷载啊.
一个四边形壳单元与四根3D梁单元的铰接问题
先在同一位置建立2个不同节点,然后耦合他们的ux,uy,uz
在板的中间施加一x方向的力,计算不了啊.显示壳单元发生了
不受约束的刚体位移!
好像这个约束只能传递竖向的力（z方向）,x,y方向的力不能通过耦合节点传到下面的梁上哦!
finish
/clear
/prep7
k,1,1,1,1
k,2,1,0,1
k,3,0,1,1
k,4,1
k,5,1,1
k,6,1,0
k,7,0,1
k,8,
k,9,1,1,1
k,10,1,0,1
k,11,0,1,1
k,12,1
et,1,shell63
et,2,beam44
mp,ex,1,2.1E11
mp,prxy,1,0.25
sectype,1,beam,csolid
secdata,0.2
l,9,5
l,10,6
l,11,7
l,12,8
latt,1,2,1
lsel,s,length,1
esize,0.2
mshape,0
mshkey,1
lmesh,all
r,1,0.05
a,4,3,1,2
aatt,1,1,1
asel,all
esize,0.1
mshape,0
mshkey,1
amesh,all
cp,1,ux,13,26
cp,2,uy,13,26
cp,3,uz,13,26
cp,4,ux,1,36
cp,5,uy,1,36
cp,6,uz,1,36
cp,7,ux,19,25
cp,8,uy,19,25
cp,9,uz,19,25
cp,10,ux,7,46
cp,11,uy,7,46
cp,12,uz,7,46
/sol
dk,5,all
dk,6,all
dk,7,all
dk,8,all
f,105,fx,-2
solve
/post1
pldisp,1
prrsol
壳单元与梁单元的铰接问题啊!
这个是我做的一个实现壳与梁单元连接的例子哦.
是可以运行的,可是变形图不对啊!支座反力没有,铰接没有传递荷载啊.
问题解答：
我来补答
首先该命令流有问题,不能直接运行,会提示1号 2号单元的惯性矩(转动惯量)是0或者负值,无法进行运算.
ansys提示是实常数的问题
我将你原先的命令改成如下：
r,2,0.05
a,4,3,1,2
aatt,1,2,1
把壳单元的实常数(ansys将实常数作为壳单元的截面厚度进行处理)的编号改成2即可
原因是：梁壳单元都会涉及到截面属性的问题,ansys定义壳单元的截面属性可以通过将实常数和截面属性进行混合,这样就不需要指定壳的截面属性了,这样实常数和截面的数据进行的传递.而你之前定义了梁的截面属性,你又定义了壳单元的实常数编号为1,这样ansys对实常数的调用就出现问题了,就会出现两种单元同事调用相同的实常数的错误.
其次,你的原先的耦合设置是对的,下面是我的运算结果,不过我个人认为,耦合的时候也可以讲旋转自由度同时耦合.
耦合的命令流如下：
cp,2001,all,13,26
cp,next,all,1,36
cp,next,all,19,25
cp,next,all,7,46
给出结果的变形图
这里面都能看到壳单元的变形影响到梁单元的变形
支座处的粉色单箭头的力就是反力
ansys提示是实常数的问题
我将你原先的命令改成如下：
r,2,0.05
a,4,3,1,2
aatt,1,2,1
把壳单元的实常数(ansys将实常数作为壳单元的截面厚度进行处理)的编号改成2即可
原因是：梁壳单元都会涉及到截面属性的问题,ansys定义壳单元的截面属性可以通过将实常数和截面属性进行混合,这样就不需要指定壳的截面属性了,这样实常数和截面的数据进行的传递.而你之前定义了梁的截面属性,你又定义了壳单元的实常数编号为1,这样ansys对实常数的调用就出现问题了,就会出现两种单元同事调用相同的实常数的错误.
其次,你的原先的耦合设置是对的,下面是我的运算结果,不过我个人认为,耦合的时候也可以讲旋转自由度同时耦合.
耦合的命令流如下：
cp,2001,all,13,26
cp,next,all,1,36
cp,next,all,19,25
cp,next,all,7,46
给出结果的变形图
这里面都能看到壳单元的变形影响到梁单元的变形
支座处的粉色单箭头的力就是反力
展开全文阅读