现在施加模型边界条件遇到了问题 香皂模型边界同时施加 对流边界 和热流边界,便采用表面效应单元,想将对流条件施加在实体模型上,热流边界施加在表面效应单元上,不知道自己的命令流是否正确?总感觉计算出来的结果 和单独施加对流边界差不多···
!!!!路面结构几何尺寸!!!!
LONG=20
WIDTH=18
CDK=3.75 !!!车道宽!!!
D=0.2 !!!轮迹宽,正方形!!!
ZJ=1.8 !!!车轮轴距!!!
LJ=CDK/2-ZJ/2-(5/4)*D !!!第一条轮迹距路面边缘的距离!!!
LJ2=LONG/2-D/2 !!!轮迹距模型行车方向的距离!!!
H1=0.18 !!!!沥青混凝土面层!!!!
H2=0.20 !!!!水泥稳定砂砾基层!!!!
H3=0.3 !!!!二灰土底基层!!!!
H4=1.0 !!!!土基!!!!
!!!划分数设置!!!
ZHF1=8
ZHF2=2
ZHF3=1
ZHF4=2
ZHF5=12
ZHF6=2
ZHF7=1
ZHF8=2
ZHF9=40
XHF1=20
XHF2=2
TUF=10
LUJF=5
DJCF=10
JCF=10
XMCF=2
ZMCF=2
SMCF=2
!!!!路面各层材料属性!!!!
!!!!上面层材料属性!!!!
C1=900
KXX1=1.0
DENS1=2200
!!!!中面层材料属性!!!!
C2=890
KXX2=1.2
DENS2=2200
!!!!下面层材料属性!!!!
C3=920
KXX3=1.1
DENS3=1800
!!!!基层材料属性!!!!
C4=860
KXX4=1.2
DENS4=1700
/FILNAME,LUJIEJIEGOUwenduchang
/TITLE,WANG YU LIAN XI
/config,fsplit,750
/PREP7
!!!!单元类型,材料属性!!!!
ET,1,SOLID70
ET,2,SURF152
KEYOPT,2,8,1
*DIM,WY2,TABLE,25,1,1
WY2(1,0)=0,3600,7200,10800,14400,18000,21600,25200,28800,32400
WY2(11,0)=36000,39600,43200,46800,50400,54000,57600,61200,64800,68400
WY2(21,0)=72000,75600,79200,82800,86400
WY2(1,1)=0,0,0,0,0,0,0,0,3.353785127,11.02772282
WY2(11,1)=17.53647244,21.73535221,23.17944591,21.73535221,17.53647244,11.02772282,3.353785127,0,0,0
WY2(21,1)=0,0,0,0,0
MP,DENS,1,DENS1
MP,KXX,1,KXX1
MP,C,1,C1
MP,DENS,2,DENS2
MP,KXX,2,KXX2
MP,C,2,C2
MP,DENS,3,DENS3
MP,KXX,3,KXX3
MP,C,3,C3
MP,DENS,4,DENS4
MP,KXX,4,KXX4
MP,C,4,C4
!!!!实体建模!!!!
BLOCK,0,LONG,0,H4,0,WIDTH
BLOCK,0,LONG,H4,H4+H3,0,WIDTH
BLOCK,0,LONG,H4+H3,H4+H3+H2,0,WIDTH
BLOCK,0,LONG,H4+H3+H2,H4+H3+H2+H1,0,WIDTH
!!!布尔粘结元素,完全连续!!!
VSEL,ALL
VGlUE,ALL
!!!压缩编号!!!
NUMCMP,VOLU
NUMCMP,AREA
NUMCMP,LINE
!!!为各层赋予属性!!!
!!!土基!!!
VSEL,S,LOC,Y,0,H4
VPLOT,ALL
VATT,4,,1
!!!路基!!!
VSEL,S,LOC,Y,H4,H4+H3
VPLOT,ALL
VATT,3,,1
!!!底基层!!!
VSEL,S,LOC,Y,H4+H3,H4+H3+H2
VPLOT,ALL
VATT,2,,1
!!!基层!!!
VSEL,S,LOC,Y,H4+H3+H2,H4+H3+H2+H1
VPLOT,ALL
VATT,1,,1
!!划分网格!!
LSEL,S,LOC,Y,0,H4
LSEL,U,LOC,Y,0
LSEL,U,LOC,Y,H4
LPLOT,ALL
LESIZE,ALL,,,5
LSEL,S,LOC,Y,H4,H4+H3
LSEL,U,LOC,Y,H4+H3
LSEL,U,LOC,Y,H4
LPLOT,ALL
LESIZE,ALL,,,3
LSEL,S,LOC,Y,H4+H3+H2,H4+H3
LSEL,U,LOC,Y,H4+H3
LSEL,U,LOC,Y,H4+H3+H2
LPLOT,ALL
LESIZE,ALL,,,2
LSEL,S,LOC,Y,H4+H3+H2+H1,H4+H3+H2
LSEL,U,LOC,Y,H4+H3+H2+H1
LSEL,U,LOC,Y,H4+H3+H2
LPLOT,ALL
LESIZE,ALL,,,2
ALLSEL,ALL
VSEL,S,LOC,Y,0,H4+H3+H2+H1
VPLOT,ALL
ESIZE,1
VMESH,ALL
!!建立表面效应单元!!
ASEL,S,LOC,Y,H4+H3+H2+H1
APLOT,ALL
NSLA,S,1
TYPE,2
ESURF,ALL
ALLSEL,ALL
SAVE
FINISH
!!!施加荷载!!!
!!!底面固定约束!!!
!/config,nproc,2!
/SOLU
ANTYPE,TRANS
TIMINT,OFF
TIME,0.01
DELTIM,0.01
ASEL,S,LOC,Y,H4+H3+H2+H1
APLOT,ALL
NSLA,S,1
D,ALL,TEMP,-8.2
ASEL,S,LOC,Y,0
APLOT,ALL
NSLA,S,1
D,ALL,TEMP,15
ALLSEL,ALL
SOLVE
EQSLV,JCG,1.0E-5
TIMINT,ON
TIME,24*3600
DELTIM,3600,3600,3600
AUTOTS,ON
DDELE,ALL,TEMP
ESEL,S,TYPE,,1
EPLOT
NSLE,S
NSEL,R,LOC,Y,H4+H3+H2+H1
NPLOT,ALL
SF,ALL,CONV,27.3,%WY1%
ESEL,S,TYPE,,2
EPLOT
NSLE,S
SF,ALL,HFLUX,%WY2%
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE
SAVE
FINISH
/POST26
NSOL,2,4669,TEMP,,T1 !!路表面!!
NSOL,3,3871,TEMP,,T2 !!面层底面!!
NSOL,4,2716,TEMP,,T3 !!基层底面!!
NSOL,5,805,TEMP,,T4 !!底基层底面!!
NSOL,6,448,TEMP,,T5 !!土基底面!!
PLVAR,2,3,4,5,6 |