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

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

twg2012 发表于 2010-11-3 12:52:12

其中一种解决方法：用体扫掠划分方法，网格比较有规律，皆为六面体，要用amap,vsweep命令……

zhanxiao86 发表于 2010-11-18 10:09:11

谢谢lz的分享{:4_117:} 资料很好先留着

lzd20011244 发表于 2010-12-6 18:04:25

1、自由划分是一种自动化程度较高的划分方法，一般在体不规则，通过切分仍然不容易切出规则形状的情况下使用，精度一般，容易造成奇异阵。
2、实体划分单元采用体扫掠的方法划分出的单元形状规则，计算精度有保障，实体各部分均能人为干预。一般有条件的话此为推荐方案。
页: [1]
查看完整版本: 铰接空心斜板桥网格划分方法