| 查看: 1020 | 回复: 2 | |||
11050209银虫 (正式写手)
|
[求助]
生死单元求助
|
|
生死单元求助！本人在做一个连接件的拉伸有限元分析，想在模拟的过程中将超过材料抗拉强度的单元杀死，可是现在遇到问题了，只能计算一个子步，感觉杀死单元这块出了问题，但又不知如何解决，特来求助。(不杀死单元可以运行通过) /filname,sbolt-16-hb6319-bu /prep7 /TITLE, lashen !定义单元类型及材料属性 et,1,solid185 et,2,plane182 mptemp,1,25,100,200,300,400,500 mpdata,ex,1,1,114e9,110e9,105e9,99.9e9,95e9,90.1e9 mpdata,prxy,1,1,0.32,0.33,0.33,0.34,0.35,0.35 mpdata,alpx,1,1,7e-6,7.71e-6,8.62e-6,9.25e-6,9.91e-6,10.2e-6 MPTEMP,,,,,,,, MPTEMP,1,25 MPTEMP,2,400 MPTEMP,3,500 MPTEMP,4,600 MPDATA,EX,2,,1.9E+011 MPDATA,EX,2,,1.52E+011 MPDATA,EX,2,,1.35E+011 MPDATA,EX,2,,1.31E+011 MPDATA,PRXY,2,,0.3 MPDATA,PRXY,2,,0.3 MPDATA,PRXY,2,,0.3 MPDATA,PRXY,2,,0.3 MPTEMP,,,,,,,, MPTEMP,1,100 MPTEMP,2,200 MPTEMP,3,300 MPTEMP,4,400 MPTEMP,5,500 MPTEMP,6,600 MPDATA,ALPX,2,,15.37e-6 MPDATA,ALPX,2,,16.09e-6 MPDATA,ALPX,2,,16.31e-6 MPDATA,ALPX,2,,16.84e-6 MPDATA,ALPX,2,,17.58e-6 MPDATA,ALPX,2,,18.06e-6 MPTEMP,,,,,,,, MPTEMP,1,25 MPTEMP,2,100 MPTEMP,3,150 MPTEMP,4,175 MPDATA,ex,3,,69e9 MPDATA,eX,3,,63e9 MPDATA,eX,3,,60e9 MPDATA,eX,3,,59e9 MPDATA,prxy,3,,0.33 MPDATA,prxy,3,,0.33 MPDATA,prxy,3,,0.33 MPDATA,prxy,3,,0.33 MPTEMP,,,,,,,, MPTEMP,1,100 MPTEMP,2,200 MPTEMP,3,300 MPDATA,alpx,3,,23e-6 MPDATA,alpx,3,,24.3e-6 MPDATA,alpx,3,,25.2e-6 TB,BKIN,1,5 TBTEMP,25 TBDATA,,1.07E9,11.4E9 TBTEMP,100 TBDATA,,0.99E9,11E9 TBTEMP,200 TBDATA,,0.89E9,10.5E9 TBTEMP,300 TBDATA,,0.822E9,9.99E9 TBTEMP,500 TBDATA,,0.715E9,9.01E9 TB,BKIN,2,1 TBTEMP,25 TBDATA,,1.1E9,19E9 TB,BKIN,3,4 TBTEMP,25 TBDATA,,530E6,6.9E9 TBTEMP,100 TBDATA,,450E6,6.3E9 TBTEMP,125 TBDATA,,420E6,6.0E9 TBTEMP,175 TBDATA,,350E6,5.9E9 !定义参数 l1=80/1000 !试样长度 l2=34/1000 !夹头长度 w1=16/1000 !试样宽度 w2=30/1000 !夹头宽度 h=3/1000 !试样厚度 l3=16/1000 !试样交叠部分长度 r=7/1000 !连接圆弧半径 r1=4/2/1000 !圆台下表面半径 r2=8.75/2/1000 !圆台上表面半径 r3=4/1000 !螺母半径 pload1=5000 pload2=-pload1 !建立几何模型(TA32) block,-l3/2,0,0,w1/2,0,h block,0,l3/2,0,w1/2,0,h block,-l3/2,0,-w1/2,0,0,h block,0,l3/2,0,-w1/2,0,h cone,r1,r2,0,h vsel,s,,,1,4,1 allsel vsbv,all,5 numcmp,all block,-31/1000,-l3/2,-w1/2,w1/2,0,h k,41,-38/1000,(w1/2+r),h k,42,-72/1000,(w1/2+r),h k,43,-72/1000,-(w1/2+r),h k,44,-38/1000,-(w1/2+r),h larc,40,41,12,r l,41,42 l,42,43 l,43,44 larc,44,37,22,r al,64,69,70,71,72,73 vdrag,35,,,,,,65 block,-72/1000,-38/1000,-w2/2,w2/2,-h,0 a,44,47,50,41 vsba,6,49 allsel vglue,all numcmp,all !建立几何模型(7050) block,-l3/2,0,0,w1/2,-h,0 block,0,l3/2,0,w1/2,-h,0 block,-l3/2,0,-w1/2,0,-h,0 block,0,l3/2,-w1/2,0,-h,0 cylind,0,r1,0,-h vsel,s,,,9,12,1 allsel vsbv,all,13 block,l3/2,31/1000,-w1/2,w1/2,-h,0 k,91,38/1000,(w1/2+r),0 k,92,72/1000,(w1/2+r),0 k,93,72/1000,-(w1/2+r),0 k,94,38/1000,-(w1/2+r),0 larc,65,91,56,r l,91,92 l,92,93 l,93,94 larc,94,63,69,r al,89,106,107,111,113,114 vdrag,53,,,,,,97 block,38/1000,72/1000,-w2/2,w2/2,0,h a,87,88,86,83 vsba,10,93 numcmp,all allsel vsel,s,,,9,16,1 vglue,all allsel !建立几何模型(1s627) cone,r1,r2,0,h cylind,0,r1,0,-(h+h) cylind,r1,r3,-h,-(h+h) vsel,s,,,11,12,1 vsel,a,,,9 vglue,all allsel !划分网格(TA32) MSHAPE,0,2D MSHKEY,1 mat,1 lsel,s,,,3,4,1 lesize,all,,,w1/2*1000 lsel,s,,,28,30,2 lesize,all,,,4 allsel AMAP,9,4,23,22,5 lsel,s,,,9,72,63 lesize,all,,,w1/2*1000 AMAP,38,5,22,21,9 lsel,s,,,74,76,2 lesize,all,,,w1/2*1000 AMAP,44,9,21,24,14 lsel,s,,,13,66,53 lesize,all,,,w1/2*1000 AMAP,40,14,24,23,4 vsel,s,,,6,8,1 vsel,a,,,1 vsweep,all allsel asel,s,,,9,38,29 asel,a,,,40,44,4 aclear,all lsel,s,,,68,70,2 lesize,all,,,(l1-l3-l2-r)*1000 vsweep,5 lsel,s,,,41 lesize,all,,,3,2 lsel,s,,,45 lesize,all,,,3,0.5 lsel,s,,,51 lesize,all,,,w1*1000 vsweep,2 lsel,s,,,42,44,2 lesize,all,,,17 vsweep,4 vsweep,3 allsel !划分网格(7050) mat,3 lsel,s,,,104,188,84 lesize,all,,,w1/2*1000 lsel,s,,,135,144,9 lesize,all,,,4 allsel AMAP,109,42,75,73,60 lsel,s,,,184,192,8 lesize,all,,,w1/2*1000 AMAP,115,60,73,74,52 lsel,s,,,186,190,4 lesize,all,,,w1/2*1000 AMAP,112,52,74,76,43 lsel,s,,,77,78,1 lesize,all,,,w1/2*1000 AMAP,74,43,76,75,42 vsel,s,,,20,22,1 vsel,a,,,13 vsweep,all allsel lsel,s,,,80,82,2 lesize,all,,,(l1-l3-l2-r)*1000 vsweep,19 lsel,s,,,178,179,1 lesize,all,,,3,2 lsel,s,,,162 lesize,all,,,w1*1000 vsweep,17 lsel,s,,,166,168,2 lesize,all,,,17 vsweep,10 vsweep,18 asel,s,,,109,115,6 asel,a,,,74,112,38 aclear,all allsel !划分网格(1s627) mat,2 lsel,s,,,94,95,1 lsel,a,,,99,100,1 lesize,all,,,w1*1000 vsweep,9 allsel lsel,s,,,107,108,1 lsel,a,,,111,112,1 lesize,all,,,w1*1000 lsel,s,,,117,118,1 lesize,all,,,6 vsweep,14 lsel,s,,,147,154,1 lesize,all,,,w1*1000 vsweep,15 !建立接触对 mp,mu,4,0.3 mat,4 r,4 r, 4, 0, 0, 0.3, 0.1, 0, 0 ! contact pair 4 real,4 et,3,targe170 et,4,conta173 !generate the target surface asel,s,,,51,53,2 cm,up_bolt,area cmsel,s,up_bolt type,3 nsla,s,1 esln,s,0 esurf,all allsel !generate the contact surface asel,s,,,7,11,4 asel,a,,,13,14,1 cm,up_plate,area cmsel,s,up_plate type,4 nsla,s,1 esln,s,0 esurf,all allsel !generate the target surface mat,4 r,5 r, 5, 0, 0, 0.3, 0.1, 0, 0 ! contact pair 5 real,5 et,5,targe170 et,6,conta173 asel,s,,,77,79,2 cm,down_bolt,area cmsel,s,down_bolt type,5 nsla,s,1 esln,s,0 esurf,all allsel !generate the contact surface asel,s,,,73,78,5 asel,a,,,83,88,5 cm,down_plate,area cmsel,s,down_plate,area type,6 nsla,s,1 esln,s,0 esurf,all allsel !generate the target surface mat,4 r,6 real,6 r, 6, 0, 0, 0.3, 0.1, 0, 0 ! contact pair 6 et,7,targe170 et,8,conta173 asel,s,,,60 cm,botom_bolt,area cmsel,s,botom_bolt type,7 nsla,s,1 esln,s,0 esurf,all allsel !generate the contact surface asel,s,,,74,109,35 asel,a,,,112,115,3 cm,botom_plate,area cmsel,s,botom_plate type,8 nsla,s,1 esln,s,0 esurf,all allsel !generate the target surface mat,4 r,7 real,7 r, 7, 0, 0, 0.3, 0.1, 0, 0 ! contact pair 7 et,7,targe170 et,8,conta173 asel,s,,,75,110,35 asel,a,,,113,116,3 cm,botom_plate,area cmsel,s,botom_plate type,7 nsla,s,1 esln,s,0 esurf,all allsel !generate the contact surface asel,s,,,39,43,4 asel,a,,,8,37,29 cm,up_plate,area cmsel,s,up_plate type,8 nsla,s,1 esln,s,0 esurf,all allsel /sol antype,static nlgeom,on lnsrcch,on pred,on outres,all,all cnvtol,f,,0.05 nsel,s,loc,x,(l1-l3+l3/2),(l1-l3+l3/2) cp,1,ux,all allsel nsel,s,loc,x,(l1-l3+l3/2),(l1-l3+l3/2) nsel,r,loc,y,0 nsel,r,loc,z,0 f,all,fx,pload1 allsel nsel,s,loc,x,-(l1-l3+l3/2),-(l1-l3+l3/2) d,all,ux,0 d,all,uy,0 d,all,uz,0 allsel bf,all,temp,125 *do,i,1,10 time,i*1/10 autots,on solve allsel ESEL,s,MAT,,1 !将1号材料超过其抗拉强度的单元杀死 nsle,s nsel,r, s,eqv,1100e6,3000e13 esln,s ekill,all allsel ESEL,s,MAT,,3 nsle,s nsel,r, s,eqv,530e6,3000e13 !将3号材料超过其抗拉强度的单元杀死 esln,s ekill,all allsel *enddo finish |
回复此楼» 猜你喜欢
求助一下有机合成大神
已经有3人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有13人回复
-
存款400万可以在学校里躺平吗
已经有30人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2017-05-23 23:30:24
3楼2017-06-23 11:14:53