| 查看: 287 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
两块石头
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-21 08:48:39
3楼2016-06-27 09:13:46
