版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

yxy657828314

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 17.5
帖子: 11
在線: 4.4小時(shí)
蟲號(hào): 4195348
注冊(cè): 2015-11-04

[求助] udf編程問題已有2人參與

最近在使用udf進(jìn)行編程，主要使用兩種宏，difine和profile，程序目的主要是截取三維模型進(jìn)出口參數(shù)利用一維進(jìn)行計(jì)算，三維和一維不斷進(jìn)行相互參數(shù)互換。
在difine宏里面插入一部分C++的程序，但是fluent迭代的時(shí)候就卡機(jī)了，不知道是因?yàn)槭裁丛颍雴栆幌率欠袷且驗(yàn)镃++程序里面使用了do-while語句。
下面是我的源程序
#include "udf.h"
real NV_VEC(A);
real avg_tempout;
real avg_tempin;
real sum_T_A=0.0;
real sum_A=0.0;
real sum_V_A=0.0;
real V2;
real V1;
Thread *thread_out;
Thread *thread_in;
face_t f;
Domain *domain;
int i;
int j;
int k;
real q[200];
real F[200];
real F1[200];
real F2[200];
real Y[200];
real Fmax;
real pj[200];
real Ar;
real Aj;
real vj;
real zsxs;
real dpj[200];
real Ck[3][6];
real ph[6];
real R[6];
real G;
real Q[6];
real fjts[6];
real hfy;
real f1;
real f2;
real Ar=80;
real Aj=0.8;
real vj=30;
real zsxs=1;
real Ck[3][6]={0,0,0,0,0,0,0,1,0,1,-1,0,0,0,1,1,0,-1};
real R[6]={0,0.02,0.02,0.01,0.02,0.02};
real fjts[6]={0,6,12,0,0,0};
real Ph[6]={0,0,0,0,0,0};
real Q[6]={0,295,20,15,5,280};
DEFINE_ADJUST(inlet_T,d)
{
domain=Get_Domain(1);
thread_in=Lookup_Thread(domain,6);
begin_f_loop(f,thread_in)
{
F_AREA(A,f,thread_in);
sum_A+=NV_MAG(A);
sum_T_A+=NV_MAG(A)*F_T(f, thread_in);
}
end_f_loop(f,thread_in)
avg_tempin=sum_T_A/sum_A;

}
DEFINE_ADJUST(inlet_v,d)
{
domain=Get_Domain(1);
thread_in=Lookup_Thread(domain,6);
begin_f_loop(f,thread_in)
{
  F_AREA(A,f,thread_in);
  sum_A+=NV_MAG(A);
  sum_V_A+=NV_MAG(A)*F_W(f,thread_in);
}
end_f_loop(f,thread_in);
V1=sum_V_A/sum_A;

}
DEFINE_ADJUST(outlet_T,d)
{domain=Get_Domain(1);
thread_out=Lookup_Thread(domain,5);
begin_f_loop(f,thread_out)
{
F_AREA(A,f,thread_out);
sum_A+=NV_MAG(A);
sum_T_A+=NV_MAG(A)*F_T(f, thread_out);
}
end_f_loop(f,thread_out)
avg_tempout=sum_T_A/sum_A;

}
DEFINE_ADJUST(outlet_v,d)
{
domain=Get_Domain(1);
thread_out=Lookup_Thread(domain,5);
begin_f_loop(f,thread_out)
{
  F_AREA(A,f,thread_out);
  sum_A+=NV_MAG(A);
  sum_V_A+=NV_MAG(A)*F_W(f,thread_out);
}
end_f_loop(f,thread_out);
V2=sum_V_A/sum_A;

f1=101000*29/8.314/avg_tempin;
   f2=101000*29/8.314/avg_tempout;
Ph[1]=0.5*f1*V2*V2-0.5*f2*V1*V1+0.5*9.8*12*f1-0.5*9.8*12*f2;
for(i=1;i<=200;i++)
F=0.0;
for(i=1;i<=200;i++)
F1=0.0;
for(i=1;i<=200;i++)
F2=0.0;
for(i=1;i<=200;i++)
Y=0.0;
for(i=1;i<=200;i++)
q=10000.0;
for(i=1;i<=200;i++)
pj=0.0;
for(i=1;i<=5;i++)
{
pj=fjts*(-1.2*vj*Aj*zsxs*fabs(Q)/Ar/Ar+1.2*vj*Aj*zsxs*vj/Ar);
dpj=fjts*(-1.2*vj*Aj*zsxs/Ar/Ar);}
   do{for(i=1;i<=2;i++)
   {  F=0;
      Y=0;
for(j=1;j<=5;j++)
{
   F=F+Ck[j]*R[j]*Q[j]*fabs(Q[j])-Ck[j]*pj[j];
   Y=Y+2*Ck[j]*Ck[j]*R[j]*fabs(Q[j])-Ck[j]*Ck[j]*dpj[j];}
   q=-F/Y;
for(j=1;j<=5;j++)
{if(Ck[j]>0)
Q[j]=Q[j]+q;
   else if(Ck[j]<0)
Q[j]=Q[j]-q;
}

F1=0;
F2=0;
for(j=1;j<=5;j++)
{if(Ck[j]/Q[j]>0)
F1=F1+Ck[j]*R[j]*Q[j]*fabs(Q[j]);}
for(j=1;j<=5;j++)
{if(Ck[j]/Q[j]<0)
F2=F2+Ck[j]*R[j]*Q[j]*fabs(Q[j]);}

for(j=1;j<=5;j++)
{if(Ck[j]*Ph[6]<0)
F1=F1-Ck[j]*Ph[6];}
for(j=1;j<=5;j++)
{if(Ck[j]*Ph[6]>0)
F2=F2-Ck[j]*Ph[6];}

for(j=1;j<=5;j++)
{if(Ck[j]*pj[j]<0)
F1=F1-Ck[j]*pj[j];}
for(j=1;j<=5;j++)
{if(Ck[j]*pj[j]>0)
F2=F2-Ck[j]*pj[j];}
for(j=1;j<=5;j++)
{ pj[j]=fjts[j]*(-1.2*vj*Aj*zsxs*fabs(Q[j])/Ar/Ar+1.2*vj*Aj*zsxs*vj/Ar);
      dpj[j]=fjts[j]*(-1.2*vj*Aj*zsxs/Ar/Ar);
      }
   }
Fmax=0;
for(k=1;k<=2;k++)
   {if(fabs(fabs(F1[k]/F2[k])-1)>Fmax)
Fmax=fabs(fabs(F1[k]/F2[k])-1);
}
} while(Fmax>0.0001);
DEFINE_PROFILE(in,thread,position)
{

begin_f_loop(f,thread)
{
  if(fabs((Q[1]-295)/295)>0.1)
  F_PROFILE(f,thread,position)=Q[1];
  else
  F_PROFILE(f,thread,position)=Q[1]/50;
  Message("\\nQ = %g\\n",Q[1]);
}
end_f_loop(f,thread)

}
DEFINE_PROFILE(out,thread,position)
{

begin_f_loop(f,thread)
{
   if(fabs((Q[1]-295)/295)>0.1)
  F_PROFILE(f,thread,position)=Q[1];
  else
  F_PROFILE(f,thread,position)=Q[1]*Q[1]/80/80;
   Message("\\nQ = %g\\n",Q[1]);
}
end_f_loop(f,thread)

}

麻煩大神幫忙看看是因?yàn)槭裁丛虺绦虿荒苷＿\(yùn)行

