本文探討了使用ANSYS APDL和C語言實現混合編程的相關內容。

　　本文討論的不是利用C語言為ANSYS寫擴展(或者說是用戶子程序)，而是討論結合APDL對ANSYS的現有的計算程序進行復用，例如對ANSYS強大的前后處理程序，包括網格剖分，計算結果顯示等，或者利用有限元計算部分的模塊，使用ANSYS方便的邊界條件加載，方程組求解等。

　　0. 為什么進行混合編程

　　既然ANSYS已經提供了APDL這么強大的編程語言，為什么還要用C語言來與APDL混合編程呢?下面的幾個原因，首先是前處理建模的要求，APDL并不能應對實際中非常復雜的模型的構建，本文的出發點之一就是復合材料分析中的幾何模型構建問題，由于材料中的增強體為隨機的顆粒，每一次模擬中，顆粒都是隨機生成的，使用APDL并不方便;其次就是利用ANSYS的計算資源的考慮，比如上述問題，就可以利用C語言來生成模型的數據，利用APDL來構建幾何模型，以及剖分有限元網格，用于自己的計算程序;還有在一些問題中，其邊界條件，載荷等也是通過程序算出來的話，也可以通過混合編程對問題求解;另外，可以發揮C語言的優勢，比如說計算效率，現有的庫等，總而言之，混合編程能更充分利用ANSYS的計算程序以及C語言的優勢，并且可以實現計算機上的自動化計算。

　　1. 混合編程的原理

　　ANSYS與C語言混合編程的原理很簡單，實際上就是用C語言生成ANSYS的apdl命令流文件，然后通過命令行調用ANSYS進行計算。這個過程同樣可以使用其他的編程語言(比如FORTRAN，Python等)來進行，可以發揮各種語言的優勢以及對現有的程序庫的使用。

　　2. 如何使用C語言的實現

　　知道了混合編程的原理，那么就很簡單了，其過程無非就是利用C語言打開一個文件，寫入命令流，關閉文件，調用ANSYS進行計算。

　　文件操作需要使用的函數有

　　1 fp=fopen("mac.inp","r"); //打開文件

　　2 fprintf(fp, "APDL_Commands,%f",SomeData); //寫入文件

　　3 fclose(fp); //關閉文件

　　調用ANSYS的命令行計算則使用

　　system("D:\Progra~1\AnsysI~1\v100\ANSYS\bin\intel\ansys100 -b -i temp.inp -o mesh.log");

　　3. 實例

　　下面給幾個APDL和C語言混合編程的實例

　　3.1 復合材料中的圓性顆粒網格剖分

　　圖中的顆粒是隨機C語言中隨機生成的圓形，結果如圖所示：

　　相應的生成上述網格的函數為

　　1 /*調用ANSYS生成網格*/

　　2 void GenMesh(struct circle0 *circle,double xmax,double xmin,double ymax,double ymin)

