| 查看: 764 | 回复: 10 | |||
11050209银虫 (正式写手)
|
[求助]
新手求助,帮看看这个加载热源时定义的表格哪错了已有1人参与
|
||
|
!删除以下参数变量 *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNCCSYS !定义表格各行 *SET,_FNCNAME,'YUANZHUI' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 !表格各行赋值 *SET,_FNC_C1(1),W *SET,_FNC_C2(1),PI *SET,_FNC_C3(1),PI/6 *SET,_FNC_C4(1),V *SET,_FNC_C5(1),R *SET,_FNCCSYS,0 !表格各数据赋值 *DIM,%_FNCNAME%,TABLE,7,21,1,,,,%_FNCCSYS% *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), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(0,1,1),1.0, -1, 0, 1, 18, 3, 4 *SET,%_FNCNAME%(0,2,1),0.0, -2, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,3,1),0, -3, 0, 3, 0, 0, -2 *SET,%_FNCNAME%(0,4,1),0.0, -4, 0, 1, -2, 4, -3 *SET,%_FNCNAME%(0,5,1),0.0, -2, 11, 1, 19, 0, 0 *SET,%_FNCNAME%(0,6,1),0.0, -3, 0, 2, 0, 0, -2 *SET,%_FNCNAME%(0,7,1),0.0, -5, 0, 1, -2, 17, -3 *SET,%_FNCNAME%(0,8,1),0.0, -2, 0, 1, -4, 3, -5 *SET,%_FNCNAME%(0,9,1),0.0, -3, 0, 2, 0, 0, 4 *SET,%_FNCNAME%(0,10,1),0.0, -4, 0, 1, 4, 17, -3 *SET,%_FNCNAME%(0,11,1),0.0, -3, 0, 1, -2, 3, -4 *SET,%_FNCNAME%(0,12,1),0.0, -2, 11, 1, 19, 0, 0 *SET,%_FNCNAME%(0,13,1),0.0, -4, 0, 1, -2, 3, 20 *SET,%_FNCNAME%(0,14,1),0.0, -2, 0, 1, -4, 3, 1 *SET,%_FNCNAME%(0,15,1),0.0, -4, 0, 1, -3, 1, -2 *SET,%_FNCNAME%(0,16,1),0.0, -2, 0, 2, 0, 0, 21 *SET,%_FNCNAME%(0,17,1),0.0, -3, 0, 1, 21, 17, -2 *SET,%_FNCNAME%(0,18,1),0.0, -2, 0, 1, -4, 1, -3 *SET,%_FNCNAME%(0,19,1),0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,20,1),0.0, -1, 0, 1, 17, 4, -3 *SET,%_FNCNAME%(0,21,1),0.0, 99, 0, 1, -1, 0, 0 |
» 猜你喜欢
投稿精细化工
已经有6人回复
博士读完未来一定会好吗
已经有36人回复
之前让一硕士生水了7个发明专利,现在这7个获批发明专利的维护费可从哪儿支出哈?
已经有10人回复
博士申请都是内定的吗?
已经有9人回复
心脉受损
已经有8人回复
读博
已经有5人回复
» 本主题相关价值贴推荐,对您同样有帮助:
新手,求助一下
已经有1人回复
新手求助一下,导师的回复,请帮我分析下
已经有4人回复
求助XPS分峰拟合,α-氧化铁的,想看看反应前后的二价铁三价铁的比例。
已经有13人回复
请大神帮我看看这个XPS的数据,新手求助
已经有0人回复
新手求助!elsevier在线投稿时通讯作者搞错了~~怎么办啊?
已经有4人回复
大神们,帮看一下什么意思,已接受,要怎么做!
已经有2人回复
化学电源方面新手求助
已经有11人回复
找过渡态,虚频太小,结构接下来怎么调才能找到过渡态?
已经有7人回复
新手求助关于Biotechnology Letters投稿格式方面的问题
已经有5人回复
新手投稿求助,请大家帮我看看
已经有3人回复
新手求助
已经有0人回复
求助,电化学活性比表面积的计算!!!
已经有2人回复
求助,抑郁和考博
已经有13人回复
【求助】我是新手~想公派读博~传播学或广告学,对很多东西都不了解,哪位好心人能帮帮
已经有1人回复
【求助】分析阻抗,小弟是新手,上次发个贴子,没说明白,这次大家帮一下啊
已经有1人回复
【求助】紧急求助——伯胺的N-烃基化,缚酸剂的选择!
已经有2人回复
【求助】关于聚合物胶束的几个问题,求高人解答!
已经有6人回复
【求助】大家来帮忙看看这个反应路线,
已经有2人回复
【求助】大家帮我看看这个问题是否属于流固耦合问题
已经有1人回复
shgao20
专家顾问 (职业作家)
-

