混凝土网格划分不了，请高手帮忙，灰常着急

蓝色雪片

finish
/clear
/prep7
et,1,solid65
et,2,link8
mp,ex,1,21945    !定义混凝土的弹性模量、泊松比和密度
mp,prxy,1,0.2
mp,dens,1,2300e-12
fc=23.1           !混凝土单轴抗压强度和单轴抗拉强度
ft=1.89
tb,concr,1        !定义混凝土材料及相关参数，关闭压碎
tbdata,,0.4,0.90,ft,-1
tb,miso,1,,11     !定义混凝土应力应变曲线，用miso模型
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.001,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
tbpt,,0.0033,fc*0.85
mp,ex,2,2.0e5    !钢筋材料的弹性模量和泊松比
mp,prxy,2,0.3
mp,ex,3,1.95e5   !预应力钢筋的弹性模量、泊松比和密度
mp,prxy,3,0.3
mp,dens,3,7921e-12
tb,biso,2        !定义钢筋应力应变曲线，用biso模型
tbdata,,300,0
r,1,140    !预应力钢绞线的公称截面面积
r,2,490.9  !直径25mm的钢筋截面面积490.9mm2
r,3,113.1  !直径12mm的钢筋截面面积113.1mm2
r,4,50.3   !直径8mm的钢筋截面面积50.3mm2
r,5,
!2.创建几何模型
k,1  !混凝土
k,2,,100.0010
k,3,100.9702,137.4000
k,4,46.1538,850.0000
k,5,943.8450,850.0000
k,6,889.0319,137.4000
k,7,990,100.0010
k,8,990,0
a,1,2,3,4,5,6,7,8
/view,1,1,1
/ang,1
vext,all,,,,,20000
rpr4,8,495.000,450.000,312.6000,,20000
vsbv,1,2
wpcsys,-1    !预应力筋
lsel,none
k,33,95,36.4000,3400.0000
k,34,95,36.4000,16600.0000
l,33,34
k,35,145,36.4000,20.0000
k,36,145,36.4000,19980.0000
l,35,36
k,37,195,36.4000,6100.0000
k,38,195,36.4000,13900.0000
l,37,38
k,39,245,36.4000,1600.0000
k,40,245,36.4000,18400.0000
l,39,40
k,41,295,36.4000,2400.0000
k,42,295,36.4000,17600.0000
l,41,42
k,43,345,36.4000,20.0000
k,44,345,36.4000,19980.0000
l,43,44
k,45,395,36.4000,4500.0000
k,46,395,36.4000,15500.0000
l,45,46
k,47,495,36.4000,6100.0000
k,48,495,36.4000,13900.0000
l,47,48
k,49,595,36.4000,4500.0000
k,50,595,36.4000,15500.0000
l,49,50
k,51,645,36.4000,20.0000
k,52,645,36.4000,19980.0000
l,51,52
k,53,695,36.4000,2400.0000
k,54,695,36.4000,17600.0000
l,53,54
k,55,745,36.4000,1600.0000
k,56,745,36.4000,18400.0000
l,55,56
k,57,795,36.4000,6100.0000
k,58,795,36.4000,13900.0000
l,57,58
k,59,845,36.4000,20.0000
k,60,845,36.4000,19980.0000
l,59,60
k,61,895,36.4000,3400.0000
k,62,895,36.4000,16600.0000
l,61,62
cm,sline,line
lsel,all
wpcsys,-1
wpcsys,-1
wprota,,,-90
wpoff,,,-139.0012
vsbw,all
wpoff,,,-134.9572
vsbw,all
wpoff,,,-71.0435
vsbw,all
wpoff,,,-149.9981
vsbw,all
wpoff,,,-150.0031
vsbw,all
wpoff,,,-71.0384
vsbw,all
wpoff,,,-134.9623
vsbw,all
wpcsys,-1
vsel,all  
wprota,,90
wpoff,,,-36.4000
vsbw,all
wpoff,,,-41.1304
vsbw,all
wpoff,,,-22.4976
vsbw,all
wpoff,,,-27.5039
vsbw,all
wpoff,,,-9.8960
vsbw,all
wpoff,,,-91.5584
vsbw,all
wpoff,,,-221.0416
vsbw,all
wpoff,,,-221.0415
vsbw,all
wpoff,,,-91.5585
vsbw,all
wpoff,,,-56.2029
vsbw,all

