箱涵问题

xugangseu

我正在做一个箱涵，主要是求出弯距，再用弯距来配筋，首先考虑的是简单情况，考虑平面应变问题，可是plane4没有弯距输出，beam又不能网格划分，郁闷了 请问该怎么办那？

xugangseu

现在是用solid65单元，首先建立面单元，mesh200，划分网格比较麻烦，原来是打算把所有面合并后再划分网格，不过那样网格就划分不了了，这能每个面单独画好网格再合并，但是vdrag成体的时候，出现了很多面，郁闷，附命令流

xugangseu

/prep7
et,1,solid65
mp,ex,1,2.8e10
mp,nuxy,1,0.2
mp,dens,1,2500
et,2,combin14
real,2
r,1,60000
et,3,mesh200
KEYOPT,3,1,6
KEYOPT,3,2,0
k,1,0,0
k,2,0.8,0
k,3,0.8,0.3
k,4,0.55,0.55
k,5,1.55,0
k,6,1.55,0.3
k,7,0.3,0.8
k,8,0,0.8
k,9,0.3,1.8
k,10,0,1.8
a,1,2,3,4
a,2,5,6,3
a,1,4,7,8
a,7,9,10,8
/pnum,area,1  ！实体的颜色和编号
/pnum,line,1
lplot
aplot
wpstyl,,,,,,,2,1 !控制工作平面的显示
csys,4
kwpave,5
arsym,x,all,,,,0
kwpave,10
arsym,y,all,,,,0
csys,0
asel,s,loc,z,-0.1,0.1,,0  ！仅选择面
aadd,all ！布尔运算，面积相加
asel,all
k,39,0,0,10
l,1,39
asel,s,loc,z,-0.1,0.1,,0
vdrag,all,,,,,,2
asel,all
aplot
kwpave,1
asel,s,loc,z,-0.1,0.1,,0
mshape,0,2d
mshkey,0
aesize,all,0.1
amesh,all

hanty

你的做法也算对，较多的面是由于你开始的时候，将边分成了多份。可以按下面的命令流来完成网格。

/prep7
et,1,solid65
mp,ex,1,2.8e10
mp,nuxy,1,0.2
mp,dens,1,2500
et,2,combin14
real,2
r,1,60000
et,3,mesh200
KEYOPT,3,1,6
KEYOPT,3,2,0
k,1,0,0
k,2,0.8,0
k,3,0.8,0.3
k,4,0.55,0.55
k,5,1.55,0
k,6,1.55,0.3
k,7,0.3,0.8
k,8,0,0.8
k,9,0.3,1.8
k,10,0,1.8
a,1,2,3,4
a,2,5,6,3
a,1,4,7,8
a,7,9,10,8
/pnum,area,1  ！实体的颜色和编号
/pnum,line,1
lplot
aplot
wpstyl,,,,,,,2,1 !控制工作平面的显示
csys,4
kwpave,5
arsym,x,all,,,,0
kwpave,10
arsym,y,all,,,,0
csys,0
asel,s,loc,z,-0.1,0.1,,0  ！仅选择面
aadd,all ！布尔运算，面积相加
asel,all
k,39,0,0,10
l,1,39
asel,s,loc,z,-0.1,0.1,,0
vdrag,all,,,,,,2
asel,all
aplot
kwpave,1
asel,s,loc,z,-0.1,0.1,,0
mshape,0,2d
mshkey,0
aesize,all,0.1
amesh,all
EXTO,VSWE,AUTO,OFF  
VSWEEP,1,17,34  
ALLSEL,ALL  
EPLOT

xugangseu

楼上的能否解释下exto的具体用法？

xugangseu

EXTOPT，Lab，Val1，Val2，Val3（扫掠体生成控制选项）
EXTOPT，ACLEAR，Val1（指定在体扫掠完成后是否删除源面上的网格）
EXTOPT，ATTR，Val1，Val2，Val3（定义扫掠体单元属性）
EXTOPT，ESIZE，Val1，Val2（定义扫掠方向的单元尺寸）
EXTOPT，VSWE，AUTO，Val2（指定在体扫掠中是否自动选择源面和目标面）
EXTOPT，VSWE，TETS，Val2（指定在体扫掠中对无法扫掠的体用四面体划分网格）

