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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

焊接網(wǎng)格劃分

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 906.1
帖子: 24
在線: 39.2小時
蟲號: 2111966

[交流] 激光反射吸收熱源模型�。。。〈a重量爆料�。。。。∏蠼獯穑。。。。。。。。�

代碼有點錯誤，我現(xiàn)在是想用最小二乘法擬合自由表面的曲線，首先輪循所有自由表面網(wǎng)格，計算自由表面網(wǎng)格數(shù)量N，開辟數(shù)值存儲自由表面上網(wǎng)格的重心（優(yōu)良擬合三次曲線方程），gt,GS,max三個函數(shù)是用來求三次曲線方程y=a+bx+cx2+dx*x*x四個系數(shù)。然后把求出的四個系數(shù)保存到C_UDMI中，用于與離散的激光光線求交。代碼編譯之后有錯誤，可能思路上都已經(jīng)與udf編程思想向偏離了，求各路大神指點，重謝！�。。。。。。。。。≈辽�400金幣�。。。『呛�
#include "udf.h"

#include "sg.h"

#include "sg_mphase.h"

#include "flow.h"

#include <math.h>

double gt(double x[],double y[]);
double GS(double a[4][4],double b[4],double x[4]);
double max(double array[4]);

DEFINE_ADJUST(store_gradient,domain)
{
  int phase_domain_index = 0;
  real xc[ND_ND],xd[ND_ND],x[4],n1,n2,thita_l,thita_g;
  Thread *t;
  Thread **pt;
  cell_t c;

  Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,phase_domain_index);
  {
Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,Vof_Deriv_Accumulate);
  }
  mp_thread_loop_c(t,domain,pt)
{
  if (FLUID_THREAD_P(t))
  {
int N=0;
Thread *ppt = pt[phase_domain_index];

begin_c_loop (c,t)
{
   C_CENTROID(xd,c,t);
   if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xd[0]-0)>0.0005)
      N++;
}
end_c_loop (c,t)

double a[N],b[N],c[N],d[N],e[N],f[N],g[N];

int i;
for(i=0; i<N; i++)
{
   a=1;
}

int j=0;
begin_c_loop (c,t)
{
   C_CENTROID(xc,c,t);
   if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xc[0]-0)>0.0005)
   {
      b[j]=xc[0];
      e[j]=xc[1];
      c[j]=b[j]*b[j];
      d[j]=b[j]*b[j]*b[j];
j++;
   }
}
end_c_loop (c,t)
f[0]=gt(a,a);f[1]=gt(a,b);f[2]=gt(a,c);f[3]=gt(a,d);f[4]=gt(b,b);f[5]=gt(b,c);f[6]=gt(b,d);
f[7]=gt(c,c);f[8]=gt(c,d);f[9]=gt(d,d);g[0]=gt(a,e);g[1]=gt(b,e);g[2]=gt(c,e);g[3]=gt(d,e);
double array[4][4]={{f[0],f[1],f[2],f[3]},{f[1],f[4],f[5],f[6]},{f[2],f[5],f[7],f[8]},{f[3],f[6],f[8],f[9]}};
GS(array,g,x);
C_UDMI(c,t,0) = x[0];
C_UDMI(c,t,1) = x[1];
C_UDMI(c,t,2) = x[2];
C_UDMI(c,t,3) = x[3];
  }
}
  Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
}

double gt(double x[],double y[])
{
double sum=0;
int i;
for(i=0;i<N;i++)
   sum=x*y+sum;
return sum;
}

double GS(double a[4][4],double b[4],double x[4])
{
double c[4]={0};
double x0[4]={0};
int i,j,k;
double sum=0;
for(k=1;;k++)
{
      for(i=0;i<4;i++)
      {
         for(j=0;j<4;j++)
         {
            sum=a[j]*x0[j]+sum;
         }
         x=x0+(b-sum)/a;
      c=fabs(x-x0);
      x0=x;
      sum=0;
}
r=max(c);
if(r<0.0001)
{
         break;
   }
}
return 0;
}

double max(double array[4])
{
double a=array[0];
int i;
for(i=0;i<4;i++)
{
if(a<array)
   a=array;
}
return a;
}

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有191人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2013-12-02 10:11:07

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

qingxiang827

鐵蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 6323.4
帖子: 7245
在線: 49.4小時
蟲號: 2608424

★
焊接網(wǎng)格劃分(金幣+3): 謝謝參與

愿樓主心想事成

回復(fù)此樓

11樓2013-12-02 11:06:04

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

gitme3388

鐵蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 6337.6
帖子: 7309
在線: 49小時
蟲號: 2607600

★
焊接網(wǎng)格劃分(金幣+3): 謝謝參與

Good luck

回復(fù)此樓

12樓2013-12-02 11:06:09

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

我愛喝奶

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 260.4
帖子: 199
在線: 13.8小時
蟲號: 2828934

★
焊接網(wǎng)格劃分(金幣+3): 謝謝參與

不太了解祝福樓主吉星高照

贊一下

回復(fù)此樓

20樓2013-12-02 11:26:38

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

