24小时热门版块排行榜

返回列表

桂~~~

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 26.1
散金: 10
红花: 1
帖子: 51
在线: 16.3小时
虫号: 10326881
注册: 2018-09-17
专业: 零件成形制造

[求助] 求助高斯光源加载不上去，初学者。已有2人参与

我做的是选择性激光烧结温度场的仿真，材料为聚苯乙烯，改了一下别人的程序后死活热源显示加载了，但是求解中就没有了。
/PREP7
!*
!======================================================== =======================================!指定单元
ET,1,SOLID70
!*
!*
!======================================================== =======================================!材料属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,25
MPTEMP,3,100
MPTEMP,4,165
MPTEMP,5,227
MPTEMP,6,327
MPDATA,KXX,1,,0.041
MPDATA,KXX,1,,0.047
MPDATA,KXX,1,,0.048
MPDATA,KXX,1,,0.049
MPDATA,KXX,1,,0.05
MPDATA,KXX,1,,0.05
MPTEMP,,
MPTEMP,1,0
MPTEMP,2,25
MPTEMP,3,100
MPTEMP,4,165
MPTEMP,5,227
MPTEMP,6,327
MPDATA,C,1,,2810
MPDATA,C,1,,2904
MPDATA,C,1,,3046
MPDATA,C,1,,3612
MPDATA,C,1,,3821
MPDATA,C,1,,3912
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,472
*dim,conve,table,8,1,1,temp       !总换热系数(W/m^2*C)
conve(0,1,1)=1
conve(1,0,1)=0
conve(1,1,1)=12.12e-3
conve(2,0,1)=200
conve(2,1,1)=24.23e-3
conve(3,0,1)=400
conve(3,1,1)=40.38e-3
conve(4,0,1)=600
conve(4,1,1)=72.69e-3
conve(5,0,1)=800
conve(5,1,1)=125e-3
conve(6,0,1)=1000
conve(6,1,1)=150e-3
conve(7,0,1)=1500
conve(7,1,1)=200e-3
conve(8,0,1)=2000
conve(8,1,1)=250e-3
!======================================================== =======================================!定义常量
Width_Base=0.025  !基底宽度
Height_Base=0.002  !基底高度
Length_Base=0.05 !基底长度
Length_Clad=0.05 !熔覆层长度
Height_Clad=0.0025 !熔覆层高度
Width_Clad=0.025  !覆层总宽度
!======================================================== =======================================!定义常量
BLOCK,0,Length_Clad,-Width_Clad/2,Width_Clad/2,0,Height_Clad, !覆层
BLOCK,0,Length_Base,-Width_Base/2,Width_Base/2,-Height_Base,0, !基底
esize,0.001
vmesh,all
SFA,3,1,CONV,%conve%,90
SFA,4,1,CONV,%conve%,90
SFA,5,1,CONV,%conve%,90
SFA,6,1,CONV,%conve%,90
SFA,7,1,CONV,%conve%,90
SFA,8,1,CONV,%conve%,90
SFA,9,1,CONV,%conve%,90
SFA,10,1,CONV,%conve%,90
SFA,11,1,CONV,%conve%,90
V=2  !扫描速度
P=20  !激光功率
r=0.00015  !激光光斑半径
Area=3.14159265*(r**2)  !激光光斑面积

Qm=2*P/Area
StepDis=0.001  !每个载荷步位移
StepTime=StepDis/V  !每个载荷步时间
StepNum=(Length_Clad)/StepDis  !载荷步数
!************************************非线性求解选项****************************************
/solu
antype,4
trnopt,full          !完全瞬态积分法
pred,on             !打开预测校正
solcontrol,on       !自动控制
nropt,full,,on       !完全牛顿-拉普森方法
!cnvtol,u,,0.5       !设置收敛值，位移
!cnvtol,f,,0.05       !设置收敛值，力
cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少
!eqslv,sparse,,2       !求解器类型
lnsrch,on             !线性搜索
autots,on             !自动步长
tintp,,,,1,,          !向后积分
timint,on             !瞬态效果
tref,90             !参考温度
kbc,0                !连续性载荷
neqit,50             !规定每个子步中最大迭代次数,默认为25
!======================================================== =======================================!高斯单元
!**** 高斯光束能量公式定义（GUI方式自动生成）****
!****GAUSS=2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)****
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS
*SET,_FNCNAME,'GAOSI'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),P  !
*SET,_FNC_C2(1),R !
*SET,_FNCCSYS,11  !将方程坐标系设定为自己定义的局部坐标系
! /INPUT,GAUSS.func,,,1
*DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS%
!
! Begin of equation: 2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17
*SET,%_FNCNAME%(0,3,1), 0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 18, 17, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 4, -4
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -2
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 2, 17, -2
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, 3, 17, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -3, 3, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 18, 17, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -3
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation:2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)
!-->
*DO, k, 1, StepNum, 1

SFA,2,1,HFLUX, %GAOSI%       !激活单元的上表面,加热流密度

OUTRES,ALL,ALL,

TIME,StepTime*k  !载荷步结束时间
NSUBST, 1  !指定载荷子步数目
DELTIM,0.1,0.1,0.1 !载荷子步

solve

SFADELE,all,HFLUX          !删除热流密度载荷

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2019-07-03 09:12:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

klzg77541812

