ansys焊接热分析直接提示不收敛。大佬看下哪里有错。到底哪里有错啊?

浏览:1317 回答:3

/prep7

/units,si

mptemp,1,20,250,500,750,1000,1500

mptemp,7,1700,2500

mpdata,kxx,1,1,50,47,40,27,30,35

mpdata,kxx,1,7,140,142

mpdata,ex,1,1,2.05e+11,1.87e+11,1.50e+11,7.00e+10,2.00e+10,1.90e+10

mpdata,ex,1,7,1.80e+10,1.20e+10

mpdata,dens,1,1,7820,7700,7610,7550,7490,7350

mpdata,dens,1,7,7300,7090

mpdata,c,1,1,460,480,530,675,670,660

mpdata,c,1,7,780,820

mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45

mpdata,prxy,1,7,0.48,0.5

mpdata,alpx,1,1,1.1e-05,1.22e-05,1.39e-05,1.48e-05,1.34e-05,1.33e-05,1.32e-05,1.31e-05

mpdata,alpx,1,7,1.32e-05,1.31e-05

mptemp,,,,,,,,

mptemp,1,20,1533,1590,1670

mpdata,enth,1,1,0,7.5e+009,9.6e+009,1.1e+010

rb=0.002                                        !单元尺寸     

length=0.12                                    !长方体几何参数

width=0.06

height=0.006

et,1,solid70

block,0,length,0,width,0,height

Block,0,length,0,3*Rb,0,height

Block,0,length,0,6*Rb,0,height

Block,0,length,0,12*Rb,0,height

VOVLAP,ALL

Lesize,4,rb/2                        !mapped划分单元尺寸

lesize,5,rb/2

lesize,19,rb/2

lesize,14,rb/2

lesize,50,rb/2

lesize,10,rb/2

lesize,52,rb/2

lesize,23,rb/2

lesize,51,rb/2

lesize,9,rb/2

lesize,49,rb/2

lesize,24,rb/2

lesize,31,rb

lesize,26,rb

lesize,43,4*rb

lesize,38,4*rb

lesize,2,4*rb

lesize,7,4*rb

lesize,53,4*rb

lesize,54,4*rb

lesize,55,4*rb

lesize,56,4*rb

lesize,11,0.003

lesize,12,0.003

lesize,47,0.003

lesize,48,0.003

lesize,35,0.003

lesize,36,0.003

lesize,64,0.004

lesize,61,0.004

lesize,62,0.004

lesize,63,0.004

lesize,60,0.002

lesize,57,0.002

lesize,58,0.002

lesize,59,0.002

type,1

mshkey,1

mshape,0

vmesh,5

vmesh,6

type,1

mshkey,0

mshape,1

vmesh,7

vmesh,8

/solu                !求解设置

antype,trans

timint,on,therm

tref,25

nropt,full,

kbc,1

autots,on

*do,i,1,60,1   !高斯移动热源加载

b=0.2*i            !时间步

time,b

deltim,0.2,0.01,0.2

xcenter=i*0.002 !加载中心点

*DEL,_FNCNAME   

*DEL,_FNCMTID   

*DEL,_FNCCSYS   

*SET,_FNCNAME,'flux2'   

*SET,_FNCCSYS,0 

! /INPUT,C:\Users\Administrator\Desktop\gaosi.func,,,1  

*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%  

!   

! Begin of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+   

! ({z}-0.006)^2)

*SET,%_FNCNAME%(0,0,1), 0.0, -999   

*SET,%_FNCNAME%(2,0,1), 0.0 

*SET,%_FNCNAME%(3,0,1), 0.0 

*SET,%_FNCNAME%(4,0,1), 0.0 

*SET,%_FNCNAME%(5,0,1), 0.0 

*SET,%_FNCNAME%(6,0,1), 0.0 

*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.01, 0, 0, 1   

*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 1 

*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 2, 2, -2 

*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 2, 0, 0, -1 

*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 17, -2   

*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 222222, 0, 0, -3

*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, -1, 3, -3

*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, 3  

*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 3, 17, -1

*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -2, 1, -3   

*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 0.006, 0, 0, 4 

*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, 4, 2, -2

*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, -3

*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -3, 17, -2  

*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, -1, 1, -4   

*SET,%_FNCNAME%(0,16,1), 0.0, -1, 7, 1, -2, 0, 0

*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 89524655, 0, 0, -1 

*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -2, 4, -1   

*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -3, 0, 0

! End of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+({z}-0.006)^ 

! 2)

allsel,all                            !选取加载中心点半径6mm内的点施加热流密度

nsel,s,loc,z,0.006

nsel,r,loc,x,xcenter-0.006,xcenter+0.006

nsel,r,loc,y,0,0.006

local,11,1,xcenter,0,0.006

csys,11

nsel,r,loc,x,0,0.006

csys,0

sf,all,hflux,%flux2%allsel,all            !高度方向4个面施加对流载荷

nsel,s,loc,x,0

sf,all,conv,10,20

allsel,all

nsel,s,loc,x,0.12

sf,all,conv,10,20

allsel,all

nsel,s,loc,y,0.06

sf,all,conv,10,20

allsel,all

nsel,s,loc,y,0

sf,all,conv,10,20

allsel,all

solve

*enddo

ansys焊接热分析直接提示不收敛。大佬看下哪里有错。到底哪里有错啊?的图1

邀请回答 我来回答

全部回答

(2)
默认 最新
邹正刚

  你的命令:

   sf,all,hflux,%flux2%allsel,all            !高度方向4个面施加对流载荷

  是在节点上施加热流载荷,不是在面上施加对流载荷,不知是否你需要的。检查一下载荷值是否太大了。否则,温度结果太大,难以收敛。

  而且该命令最后的 allsel,all 也应该另起一行,放在这里不知会造成什么后果。

  另外,你的求解设置过于简单。可以在求解前增加以下命令试试:

    nsubst, 10, 100,2

    autots,on

    kbc,0

    outres,all,all

    /gst,1

    ncnv,2

2019年9月27日
评论 1 点赞
小杜_6401
2019年9月26日
评论 点赞

没解决?试试专家一对一服务

换一批