铰接空心斜板桥网格划分方法

twg2012 发表于 2010-10-27 21:52:38 | 显示全部楼层 | 阅读模式
我采用的是四面体自由划分的方法,但是觉得划分得并不是很规则,我想用六面体划分,但是提示不行,请问怎样划分较好?
命令流如下(如有不当之处,请指正):
finish
/clear
/prep7
/title External Transverse Prestress for Sressgthen Girder Bridges
et,1,solid45
et,2,solid45
et,3,mesh200
mp,dens,1,2400
mp,ex,1,3e10
mp,prxy,1,0.1667
mp,dens,2,2600
mp,ex,2,3.3e10
mp,prxy,2,0.1667
k,1,0.99,0,0
k,2,1,0,0
k,3,1,0.6,0
k,4,0.99,0.6,0
a,1,2,3,4
rectng,0,0.99,0,0.6
cyl4,0.27,0.26,0.19,-180
cyl4,0.27,0.34,0.19,180
rectng,0.08,0.46,0.26,0.34
agen,2,3,5,1,0.45,0,0,100
asel,s,,,3,8,1
cm,ar,area
allsel
asba,2,ar
allsel
vext,1,9,8,0,0,15.0
vsel,s,loc,x,0,0.99
vatt,2,2,2
allsel
vsel,s,loc,x,0.99,1
vatt,1,1,1
k,1000,2.641873,0,4.575858
k,1001,2.641873,0.6,4.575858
a,5,8,1001,1000
agen,2,21,0,0,0,0,12.6
cskp,100,0,5,1000,8
csys,0
vsel,s,,,1,2,1
cm,vv,volu
asel,s,,,21,22,1
cm,aa,area
vsba,vv,aa
csys,100
vsel,s,loc,z,-0.866,0
vsel,r,loc,y,0,0.6
cm,vvv,volu
vdele,vvv
asel,s,loc,z,-0.866,0
asel,r,loc,y,0,0.6
cm,aaa,area
adele,aaa
csys,0
wpave,0,0,12.6
wprota,0,0,-60
cswpla,101,0
csys,101
vsel,s,loc,z,0,1.20
vsel,r,loc,y,0,0.6
cm,vvvv,volu
vdele,vvvv
asel,s,loc,z,0,1.20
asel,r,loc,y,0,0.6
cm,aaaa,area
adele,aaaa
vsel,all
vglue,all
esize,0.05
mshkey,0
mshape,1,3d
vmesh,all
vsel,all
vgen,9,all,,,1,0,0,1000
横截面.JPG
精彩评论倒序浏览

2355查看4评论