　　3 {

　　4

　　5 FILE* ansysfile;

　　6 ansysfile=fopen("temp.inp","w");

　　7

　　8 fprintf(ansysfile,"/CLEAR,NOSTART");

　　9 fprintf(ansysfile,"/PREP7");

　　10

　　11 //創建基體

　　12 fprintf(ansysfile,"RECTNG,%f,%f,%f,%f,",xmin,xmax,ymin,ymax);

　　13

　　14 //創建顆粒

　　15 for (int i=0; i

　　16 {

　　17 if (circle[i].r<0.00000001)

　　18 {

　　19 continue;

　　20 }

　　21 fprintf(ansysfile,"K,%d,%f,%f",i+1000,

　　22 circle[i].centre_x,circle[i].centre_y);

　　23 fprintf(ansysfile,"CIRCLE,%d,%f",i+1000,circle[i].r);

　　24 }

　　25 fprintf(ansysfile,"ASBL,1,ALL");

　　26 fprintf(ansysfile,"AGLUE,All");

　　27

　　28 //創建單元類型

　　29 fprintf(ansysfile,"ET,1,PLANE42");

　　30 //創建材料

　　31 fprintf(ansysfile,"MP,DENS,1,0.1");//材料一

　　32 fprintf(ansysfile,"MP,DENS,2,10");//材料二

　　33 //為基體賦一號材料值

　　34 fprintf(ansysfile,"ASEL,S,AREA,,ARINQR(0,14)");

　　35 fprintf(ansysfile,"AATT,1");

　　36

　　37 //為顆粒賦二號材料值

　　38 fprintf(ansysfile,"ASEL,ALL");

　　39 fprintf(ansysfile,"ASEL,U,AREA,,ARINQR(0,14)");

　　40 fprintf(ansysfile,"AATT,2");

　　41 fprintf(ansysfile,"ALLSEL,ALL");

　　42

　　43 //生成網格

　　44 //網格尺寸

　　45 fprintf(ansysfile,"SMRT,3");

　　46 //fprintf(ansysfile,"AESIZE,ALL,0.6,

　　");

　　47

　　48 fprintf(ansysfile,"MSHAPE,1,2D");

　　49 fprintf(ansysfile,"MSHKEY,0");

　　50 //開始劃分網格

　　51 fprintf(ansysfile,"AMESH,ALL");

　　52

　　53 //輸出網格

　　54 // fprintf(ansysfile,"/INPUT,ExportMesh,inp

　　");

　　55

　　56 fclose(ansysfile);

　　57 //---------------調用ANSYS生成網格

　　58 system("D:\Progra~1\AnsysI~1\v100\ANSYS\bin\intel\ansys100 -b -i temp.inp -o mesh.log");

　　59 // system("tec360 tecplot.plt");

　　60 }

　　3.2 復合材料中的圓性顆粒網格剖分

　　對于橢圓形顆粒，網格剖分結果如圖所示

　　生成上述網格的代碼為

　　1

　　2 //生成顆粒是橢圓時的網格

　　3 //橢圓的數據結構是下面的struct

　　4 //Genmesh函數是調用ANSYS生成網格的函數

　　5 void GenMesh(struct elliptic0 *ellipse,int ellipsenum,double xmax,double xmin,double ymax,double ymin)

　　6 {

　　7

　　8 FILE* ansysfile;

　　9 ansysfile=fopen("temp.inp","w");

　　10 fprintf(ansysfile,"FINISH");

　　11 fprintf(ansysfile,"/CLEAR,NOSTART");

　　12 fprintf(ansysfile,"/PREP7");

　　13

　　14 //創建基體

　　15 fprintf(ansysfile,"RECTNG,%f,%f,%f,%f,",xmin,xmax,ymin,ymax);

　　16

　　17 //創建橢圓顆粒

　　18 for (int i=0; i

　　19 {

　　20 //把工作平面移到橢圓圓心。

　　21 fprintf(ansysfile,"wpave,%f,%f",ellipse[i].center_x,ellipse[i].center_y);

　　22 //旋轉工作平面

　　23 fprintf(ansysfile,"wprota,%f",180.0/pi*ellipse[i].theta);

　　24 //激活當前活動坐標系為工作平面

　　25 fprintf(ansysfile,"csys,4");

　　26 //生成上半個橢圓

　　27 fprintf(ansysfile,"k,%d,%f",2*i+1000,ellipse[i].a);

　　28 fprintf(ansysfile,"k,%d,%f",2*i+1001,-ellipse[i].a);

　　29 fprintf(ansysfile,"cswpla,%d,2,%f",2000+2*i,(ellipse[i].b/ellipse[i].a));

　　30 fprintf(ansysfile,"l,%d,%d",2*i+1000,2*i+1001);

　　31 fprintf(ansysfile,"reflect_line=_return");//記錄線的編號

　　32 fprintf(ansysfile,"csys,4"); //切換活動坐標系為工作平面

　　33 fprintf(ansysfile,"lsymm,y,reflect_line");//反射生成橢圓

　　34 //調整工作平面使其與整體坐標系重合

　　35 fprintf(ansysfile,"wprota,%f",-180.0/pi*ellipse[i].theta);

　　36 fprintf(ansysfile,"CSYS,0");

　　37 fprintf(ansysfile,"WPAVE,0,0,0");

　　38 }

　　39

　　40 fprintf(ansysfile,"ASBL,1,ALL");

　　41 fprintf(ansysfile,"AGLUE,All");

　　42

　　43 //創建單元類型

　　44 fprintf(ansysfile,"ET,1,PLANE42");

　　45 //創建材料

　　46 fprintf(ansysfile,"MP,DENS,1,0.1");//材料一

　　47 fprintf(ansysfile,"MP,DENS,2,10");//材料二

　　48 //為基體賦一號材料值

　　49 fprintf(ansysfile,"ASEL,S,AREA,,ARINQR(0,14)");

　　50 fprintf(ansysfile,"AATT,1");

　　51 //為顆粒賦二號材料值

　　52 fprintf(ansysfile,"ASEL,ALL");

　　53 fprintf(ansysfile,"ASEL,U,AREA,,ARINQR(0,14)");

　　54 fprintf(ansysfile,"AATT,2");

　　55 fprintf(ansysfile,"ALLSEL,ALL");

　　56 //生成網格

　　57 //網格尺寸

　　58 fprintf(ansysfile,"SMRT,3");

　　59 //fprintf(ansysfile,"AESIZE,ALL,0.6,

　　");

　　60

　　61 fprintf(ansysfile,"MSHAPE,1,2D");

　　62 fprintf(ansysfile,"MSHKEY,0");

　　63 //劃分網格

　　64 fprintf(ansysfile,"AMESH,ALL");

　　65

　　66 //輸出網格

　　67 fprintf(ansysfile,"/INPUT,ExportMesh,inp");//輸出網格到Tecplot

　　68 fprintf(ansysfile,"EWRITE,elemlist1,txt");//輸出單元

　　69 fprintf(ansysfile,"NWRITE,nodelist1,txt");//輸出節點

　　70

　　71 fclose(ansysfile);

　　72 //---------------調用ANSYS生成網格

　　73 system("D:\Progra~1\AnsysI~1\v100\ANSYS\bin\intel\ansys100 -b -i temp.inp -o mesh.log");

　　74 system("tec360 tecplot.plt");

　　75 }