专家顾问 (正式写手)

NVH-FEA工程师

应助: 117 (高中生)
金币: 2141.8
红花: 28
帖子: 621
在线: 606.2小时
虫号: 1479383
注册: 2011-11-06
性别: GG
专业: 固体力学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

!1定义局部坐标系 2 加载命令放在循环语句外面，因为加载的是函数，没有必要重复加 3 热源处网格太少了，加密网格 4 计算总时间太短，可能你还没调整 5 热流密度的大小不知道对不对，显示加载数值很小0.12

/PREP7
!*
!======================================================== =======================================!指定单元
ET,1,SOLID70
!*
!*
!======================================================== =======================================!材料属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,25
MPTEMP,3,100
MPTEMP,4,165
MPTEMP,5,227
MPTEMP,6,327
MPDATA,KXX,1,,0.041
MPDATA,KXX,1,,0.047
MPDATA,KXX,1,,0.048
MPDATA,KXX,1,,0.049
MPDATA,KXX,1,,0.05
MPDATA,KXX,1,,0.05
MPTEMP,,
MPTEMP,1,0
MPTEMP,2,25
MPTEMP,3,100
MPTEMP,4,165
MPTEMP,5,227
MPTEMP,6,327
MPDATA,C,1,,2810
MPDATA,C,1,,2904
MPDATA,C,1,,3046
MPDATA,C,1,,3612
MPDATA,C,1,,3821
MPDATA,C,1,,3912
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,472
*dim,conve,table,8,1,1,temp       !总换热系数(W/m^2*C)
conve(0,1,1)=1
conve(1,0,1)=0
conve(1,1,1)=12.12e-3
conve(2,0,1)=200
conve(2,1,1)=24.23e-3
conve(3,0,1)=400
conve(3,1,1)=40.38e-3
conve(4,0,1)=600
conve(4,1,1)=72.69e-3
conve(5,0,1)=800
conve(5,1,1)=125e-3
conve(6,0,1)=1000
conve(6,1,1)=150e-3
conve(7,0,1)=1500
conve(7,1,1)=200e-3
conve(8,0,1)=2000
conve(8,1,1)=250e-3
!======================================================== =======================================!定义常量
Width_Base=0.025  !基底宽度
Height_Base=0.002  !基底高度
Length_Base=0.05 !基底长度
Length_Clad=0.05 !熔覆层长度
Height_Clad=0.0025 !熔覆层高度
Width_Clad=0.025  !覆层总宽度
!======================================================== =======================================!定义常量
BLOCK,0,Length_Clad,-Width_Clad/2,Width_Clad/2,0,Height_Clad, !覆层
BLOCK,0,Length_Base,-Width_Base/2,Width_Base/2,-Height_Base,0, !基底
esize,0.001
vmesh,all
SFA,3,1,CONV,%conve%,90
SFA,4,1,CONV,%conve%,90
SFA,5,1,CONV,%conve%,90
SFA,6,1,CONV,%conve%,90
SFA,7,1,CONV,%conve%,90
SFA,8,1,CONV,%conve%,90
SFA,9,1,CONV,%conve%,90
SFA,10,1,CONV,%conve%,90
SFA,11,1,CONV,%conve%,90
V=2  !扫描速度
P=20  !激光功率
r=0.00015  !激光光斑半径
Area=3.14159265*(r**2)  !激光光斑面积

Qm=2*P/Area
StepDis=0.001  !每个载荷步位移
StepTime=StepDis/V  !每个载荷步时间
StepNum=(Length_Clad)/StepDis  !载荷步数
!************************************非线性求解选项****************************************
/solu
antype,4
trnopt,full          !完全瞬态积分法
pred,on             !打开预测校正
solcontrol,on       !自动控制
nropt,full,,on       !完全牛顿-拉普森方法
!cnvtol,u,,0.5       !设置收敛值，位移
!cnvtol,f,,0.05       !设置收敛值，力
cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少
!eqslv,sparse,,2       !求解器类型
lnsrch,on             !线性搜索
autots,on             !自动步长
tintp,,,,1,,          !向后积分
timint,on             !瞬态效果
tref,90             !参考温度
kbc,0                !连续性载荷
neqit,50             !规定每个子步中最大迭代次数,默认为25

CSWPLA,11,0,1,1,  !定义局部坐标系11

!======================================================== =======================================!高斯单元
!**** 高斯光束能量公式定义（GUI方式自动生成）****
!****GAUSS=2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)****
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS
*SET,_FNCNAME,'GAOSI'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),P  !
*SET,_FNC_C2(1),R !
*SET,_FNCCSYS,11  !将方程坐标系设定为自己定义的局部坐标系
! /INPUT,GAUSS.func,,,1
*DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS%
!
! Begin of equation: 2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17
*SET,%_FNCNAME%(0,3,1), 0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 18, 17, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 4, -4
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -2
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 2, 17, -2
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, 3, 17, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -3, 3, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 18, 17, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -3
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation:2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)
!-->

  SFA,2,1,HFLUX, %GAOSI%       !激活单元的上表面,加热流密度

*DO, k, 1, StepNum, 1

OUTRES,ALL,ALL,

TIME,StepTime*k  !载荷步结束时间
NSUBST, 1  !指定载荷子步数目
DELTIM,0.1,0.1,0.1 !载荷子步

solve