| 5 | 1/1 | 返回列表 |
| 查看: 1165 | 回復(fù): 4 | |||
jiangge2585銅蟲 (小有名氣)
|
[求助]
udf 熱源 初始化
|
|
求大神幫我看看我的程序哪里出問題了?這是一個激光高斯面熱源程序,我的意圖是:激光脈寬0.008s,頻率是10Hz,就是1秒激光打10次,每次持續(xù)0.008s。但是 模擬的結(jié)果卻是不能初始化,初始化時就卡主了,一直停在那。求大神救我! 這是程序: DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a,time_step; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; a=0; n=0; if(a<=1) { if(a>=0.108*n&&a<=0.108*n+0.08) { if(sqrt(x*x+y*y)<=0.001) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); } else { F_PROFILE(f,t,i)=0; } RP_Get_Real("flow-time" ;a=RP_Get_Real("flow-time" ;n=(int)a/0.108; } else { F_PROFILE(f,t,i)=0; } } else { Error("error reading file" ;} } end_f_loop(f,t) } |
銅蟲 (小有名氣)
|
這些天想了想,好像是begin_f_loop(f,t) 和循環(huán)程序沖突了,所以才會卡主,begin_f_loop(f,t) 掃描出一個面網(wǎng)格中心P,然后把符合(sqrt(x*x+y*y)<=0.001) 的P丟進循環(huán)語句執(zhí)行,執(zhí)行過一次之后,a和n的值就已經(jīng)超出了循環(huán)條件所規(guī)定的范圍,不能再進行循環(huán)了。 我就改了一下,但下面這個程序卻沒有發(fā)熱,不知道哪有問題了?求大神救我! #include "udf.h" DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; if(sqrt(x*x+y*y)<=0.001) for(a=n=0;a>=0.1*n&&a<=0.1*n+0.008&&a<=1;n=(int)a/0.1) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); a=RP_Get_Real("flow-time" ;} } end_f_loop(f,t) } |
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
新蟲 (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|
[考研] 環(huán)境學碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 南昌大學材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 材料學碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|