| 查看: 288 | 回复: 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 |
回复此楼» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有8人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
Ei源刊怎么投
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
3楼2016-06-27 09:13:46
2楼2016-06-21 08:48:39