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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

殺人者武三郎

新蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 151.2
散金: 16
紅花: 3
帖子: 154
在線: 40.5小時
蟲號: 2972508
注冊: 2014-02-17
專業(yè): 機構(gòu)學(xué)與機器人

[求助] 并行計算udf

昨天下午開始著手并行計算，遇到很多問題，從64位系統(tǒng)編譯udf，到單機多核計算的設(shè)置，晚上終于能算了。但是上午一看結(jié)果，動網(wǎng)格在動了一個周期之后就沒有運動了。這個udf是單核計算的udf，并行計算下需要修改嗎？udf如下：

#include <stdio.h>
#include "udf.h"
//調(diào)用fluent udf頭文件
//udf.h中的運動控制子程序
DEFINE_CG_MOTION(finright, dt, vel, omega, time, dtime)
// finright在udf里的名字
{
//# define PI 3.1415926

real freq=0.3;
real T = 1.0/freq;

//定義上下?lián)u翼前后平均值(初值)
real faiLC=0.0*M_PI/180;
real faiFEC=0.0*M_PI/180;
real faiFLC=0.0*M_PI/180;
//定義上下?lián)u翼前后幅值
real faiLA=20.0*M_PI/180;
real faiFEA=32.4406*M_PI/180;
real faiFLA=0.0*M_PI/180;
//上下相對搖翼的相位差
real dfaiL=0.0*M_PI/180;
//前后相對搖翼的相位差
real dfaiFL=0.0*M_PI/180;
//前后相對搖翼的相位差
real dfaiFE=90.0*M_PI/180;
//35度的角fai,翼型初始時y*軸與固定坐標系y軸的夾角
real fai=0.0*M_PI/180;
//faiL=faiLC-faiLAcos(w*t+dfaiL)當t=0時的faiL用faiLini表示
real faiLini;
//faiFL=faiFLC-faiFLAcos(w*t+dfaiFL)當t=0時的faiFL用faiFLini表示
real faiFLini;
//搖翼角速度
real omegaFE;
//上下拍翼角速度
real omegaL;
//前后拍翼角速度
real omegaFL;
//上下拍翼角度
// real faiL;
//搖翼角度
real faiFE;
//前后拍翼角度
real faiFL;



/* reset velocities */
//初始化速度和角速度
  NV_S (vel, =, 0.0);
  NV_S (omega, =, 0.0);
  // 定義
/*faiFE=faiFEC-faiFEA*cos(2*M_PI*freq*time+dfaiFE);
omegaFE=faiFEA*2*M_PI*freq*cos(2*M_PI*freq*time+dfaiFE);
omegaL=faiLA*2*M_PI*freq*cos(2*M_PI*freq*time+ dfaiL);*/
//
//faiLini=faiLC-faiLA*cos(dfaiL);

  if (!Data_Valid_P ())
return;

// 在0~T/4時間內(nèi)將faiL轉(zhuǎn)到faiLini，即先將翼轉(zhuǎn)到初始的位置
// 這里實際就是繞x星轉(zhuǎn)faiLini角，下面的程序就是要找出在固定坐標系下怎么轉(zhuǎn)得到這個結(jié)果
// x星=cosfai*x+sinfai*y
// x星的單位矢量=cosfai*i+sinfai*j
// w=wx星*x星的單位矢量
// 現(xiàn)在未知的是wx星，需要wx星在0~T/4內(nèi)積分=faiLini
// 試湊出wx星=faiLini*w*coswt
// 所以w=faiLini*w*coswt*cosfai*i+faiLini*w*coswt*sinfai*j
// 實際上認為t=0時,faiL=faiLini，但模型初始不是的
// T周期w=2pi/T
if (time >=0 && time<0.25*T)
{


         //x軸omega[0]，y軸omega[1]，z軸omega[2]
         faiFEini= faiFEC+faiFEA*sin(dfaiFE);
         faiFLini=faiFLC+faiFLA*sin(dfaiFL);

omega[0]= faiFLini*2.0*M_PI*freq*cos(2.0*M_PI*freq*time)*sin(faiFEini);
omega[1]= faiFEini*2.0*M_PI*freq*cos(2.0*M_PI*freq*time);
omega[2]= faiFLini*2.0*M_PI*freq*cos(2.0*M_PI*freq*time)*cos(faiFEini);


}
// 運動控制
if (time >=0.25*T)


         {
      faiFE=faiFEC+faiFEA*sin(2*M_PI*freq*(time-0.25*T)+ dfaiFE);
faiFL=faiFLC-faiFLA*sin(2*M_PI*freq*(time-0.25*T)+dfaiFL);
   omegaFE=faiFEA*2*M_PI*freq*cos(2*M_PI*freq*(time-0.25*T) + dfaiFE);
omegaFL=faiFLA*2*M_PI*freq*cos(2*M_PI*freq*(time-0.25*T)+dfaiFL);
omegaL=faiLA*2*M_PI*freq*cos(2*M_PI*freq*(time-0.25*T));

         omega[0]=-omegaFE*sin(fai)+omegaFL*cos(fai)*sin(faiFE)+omegaL*cos(faiFL)*cos(faiFE)*cos(fai)-omegaL*sin(faiFL)*sin(fai);
omega[1]=omegaFE*cos(fai)+omegaFL*sin(faiFE)*sin(fai)+omegaL*cos(faiFL)*cos(faiFE)*sin(fai)+omegaL*sin(faiFL)*cos(fai);
omega[2]=omegaFL*cos(faiFE)-omegaL*cos(faiFL)*sin(faiFE);

         }

}

回復(fù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

1樓 2014-03-26 15:27:04

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

fluent12.0單機并行計算已經(jīng)有10人回復(fù)
UDF并行編譯求助-附程序段已經(jīng)有3人回復(fù)
動網(wǎng)格的UDF問題已經(jīng)有21人回復(fù)
udf在非host節(jié)點下并行計算message 已經(jīng)有11人回復(fù)
open_udf_library:the system cannot find the file specified.如何解決？已經(jīng)有6人回復(fù)
關(guān)于fluentUDF的一個小問題！已經(jīng)有5人回復(fù)
單機多核并行計算下UDF的問題已經(jīng)有9人回復(fù)
一個UDF并行報錯已經(jīng)有17人回復(fù)
64位ANSYS FLUENT不能compile UDF的問題已經(jīng)有20人回復(fù)
UDF 改為并行計算已經(jīng)有9人回復(fù)
并行計算沒有輸出？怎么回事已經(jīng)有16人回復(fù)
單機多核一個cpu的串行并行求助已經(jīng)有26人回復(fù)
FLUENT讀取并行case自動刪除文件已經(jīng)有3人回復(fù)
【求助】ms5.5單機多核并行問題已經(jīng)有9人回復(fù)
【求助】單機多核UDF怎么并行計算已經(jīng)有15人回復(fù)
【討論】fluent12.0和12.1的并行速度已經(jīng)有9人回復(fù)
【求助】udf出錯已經(jīng)有13人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主學(xué)員ev57BJ 的主題更新