| 查看: 1531 | 回复: 3 | ||
[求助]
求助高斯光源加载不上去，初学者。 已有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 !删除热流密度载荷 |
回复此楼» 猜你喜欢
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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 |
2楼2019-07-03 11:51:17
klzg77541812
专家顾问 (正式写手)
NVH-FEA工程师
- 应助: 117 (高中生)
- 金币: 2141.8
- 红花: 28
- 帖子: 621
- 在线: 606.2小时
- 虫号: 1479383
- 注册: 2011-11-06
- 性别: GG
- 专业: 固体力学
- 管辖: 计算模拟区
3楼2019-07-03 13:01:52
4楼2019-07-03 22:37:03
