反射裂缝路面应力分析求救！！！

moonhu371

我写了一段求路面应力分析的命令流（见下面），运行中总是出现错误“大致意思是硬盘不够” 可是查看下硬盘还有十几G的空间呢，不知道是不是命令流出了问题，请指教。以前曾经算出一次结果，运行了十四五个小时，这次只是去掉了面层的边界约束，约束的改变会对输出结果大小有这么大的影响嘛？

/title,suanli
/units,si

!假设各层连续接触，接缝间无传力效应
!设置各层参数h为厚度、e为弹模、v为波松比、l为长度、，共四层
/prep7
crack=0.01
pave_l=10+crack
pave_w=4
base_l=16+crack
base_w=10

h1=0.05
e1=1.2e9
v1=0.25

h2=0.10
e2=3.0e10
v2=0.15

h3=0.20
e3=3.0e10
v3=0.15

h4=0.30
e4=3.0e10
v4=0.15

h5=9.0
e5=8.0e7
v5=0.30

p_load=7.0e5
load_l=0.189
!定义单元类型及材料属性号
et,1,solid45
et,2,solid64
mp,ex,1,e1
mp,prxy,1,v1
r,1
mp,ex,2,e2
mp,prxy,2,v2
!mp,ex,3,e3
!mp,prxy,3,v3
!mp,ex,4,e4
!mp,prxy,4,v4
!mp,ex,5,e5
!mp,prxy,5,v5

blc4,,,pave_w,pave_l,-h1-h2-h3-h4                            !v1
wpoffs,-(base_w-pave_w)/2,-(base_l-pave_l)/2,-h1-h2-h3-h4
blc4,,,base_w,base_l,-h5                                 !v2
vglue,all
nummrg,all
numcmp,all

wpoffs,(base_w-pave_w)/2,(base_l-crack)/2,h2+h3+h4
vsbw,1
!vsbw,all  
wpoffs,,,-h2
vsbw,3
!vsbw,all                 
blc4,,,pave_w,crack,-h3                                      !
vsbv,1,3
wpoffs,,,-h3
vsbw,all

wpoffs,,,h3+h2+h1
wprota,,-90
vsbw,all
wpoffs,,,crack
vsbw,all
wpoffs,,,load_l
vsbw,all

wprota,,,90
wpoffs,,,load_l
vsbw,all
wpoffs,,,0.131
vsbw,all
wpoffs,,,load_l
vsbw,all
wpoffs,,,1.311
vsbw,all
wpoffs,,,load_l
vsbw,all
wpoffs,,,0.131
vsbw,all
wpoffs,,,load_l
vsbw,all

wpoffs,,,-(4*load_l+1.573)
vsbw,all
wpoffs,,,pave_w
vsbw,all
wpoffs,,,-pave_w
wprota,,,-90
wpoffs,,,(pave_l-crack)/2-load_l
vsbw,all
wpoffs,,,-pave_l
vsbw,all
wprota,,90
allsel
nummrg,all
numcmp,all
vplot
fini

/prep7
!划分网格
!!划线
allsel
lsel,s,loc,x,0.001,load_l-0.001
lsel,a,loc,x,load_l+0.001,load_l+0.130
lsel,a,loc,x,load_l+0.132,2*load_l+0.130
lsel,a,loc,x,2*load_l+1.443,3*load_l+1.441
lsel,a,loc,x,3*load_l+1.443,3*load_l+1.572
lsel,a,loc,x,3*load_l+1.574,4*load_l+1.572
lesize,all,,,2
lsel,a,loc,x,2*load_l+0.132,2*load_l+1.441
lsel,a,loc,x,-2.9,-0.001
lsel,a,loc,x,pave_w+0.001,pave_w+2.99
lesize,all,,,12
lsel,a,loc,x,4*load_l+1.574,pave_w-0.001
lesize,all,,,13

lsel,s,loc,y,pave_l+0.001,pave_l+2.9
lsel,a,loc,y,-2.99,-0.001
lesize,all,,,12
lsel,s,loc,y,5.001,5.009
lesize,all,,,1
lsel,s,loc,y,5.011,5.198
lesize,all,,,2

lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2
lsel,r,loc,x,0.001,pave_w-0.001
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2
lsel,r,loc,x,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2
lsel,r,loc,x,pave_w
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,0.001,pave_w-0.001
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,pave_w
lesize,all,,,25,0.4,,,,1

lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,-(base_w-pave_w)/2
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,(base_w+pave_w)/2
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,0.001,4.99
lsel,r,loc,z,-h1-h2-h3-h4-h5
lesize,all,,,25,2.5,,,,1

lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2
lesize,all,,,25,0.4,,,,1

lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3
lsel,r,loc,x,0.001,pave_w-0.001
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3
lsel,r,loc,x,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3
lsel,r,loc,x,pave_w
lesize,all,,,25,0.4,,,,1

lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3-h4
lsel,u,loc,x,0
lsel,u,loc,x,pave_w
lesize,all,,,25,2.5,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,0
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3-h4
lsel,r,loc,x,pave_w
lesize,all,,,25,0.4,,,,1
lsel,s,loc,y,pave_l/2+crack+load_l+0.001,pave_l-0.001
lsel,r,loc,z,-h1-h2-h3-h4-h5
lesize,all,,,25,2.5,,,,1


lsel,s,loc,z,-0.001,-h1+0.001
lesize,all,,,2
lsel,s,loc,z,-h1-0.001,-h1-h2+0.001
lesize,all,,,3
lsel,s,loc,z,-h1-h2-0.001,-h1-h2-h3+0.001
lesize,all,,,4
lsel,s,loc,z,-h1-h2-h3-0.001,-h1-h2-h3-h4+0.001
lesize,all,,,5
lsel,s,loc,z,-h1-h2-h3-h4+0.001,-h1-h2-h3-h4-h5+0.001
lesize,all,,,25

mp,ex,3,e3
mp,prxy,3,v3
mp,ex,4,e4
mp,prxy,4,v4
mp,ex,5,e5
mp,prxy,5,v5

!!!!!!!!!!
/prep7
!划体
!第一层
allsel
mshape,0,3d
mshkey,1
vsel,s,loc,z,-h1,0
vmesh,all
!第二层
mat,2
type,2
vsel,s,loc,z,-h1-h2,-h1
vmesh,all
!第三层
mat,3
type,1
vsel,s,loc,z,-h1-h2-h3,-h1-h2
vmesh,all
!第四层
mat,4
vsel,s,loc,z,-h1-h2-h3-h4,-h1-h2-h3
vmesh,all
!第五层
mat,5
vsel,s,loc,z,-h1-h2-h3-h4-h5,-h1-h2-h3-h4
vmesh,all
allsel
/rep,fast
fini


!!!!!加载
/solu
!!!!!!!dof
asel,s,loc,z,-h1-h2-h3-h4-h5
da,all,all

asel,s,loc,y,-(base_l-pave_l)/2
asel,a,loc,y,(base_l+pave_l)/2
da,all,uy
asel,s,loc,x,-(base_w-pave_w)/2
asel,a,loc,x,(base_w+pave_w)/2
da,all,ux

allsel


!!!加载
asel,s,loc,z,0
asel,r,loc,x,0,load_l
asel,r,loc,y,pave_l/2+0.01,pave_l/2+load_l
sfa,all,,pres,p_load
asel,s,loc,z,0
asel,r,loc,x,load_l+0.131,2*load_l+0.131
asel,r,loc,y,pave_l/2+0.01,pave_l/2+load_l
sfa,all,,pres,p_load
asel,s,loc,z,0
asel,r,loc,x,2*load_l+1.443,3*load_l+1.441
asel,r,loc,y,pave_l/2+0.01,pave_l/2+load_l
sfa,all,,pres,p_load
asel,s,loc,z,0
asel,r,loc,x,3*load_l+1.574,4*load_l+1.572
asel,r,loc,y,pave_l/2+0.01,pave_l/2+load_l
sfa,all,,pres,p_load
allsel
/rep,fast
fini

!!求解
/solu
antype,0
nsubst,30
outres,all,all
save
solve

隐身

可能是你的硬盘格式的问题
  fat32格式允许的文件容量好像是4个G左右
而ntfs格式就能够容纳更多的容量

xiaoxi_321

楼上说的不错，你试试

陈一彬

楼主好强大~我是初学者开始做路面分析~~