| 查看: 244 | 回复: 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 |
回复此楼» 猜你喜欢
招博士
已经有4人回复
-
救命帖
已经有7人回复
-
限项规定
已经有6人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
最失望的一年
已经有18人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-21 08:48:39
3楼2016-06-27 09:13:46