crysmat

鐵桿木蟲 (知名作家)

應(yīng)助: 33 (小學生)
金幣: 10713.8
帖子: 8155
在線: 223.1小時
蟲號: 596253

★
焊接網(wǎng)格劃分(金幣+3): 謝謝參與

頂一個，祝樓主一切順利

贊一下

回復(fù)此樓

22樓2013-12-02 11:37:53

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

焊接網(wǎng)格劃分

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 906.1
帖子: 24
在線: 39.2小時
蟲號: 2111966

蟲友們！給大家提個方向去修改，可以用FORTRAN編寫函數(shù)gt,GS,max,先預(yù)編譯成.obj文件然后在ADJUST宏中調(diào)用。希望高手們多多指教�。。。。。。�！

贊一下

回復(fù)此樓

37樓2013-12-03 10:27:32

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

ttmany7788

鐵蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 6408.5
帖子: 7308
在線: 50.2小時
蟲號: 2607370

★
焊接網(wǎng)格劃分(金幣+1): 謝謝參與

愿樓主心想事成

回復(fù)此樓

42樓2013-12-03 11:03:56

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

freedom767

鐵蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 6425.3
帖子: 7369
在線: 50.3小時
蟲號: 2606677

★
焊接網(wǎng)格劃分(金幣+1): 謝謝參與

愿樓主心想事成

回復(fù)此樓

43樓2013-12-03 11:04:57

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

zhangwei19

榮譽版主 (文學泰斗)

★
焊接網(wǎng)格劃分(金幣+1): 謝謝參與

不會

回復(fù)此樓

48樓2013-12-03 11:23:48

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

簡單回復(fù)

湮沒黑洞2樓

2013-12-02 10:21 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

051031113樓

2013-12-02 10:25 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

風雨飄然4樓

2013-12-02 10:29 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

virologik5樓

2013-12-02 10:30 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

haixiawu6樓

2013-12-02 10:30 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

tzynew7樓

2013-12-02 10:33 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

slz20088樓

2013-12-02 10:34 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

yaoyao1459樓

2013-12-02 10:48 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

wlxiao10樓

2013-12-02 10:57 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

峰鋒風逢烽13樓

2013-12-02 11:08 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

doctorbang14樓

2013-12-02 11:14 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

wullww15樓

2013-12-02 11:18 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

aboluo19833116樓

2013-12-02 11:23 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

pzchen17樓

2013-12-02 11:24 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

xianzhi201018樓

2013-12-02 11:24 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

fujiating19樓

2013-12-02 11:24 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

wendytxy21樓

2013-12-02 11:33 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

abczhang23樓

2013-12-02 11:41 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

Spinage24樓

2013-12-02 11:41 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

頂

chlq314101a25樓

2013-12-02 11:43 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

宋全偉26樓

2013-12-02 11:44 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

， [ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

會飛的烏龜--27樓

2013-12-02 11:44 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

wander051528樓

2013-12-02 11:45 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

yuzhichao29樓

2013-12-02 11:50 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

祝福

Sammy230樓

2013-12-02 11:52 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

doublehaohao31樓

2013-12-02 11:58 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

塵軒無夢32樓

2013-12-02 11:59 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

、 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

lgward33樓

2013-12-02 12:01 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

huiminchlq34樓

2013-12-02 12:01 回復(fù)

焊接網(wǎng)格劃分(金幣+3): 謝謝參與

xiejf35樓

2013-12-02 12:09 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

silence136樓

2013-12-02 12:25 回復(fù)

ts_chenchao38樓

2013-12-03 11:00 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

天涯201239樓

2013-12-03 11:02 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

neu23440樓

2013-12-03 11:02 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

miaojiabing41樓

2013-12-03 11:03 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

gyfwon44樓

2013-12-03 11:09 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

gander45樓

2013-12-03 11:09 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

liu12333846樓

2013-12-03 11:11 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

lefantian12247樓

2013-12-03 11:18 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

zhoupeng8749樓

2013-12-03 11:32 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

base0803030850樓

2013-12-03 11:34 回復(fù)

焊接網(wǎng)格劃分(金幣+1): 謝謝參與

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主焊接網(wǎng)格劃分的主題更新

返回列表

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

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

[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 284求調(diào)劑 +5	Zhao anqi 2026-03-22	5/250	2026-03-22 17:38 by barlinike
[考研] 293求調(diào)劑 +12	zjl的號 2026-03-16	17/850	2026-03-22 16:51 by i_cooler
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 一志愿北京化工大學070300 學碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 一志愿深大，0703化學，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學員8dgXkO
[基金申請] 學校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 083200學碩321分一志愿暨南大學求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 南京大學化學376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂

24小時熱門版塊排行榜

[交流] 激光反射吸收熱源模型�。。。〈a重量爆料�。。。。∏蠼獯穑。。。。。。。。�

» 猜你喜歡

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

[交流] 激光反射吸收熱源模型�。。。〈a重量爆料�。。。。∏蠼獯穑。。。。。。。。�

» 搶金幣啦！回帖就可以得到:

查看全部散金貼