ansys-SHELL单元

用板壳元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