我的汽车移动荷载命令流及请教 0 汽车驶过路面的命令流，路面中有个检查井，请高手看看是否有错？
现在我想看看整个过程中检查井的竖向位移、y向应力最大，也就是最容易破坏的地方该怎么操作啊。
另外这个阻尼系数是否有必要？
谢谢了！
Quote ini
/clear
/title,"transient"
/filename,transient
/config,nres,3000
!---------------------------------------------
*dim,nodenum1,array,500
*dim,nodenum2,array,500
*dim,nod,array,500
*dim,nodez,array,500
*dim,nodey,array,500
*dim,nodex,array,500
*dim,nn,array,500
!定义关键尺寸，用数组表示
*dim,h,array,10
h(1)=0 !模型底部
h(2)=0.9 !土基顶面高度
h(3)=h(2)+0.25 !底基层顶面高度
h(4)=h(3)+0.2 !基层顶面的高度
h(5)=h(4)+0.15 !面层顶面高度
a=12.8 !模型长
b=9.02 !模型宽
jcja=0.64 !检查井尺寸
jcjb=0.66 !检查井尺寸
zc=0.115 !砖的宽度
jgh=0.02 !井盖厚度
d=0.1599 !轮子中心距轮隙中心
l=0.32 !轮载当量矩形长度
w=0.22 !轮载当量矩形宽度
tt=1/15.547 !瞬态分析中采用第2阶的频率
!pressure=52404.38321e3 !单位Pa,2为动荷系数
pressure=0.7e6
!------------------------------------------
!进入前处理模块
!------------------------------------------
/prep7
!第一步：定义单元类型
et,1,solid45
et,2,mesh200,6 !MESH200单元是仅用来划分网格的单元
et,3,surf154,,,,1 !定义表面作用单元
!------------------------------------------
!第二步：定义各种材料属性
!地基
mp,ex,1,45e6
mp,prxy,1,0.35
mp,dens,1,1900
mp,damp,1,0.2 !阻尼系数
!底基层
mp,ex,2,750e6
mp,prxy,2,0.25
mp,dens,2,1750
mp,damp,2,0.01
!基层
mp,ex,3,400e6
mp,prxy,3,0.25
mp,dens,3,1950
mp,damp,3,0.01
!沥青混凝土面层
mp,ex,4,1.2e6
mp,prxy,4,0.25
mp,dens,4,2000
mp,damp,4,0.02
mp,dens,5, !表面单元
!井盖
mp,ex,6,154e9
mp,prxy,6,0.25
mp,dens,6,7300
mp,damp,6,0.03
!井筒砖
mp,ex,7,1200e6
mp,prxy,7,0.2
mp,dens,7,1600
mp,damp,7,0.02
!------------------------------------------
!第三步：建立实体单元，并用工作平面进行分割，以便于生成规则的有限元网格
!------------------------------------------
block,,a,,b,h(1),h(5)
*do,i,2,4,
wpoffs,,,h(i)
vsbw,all
wpcsys,-1,0
*enddo
!切割出检查井来
wpoffs,,b/2-jcjb/2,
wprota,,-90,
vsbw,all
wpcsys,-1,0
wpoffs,,b/2+jcjb/2,
wprota,,-90,
vsbw,all
wpcsys,-1,0
wpoffs,a/2-jcja/2,,
wprota,,,-90
vsbw,all
wpcsys,-1,0
wpoffs,a/2+jcja/2,,
wprota,,,-90
vsbw,all
wpcsys,-1,0
wpoffs,,,h(5)-jgh !用工作面切割出井盖
vsbw,all
wpcsys,-1,0
!选择中间那些体
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,r,loc,z,h(1),h(5)-jgh !选择井盖以下部分
wpoffs,,b/2-jcjb/2+zc, !划分出井筒开始
wprota,,-90,
vsbw,all
wpcsys,-1,0
wpoffs,,b/2+jcjb/2-zc,
wprota,,-90,
vsbw,all
wpcsys,-1,0
wpoffs,a/2-jcja/2+zc,,
wprota,,,-90
vsbw,all
wpcsys,-1,0
wpoffs,a/2+jcja/2-zc,,
wprota,,,-90
vsbw,all
wpcsys,-1,0 !划分出井筒结束
!再选择中间的中间那些体
vsel,r,loc,x,a/2-jcja/2+zc,a/2+jcja/2-zc
vsel,r,loc,y,b/2-jcjb/2+zc,b/2+jcjb/2-zc
VDELE,all,,,1! 删除中间的体
alls
wpoffs,,b/2-jcjb/2+w, !建立轮迹带开始
wprota,,-90,
vsbw,all
wpcsys,-1,0
wpoffs,,b/2-jcjb/2+2*w,
wprota,,-90,
vsbw,all
wpcsys,-1,0 !建立轮迹带结束
alls
save,transient,db
!---------------------------------
!实体模型建立完成，下一步生成有限元网格
!---------------------------------
!************************************
!第四步：生成有限元网格，由面单元生成体单元
!************************************
!赋予材料参数
!1.土基：选择步骤：先选择中间那些体，然后反选，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,invert,loc,z,h(1),h(2)
vsel,r,loc,z,h(1),h(2)
type,1 !选择要划分单元的类型
mat,1 !选择要划分单元的材料模型
!2.底基层：选择步骤：先选择中间那些体，然后反选，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,invert,loc,z,h(1),h(2)
vsel,r,loc,z,h(2),h(3)
type,1 !选择要划分单元的类型
mat,2 !选择要划分单元的材料模型
!3.基层：选择步骤：先选择中间那些体，然后反选，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,invert,loc,z,h(1),h(2)
vsel,r,loc,z,h(3),h(4)
type,1 !选择要划分单元的类型
mat,3 !选择要划分单元的材料模型
!4.面层：选择步骤：先选择中间那些体，然后反选，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,invert,loc,z,h(1),h(2)
vsel,r,loc,z,h(4),h(5)
type,1 !选择要划分单元的类型
mat,4 !选择要划分单元的材料模型
!4.井盖：选择步骤：先选择中间那些体，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,r,loc,z,h(5)-0.0195,h(5)
type,1 !选择要划分单元的类型
mat,6 !选择要划分单元的材料模型
!4.井筒：选择步骤：先选择中间那些体，再选择相应的层
vsel,s,loc,x,a/2-jcja/2,a/2+jcja/2
vsel,r,loc,y,b/2-jcjb/2,b/2+jcjb/2
vsel,r,loc,z,h(1),h(5)-0.02
type,1 !选择要划分单元的类型
mat,7 !选择要划分单元的材料模型
! 划分网格设置线段参数
lsel,s,loc,x,0.1,a/2-jcja/2-0.1 !这里的0.1是选择线的需要
lesize,all,,,19 !分成19份
lsel,s,loc,x,a/2+jcja/2+0.1,a-0.1
lesize,all,,,19 !分成19份
lsel,s,loc,y,0.1,b/2-jcjb/2-0.1 !这里的0.1是选择线的需要
lesize,all,,,19 !分成19份
lsel,s,loc,y,b/2+jcjb/2+0.1,b-0.1
lesize,all,,,19 !分成19份
lsel,s,loc,x,a/2-jcja/2+0.12,a/2+jcja/2-0.12 !这里的0.12是选择线的需要
lesize,all,,,2 !分成2份
lsel,s,loc,z,h(1)+0.12,h(2)-0.12
lesize,all,0.22 !单元边长0.22m
allsel
esize,1 !指定线的缺省划分份数(已直接定义的线，关键点网格划分设置不受影响)
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
vsweep,all !划分网格
allsel
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
numstr,node
numstr,elem
save,transient_mesh,db
!设置表面单元----------------------------
nsel,s,loc,z,h(5)-0.001, h(5)+0.001
nsel,r,loc,y,b/2-jcjb/2+w,b/2-jcjb/2+2*w
type,3
esurf
!--------------------------
!第六步：加约束条件
!--------------------------
nsel,s,loc,z,0 !底部全约束
d,all,all,0
nsel,s,loc,y,0 !前后面y约束
nsel,a,loc,y,b !再选择某些节点加入Active节点中
d,all,uy,0
alls !选择所有元素
!acel,0,0,9.8 !施加重力加速度
alls !选择所有元素
!--------------------------------
!有限元建模完成，保存数据文件
!--------------------------------
save,transient,db
fini !必须项
!---------------------------------
!/////////////////////////////////
!-------------------------------------
!*****Next Step: Solution Module*****
!--------------------------------------
!第七步：求解计算：
/solu
antype,transient
nsubst,10
auto,on
kbc,0
!outres,all,none
!outres,nsol,last
!outres,esol,last
outres,all,last
*do,i,1,40,1
*do,j,1,6
nsel,s,loc,z,h(5)-0.001, h(5)+0.001
nsel,r,loc |
/2 
发表于 2009-1-2 18:16:00