twg2012 发表于 2010-11-11 13:54:44
现在上传我的命令流和划分示意图,如有不当之处,请不吝指正!
finish  
/clear  
/prep7  
/title External Transverse Prestress for Sressgthen Girder Bridges  
*SET,qmpz,2*0.8625*3.5*23000
*SET,rxdlg,8000
*SET,g1,(2*rxdlg+qmpz)/9
*SET,jfarea,0.01*0.6
*SET,g2,jfarea*26000
*SET,kxbarea,0.3123
*SET,g3,kxbarea*26000   
*SET,g,g1+g2+g3
*SET,B,1
*SET,H,1
*SET,A1,B*H
*SET,I1,B*H**3/12
*SET,I2,H*B**3/12
et,1,solid45
et,2,solid45
mp,dens,1,2600  
mp,ex,1,3.45e10
mp,prxy,1,0.1667
mp,dens,2,0.1*g/(kxbarea+jfarea)
mp,ex,2,3.45e10
mp,prxy,2,0.1667
k,1,0.99,0,0
k,2,1,0,0   
k,3,1,0.6,0
k,4,0.99,0.6,0  
a,1,2,3,4   
rectng,0,0.99,0,0.6
cyl4,0.27,0.26,0.19,-180
cyl4,0.27,0.34,0.19,180
rectng,0.08,0.46,0.26,0.34  
agen,2,3,5,1,0.45,0,0,100   
asel,s,,,3,8,1  
cm,ar,area
allsel  
asba,2,ar   
vext,1,9,8,0,0,17.0
vsel,s,loc,x,0,0.99
vatt,2,2,2  
allsel  
vsel,s,loc,x,0.99,1
vatt,1,1,1  
allsel  
vplot   
wprota,0,0,-30
cswpla,100,0
vsbw,all
asbw,all
lsbw,all
FLST,2,2,6,ORDE,2   !从log文件中复制出来的
FITEM,2,3   
FITEM,2,5   
VDELE,P51X  
FLST,2,4,5,ORDE,4   
FITEM,2,22  
FITEM,2,-23
FITEM,2,31  
FITEM,2,40  
ADELE,P51X  
FLST,2,4,5,ORDE,4   
FITEM,2,9   
FITEM,2,24  
FITEM,2,34  
FITEM,2,36  
ADELE,P51X  
FLST,2,2,5,ORDE,2   
FITEM,2,32  
FITEM,2,38  
ADELE,P51X  
FLST,2,4,5,ORDE,4   
FITEM,2,1   
FITEM,2,21  
FITEM,2,30  
FITEM,2,35  
ADELE,P51X  
ADELE,      39  
ADELE,      33  
FLST,2,5,4,ORDE,5   
FITEM,2,7   
FITEM,2,22  
FITEM,2,52  
FITEM,2,54  
FITEM,2,63  
LDELE,P51X  
FLST,2,8,4,ORDE,8   
FITEM,2,2   
FITEM,2,4   
FITEM,2,-5  
FITEM,2,19  
FITEM,2,49  
FITEM,2,62  
FITEM,2,73  
FITEM,2,75  
LDELE,P51X  
FLST,2,8,4,ORDE,8   
FITEM,2,3   
FITEM,2,6   
FITEM,2,9   
FITEM,2,12  
FITEM,2,31  
FITEM,2,50  
FITEM,2,67  
FITEM,2,72  
LDELE,P51X  
FLST,2,4,4,ORDE,4   
FITEM,2,30  
FITEM,2,32  
FITEM,2,69  
FITEM,2,76  
LDELE,P51X  
LDELE,       1  
FLST,2,4,3,ORDE,4   
FITEM,2,1   
FITEM,2,3   
FITEM,2,110
FITEM,2,112
KDELE,P51X
FLST,2,2,3,ORDE,2   
FITEM,2,2   
FITEM,2,4   
KDELE,P51X  
FLST,2,2,3,ORDE,2   
FITEM,2,6   
FITEM,2,-7  
KDELE,P51X  
allsel  
vplot
csys,0
wpcsys
wpoff,0,0,15.6
wprota,0,0,-30  
vsbw,all
asbw,all
lsbw,all
wpcsys  
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,3   
VDELE,P51X  
FLST,2,3,5,ORDE,3   
FITEM,2,15  
FITEM,2,-16
FITEM,2,24  
ADELE,P51X  
FLST,2,4,5,ORDE,4   
FITEM,2,7   
FITEM,2,17  
FITEM,2,63  
FITEM,2,72  
ADELE,P51X  
FLST,2,7,5,ORDE,7   
FITEM,2,2   
FITEM,2,13  
FITEM,2,19  
FITEM,2,21  
FITEM,2,23  
FITEM,2,64  
FITEM,2,-65
ADELE,P51X  
FLST,2,3,5,ORDE,3   
FITEM,2,18  
FITEM,2,20  
FITEM,2,61  
ADELE,P51X  
FLST,2,6,5,ORDE,4   
FITEM,2,5   
FITEM,2,66  
FITEM,2,-69
FITEM,2,71  
ADELE,P51X  
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,3   
FITEM,2,22  
ADELE,P51X  
FLST,2,3,5,ORDE,3   
FITEM,2,14  
FITEM,2,57  
FITEM,2,70  
ADELE,P51X  
ADELE,      60  
ADELE,      62  
ADELE,      58  
ADELE,       4  
ADELE,      59  
allsel  
vplot   
FLST,2,11,4,ORDE,11
FITEM,2,18  
FITEM,2,20  
FITEM,2,23  
FITEM,2,26  
FITEM,2,28  
FITEM,2,35  
FITEM,2,-36
FITEM,2,41  
FITEM,2,-42
FITEM,2,118
FITEM,2,-119
LDELE,P51X  
FLST,2,11,4,ORDE,11
FITEM,2,1   
FITEM,2,5   
FITEM,2,11  
FITEM,2,13  
FITEM,2,24  
FITEM,2,-25
FITEM,2,32  
FITEM,2,34  
FITEM,2,49  
FITEM,2,63  
FITEM,2,111
LDELE,P51X  
FLST,2,9,4,ORDE,9   
FITEM,2,7   
FITEM,2,31  
FITEM,2,44  
FITEM,2,46  
FITEM,2,50  
FITEM,2,116
FITEM,2,121
FITEM,2,126
FITEM,2,-127
LDELE,P51X  
FLST,2,6,4,ORDE,6   
FITEM,2,14  
FITEM,2,19  
FITEM,2,54  
FITEM,2,106
FITEM,2,123
FITEM,2,131
LDELE,P51X  
FLST,2,2,4,ORDE,2   
FITEM,2,21  
FITEM,2,108
LDELE,P51X  
FLST,2,3,4,ORDE,3   
FITEM,2,2   
FITEM,2,114
FITEM,2,129
LDELE,P51X  
FLST,2,3,4,ORDE,3   
FITEM,2,27  
FITEM,2,33  
FITEM,2,103
LDELE,P51X  
LDELE,     104  
LDELE,     112  
LDELE,      10  
FLST,2,11,3,ORDE,4  
FITEM,2,15  
FITEM,2,17  
FITEM,2,20  
FITEM,2,-28
KDELE,P51X  
FLST,2,2,3,ORDE,2   
FITEM,2,14  
FITEM,2,16  
KDELE,P51X  
FLST,2,2,3,ORDE,2   
FITEM,2,11  
FITEM,2,19  
KDELE,P51X  
KDELE,      18  
allsel  
vplot  
save
vsel,all
nummrg,all  
numcmp,all         
wpoffs,0.27,0,0
wprota,0,0,90   
vsbw,all
wpoffs,0,0,0.45
vsbw,all
wpoffs,0,0.26,0  
wprota,0,90,0   
vsbw,all
wpoffs,0,0,-0.08
vsbw,all
wpcsys  
wpoffs,0.495,0,0
wprota,0,0,90   
vsbw,all
vsel,all
vglue,all
wpcsys
et,5,mesh200,6
esize,0.05  
amap,51,16,39,35,55
amap,85,39,17,65,35
amap,79,20,49,44,65
amap,48,49,21,53,44
amap,30,53,54,13,12
amap,28,62,19,16,55
amap,43,18,66,65,17
amap,55,66,23,20,65
amap,65,22,60,53,21
amap,70,60,59,53,54
amap,74,19,41,37,62
amap,88,41,18,66,37
amap,83,23,47,46,66
amap,72,22,47,46,60
amap,6,15,14,60,59  
lesize,33,,,156
vsel,all
vsweep,all  
体扫掠划分示意图.JPG
横截面示意图.JPG
举报 回复
twg2012 发表于 2010-11-3 12:52:12
其中一种解决方法:用体扫掠划分方法,网格比较有规律,皆为六面体,要用amap,vsweep命令……
举报 回复
zhanxiao86 发表于 2010-11-18 10:09:11
谢谢lz的分享 资料很好先留着
举报 回复
lzd2001124... 发表于 2010-12-6 18:04:25
1、自由划分是一种自动化程度较高的划分方法,一般在体不规则,通过切分仍然不容易切出规则形状的情况下使用,精度一般,容易造成奇异阵。
2、实体划分单元采用体扫掠的方法划分出的单元形状规则,计算精度有保障,实体各部分均能人为干预。一般有条件的话此为推荐方案。
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00