| 查看: 248 | 回复: 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 |
回复此楼» 猜你喜欢
(出分啦!!!)26 考研调剂 20w+调剂信息火速更新中！【捡漏】qq 群：1087823172
已经有133人回复
-
面上模板改不了页边距吧？
已经有3人回复
-
0856材料求调剂
已经有10人回复
-
高分子化学与物理调剂
已经有5人回复
-
版面费该交吗
已经有18人回复
-
285求调剂
已经有5人回复
-
材料学调剂
已经有3人回复
-
材料类求调剂
已经有4人回复
-
298求调剂
已经有4人回复
-
材料调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2016-06-21 08:48:39
3楼2016-06-27 09:13:46
