应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 焊接应力场仿真不收敛问题
51/1返回列表
查看: 2096  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Fanph

铁虫 (小有名气)

[求助] 焊接应力场仿真不收敛问题 已有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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2017-12-15 20:49:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
你的命令流中，没有进行热分析，直接就转换为应力分析，为什么这样做？
赞一下 回复此楼
2楼2017-12-16 08:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Fanph

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by shgao20 at 2017-12-16 08:22:17
你的命令流中，没有进行热分析，直接就转换为应力分析，为什么这样做？

我已经做完热分析了，这是应力分析的部分，它的加载就是调取热分析的结果
赞一下 回复此楼
3楼2017-12-22 09:25:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zndxwdm

新虫 (初入文坛)

【答案】应助回帖

看到和我毕业论文一样的了，忘得比较多，不收敛的原因太多，可以逐步查查约束是否合理，此外，结构分析的时候可以用181单元，这样收敛会好很多。
赞一下 回复此楼
4楼2017-12-23 18:55:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奥利奥麦旋风

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by shgao20 at 2017-12-16 08:22:17
你的命令流中，没有进行热分析，直接就转换为应力分析，为什么这样做？

大神们 可以帮我看看这个程序流吗
/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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Fanph 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
关闭
请填处理意见
关闭