xugangseu

这个是完整的命令流，可惜我得截面弯距还在摸索中，不会做
最近在做一个关于箱涵的分析，因为我想考虑纵向的影响，并且想真实的模拟出车轮的作用，所以是个空间问题，遇到一下几个问题：  
1.在面积上施加梯形面荷载（用以模拟土压力）,侧面是否要加约束？  
2.底面约束加上一个弹性系数（用以模拟基床系数） ，弹簧单元底面完全约束
3.关于sweep网格，因为我是先在箱涵横断面（mesh200）上面划分平面网格，然后sweep，但是遇到的问题是，我是不是要定义平面单元和实体单元两种，平面单元实际上只是起了过渡作用呢？在我实际分析中用不到平面网格，那么这种平面单元的特性定义有何要求呢（我用的是mesh200过渡单元）？


fini
/cle
/prep7!定义三种单元类型，用以模拟砼、弹性地基、过渡单元，所有单位采用国际单位
et,1,solid65
mp,ex,1,2.8e10
mp,nuxy,1,0.2
mp,dens,1,2500
et,2,combin14
R,1,113636, , ,  
et,3,mesh200
KEYOPT,3,1,6
KEYOPT,3,2,0
k,1,0,0
k,2,0.8,0
k,3,0.8,0.3
k,4,0.55,0.55
k,5,1.55,0
k,6,1.55,0.3
k,7,0.3,0.8
k,8,0,0.8
k,9,0.3,1.8
k,10,0,1.8
a,1,2,3,4
a,2,5,6,3
a,1,4,7,8
a,7,9,10,8
/pnum,area,1  ！实体的颜色和编号
/pnum,line,1
lplot
aplot
wpstyl,,,,,,,2,1 !控制工作平面的显示
csys,4
kwpave,5
arsym,x,all,,,,0
kwpave,10
arsym,y,all,,,,0
csys,0
asel,s,loc,z,-0.1,0.1,,0  ！仅选择面
aadd,all ！布尔运算，面积相加
asel,all
k,39,0,0,10
l,1,39
asel,s,loc,z,-0.1,0.1,,0
vdrag,all,,,,,,2
asel,all
aplot
kwpave,1
asel,s,loc,z,-0.1,0.1,,0
aatt,1,,3
mshape,0,2d
mshkey,0
aesize,all,0.1
amesh,all

allsel,all
gplot
type,1

!!!


vsweep,1,17,34  !荷载

numcmp,all

nsel,s,loc,y,0  !建立弹簧单元
nplot!!
FLST,4,528,1,ORDE,5
FITEM,4,80  
FITEM,4,88  
FITEM,4,-119
FITEM,4,2177
FITEM,4,-2671   
NGEN,2,9008,P51X, , ,0,-1,0,1  !!
type,2

e,80,9088
*do,i,88,119,1
e,i,i+9008
*enddo
*do,i,2177,2671,1
e,i,i+9008
*enddo
allsel,all
gplot
finish
/solu
nsel,s,loc,y,-1
d,all,all  !!fini
allsel,all
gplot!!

sfgrad,pres,0,y,0,-5166
esel,s,type,,1
nsel,r,loc,x,0
sf,all,pres,23400
allsel,all
sfgrad,pres,0,y,0,-5166
esel,s,type,,1
nsel,r,loc,x,3.1
sf,all,pres,23400
sfgrad,pres,0,y,0,0
esel,s,type,,1
nsel,r,loc,x,0
sf,all,pres,13720
!!
allsel,all

nsel,s,loc,y,3.6
sf,all,pres,63060
allsel,all
gplot
/psf,pres,norm,2,0,1
acel,0,9.8,0
antype,0
/status,solu
solve
finish
/post1
pldisp,1
plesol,s,eqv,0,1

