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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

[交流] UDF的問題已有2人參與

UDF的內(nèi)容如下:
#include "udf.h"
#include <stdio.h>
#include "dynamesh_tools.h"

static int nozzle_tid = 10;
static real g_c = 9.81;
static real loc = 0.0;

static real initial_weight = 10.0;
static real burn_rate = 0.0;

static real current_time = -1.0;
static real thrust_prev = 0.0;
static real missile_velocity = 0.0;

static real
missile_weight (real time)
{
  return (initial_weight - burn_rate * time);
}

DEFINE_ON_DEMAND(reset_velocity)
{
  missile_velocity = 0.0;
  thrust_prev    = 0.0;
  current_time    = -1.0;
}

DEFINE_CG_MOTION(missile_1dof, dt, cg_vel, cg_omega, time, dtime)
{
#if 1
  FILE *fp;
#endif

  /* reset velocities */
  NV_S (cg_vel, =, 0.0);
  NV_S (cg_omega, =, 0.0);

  if (!Data_Valid_P ())
return;

/* Give rocket time to establish thrust for specified time */
  if (time < 0.1)
Message (" NO MOTION YET: TIME LESS THAN 0.1 SECONDS "

;

  if (time < 0.1)
return;

  /* update missile velocity only if we are at
   the next new time level */
  if ((time - 0.1*dtime) > current_time)
{
   Domain *domain;
   Thread *t;
   face_t f;
   real force0, force1, force;
   real area, total, v_avg, dv;
   real w0, w1, mass_flow;

   /* update time stamp */
   current_time = time;

   /* get nozzle exit thread (from predefined thread id) */
   domain = THREAD_DOMAIN (DT_THREAD ((Dynamic_Thread *)dt));
   if (NULLP (t = Lookup_Thread (domain, nozzle_tid)))
return;

   /* compute average exit velocity (weighted by area) */
   v_avg = 0.0;
   total = 0.0;
   begin_f_loop (f, t)
{
   area = NV_MAG (F_AREA_CACHE (f, t));

#if RP_3D
   v_avg += area * ND_MAG (F_U (f, t), F_V (f, t), F_W (f, t));
#else
   v_avg += area * ND_MAG (F_U (f, t), F_V (f, t), 0.0);
#endif
   total += area;
}

   end_f_loop (f, t)
   v_avg /= total;

/* subtract off missile_velocity to get relative velocity */
/* ???    v_avg = v_avg - missile_velocity; */

   /* compute thrust from mdot and v_e
   note: assume nozzle_exit is type mass-flow-exit  */
   if (THREAD_VAR(t).mfi.flow_spec == MASS_FLOW_TYPE)
mass_flow = THREAD_VAR(t).mfi.mass_flow;
   else
{
   real mass_flow;

   /* if mass_flux given as profile, then sum up
      area * mass_flux over face thread, else, mass_flux
      is constant */
   if (IS_PROFILE (THREAD_VAR (t).mfi.mass_flux))
      {
      begin_f_loop (f, t)
mass_flow = NV_MAG (F_AREA_CACHE (f, t)) *
            F_VAR (f, t, THREAD_VAR (t).mfi.mass_flux);
      end_f_loop (f, t)
      }
   else
      {
      mass_flow = THREAD_VAR(t).mfi.mass_flux.constant * total;
      }
#if RP_2D
   if (rp_axi)
      mass_flow *= 2.0 * M_PI;
#endif
}
   force = v_avg * F_VAR (0, t, THREAD_VAR(t).mfi.mass_flux) / total;

   /* compute change in velocity (use trapezoidal rule) */
   w0 = missile_weight (time - dtime);
   force0 = (thrust_prev - w0 * g_c) / w0;
   w1 = missile_weight (time);
   force1 = (force - w1 * g_c) / w1;
   dv = MAX (0.0, 0.5 * dtime * (force0 + force1));

   missile_velocity += dv;

   loc += missile_velocity * dtime;
   thrust_prev = force;

#if 1
  if( (fp=fopen("silo.dat","a"

) !=NULL )
{
   float check = F_FLUX(f,t)+F_GRID_FLUX(f,t);
   fprintf(fp,"%f %f %f    %f %f %f \n",
            time,missile_velocity,loc,v_avg,F_FLUX(f,t),F_GRID_FLUX(f,t));
   fclose(fp);
}
#endif

   Message ("time = %12.5e, x_vel = %12.5e, force = %12.5e, loc(m)= %1.5e\n",
         time, missile_velocity, 0.5*(force0 + force1), loc);
}

  /* set missile velocity */
  cg_vel[0] = missile_velocity;
}
編譯型UDF，build后的提示如下：
silo.c
..\..\src\silo.c(94) : error C2440: “=”: 無法從“Profile”轉(zhuǎn)換為“real”
連著出了好幾個問題，現(xiàn)在都已經(jīng)暈乎了，求助�。。。�！

回復(fù)此樓

» 猜你喜歡

283求調(diào)劑已經(jīng)有5人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有7人回復(fù)
數(shù)一英一271專碩（085401）求調(diào)劑，可跨已經(jīng)有4人回復(fù)
0856，材料與化工321分求調(diào)劑已經(jīng)有11人回復(fù)
0703化學(xué)求調(diào)劑已經(jīng)有8人回復(fù)
277跪求調(diào)劑已經(jīng)有9人回復(fù)
材料277求調(diào)劑已經(jīng)有6人回復(fù)
085602 307分求調(diào)劑已經(jīng)有7人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有7人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有6人回復(fù)

心情不爽！

1樓 2015-11-03 15:56:11

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

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

這么快就沉了，趕緊自己頂一下

贊一下

回復(fù)此樓

心情不爽！

2樓2015-11-03 16:48:11

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

ls315920

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 140.5
帖子: 64
在線: 20.9小時
蟲號: 4120687
注冊: 2015-10-05
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

貌似遇到大神，趕緊拜拜。。。

贊一下

回復(fù)此樓

虛心學(xué)習(xí)

3樓2015-11-03 18:44:31

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

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

大家給點建議吧，應(yīng)該主要是 if (THREAD_VAR(t).mfi.flow_spec == MASS_FLOW_TYPE)
mass_flow = THREAD_VAR(t).mfi.mass_flow;
else
這一句中的=兩邊有問題，但是不知道該怎么解決。

贊一下

回復(fù)此樓

心情不爽！

4樓2015-11-04 11:42:11

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

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

自己繼續(xù)頂吧，實在頭大。

贊一下

回復(fù)此樓

心情不爽！

5樓2015-11-04 15:55:15

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

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

遇到問題還是靠自己解決，靠論壇什么的實在是不靠譜。

贊一下

回復(fù)此樓

心情不爽！

6樓2015-11-05 08:53:32

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

死天

專家顧問 (著名寫手)

專家經(jīng)驗: +60
應(yīng)助: 399 (碩士)
金幣: 3648.9
紅花: 89
帖子: 1336
在線: 443.8小時
蟲號: 2643057
注冊: 2013-09-09
專業(yè): 流體力學(xué)
管轄: 仿真模擬

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

mass_flow = THREAD_VAR(t).mfi.mass_flow;這個有問題

THREAD_VAR(t).mfi.mass_flow的意思是，通俗的說THREAD_VAR(t)列表下的mfi列表下的mass_flow的變量取值，如果我沒有記錯的話（手頭沒有UDF，沒法查詢），這個變量的返回值不是可用于計算的值，而是fluent自定義的變量類型。

贊一下

回復(fù)此樓

7樓2015-11-05 12:21:16

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

木風(fēng)10086

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 117
紅花: 5
帖子: 53
在線: 40.6小時
蟲號: 1266670
注冊: 2011-04-15
性別: GG
專業(yè): 流體力學(xué)

引用回帖:

7樓: Originally posted by 死天 at 2015-11-05 12:21:16
mass_flow = THREAD_VAR(t).mfi.mass_flow;這個有問題

THREAD_VAR(t).mfi.mass_flow的意思是，通俗的說THREAD_VAR(t)列表下的mfi列表下的mass_flow的變量取值，如果我沒有記錯的話（手頭沒有UDF，沒法查詢），這 ...

對，確實是這個，我也是找了好久，在文件里一個個找，最后發(fā)現(xiàn)就是這個的問題。在老版本的時候能用，現(xiàn)在的不行了。我不知道想用這個量的話該怎么辦，我現(xiàn)在是用別量通過公式反算這個量了，直接提取還是不會。

贊一下

回復(fù)此樓

心情不爽！

8樓2015-11-11 15:10:00

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

死天

專家顧問 (著名寫手)

專家經(jīng)驗: +60
應(yīng)助: 399 (碩士)
金幣: 3648.9
紅花: 89
帖子: 1336
在線: 443.8小時
蟲號: 2643057
注冊: 2013-09-09
專業(yè): 流體力學(xué)
管轄: 仿真模擬

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

8樓: Originally posted by 木風(fēng)10086 at 2015-11-11 15:10:00
對，確實是這個，我也是找了好久，在文件里一個個找，最后發(fā)現(xiàn)就是這個的問題。在老版本的時候能用，現(xiàn)在的不行了。我不知道想用這個量的話該怎么辦，我現(xiàn)在是用別量通過公式反算這個量了，直接提取還是不會。...

THREAD_VAR(t).mfi.mass_flow只是個返回值，調(diào)用該值，下面應(yīng)該有個函數(shù)來通過該值直接調(diào)用庫

贊一下

回復(fù)此樓

9樓2015-11-11 18:42:32

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主木風(fēng)10086 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

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

[考研] 材料與化工考研調(diào)劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	13/650	2026-03-27 22:42 by Micky11223
[考研] 086000調(diào)劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 求調(diào)劑 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6	小小陳小小 2026-03-22	9/450	2026-03-26 15:58 by dick_runner
[考研] 309求調(diào)劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7	niko- 2026-03-22	7/350	2026-03-25 20:14 by qingfeng258
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 296求調(diào)劑 +4	汪！？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 求調(diào)劑 +6	研研，接電話 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	6/300	2026-03-24 15:03 by 余意卿
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike

24小時熱門版塊排行榜

[交流] UDF的問題 已有2人參與

» 猜你喜歡

[交流] UDF的問題已有2人參與