| 查看: 667 | 回复: 1 | ||
[求助]
我要做一个接触分析的，算不出结果，求高手帮我看看
|
|
我写了两段命令，一个算不出结果，另一个不收敛，求高手帮我看看。下面是命令： h=6 !.....结构层厚度 h1=0.30 !.....面层厚度 h2=0.30 !.....级配碎石层厚度 h3=h-h1-h2 !.....土基厚度 b=0.189 !..........轮胎尺寸 d=0.32 !..........双轮间距 w=5 !..........模型宽度 FINISH /prep7 *afun,deg et,1,plane42 blc4,0,0,w,h3,0 blc4,0,h3,w,h2,0 blc4,0,h2+h3,w,h1,0 wprota,0,0,90 wpave,w/2-d/2-b/2,h,0 asbw,all wpave,w/2-d/2+b/2,h,0 asbw,all wpave,w/2+d/2-b/2,h,0 asbw,all wpave,w/2+d/2+b/2,h,0 asbw,all asel,s,loc,y,0,h3+h2-0.0001 aglue,all asel,s,loc,y,h2+h3+0.0001,h aglue,all MP,ex,1,1200e6 mp,prxy,1,0.35 mp,ex,2,300e6 mp,prxy,2,0.3 tb,dp,2 TBDATA,1,54000 TBDATA,2,47 mp,ex,4,40e6 mp,prxy,4,0.35 asel,s,loc,y,0,h3 AATT,4,,,,0 asel,s,loc,y,h3,h3+h2 AATT,2,,,,0 asel,s,loc,y,h3+h2,h AATT,1,,,,0 lsel,s,length,,h1 lesize,all,,,5 lsel,s,length,,h2 LESIZE,all,,,5 lsel,s,length,,h3 lesize,all,,,50 lsel,s,length,,b lesize,all,,,3 lsel,s,length,,w/2-d/2-b/2 lesize,all,,,20 lsel,s,length,,d-b lesize,all,,,2 allsel,all amesh,all /COM, CONTACT PAIR CREATION - START MP,MU,1,0.5 MAT,1 !MP,EMIS,1,7.88860905221e-031 R,3 REAL,3 ET,2,169 ET,3,172 R,3,,,1.0,0.1,0, RMORE,,,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,,1.0,0.5 RMORE,0,1.0,1.0,0.0,,1.0 RMORE,10.0 KEYOPT,3,3,0 KEYOPT,3,4,2 KEYOPT,3,5,0 KEYOPT,3,7,0 KEYOPT,3,8,0 KEYOPT,3,9,0 KEYOPT,3,10,2 KEYOPT,3,11,0 KEYOPT,3,12,6 KEYOPT,3,2,2 ! 生成目标面 asel,s,mat,,2 lsla,r lsel,r,loc,y,h2+h3 CM,_TARGET,LINE TYPE,2 NSLL,S,1 ESLN,S,0 ESURF CMSEL,S,_ELEMCM ! 生成接触面 lsel,all asel,s,mat,,1 lsla,r lsel,r,loc,y,h2+h3 CM,_CONTACT,LINE TYPE,3 NSLL, ESLN,S,0 ESURF,all ALLSEL,all ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /PSYMB,ESYS,1 /PNUM,TYPE,1 /NUM,1 EPLOT ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /COM, CONTACT PAIR CREATION - END /MREP,EPLOT /solu !ANTYPE,0 !NLGEOM,ON !NROPT,FULL,,off !pred,on !AUTOTS,1 !TIME,1 !NSUBST,50,,,1 !TIME,0.1 !NSUBST,10,20,5 antype,0 nlgeom,on autots,on lsel,s,loc,y,0 d,all,all,0 allsel,all lsel,s,loc,x,0 lsel,a,loc,x,w d,all,,ux,0 allsel,all asel,s,mat,,1 lsla,r lsel,r,loc,y,h lsel,r,length,,b SFL,all,pres,0.7e6 allsel,all solve 一上命令算不出结果，好像是荷载没加上一样。 第二个 h=6 !.....结构层厚度 h1=0.30 !.....面层厚度 h2=0.30 !.....级配碎石层厚度 h3= 0.30 !.....半刚性基层厚度 h4=h-h1-h2-h3 !.....土基厚度 b=0.189 !..........轮胎尺寸 d=0.32 !..........双轮间距 w=5 !..........模型宽度 FINISH /prep7 *afun,deg et,1,plane42 blc4,0,0,w,h4,0 blc4,0,h4,w,h3,0 blc4,0,h4+h3,w,h2,0 blc4,0,h4+h3+h2,w,h1,0 wprota,0,0,90 wpave,w/2-d/2-b/2,h,0 asbw,all wpave,w/2-d/2+b/2,h,0 asbw,all wpave,w/2+d/2-b/2,h,0 asbw,all wpave,w/2+d/2+b/2,h,0 asbw,all asel,s,loc,y,0,h4+h3+h2-0.0001 aglue,all asel,s,loc,y,h2+h3+h4+0.0001,h aglue,all MP,ex,1,1200e6 mp,prxy,1,0.35 mp,ex,2,300e6 mp,prxy,2,0.3 tb,dp,2 TBDATA,1,5400 TBDATA,2,47 mp,ex,3,1500e6 mp,prxy,3,0.25 mp,ex,4,40e6 mp,prxy,4,0.35 asel,s,loc,y,0,h4 AATT,4,,,,0 asel,s,loc,y,h4,h4+h3 AATT,3,,,,0 asel,s,loc,y,h4+h3,h4+h3+h2 AATT,2,,,,0 asel,s,loc,y,h4+h3+h2,h AATT,1,,,,0 lsel,s,length,,h1 lesize,all,,,5 lsel,s,length,,h2 LESIZE,all,,,5 lsel,s,length,,h3 lesize,all,,,5 lsel,s,length,,h4 lesize,all,,,50 lsel,s,length,,b lesize,all,,,3 lsel,s,length,,w/2-d/2-b/2 lesize,all,,,20 lsel,s,length,,d-b lesize,all,,,2 allsel,all amesh,all /COM, CONTACT PAIR CREATION - START MP,MU,1,0.5 MAT,1 MP,EMIS,1,7.88860905221e-031 R,3 REAL,3 ET,2,169 ET,3,172 R,3,,,1.0,0.1,, NROPT,UNSYM RMORE,,,1.0E20,0.0,1.0, RMORE,0.0,0,1.0,,1.0,0.5 RMORE,0,1.0,1.0,0.0,,1.0 RMORE,10.0 KEYOPT,3,3,0 KEYOPT,3,4,2 KEYOPT,3,5,0 KEYOPT,3,7,0 KEYOPT,3,8,0 KEYOPT,3,9,0 KEYOPT,3,10,2 KEYOPT,3,11,0 KEYOPT,3,12,6 KEYOPT,3,2,2 ! 生成目标面 asel,s,mat,,2 lsla,r lsel,r,loc,y,h2+h3+h4 CM,_TARGET,LINE TYPE,2 NSLL,S,1 ESLN,S,0 ESURF,all CMSEL,S,_ELEMCM ! 生成接触面 lsel,all asel,s,mat,,1 lsla,r lsel,r,loc,y,h2+h3+h4 CM,_CONTACT,LINE TYPE,3 NSLL, ESLN,S,0 ESURF,all ALLSEL,all ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /PSYMB,ESYS,1 /PNUM,TYPE,1 /NUM,1 EPLOT ESEL,ALL ESEL,S,TYPE,,2 ESEL,A,TYPE,,3 ESEL,R,REAL,,3 /COM, CONTACT PAIR CREATION - END /MREP,EPLOT /solu !ANTYPE,0 !NLGEOM,1 !TIME,1 !AUTOTS,1 !NSUBST,20, , ,1 !TIME,0.1 !NSUBST,10,20,5 antype,0 nlgeom,on autots,on nsel,s,loc,y,0 d,all,all,0 allsel,all nsel,s,loc,x,0 nsel,a,loc,x,w d,all,ux,0 allsel,all asel,s,mat,,1 lsla,r lsel,r,loc,y,h lsel,r,length,,b SFL,all,pres,0.7e6 solve 这段命令好像不收敛，我对比两段命令，找不出错误，求高手帮我看看，谢谢了。 我是新手，没有太多金币，谢谢了。 |
回复此楼» 猜你喜欢
298求调剂
已经有14人回复
-
成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。
已经有8人回复
-
课题组接收材料类调剂研究生
已经有6人回复
-
338求调剂
已经有5人回复
-
11408，学硕276求调剂
已经有5人回复
-
材料化工调剂
已经有16人回复
-
267求调剂
已经有3人回复
-
调剂
已经有6人回复
-
290分材料工程085601求调剂 数二英一
已经有9人回复
-
085600 英一数二272求调剂
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求高手帮看看我的笔记本问题出在哪？
已经有6人回复
- 求高手帮忙解析这个氢谱，看看有没有必要做升温核磁？
已经有12人回复
- 高手帮我看看异构体NMR如何解析
已经有8人回复
- DSC作图，做不出来了，求哪位高手帮忙做一下图，非常感谢~
已经有7人回复
- 大家来帮我看看这两个INCAR 文件有什么区别哦。。为啥计算结果不一样呢
已经有4人回复
- 求高手帮我看看DGGE图片，为什么跑不开，非常感谢
已经有32人回复
- 跟各位高手交流一下：大家是什么时候开始接触量化的？
已经有25人回复
- （BiFeO3）GGA和LDA算出来的V0都比实验值小约10~20个立方埃？请各位高手帮我检查一下
已经有27人回复
- 求高手帮忙分析一下NMR和IR谱图，第一次做，无从下手，哪怕进来给点建议也好?
已经有43人回复
- 【求助】那位高手有没有接触过：二氧化钛掺杂元素，使其禁带宽度变宽的
已经有13人回复
- 【求助】求分析一个XRD物相组成，请高手帮忙，谢谢
已经有8人回复
- 【求助】最近做机器人运动仿真，希望高手帮我画出它的第3关节运动轨迹！
已经有13人回复
- 【求助】求助高手帮忙看看这张怪异的谱图
已经有10人回复
hudie012
荣誉版主 (文坛精英)
小兔爱吃萝卜
- 应助: 123 (高中生)
- 贵宾: 11.438
- 金币: 18381.5
- 散金: 29063
- 红花: 143
- 沙发: 43
- 帖子: 12339
- 在线: 1962.1小时
- 虫号: 1296312
- 注册: 2011-05-15
- 性别: GG
- 专业: 工程热物理相关交叉领域
- 管辖: 休闲灌水
2楼2012-05-05 12:52:56
10