[如是我闻]壳与实体的连接

pfsr

壳与实体的连接一向以来是一个问题。现在总结了几种方法，抛块砖头，希望能有高手前来指点一二。

pfsr

一个管子，分为两段，一段用壳单元，另一段用实体单元。考虑对称。

pfsr

先用bonded contact

pfsr

再试试用SHSD。 由于shsd是根据target建立虚拟单元，而WB是根据contact的ID确认接触对，那只能借助于commands建立接触对在利用shsd。在commands中NamedSelection不可少。

pfsr

所用到的commands

/prep7

ET,5,TARGE170
keyopt,5,5,1
r,5
ET,6,Conta175
KEYOPT,6,2,2
KEYOPT,6,12,5
keyopt,6,4,0
R,5

cmsel,s,solid
TYPE,5
REAL,5
ESURF

cmsel,s,skal
TYPE,6
real,5
ESURF

ALLSEL,ALL
cncheck,adjust,5
SHSD,5,CREATE

fini
/solu

pfsr

第三种方法，在壳下面再建立一层很薄的壳。两个壳间共用节点，薄壳跟实体bonded contact

zgp0816

对于第二种方法,实体单位厚度方向也分为一层,这个是否考虑正确呢,精度能否保证呢

pfsr

最后用CP把壳上节点跟实体上的节点耦合自由度。这是用到的commands
cmsel,s,skal
*get,nnode,node,,count

*dim,nodlist,array,nnode
*vget,nodlist,node,,nlist

allsel

*do,i,1,nnode
cmsel,s,solid

nskal=nodlist(i)
nsel,r,loc,x,nx(nskal)
nsel,r,loc,y,ny(nskal)
nsel,r,loc,z,nz(nskal)
*get,nsolid,node,,num,min
nsel,a,,,nskal
cp,next,ux,nskal,nsolid
cp,next,uy,nskal,nsolid
cp,next,uz,nskal,nsolid

*enddo

allsel

[ 本帖最后由 pfsr 于 2008-3-10 20:20 编辑 ]

pfsr

到Classic检查一下节点耦合

pfsr

这是dsdb文件，在V11 SP1中完成。感谢各位指点。

pfsr

原帖由 zgp0816 于 2008-3-8 20:48 发表
对于第二种方法,实体单位厚度方向也分为一层,这个是否考虑正确呢,精度能否保证呢

可能是图小，看不清。在实体上我一共分了四层单元。

sh_lin30

原帖由 zgp0816 于 2008-3-8 20:48 发表
对于第二种方法,实体单位厚度方向也分为一层,这个是否考虑正确呢,精度能否保证呢

对solid186单元，单元算法取完全积分的话，可以只划分为一层

曾经做过测试，完全积分的一层单元和缩减积分的三层单元，应力结果相差不大

jania

对commond还不是很熟，但还是非常感谢版主pfsr

linuaries

多谢psfr的热心分享：）

whatinrain

[quote]原帖由 pfsr 于 2008-3-8 20:54 发表
最后用CP把壳上节点跟实体上的节点耦合自由度。这是用到的commands
msel,s,skal
*get,nnode,node,,count

非常谢谢老大分享。这个应该是cmsel吧。

pfsr

已更正，多谢whatinrain兄

a3050311118

谢谢，顶下

swan321

我曾经做过壳体和实体的接触与纯实体接触的比较，发现有很大的不同，一直没有解决

swan321

部门的工程师说：可能是自由度的传递有问题。建议带有接触的地方最好用实体。
或者分析的时候最好用一种类型

tjknight11

实体与壳连接部位的应力不可信，因为实体与壳每个节点代表的意义不同。
但根据圣维南原理，较远处的应力是可信的