| 查看: 1947 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
焊接应力场仿真不收敛问题 已有2人参与
|
|||
|
大神们好，刚开始做焊接温度场和应力场仿真，各种问题不断。 下面是我应力场的APDL（刚开始，写的不好请见谅），麻烦大神们看下有什么问题？ 我在计算的时候，第1个载荷步就不收敛了，不知道问题所在。 金币也没多少，还望不吝赐教！谢谢！ /prep7 /UNITS,SI et,1,solid70 et,2,plane55 !材料参数 mptemp,1,0,100,200,400,600,800 !MPTEMP, STLOC, T1, T2, T3, T4, T5, T6 mptemp,7,1000,1200,1400,1550,1800,2000 mpdata,kxx,1,1,50,45,40,36,34,25 !MPDATA, Lab, MAT, STLOC, C1, C2, C3, C4, C5, C6 mpdata,kxx,1,7,26,28,57,67,215,280 mpdata,c,1,1,480,500,520,650,750,1000 mpdata,c,1,7,1200,1400,1600,1700,1770,1800 mpdata,dens,1,1,7800,7800,7800,7760,7760,7520 mpdata,dens,1,7,7390,7300,7250,7180,7180,7180 mpdata,ex,1,1,2.1e11,2e11,2e11,1.7e11,0.8e11,0.35e11 mpdata,ex,1,7,0.2e11,0.15e11,0.1e11,0.1e11,0.1e11,0.1e11 mpdata,prxy,1,1,0.27,0.27,0.28,0.28,0.29,0.3 mpdata,prxy,1,7,0.31,0.32,0.32,0.32,0.32,0.32 mpdata,alpx,1,1,1.15e-5,1.2e-5,1.3e-5,1.42e-5,1.45e-5,1.45e-5 mpdata,alpx,1,7,1.45e-5,1.45e-5,1.45e-5,1.45e-5,1.45e-5,1.45e-5 tb,biso,1,6 tbtemp,20 tbdata,,2.2e8,2.05e10 tbtemp,500 tbdata,,1.3e8,1.5e10 tbtemp,750 tbdata,,4e7,7e9 tbtemp,1000 tbdata,,2.5e7,2e9 tbtemp,1700 tbdata,,1e6,1.8e9 tbtemp,2500 tbdata,,1e5,1.2e9 !建立几何模型 k,1,-0.06,0,0 k,2,0.06,0,0 k,3,0.06,0.01,0 k,4,-0.06,0.01,0 k,5,-0.003,0.01,0 k,6,0.003,0.01,0 k,7,0.003,0.13,0 k,8,-0.003,0.13,0 k,9,-0.06,0.13,0 k,10,0.06,0.13,0 k,11,0.06,0.14,0 k,12,-0.06,0.14,0 k,13,-0.0109,0.01,0 k,14,-0.003,0.0179,0 k,15,0.0109,0.01,0 k,16,0.003,0.0179,0 k,17,-0.0109,0.13,0 k,18,-0.003,0.1221,0 k,19,0.0109,0.13,0 k,20,0.003,0.1221,0 k,21,-0.0109,0,0 k,22,-0.003,0,0 k,23,0.003,0,0 k,24,0.0109,0,0 k,25,-0.0109,0.14,0 k,26,-0.003,0.14,0 k,27,0.003,0.14,0 k,28,0.0109,0.14,0 k,29,-0.06,0,0.1 a,5,13,14 a,6,15,16 a,7,20,19 a,8,17,18 a,1,21,13,4 a,24,2,3,15 a,19,10,11,28 a,9,17,25,12 a,21,22,5,13 a,22,23,6,5 a,23,24,15,6 a,5,6,16,14 a,14,16,20,18 a,18,20,7,8 a,17,8,26,25 a,8,7,27,26 a,7,19,28,27 l,1,29 !划分网格 lsel,s,,,16 lsel,a,,,18 lsel,a,,,28 lsel,a,,,22 lsel,a,,,14 lsel,a,,,30 lsel,a,,,32 lsel,a,,,20 lsel,a,,,26 lsel,a,,,40 lsel,a,,,42 lsel,a,,,24 lesize,all,,,4,,1 allsel lsel,s,,,31 lsel,a,,,33 lsel,a,,,35 lsel,a,,,37 lsel,a,,,39 lsel,a,,,43 lesize,all,,,3,,1 allsel lsel,s,,,1 lsel,a,,,3 lsel,a,,,6 lsel,a,,,4 lsel,a,,,10 lsel,a,,,12 lsel,a,,,7 lsel,a,,,9 lsel,a,,,29 lsel,a,,,34 lsel,a,,,41 lsel,a,,,44 lesize,all,,,4,,1 allsel lesize,36,,,24,-4 lesize,38,,,24,-4 lesize,13,,,10,0.5,,1 lesize,15,,,10,2,,1 lesize,17,,,10,2,,1 lesize,19,,,10,0.5,,1 lesize,25,,,10,0.5,,1 lesize,27,,,10,2,,1 lesize,21,,,10,2,,1 lesize,23,,,10,0.5,,1 type,2 amesh,all !拉伸平面网格为体网格 type,1 extopt,esize,50 extopt,aclear,1 vdrag,all,,,,,,45 !将solid70转换为对应结构分析单元 allsel etchg,tts keyopt,1,2,2 !定义焊缝区的单元组件和节点组件 vsel,s, , ,1 eslv,s nslv,s,1 cm,e1,elem cm,n1,node allsel vsel,s, , ,2 eslv,s nslv,s,1 cm,e2,elem cm,n2,node allsel vsel,s, , ,3 eslv,s nslv,s,1 cm,e3,elem cm,n3,node allsel vsel,s, , ,4 eslv,s nslv,s,1 cm,e4,elem cm,n4,node allsel !应力应变求解 /solu lsclear,all antype,4 nropt,full timint,on pred,on nlgeom,on autots,on lnsrch,on outres,basic,last tref,20 !防止刚体位移 d,node(0.06,0,0),all !杀死焊缝单元 cmsel,s,e1 cmsel,a,e2 cmsel,a,e3 cmsel,a,e4 ekill,all allsel !开始循环加载 j=0 *do,i,1,50 cmsel,s,n2 nsel,r,loc,z,0.002*(i-1),0.002*i esln,s,1 ealive,all !临时约束未激活节点 !allsel !esel,s,live !nsle,s,1 !nsel,inve !d,all,all allsel !读入温度场结果 j=j+0.4 ldread,temp,,,j,,T_field,rth time,j nsubst,4,,,1 kbc,0 solve *enddo !冷却阶段 ldread,temp,,,500,,T_field,rth time,500 nsubst,4,,,1 solve |
回复此楼» 猜你喜欢
26申博
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有22人回复
-
最失望的一年
已经有4人回复
-
拟解决的关键科学问题还要不要写
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有19人回复
-
请教限项目规定
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
疑惑？
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
|
大神们 可以帮我看看这个程序流吗 /PREP7 ET, 1, SOLID70 MP, DENS, 1, 7.839e3 MP, EMIS, 1, 0.8 MPTEMP, 1, 20,100,200,300,400,500 MPTEMP, 7, 600, 700, 765, 800, 900, 1000 MPDATA, KXX, 1, 1, 41.84, 40.88, 39.58, 37.69, 35.35, 32.15 MPDATA, KXX, 1, 7, 28.44, 24.61, 26.77, 17.38, 22.07, 24.88 MPTEMP, 1, 20,100,200,300,400,500 MPTEMP,7,600,700,765,800,900,1000 MPDATA,C,1,1,4.67e2,4.91e2,5.27e2,5.65e2,6.11e2,6.78e2 MPDATA,C,1,7,7.87e2,10.50e2,14.23e2,6.72e2,6.34e2,6.49e2 MPTEMP,1,40,60,80,100,120,140 MPTEMP,7,160,180,220,260,300,340 MPTEMP,13,380,480,580,680,780,980 MPTEMP,19,1180,1380,1580 MPDATA,HF,1,1,4.85,6.0,6.71,7.21,7.57,7.92 MPDATA,HF,1,7,8.21,8.45,8.73,8.89,9.91,9.34 MPDATA,HF,1,13,9.49,9.76,9.68,9.7,9.65,9.46 MPDATA,HF,1,19,9.23,9.05,8.87 TOFFST,273 /VIEW,1,1,1,1 BLOCK,0,0.03875,0,0.004,0,0.015 SAVE Wpoff,0,0,0.012 VSBW, 1 KWPAVE, 2 /DIST, 1, 1.08222638492,1 /REP, FAST FLST, 5, 6, 4, ORDE, 4 FITEM, 5, 1 FITEM,5,-4 FITEM,5,13 FITEM,5,-14 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1,0.000775/2, , , , , , ,1 !* FLST, 5, 8, 4, ORDE, 4 FITEM, 5, 5 FITEM, 5,-8 FITEM, 5, 21 FITEM, 5,-24 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1,0.003875/2, , , , , , ,1 !* FLST, 5, 2, 6, ORDE, 2 FITEM, 5, 2 FITEM, 5,-3 CM,_Y,VOLU VSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* MSHAPE,0,3d MSHKEY,1 VMESH,_Y1 MSHKEY,0 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* /AUTO, 1 /REP, FAST SAVE /SOLU Q=12984120 L=0.00775 LSIZE=0.1*L V=1/150 TINC=LSIZE/V MAX_X=50 MAX_Y=6 MAX_TIME=50 *DIM, FLUX2, TABLE, MAX_X, MAX_Y, MAX_TIME, X, Y, TIME *DO, K, 1, MAX_TIME, 1 *DO, I, 1, MAX_X, 1 FLUX2 (I, 0, K) = (I-1)*LSIZE *ENDDO *DO, J, 1, MAX_Y, 1 FLUX2 (0, J, K) = (J-1)*LSIZE *ENDDO *DO, I, 1, MAX_X, 1 *DO, J, 1, MAX_Y, 1 XCENTER=V*(K-1)*TINC YCENTER=0.002 DISTANCE=(I-1)*LSIZE-XCENTER *IF,ABS(DISTANCE),GT,L/2,THEN FLUX2(I,J,K)=0 *ELSE FLUX2(I,J,K)=Q*(1+2*DISTANCE/L) *ENDIF *ENDDO *ENDDO FLUX2(0,0,K)=(K-1)*TINC *ENDDO !* ANTYPE,TRANS TUNIF,20 KBC,1 *SET,DELT,1 TIMINT, ON AUTOTS, OFF TT=4.650001 *DO, TM, 1E-6, TT, TINC FLST, 2, 1, 5, ORDE, 1 FITEM, 2, 1 /GO !* TIME, TM AUTOTS, 0 KBC,0 !* TSRES,ERASE FLST,2,1,5,ORDE,1 FITEM, 2, 1 /GO !* !* SFA, P51X, 1, HFLUX, %FLUX2% /STATUS, SOLU SOLVE OUTRES, ALL, ALL, /PSF, HFLUX, 2 /REPLOT *ENDDO OUTRES, ALL, ALL, FINISH /PREP7 ETCHG,TTS ET, 1, 45 !* MPTEMP,,,,,,,, MPTEMP, 1,293,373,473,573,673,773 MPDATA, EX, 1, 1,209E9, 205E9, 200E9, 190E9, 183E9, 174E9 MPTEMP,,,,,,,, MPTEMP, 1,293,373,473,573,673,773 MPDATA,PRXY,1,1,0.295,0.291,0.295,0.277,0.319,0.336 MPTEMP,,,,,,,, MPTEMP, 1,373,473,573,673,773,873 MPTEMP, 7,973 MPDATA,ALpX,1,1,1.28E-5,1.34E-5,1.46E-5,1.46E-5,1.47E-5,1.47E-5 MPDATA, ALpX, 1, 7, 1.47E-5 TB, BKIN, 1, 6,,, TBTEMP, 293 TBDATA, 1, 9.5E8, 8.1E10 TEMP, 373 TBDATA, 1, 7.95E8, 8E10 TBTEMP, 473 TBDATA, 1, 8.35E8, 7.7E10 TBTEMP, 573 TBDATA, 1, 8.15E8, 7.5E10 TBTEMP,673 TBDATA,1,7.65E8,7E10 TEMP,773 TBDATA,1,6.75E8,6.5E10 FINI /SOLU LSCLEAR,SOLID ANTYPE,STATIC,NEW KBC,1 NSUBST,5,20,3 DELTIM, , , ,3 TIMINT，ON AUTOTS, ON L=0.00775 LSIZE=0.1*L V=1/150 TINC=LSIZE/V TT=2.90625 *DO, TM, 1E-6, TT, TINC ASEL, S, LOC, Z, 0.015 LSEL, S, EXT NSLL, S, L D, ALL, ALL ALLSEL TREF, 293 LDREAD, TEMP,,, TM,,TEMP2,RTH FINISH /SOLU ANTYPE,, REST TIME, TM SOLVE SAVE *ENDDO OUTRES, ALL, ALL, FINISH |
5楼2018-02-06 18:03:12
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2017-12-16 08:22:17
3楼2017-12-22 09:25:55
4楼2017-12-23 18:55:23