wpcsys,-1               !切出拉区纵筋竖向位置
vsel,s,loc,y,0,77.5034   
wprota,,,-90
wpoff,,,-45
vsbw,all
wpoff,,,-400
vsbw,all
wpoff,,,-100
vsbw,all
wpoff,,,-400
vsbw,all
vsel,all
vplot
wpcsys,-1       !切出加载位置和加载面
wpoff,,,525
vsbw,all
wpoff,,,8475
vsbw,all
wpoff,,,2000
vsbw,all
wpoff,,,8475
vsbw,all
nummrg,all,,,,low !整理
numcmp,all !压缩编号
!3.划分钢筋网格
elemsiz=50
lsel,s,loc,x,45  !拉区普通钢筋，定义组件，设置属性等
lsel,a,loc,x,445
lsel,a,loc,x,545
lsel,a,loc,x,945
lsel,r,loc,y,36.4000
cm,lqp,line
latt,2,3,2
lesize,all,elemsiz
lsel,s,loc,x,45   !洞口下方普通钢筋，定义组件，设置属性等
lsel,a,loc,x,945
lsel,r,loc,y,77.5034
cm,lqpt,line
latt,2,3,2
lesize,all,elemsiz
lsel,s,loc,x,139.0012  
lsel,a,loc,x,851.0038
lsel,r,loc,y,127.504
cm,lqpto,line
latt,2,3,2
lesize,all,elemsiz
lsel,s,loc,x,139.0012 !压区普通钢筋，定义组件，设置属性等
lsel,a,loc,x,345.0019
lsel,a,loc,x,645.0031
lsel,a,loc,x,851.0038
lsel,r,loc,y,818.8029
cm,yqp,line
latt,2,3,2
lesize,all,elemsiz
nummrg,all,,,,low !整理
numcmp,all !压缩编号
cmsel,s,lqp
cmsel,a,lqpt
cmsel,a,lqpto
cmsel,a,yqp
cm,gj,line
lmesh,all
lsel,all
!/eshape,1
eplot
vsel,all
vplot
钢筋网格已经划分，solid65单元建的实体网格划分不了，请高手帮忙，看看有没有什么好的划分办法，谢谢

蓝色雪片

采用分离式模型加的普通钢筋，在我能力范围内，尝试了不少次，结果混凝土网格划分不了，纠结了好久了，所以非常着急

senja

试了一下，发现我横向切割（楼主指加载位置的切割）的时候，有两个体切不了（可能是我软件问题，具体还不知道怎么回事），即为程序最后的561和562号体，这两体最后使用了vsweep划分，其他的mapped。







命令流在附件中，少部分修改过。

蓝色雪片

senja 发表于 2012-5-8 20:47
试了一下，发现我横向切割（楼主指加载位置的切割）的时候，有两个体切不了（可能是我软件问题，具体还不知 ...

很很感谢你，谢谢哦！

蓝色雪片

蓝色雪片 发表于 2012-5-8 21:23
很很感谢你，谢谢哦！

我自己在用自由网格划分时也是遇到了同样的问题，估计是我模型的问题。这样有自由网格划分的精度是不是比较差啊？真的十分感谢你。

senja

蓝色雪片 发表于 2012-5-8 21:34
我自己在用自由网格划分时也是遇到了同样的问题，估计是我模型的问题。这样有自由网格划分的精度是不是比 ...

自由网格划分，适应所有复杂难度的模型。
不过，如果单元为六面体单元，不推荐使用其退化单元，而自由网格划分只能生成四面体网格，也是不推荐使用的。也就是说，使用六面体单元，尽量划分六面体网格，如果操作难度很大，那么局部少量使用四面体单元。另外，自由网格划分，同等网格尺寸下网格的数量会成倍的增加，也就增加了计算量。

蓝色雪片

senja 发表于 2012-5-8 23:14
自由网格划分，适应所有复杂难度的模型。
不过，如果单元为六面体单元，不推荐使用其退化单元，而自由网 ...

嗯，谢谢你