专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
8楼2016-05-17 07:47:54
shgao20
专家顾问 (职业作家)
-

专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2016-05-13 07:32:35
11050209
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 938.2
- 散金: 100
- 帖子: 321
- 在线: 68.3小时
- 虫号: 3896498
- 注册: 2015-05-28
- 专业: 金属材料的凝固与结晶学
3楼2016-05-13 08:34:58
shgao20
专家顾问 (职业作家)
-

专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2016-05-14 08:54:54
11050209
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 938.2
- 散金: 100
- 帖子: 321
- 在线: 68.3小时
- 虫号: 3896498
- 注册: 2015-05-28
- 专业: 金属材料的凝固与结晶学
5楼2016-05-14 14:56:10
shgao20
专家顾问 (职业作家)
-

专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
6楼2016-05-15 07:58:31
11050209
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 938.2
- 散金: 100
- 帖子: 321
- 在线: 68.3小时
- 虫号: 3896498
- 注册: 2015-05-28
- 专业: 金属材料的凝固与结晶学
|
谢谢大神,下面是我在某个论坛看到的一个案例,楼主说是热源加不上,我也看不出问题在哪,不知道大神能不能给看看 /units,si /prep7 lx=0.03 ly=0.03 lz=0.004 v=0.01 p=2000 ef=0.3 lsize=0.0002 pi=3.1415926 h=0.003 ri=0.0002 re=0.001 n1_fz=h/lsize+1 mptemp MPTEMP,1,20 MPTEMP,2,100 MPTEMP,3,200 MPTEMP,4,300 MPTEMP,5,400 MPTEMP,6,500 MPTEMP,7,600 MPTEMP,8,700 MPTEMP,9,730 MPTEMP,10,750 MPTEMP,11,800 mptemp,12,1000 mptemp,13,1500 mptemp,14,1840 mptemp,15,1890 mptemp,16,2860 MPDATA,KXX,1,,51.9 MPDATA,KXX,1,,49.4 MPDATA,KXX,1,,49 MPDATA,KXX,1,,46.1 MPDATA,KXX,1,,42.7 MPDATA,KXX,1,,39.4 MPDATA,KXX,1,,35.6 MPDATA,KXX,1,,31.8 MPDATA,KXX,1,,30.1 MPDATA,KXX,1,,28.9 MPDATA,KXX,1,,26 mpdate,kxx,1,,27.2 mpdata,kxx,1,,29.7 mpdata,kxx,1,,29.7 mpdata,kxx,1,,29.7 mpdata,kxx,1,,29.7 MPDATA,C,1,,450 MPDATA,C,1,,494 MPDATA,C,1,,526 MPDATA,C,1,,566 MPDATA,C,1,,615 MPDATA,C,1,,684 MPDATA,C,1,,1139 MPDATA,C,1,,1384 MPDATA,C,1,,1191 MPDATA,C,1,,931 MPDATA,C,1,,779 MPDATA,C,1,,400 MPDATA,C,1,,847 MPDATA,C,1,,400 MPDATA,C,1,,400 MPDATA,DENS,1,,7872 MPDATA,DENS,1,,7845 MPDATA,DENS,1,,7816 MPDATA,DENS,1,,7740 MPDATA,DENS,1,,7733 MPDATA,DENS,1,,7711 MPDATA,DENS,1,,7669 MPDATA,DENS,1,,7625 MPDATA,DENS,1,,7612 MPDATA,DENS,1,,7602 MPDATA,DENS,1,,7578 MPDATA,DENS,1,,7552 MPDATA,DENS,1,,7268 MPDATA,DENS,1,,6757 MPDATA,DENS,1,,6715 MPDATA,DENS,1,,5902 mptemp MPTEMP,1,20 MPTEMP,2,100 MPTEMP,3,300 MPTEMP,4,600 MPTEMP,5,900 MPTEMP,6,1000 MPTEMP,7,1300 MPTEMP,8,1600 MPTEMP,9,2000 MPTEMP,10,2300 MPTEMP,11,2600 mptemp,12,3000 mptemp,13,3100 mptemp,14,3300 mptemp,15,3500 MPDATA,HF,1,,5 MPDATA,HF,1,,12 MPDATA,HF,1,,32 MPDATA,HF,1,,73 MPDATA,HF,1,,123.8 MPDATA,HF,1,,141.2 MPDATA,HF,1,,198.5 MPDATA,HF,1,,290.8 MPDATA,HF,1,,364 MPDATA,HF,1,,446 MPDATA,HF,1,,532.4 MPDATA,HF,1,,724 MPDATA,HF,1,,819 MPDATA,HF,1,,963 et,1,solid70 et,2,surf152 keyopt,2,4,1 keyopt,2,5,0 keyopt,2,8,3 block,0,lx,0,ly,0,-lz block,0,lx,0,0.004,0,-lz block,0,lx,0.004,0.0044,0,-lz block,0,lx,0.0044,0.008,0,-lz block,0,lx,0.008,0.009,0,-lz vovlap,1,2,3,4,5 vglue,6,7,8,9,10 esize,0.0002 vmesh,6 esize,0.0005 vmesh,9 esize,0.002 vmesh,7 mshkey,0 mshape,1,3d CM,_Y,VOLU VSEL, , , , 8 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 mshkey,0 mshape,1,3d CM,_Y,VOLU VSEL, , , , 10 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 type,2 asel,s,,,31,50,1 asel,a,,,4 amesh,all save step=0.0002 n=lx/step *do,kk,0,n,1 xx=kk*step tm=(xx+step)/v /SOLU ANTYPE,4 TRNOPT,FULL TUNIF,20, OUTRES,NSOL,last NROPT,FULL,,ON time,tm NSUBST,5, , ,1 !!! KBC,1 lnsrch,1 !!! TSRES,ERASE TINTP,,,,1,,,,,!!! TIMINT,0,struct TIMINT,1,THERM TIMINT,0,MAG esel,s,type,,2 !第一个载荷步,有对流和热生成率 nsle,s sf,all,conv,-1,20 allsel *afun,deg !!!选择锥形热源内的节点 local,12,1,xx,0,0 esel,s,type,,1 nsle,s nsel,r,loc,z,-h,-h nsel,r,loc,x,0,ri nsel,r,loc,y,0,180 *do,k,1,n1_fz,1 hk=k*lsize rk=ri+hk*(re-ri)/h nsel,a,loc,z,-h+hk,-h+hk nsel,r,loc,x,0,rk nsel,r,loc,y,0,180 cm,node_f,node *enddo local,14,0,xx,0,0 *DEL,_FNCNAME !!!!!定义锥形热源函数 *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNCCSYS *del,'yzx' *SET,_FNCNAME,'yzx' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *SET,_FNC_C1(1),ef *SET,_FNC_C2(1),p *SET,_FNC_C3(1),re *SET,_FNC_C4(1),ri *SET,_FNC_C5(1),h *SET,_FNCCSYS,14 ! /INPUT,yzx.func,,,1 *DIM,%_FNCNAME%,TABLE,7,30,1,,,,%_FNCCSYS% ! ! Begin of equation: 2*ef*p*exp(-2*({X}^2+{Y}^2)/(re+(re-ri)*{Z}/h)^2)/{PI}*re^ ! 2*h *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), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *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, 1, -2, 3, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 2, 0, 0, -4 *SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 3, 17, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 19, 2, 20 *SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -2, 3, 4 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -3, 4, 21 *SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, 19, 1, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, -3 *SET,%_FNCNAME%(0,20,1), 0.0, -5, 0, 1, -3, 17, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -5 *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, -1, 0, 3.14159265358979310, 0, 0, -3 *SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 1, -3, 4, -1 *SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 2, 0, 0, 19 *SET,%_FNCNAME%(0,27,1), 0.0, -3, 0, 1, 19, 17, -1 *SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -2, 3, -3 *SET,%_FNCNAME%(0,29,1), 0.0, -2, 0, 1, -1, 3, 21 *SET,%_FNCNAME%(0,30,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: 2*ef*p*exp(-2*({X}^2+{Y}^2)/(re+(re-ri)*{Z}/h)^2)/{PI}*re^2* ! h !--> bf,node_f,hgen,%yzx% !!!!加载圆锥热源 cmdele,node_f allsel,all solve bfdele,all,hgen *enddo save |
7楼2016-05-16 17:31:28
11050209
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 938.2
- 散金: 100
- 帖子: 321
- 在线: 68.3小时
- 虫号: 3896498
- 注册: 2015-05-28
- 专业: 金属材料的凝固与结晶学
9楼2016-05-17 09:00:48
shgao20
专家顾问 (职业作家)
-

专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
10楼2016-05-18 07:34:09













回复此楼