xugangseu

这个是我用shell93单元做的箱涵模型 ，可惜厚度大的（0.3m）的时候模型在结合处好像不是太好，不知对结果是否有影响
/prep7
et,1,shell93
r,1,0.3
mp,ex,1,2.8e10
mp,prxy,1,0.2
mp,dens,1,2500
k,1,
k,2,3.1
k,3,3.1,3.6
k,4,0,3.6
a,1,2,3,4
VEXT,1, , ,0,0,10,,,,
vdele,1
adele,1,2,1
esize,0.5
amesh,all
fini

xugangseu

箱涵mesh

xugangseu

这个对箱涵结构分析有影响吗？

wei1012

建议你用beam单元，采用自定义截面的方式进行分析，好处能直接获得单元的弯矩！

xugangseu

按楼上的建议，参考了一本《ansys土木...》的隧道例子，现附上命令流，不过有几点疑惑：
1，隧道两侧的节点荷载是线性变化的，不知是否可以施加坡度载荷，我是在一个一个节点上施加的，要是能在线上施加就ok了，sfl好像不行
2，有必要在四周都加弹簧约束吗？在底部加一个不就ok了吗？
3，另外如果把所有的弹簧单元都杀死，是否意味这个模型没有约束？那这样还可以求解吗？

xugangseu

/nopr
keyw,pr_struc,1
/prep7
et,1,beam3
et,2,combin14
r,1,0.5,0.0104167,0.5
r,2,30e6
mp,ex,1,30e9
mp,prxy,1,0.2
mp,dens,1,2500

!!
k,1
k,2,5
k,3,5,5
k,4,,5
k,5,-5,5
k,6,-5
lstr,1,2
lstr,2,3
lstr,2,3
lstr,3,4
lstr,4,5
lstr,5,6
lstr,6,1
lstr,1,4


!!
lesize,all,0.5
type,1
mat,1
real,1
lmesh,all
/pnum,elem,1

!!
nsel,s,loc,y,0
ngen,2,100,all,,,,-1
nsel,s,loc,y,5
ngen,2,200,all,,,,1
nsel,s,loc,x,5
ngen,2,300,all,,,1
nsel,s,loc,x,-5
ngen,2,400,all,,,-1
allsel,all
nplot

ndele,102,302,200
ndele,212,312,100
ndele,232,432,200
ndele,142,442,300
nummrg,all
numcmp,all

/rep
!!
type,2
mat,1
real,2
e,1,70
e,2,117
e,12,127
e,32,129
e,42,128
*do,i,3,11,1
e,i,i+68
*enddo
*do,i,13,21,1
e,i,i+95
*enddo
*do,i,22,31,1
e,i,i+67
*enddo
*do,i,33,41,1
e,i,i+66
*enddo
*do,i,43,51,1
e,i,i+75
*enddo
*do,i,52,60,1
e,i,i+28
*enddo
eplot
!!
/sol
nsel,s,loc,y,-1
nsel,a,loc,y,6
nsel,a,loc,x,6
nsel,a,loc,x,-6
d,all,all
allsel,all
/rep
!!
acel,0,10
nsel,s,loc,y,5
f,all,fy,-60000
nsel,s,loc,y,0
f,all,fy,76850
!!
allsel,all
f,10,fx,-32700
f,11,fx,-62400
f,12,fx,-59400
f,13,fx,-56400
f,14,fx,-53400
f,15,fx,-50400
f,16,fx,-47400
f,17,fx,-44400
f,18,fx,-41400
f,19,fx,-38400
f,20,fx,-17700
f,10,fx,-32700

f,41,fx,17700
f,42,fx,38400
f,43,fx,41400
f,44,fx,44400
f,45,fx,47400
f,46,fx,50400
f,47,fx,53400
f,48,fx,56400
f,49,fx,59400
f,50,fx,62400
f,51,fx,32700

