| 查看: 275 | 回复: 2 | ||
[求助]
求大神们帮我改一下程序，用的是生死单元。是一个多孔材料的压缩，随机孔
|
|
fini /clear /filname,infil_pre /prep7 et,1,solid185 ET,2,TARGE170 ET,3,CONTA174,,,,,1 mp,ex,1,2.204042647e10 mp,nuxy,1,0.34 mp,mu,1,0.10 mp,mu,3,0.10 TB,MISO,1,1,43 TBTEMP,25 tbpt,,2.31482E-06,51019.62 tbpt,,4.62963E-06,51019.62 tbpt,,1.38889E-05,51019.62 tbpt,,7.63889E-05,51019.62 tbpt,,0.000150463,51019.62 tbpt,,0.000252315,51215.84 tbpt,,0.000326389,52589.45 tbpt,,0.000402778,135005.8 tbpt,,0.000532407,147368.2 tbpt,,0.004493056,1572778 tbpt,,0.006789352,3092573 tbpt,,0.00886574,4934774 tbpt,,0.01103241,7460834 tbpt,,0.01221759,9172345 tbpt,,0.0139375,12042980 tbpt,,0.0154213,14788230 tbpt,,0.01705093,17996190 tbpt,,0.01966204,22702750 tbpt,,0.02328935,29231100 tbpt,,0.02649537,35805370 tbpt,,0.03384491,46561090 tbpt,,0.04443982,55628650 tbpt,,0.05241898,60405070 tbpt,,0.07268055,69181420 tbpt,,0.1044213,78260160 tbpt,,0.1618287,89023140 tbpt,,0.2191782,100588900 tbpt,,0.4472755,209651000 tbpt,,0.4778287,240189200 tbpt,,0.5093819,279676000 tbpt,,0.5367338,321502600 tbpt,,0.5617916,367557500 tbpt,,0.582382,411481000 tbpt,,0.6029213,461667300 tbpt,,0.617088,500828800 tbpt,,0.6334027,550832800 tbpt,,0.6538842,622046300 tbpt,,0.6635301,659255200 tbpt,,0.684956,750819500 tbpt,,0.7010486,828383100 tbpt,,0.7141621,898207800 tbpt,,0.7219189,942393600 tbpt,,0.7242847,956303800 ！这里的点太多，大家随意啦 R,1,,,0.1 block,0,0.005,0,0.005,0,0.005 wpoff,0.005/2,0.005,0.005/2 wprota,0,90 rectng,-0.005,0.005,-0.005,0.005 lesize,11,,,20 lesize,6,,,20 lesize,7,,,20 mat,1 mshape,0 mshkey,1 vmesh,1 lesize,13,,,5 lesize,14,,,5 type,2 mshape,0 mshkey,1 amesh,7 *set,N_ele,8000 *do,i,1,N_ele,1 porosity=0.54 r_any=rand(0,1) *if,r_any,gt,porosity,then ekill,i *else *endif *enddo esel,s,live nsle,s nsel,s,loc,y,0.005 type,3 r,1 esurf,all alls finish /solu esel,s,live nsle,s nsel,inve d,all,all alls !不在任何活单元上的节点完全约束 nsel,s,loc,y,0 d,all,uy,0 nsel,s,loc,x,0.0025 d,all,ux,0 nsel,s,loc,z,0.0025 d,all,uz,0 alls antype,0 nlgeom,1 nropt,full autots,0 time,100 solve finish /solu nsubst,3000,1000000,10 autots,1 outres,all,1 time,500 kbc,0 ！加载类型 asel,s,area,7 //这句貌似是有问题，反正就是在刚性压板上加载 nsla,s d,all,vely,-0.0035/500 alls solve 这里还没有考虑后处理，所以没有将生死单元分开。请大神帮我调调程序，所有的参数都在上面了@wuming524 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-21 08:48:39
3楼2016-06-27 09:13:46
