中文幫助fluent udf第二章語言基礎(chǔ)

第二章.UDFCC用#define用#include本章介紹了C語言的一些基本信息，這些信息對(duì)處理FLUENT的UDF很有C語言的初級(jí)介紹。本章不會(huì)介紹諸如while-do循環(huán)，聯(lián)合，遞歸，結(jié)構(gòu)以及讀寫文件的基礎(chǔ)知識(shí)。如果你對(duì)C語言不熟悉可以參閱C語言的相關(guān)書籍。如：/*ThisishowIputacommentinmyC 定義real為double型。因?yàn)榻忉尯瘮?shù)自動(dòng)作如此分配，所以使用在UDF中所有的floatdouble數(shù)據(jù)變量時(shí)使用real數(shù)據(jù)類型是很好的編程習(xí)慣。C程序中用語句#define來定義。最簡單的常數(shù)是十進(jìn)制整數(shù)（如：0，1，2）e的十進(jìn)制數(shù)IDYMINYMAX如下：#defineWALL_ID5#defineYMIN#defineYMAX值。變量在使用之前必須在C程序中。這樣，計(jì)算機(jī)才會(huì)提前知道應(yīng)該如何分配給相應(yīng)變量的類型。變量的結(jié)構(gòu)如下：首先是數(shù)據(jù)類型，然后是具有相應(yīng)類型的一個(gè)或多個(gè)C所允許的合法字符，變量名字中可以有字母，數(shù)字和下劃線。需要注意的是，在C程序中，字母是區(qū)分大小寫的。下面是變量的例子：int /*變量n為整型inti1, /*變量i1和i2為整型floattmax realaverage_temp ）mu_lam和temp是局部變量。）DEFINE_PROPERTY(cell_viscosity,cell,{realrealtemp=C_T(cell,if(temp>288.)mu_lam=5.5e-elseif(temp>mu_lam=143.2135-0.49725*temp;mu_lam=return}全局變量UDF源文件中是對(duì)所有的函數(shù)都起作用（調(diào)用一UDF源方便的處理方法就是在頭文件（.h）中加上extern的定義，然后在所有的.c文如下所示。注意：extern只用于編譯過的UDF。/* externrealvolume;/* #include"udf.h"realvolume;{volume=....}/* #include"udf.h"#include"filea.h"{/*用總數(shù)來計(jì)算每個(gè)單位體積的源項(xiàng)的代碼realtotal_source=...;realreturnsource;}static對(duì)于全局變量和局部變量的影響是不一樣的。靜態(tài)局部變量在函之外對(duì)任何函數(shù)保持不可見。靜態(tài)也可以用于函數(shù)，使該函數(shù)只對(duì)定義它的.c源文件保持可見。下面是靜態(tài)全局變量的例子。注意：extern只用于編譯過的UDF。#includestaticrealabs_coeff=1.0; /*吸收系數(shù)*/realsource;DEFINE_SOURCE(energy_source,c,t,dS,{intP1 dS[eqn]=-16.*abs_coeff*SIGMA_SBC*source=-abs_coeff*(4.*SIGMA_SBC*pow(C_T(c,t),4.)-C_UDSI(c,t,P1));returnsource;}DEFINE_SOURCE(p1_source,c,t,dS,{intP1=dS[eqn]=-source=abs_coeff*(4.*SIGMA_SBC*pow(C_T(c,t),4.)-C_UDSI(c,t,P1));returnsource;}注意：typedef只用于編譯過的UDF。typedefstruct{inta;intc;} /*mylist為類型結(jié)構(gòu)列表mylist /*x,y,z為類型結(jié)構(gòu)列表中的類型包括int，float等等，如下例所示：intx=realy=intz=x+((int) /*z=4real，void等，其相應(yīng)的返回值就是該數(shù)據(jù)類型，錄A的列表。!!C函數(shù)不能改變它們的，但是可以改變這些所指向的變量數(shù)組的定義格式為：名字[數(shù)組元素個(gè)數(shù)]，C數(shù)組的下標(biāo)是從零開始的。變inta[10],reala[0]=radii[4]=b[10][10]=;/*a為一個(gè)一維數(shù)組/*radii為一個(gè)一維數(shù)組/*b為一個(gè)二維數(shù)組2.10須以*開頭。指針廣泛用于提取結(jié)構(gòu)中的數(shù)據(jù)，以及在多個(gè)函數(shù)中通過數(shù)據(jù)例如：intip=*ip=inta=1;int*ip;ip= /*&aa的地址值printf("contentofaddresspointedtobyip=%d\n",*ip= /*a= printf("nowa=%d\n",數(shù)組的起始地址，在C中指針和數(shù)組具有緊密的聯(lián)系。2.10.1CFLUENT中，線程和域指針是UDF常用的自變量。當(dāng)你在UDF中指定這些自變量時(shí)，F(xiàn)LUENT解算器會(huì)UDF，從而使你的函數(shù)可以存取解算器的數(shù)據(jù)定（DEFINE_PROFILE宏來定義的）自定義輪UDF的自變量是一個(gè)指向應(yīng)用于邊界條件的線程的指針。DEFINE_PROFILE函數(shù)會(huì)存取線程指針?biāo)赶騣fif語句是條件控制語句的一種。格式為if(邏輯表達(dá)式{語句if(q!={a=0;b=if-elseif(邏輯表達(dá)式{語句if(x<{x=-y=}Y=X/50.X=-Y=forfor(起點(diǎn);終點(diǎn);增量{語句inti,j,n=10;for(i=1;i<=n;{j=printf("%d}N=10DOI=1,10J=I*IC 取++--減 不等Cx（y）的三角函數(shù)值。函數(shù)doubleacos(doubledoubleasin(doubledoubleatan(doubledoubleatan2(doublex,doubledoublecos(doubledoublesin(doubledoubletan(doubledoublecosh(doubledoublesinh(doubledoubletanh(doubledoublesqrt(doublexdoublepow(doublex,doublexdoubleexp(doubledoublelog(doublelnxdoublelog10(doublelog10xdoublefabs(doublexdoubleceil(doubledoublefloor(doubleFILE*fopen(char*filename,charintfclose(FILEintfprintf(FILE*fd,char*format,intprintf(char*format,intfscanf(FILE*fd,char*format,fscanfprintf是一般的輸除了fopen之外所有的函數(shù)都為整數(shù)，這是因?yàn)樵摵瘮?shù)所返回的整數(shù)會(huì)告訴在下面的例子中，需要打開的數(shù)據(jù)文件的名字引號(hào)括起來。fopen中的選項(xiàng)r表明該文件是以可讀形式打開的。fscanfd所指向的文件中讀入兩個(gè)浮點(diǎn)數(shù)并將它們?yōu)閒1和f2。關(guān)于C的標(biāo)準(zhǔn)輸入輸出函數(shù)其它的信息，你可FILEfd=fopen("data.txt","r"); opensafilenameddata.txt fscanf(fd,"%f,%f'',&f1,&f2);用#defineUDFC預(yù)處理程序命令。當(dāng)你使用#define宏置換命令，C預(yù)處理程序（如，cpp）執(zhí)行了一個(gè)簡單的置換，并用替換文本替換宏中定義#definemacrorecement-#defineRAD預(yù)處理程UDF中所有的變RAD出現(xiàn)的地方RAD替換1.2345。在#defineAREA_RECTANGLE(X,Y)你的UDF中所有的AREA_RECTANGLE(X,Y)都會(huì)被替換為(X)和(Y)的乘積用#include件時(shí)，C前處理程序會(huì)