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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 求助高斯光源加载不上去，初学者。
北京石油化工学院2026年研究生招生接收调剂公告
41/1返回列表
查看: 1365  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

桂~~~

新虫 (小有名气)

[求助] 求助高斯光源加载不上去，初学者。 已有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           !删除热流密度载荷
回复此楼

» 猜你喜欢

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

klzg77541812

专家顾问 (正式写手)

NVH-FEA工程师

【答案】应助回帖

感谢参与，应助指数 +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
赞一下 回复此楼
慢慢地，慢慢地，都随风而逝了......
2楼2019-07-03 11:51:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klzg77541812

专家顾问 (正式写手)

NVH-FEA工程师
上面少了个*enddo
赞一下 回复此楼
慢慢地，慢慢地，都随风而逝了......
3楼2019-07-03 13:01:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

346507037

铁虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
楼上大神说的对，定义局部坐标系，在局部坐标系选择单元，激活一个载荷步的单元，加载高斯热源。通过局部坐标系的移动，实现光源移动
赞一下 回复此楼
4楼2019-07-03 22:37:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 桂~~~ 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 337求调剂 +4 《树》 2026-03-29 4/200 2026-03-29 20:52 by 唐沐儿
[考研] 275求调剂 +15 Micky11223 2026-03-25 20/1000 2026-03-29 20:44 by 唐沐儿
[考研] 各位老师好，我的一志愿为北京科技大学085601材料专硕 +9 Koxui 2026-03-28 9/450 2026-03-29 20:43 by 无际的草原
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 085600 材料与化工 329分求调剂 +14 Mr. Z 2026-03-25 14/700 2026-03-29 17:27 by 学员8dgXkO
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 295材料工程专硕求调剂 +7 1428151015 2026-03-27 7/350 2026-03-28 19:58 by S240
[考研] 压国家一区线，求导师收留，有恩必谢！ +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 0703一志愿9，初试成绩：338，四六级已过，有科研经历，求调剂！ +4 Zuhui0306 2026-03-25 4/200 2026-03-28 13:07 by 唐沐儿
[考研] 292求调剂 +14 鹅鹅鹅额额额额� 2026-03-25 15/750 2026-03-28 08:45 by WYUMater
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 0703化学求调剂，各位老师看看我！！！ +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考博] 26申博 +3 加油冲啊！ 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 0703化学一志愿南京师范大学303求调剂 +3 zzffylgg 2026-03-24 3/150 2026-03-27 10:42 by shangxh
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
信息提示
关闭
请填处理意见
关闭