用板壳元shell63单元建模时,如下图示,两个平面互相垂直,如何使这两个平面保持为一个整体来受力.因为实际结构中,比如一个由钢板焊成的箱梁,内有加劲隔板,如何模拟.由于这两个面没有公共边,用了粘贴及搭接都不行,试问怎样简单实现整体性这一目的.
A:
分网时控制单元长度,使面与面之间有重合节点或相邻节点,然后使用Merge Items或Adjacent Regions节点既可把面与面给“焊起来”。
Adjacent Regions命令时,必须先将一个区域的节点和相邻区域的单元选择好,再自动生成约束方程,是这样的吗,对整体计算的精度应有保证吧
Adjacent Regions命令的使用确如楼上所说,但对整体计算精度能否保证就不好妄加评论了。其实我认为SHELL63单元连接的方法是多种多样的,如何选择关键在于快捷,准确。当同位置的结点难以实现,Adjacent Regions命令又较困难时,那么还有相邻节点重合(Coincident Nodes),只要知道要求重合节点的距离容差就行了
Q:
请教ansys中梁壳两种单元一起建模时遇到的问题。
我分析了一个地下结构,楼板用shell63单元。
立柱用BEAM4单元。
结果出现的问题好奇怪,说单元1没有定义材料类型。
其实不是的。
后来我建了一个简单的模型,还是同样的问题。
我非常不解。
我把模型发了上来,哪位好心人能否帮我诊断一下。谢谢了。
A:
/prep7
et,1,shell63
mp,ex,1,3e10
mp,prxy,1,0.2
mp,dens,1,2500
keyopt,1,3,2!梁壳一起建模时必须要考虑shell中面的转动刚度
r,1,1
et,2,beam4
mp,ex,2,3e10
mp,prxy,2,0.2
mp,dens,2,2500
r,2,0.785,0.049,0.049,1,1,,,0.098,,,,
k,1,0,0,10
k,2,0,0,0
k,3,10,0,0
k,4,10,0,10
a,1,2,3,4
k,9,0,-5,10
k,10,0,-5,0
k,11,10,-5,0
k,12,10,-5,10
a,9,10,11,12
!a,1,9,12,4
wplane,100,0,0,0,10,0,0,0,10,0
csys,4
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
wplane,100,0,0,0,0,0,-10,0,10,0
csys,wp
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
nummrg,all
numcmp,all
csys,0!恢复原始的坐标系
l,,65
lsel,s,,,27
lesize,all,1
latt,2,2,2
lmesh,all
allsel
aatt,1,1,1
mshkey,0
mshape,0
esize,2
amesh,all
allsel
nummrg,all
numcmp,all
allsel
nsel,s,loc,y,-5
d,all,all
allsel
asel,s,loc,y,0
sfa,all,,pres,57e3
allsel
/solu
antype,static
solve
finish
我把模型发了上来,出现的错误真是令人匪夷所思。
我反复试验了半天,发现了问题。
把csys,0 前面的nummrg,ALL和NUMCMP,ALL删除。
这样后面建立的立柱线条的编号就没有规律了。
然后我手动划分立柱单元。
这样后就没有问题了。当然在我发的这个简单模型中只有一根
立柱,在
l,,65
lsel,s,,,27
lesize,all,1
latt,2,2,2
lmesh,all
中只需要改变线条的编号,来划分立柱单元就可以了。
我总结了一下问题:
就是在建立立柱几何模型前,不能将已经建立的板(已经用asbw,all命令细分了)进行nummrg,all和numcmp,all。
但是我kaisi不得其解。这不是很正常的操作吗?
何况出现的问题也太奇怪了。
请有兴趣的试验一下。谢谢。
这样改了下没问题
fini
/clear
/pnum,kp,1
/prep7
et,1,shell63
mp,ex,1,3e10
mp,prxy,1,0.2
mp,dens,1,2500
keyopt,1,3,2!梁壳一起建模时必须要考虑shell中面的转动刚度
r,1,1
et,2,beam4
mp,ex,2,3e10
mp,prxy,2,0.2
mp,dens,2,2500
r,2,0.785,0.049,0.049,1,1,,,0.098,,,,
k,1,0,0,10
k,2,0,0,0
k,3,10,0,0
k,4,10,0,10
a,1,2,3,4
k,9,0,-5,10
k,10,0,-5,0
k,11,10,-5,0
k,12,10,-5,10
a,9,10,11,12
!a,1,9,12,4
wplane,100,0,0,0,10,0,0,0,10,0
csys,4
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
wplane,100,0,0,0,0,0,-10,0,10,0
csys,wp
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
!nummrg,all
numcmp,all
csys,0!恢复原始的坐标系
nn_l=lsinqr(0,14)
!numstr,line,nn_l+1
l,,65
lsel,s,,,nn_l+1
lesize,all,1
latt,2,2,2
lmesh,all
allsel
aatt,1,1,1
mshkey,0
mshape,0
esize,2
amesh,all
allsel
nsel,s,loc,y,-5
d,all,all
allsel
asel,s,loc,y,0
sfa,all,,pres,57e3
allsel
/solu
antype,static
solve
finish
这样是不是更好?所有编号自己控制,嘿嘿。
fini
/clear
/pnum,kp,1
/prep7
/ESHAPE,1.0
/VIEW,1,1,1,1
/ANG,1
et,1,shell63
mp,ex,1,3e10
mp,prxy,1,0.2
mp,dens,1,2500
keyopt,1,3,2!梁壳一起建模时必须要考虑shell中面的转动刚度
r,1,1
et,2,beam4
mp,ex,2,3e10
mp,prxy,2,0.2
mp,dens,2,2500
r,2,0.785,0.049,0.049,1,1,,,0.098,,,,
*do,kk,1,6
k,kk,0,0,(kk-1)*2
*enddo
l,(1:5),(2:6)
LGEN,6,all,,,2,,,,1
*do,kk,1,6
l,(kk: (kk+5*6-1)),((kk+6): (kk+5*6))!此处为了帖子中:不变成笑脸在后面加了空格
*enddo
lll=0
*do,jj,1,21,5
*do,kk,jj,jj+4
al,kk,kk+5,kk+30+lll,kk+31+lll
*enddo
lll=lll+1
*enddo
LGEN,2,all,,,,-5,,,1
lll=0
*do,jj,1+60,21+60,5
*do,kk,jj,jj+4
al,kk,kk+5,kk+30+lll,kk+31+lll
*enddo
lll=lll+1
*enddo
nn_l=lsinqr(0,14)
l,4*6+2,4*6+2+6*6
lsel,s,,,nn_l+1
lesize,all,1
latt,2,2,2
lmesh,all
allsel
aatt,1,1,1
mshkey,0
mshape,0
esize,2
amesh,all
allsel
nsel,s,loc,y,-5
d,all,all
allsel
asel,s,loc,y,0
sfa,all,,pres,57e3
allsel
/solu
antype,static
solve
finish
/post1
PLESOL,S,X,0,1
问题找到了,你在同一点上定义了两种材料,实际上你的材料1和2是一样的。故改成以下没问题。
如果一定要两种材料,那只有自己控制定义后面的线号numstr,line,nn_l+1,不进行nummrg,all也一样可以通过。
fini
/clear
/pnum,kp,1
/prep7
et,1,shell63
mp,ex,1,3e10
mp,prxy,1,0.2
mp,dens,1,2500
keyopt,1,3,2!梁壳一起建模时必须要考虑shell中面的转动刚度
r,1,1
et,2,beam4
r,2,0.785,0.049,0.049,1,1,,,0.098,,,,
k,1,0,0,10
k,2,0,0,0
k,3,10,0,0
k,4,10,0,10
a,1,2,3,4
k,9,0,-5,10
k,10,0,-5,0
k,11,10,-5,0
k,12,10,-5,10
a,9,10,11,12
!a,1,9,12,4
wplane,100,0,0,0,10,0,0,0,10,0
csys,4
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
wplane,100,0,0,0,0,0,-10,0,10,0
csys,wp
*do,i,1,4,1
wpave,0,0,2
asbw,all
*enddo
nummrg,all
numcmp,all
csys,0!恢复原始的坐标系
nn_l=lsinqr(0,14)
numstr,line,nn_l+1
l,,65
lsel,s,,,nn_l+1
lesize,all,1
latt,1,2,2!请注意此处材料改为1
lmesh,all
allsel
aatt,1,1,1
mshkey,0
mshape,0
esize,2
amesh,all
allsel
nsel,s,loc,y,-5
d,all,all
allsel
asel,s,loc,y,0
sfa,all,,pres,57e3
allsel
/solu
antype,static
solve
finish下载本文