!!
nropt,full
allsel
outres,all,all
solve
fini
!!
/post1
pldisp,1
!!
etable,,smisc,6
etable,,smisc,12
etable,,smisc,1
etable,,smisc,7
etable,,smisc,2
etable,,smisc,8
plls,smis6,smis12,-1,0 !弯距
plls,smis1,smis7,1,0  ！轴力
plls,smis2,smis8,1   ！剪力

xugangseu

模型 见图

xugangseu

不清楚的话可以点击放大看

xugangseu

我后来又做了一个单箱涵的，主要问题出在四个角点的受力上不知如何取，上面的双箱涵感觉角点的受力取得没有道理，单箱涵的主要问题是在下面两个角点处弯距和剪力集中：
/nopr
keyw,pr_struc,1
/prep7
et,1,beam3
et,2,combin14
r,1,0.5,0.0104167,0.5
r,2,30e6  !这个单位应是N/m ，和弹簧单元数目应该有关系？
mp,ex,1,3e10  !一般数据中只有地基反应模量，单位N/m3，
mp,prxy,1,0.2
mp,dens,1,2500

!!
k,1
k,2,5
k,3,5,5
k,4,,5

lstr,1,2
lstr,2,3
lstr,3,4
lstr,4,1
!!
lesize,all,0.5
type,1
mat,1
real,1
lmesh,all
/pnum,elem,1
/pnum,node,1
/rep

!!
nsel,s,loc,y,0
ngen,2,100,all,,,,-1
nsel,s,loc,y,5
ngen,2,200,all,,,,1
nsel,s,loc,x,5
ngen,2,300,all,,,1
nsel,s,loc,x,0
ngen,2,400,all,,,-1
allsel,all
nplot
/pnum,node,1
/rep

ndele,102,302,200
ndele,212,312,100
ndele,222,422,200
ndele,101,401,300
nummrg,all
numcmp,all

/rep
!!
type,2
mat,1
real,2
e,1,78
e,2,68
e,12,79
e,22,80

*do,i,3,11,1
e,i,i+38
*enddo
*do,i,13,21,1
e,i,i+46
*enddo
*do,i,23,31,1
e,i,i+27
*enddo
*do,i,32,40,1
e,i,i+37
*enddo
eplot
!!
/sol
nsel,s,loc,y,-1
nsel,a,loc,y,6
nsel,a,loc,x,6
nsel,a,loc,x,-1
d,all,all
allsel,all
/rep
!!
acel,0,10
nsel,s,loc,y,5
nsel,u,node,,12,22,10
f,all,fy,-60000
nsel,s,loc,y,0
nsel,u,node,,1,2,1
f,all,fy,76850
!!
allsel,all
f,21,fx,-38400
f,20,fx,-41400
f,19,fx,-44400
f,18,fx,-47400
f,17,fx,-50400
f,16,fx,-53400
f,15,fx,-56400
f,14,fx,-59400
f,13,fx,-62400
!!

f,32,fx,38400
f,33,fx,41400
f,34,fx,44400
f,35,fx,47400
f,36,fx,50400
f,37,fx,53400
f,38,fx,56400
f,39,fx,59400
f,40,fx,62400
!!
f,2,fx,-32700   !四个角点受力取值估计有问题？
f,2,fy,38425     !f,1,fx,65400？
f,12,fx,-17700   !f,1,fy,76850？ 其余四角点类似
f,12,fy,-30000
f,22,fx,17700
f,22,fy,-30000
f,1,fx,32700
f,1,fy,38425

nropt,full
allsel
outres,all,all
solve
fini
!!
/post1
pldisp,1
!!
etable,,smisc,6
etable,,smisc,12
etable,,smisc,1
etable,,smisc,7
etable,,smisc,2
etable,,smisc,8
plls,smis6,smis12,-1,0 !弯距
plls,smis1,smis7,1,0  !轴力
plls,smis2,smis8,1   !剪力

xugangseu

注意角点的集中