回復(fù)此樓

» 猜你喜歡

材料與化工調(diào)劑一志愿大連海事085600，349 已經(jīng)有4人回復(fù)
甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理已經(jīng)有4人回復(fù)
理學(xué)07化學(xué) 303求調(diào)劑已經(jīng)有3人回復(fù)
一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有11人回復(fù)
生物考研337分求調(diào)劑已經(jīng)有4人回復(fù)
材料求調(diào)劑一志愿哈工大總分298分，前三科223分已經(jīng)有9人回復(fù)
哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分已經(jīng)有9人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有8人回復(fù)

1樓 2016-06-22 20:39:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mycc

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +181
應(yīng)助: 230 (大學(xué)生)
金幣: 2589.6
紅花: 35
帖子: 531
在線: 58.8小時(shí)
蟲號(hào): 881523
注冊(cè): 2009-10-23
性別: GG
專業(yè): 流體力學(xué)
管轄: 仿真模擬

【答案】應(yīng)助回帖

自己在UDF中添加一些輸出語句，看看錯(cuò)誤出現(xiàn)在什么位置，這么長(zhǎng)的UDF需要你自己調(diào)試，我?guī)湍愕脑�，要收費(fèi)的

贊一下

回復(fù)此樓

項(xiàng)目匯總：http://www.cfluid.com/forum.php?mod=viewthread&tid=114340&extra=

7樓2016-06-27 09:37:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個(gè)回答

mycc

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +181
應(yīng)助: 230 (大學(xué)生)
金幣: 2589.6
紅花: 35
帖子: 531
在線: 58.8小時(shí)
蟲號(hào): 881523
注冊(cè): 2009-10-23
性別: GG
專業(yè): 流體力學(xué)
管轄: 仿真模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

do-while語句是可以用的，UDF文件可以自己調(diào)試，看看問題出在什么地方

贊一下

回復(fù)此樓

項(xiàng)目匯總：http://www.cfluid.com/forum.php?mod=viewthread&tid=114340&extra=

2樓2016-06-23 08:07:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

匿名

用戶注銷 (小有名氣)

應(yīng)助: 31 (小學(xué)生)
金幣: 1248.2
散金: 68
紅花: 7
帖子: 220
在線: 78.3小時(shí)
蟲號(hào): 0
注冊(cè): 2015-01-18
專業(yè): 工程熱物理與能源利用

本帖僅樓主可見

贊一下

3樓2016-06-23 09:07:42

已閱申請(qǐng)仿真EPI 回復(fù)此樓編輯查看我的主頁

小七工作室1

新蟲 (著名寫手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣-100, 應(yīng)助指數(shù)-2, 屏蔽內(nèi)容, 違規(guī)存檔, 多次用完全一樣、無實(shí)質(zhì)性幫助的內(nèi)容，回復(fù)不同應(yīng)助帖！ 2017-03-24 11:15:47

本帖內(nèi)容被屏蔽

4樓2016-06-24 08:56:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求收留 +6	1943443204 2026-03-28	6/300	2026-03-30 20:37 by dick_runner
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9	hanamiko 2026-03-27	9/450	2026-03-30 20:10 by 無際的草原
[考研] 生物學(xué) 296 求調(diào)劑 +5	朵朵- 2026-03-26	7/350	2026-03-30 20:07 by 源_2020
[考研] 085701求調(diào)劑初試286分 +5	secret0328 2026-03-28	5/250	2026-03-30 12:54 by fangnagu
[考研] 材料與化工304求B區(qū)調(diào)劑 +4	邱gl 2026-03-26	7/350	2026-03-30 08:39 by 探123
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 275求調(diào)劑 +15	Micky11223 2026-03-25	20/1000	2026-03-29 20:44 by 唐沐兒
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

yxy657828314

[求助] udf編程問題 已有2人參與

» 猜你喜歡

mycc

【答案】應(yīng)助回帖

mycc

【答案】應(yīng)助回帖

匿名

小七工作室1

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] udf編程問題已有2人參與