木風10086

新蟲 (小有名氣)

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

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

UDF的內容如下:
#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”轉換為“real”
連著出了好幾個問題，現在都已經暈乎了，求助�。。。�！

回復此樓

» 猜你喜歡

292求調劑已經有15人回復
304求調劑已經有5人回復
材料與化工考研調劑已經有17人回復
321求調劑已經有5人回復
329求調劑已經有5人回復
材料求調劑一志愿哈工大324 已經有8人回復
286求調劑已經有7人回復
083000學碩274求調劑已經有7人回復
311求調劑已經有8人回復
一志愿C9材料與化工專業(yè)總分300求調劑已經有9人回復

心情不爽！

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

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

木風10086

新蟲 (小有名氣)

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

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

贊一下

回復此樓

心情不爽！

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

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

ls315920

新蟲 (小有名氣)

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

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

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

贊一下

回復此樓

虛心學習

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

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

木風10086

新蟲 (小有名氣)

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

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

贊一下

回復此樓

心情不爽！

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

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

木風10086

新蟲 (小有名氣)

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

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

贊一下

回復此樓

心情不爽！

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

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

木風10086

新蟲 (小有名氣)

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

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

贊一下

回復此樓

心情不爽！

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

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

死天

專家顧問 (著名寫手)

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

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

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

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

贊一下

回復此樓

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

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

木風10086

新蟲 (小有名氣)

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

引用回帖:

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ā)現就是這個的問題。在老版本的時候能用，現在的不行了。我不知道想用這個量的話該怎么辦，我現在是用別量通過公式反算這個量了，直接提取還是不會。

贊一下

回復此樓

心情不爽！

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

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

死天

專家顧問 (著名寫手)

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

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

引用回帖:

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

THREAD_VAR(t).mfi.mass_flow只是個返回值，調用該值，下面應該有個函數來通過該值直接調用庫

贊一下

回復此樓

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

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

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

[考研] 一志愿廈門大學化學學碩307求調劑 +9	y7czhao 2026-03-26	9/450	2026-03-28 07:53 by Iveryant
[考研] 085404求調劑，總分309，本科經歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +3	迷人的哈哈 2026-03-28	3/150	2026-03-28 07:34 by ACS Nano——
[考研] 265求調劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 復試調劑，一志愿南農083200食品科學與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當當
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經歷，求調劑！ +3	Zuhui0306 2026-03-25	3/150	2026-03-27 14:09 by shangxh
[考研] 303求調劑 +7	安憶靈 2026-03-22	8/400	2026-03-27 11:46 by sanrepian
[考研] 085600材料與化工306 +10	z1z2z3879 2026-03-21	11/550	2026-03-27 11:31 by wangjy2002
[考研] 材料求調劑 +8	@taotao 2026-03-21	8/400	2026-03-27 11:21 by wangjy2002
[考研] 材料求調劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 352求調劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 297求調劑 +6	田洪有 2026-03-26	6/300	2026-03-26 15:55 by 不吃魚的貓
[考研] 一志愿河工大 081700 276求調劑 +4	地球繞著太陽轉 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 347求調劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 0854人工智能方向招收調劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 食品專碩一志愿雙一流 328 +3	xiaom99 2026-03-21	4/200	2026-03-24 21:20 by lailaisimei
[考研] 0703化學調劑，求導師收 +7	天天好運來上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 300求調劑，材料科學英一數二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 341求調劑(一志愿湖南大學070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte