ANSYS输出的单元和节点文件的C语言读入
ANSYS中使用EWRITE来输出当前选择的单元到文件。他的参数为
其中前两个分别是文件名和后缀，KAPPND表明是否清空单元文件来进行输出（用0和1 表示），Format则是说明输出文
件的位宽，其中SHORT是默认选项，表示用六位位宽输出，而LONG是表示8为位宽输出。
同样的，ANSYS使用NWRITE来输出当前选择的节点到文件。其输入参数与EWRITE的输入参数的含义相同,只不过没有位宽
的选项。但是节点的输出文件的中第一列为节点号位宽是8位，而第二到第四列为节点坐标，各20为位宽，各列之间没
有空格。
在使用C语言读入网格数据的时候，可以使用上述的文件来读入.通常比较方便的做法是知道节点个数，和单元个数然后
在依次读入C程序中。但是上述输出的单元和节点的文件并没有提供这个信息。
不过可已通过一下方式来产生这种文件：在APDL中先读取单元数和节点数，分别写入对应的单元和节点文件，然后以追
加方式写节点和单元文件。
下面是一个简单的用C语言读取二维的三角网格的例子（没有写入单元和节点的个数，而是通过计算行数来得到的）
1 #include <stdio.h>
2 #include <stdlib.h>
3 //ANSYS 's element output format
4
5 char NodeFormat[]="%*8c%20s%20s"; //节点的行格式
6 char ELongFormat[]="%8d%8d%8d%*40c%8d";//单元文件三角网格的长格式
7 char EShortFormat[]="%6d%6d%6d%*30c%6d";//单元文件三角网格的短格式
8
9
10
11
12 int main(int argc, const char *argv[])
13 {
14 FILE* fp;
15 char buffer[200];
16 char xx[21];
17 char yy[21];
18 int ch;
19 int mat,num;
20 int linenum=0;
21 int a[3]={0};
22 char NodeFile[]="nodelist1.n";//节点文件
23 char ElemFile[]="elemlist1.e";//单元文件
24
25 fp=fopen(NodeFile,"r");
26
27 //计算节点数目
28 while((ch=fgetc(fp))!=EOF)
29 {
30 if(ch=='\n')linenum++;
31 }
32 rewind(fp);
33 printf("节点数是%d\n",linenum);
34 //读入节点的坐标
35 while(fgets(buffer,200,fp)!=NULL)
36 {
37 xx[0]='\0';
38 yy[0]='\0';
39 sscanf(buffer,NodeFormat,xx,yy);
40 printf("%f,%f\n",atof(xx),atof(yy));
41 }
42 fclose(fp);
43
44 fp=fopen(ElemFile,"r");
45 int linenum=0;//计算单元的数目
46 while((ch=fgetc(fp))!=EOF)
47 {
48 if(ch=='\n')linenum++;
49 }
50 rewind(fp);
51 printf("单元数是%d\n",linenum);
52 //读入单元
53 while(fgets(buffer,200,fp)!=NULL)
54 {
55 sscanf(buffer,EshortFormat,&a[0],&a[1],&a[2],&mat);
56 printf("%8d %8d %8d %8d\n",a[0],a[1],a[2],mat);
57 }
58 fclose(fp);
59
60
61 return 0;
62 }
63
2 #include <stdlib.h>
3 //ANSYS 's element output format
4
5 char NodeFormat[]="%*8c%20s%20s"; //节点的行格式
6 char ELongFormat[]="%8d%8d%8d%*40c%8d";//单元文件三角网格的长格式
7 char EShortFormat[]="%6d%6d%6d%*30c%6d";//单元文件三角网格的短格式
8
9
10
11
12 int main(int argc, const char *argv[])
13 {
14 FILE* fp;
15 char buffer[200];
16 char xx[21];
17 char yy[21];
18 int ch;
19 int mat,num;
20 int linenum=0;
21 int a[3]={0};
22 char NodeFile[]="nodelist1.n";//节点文件
23 char ElemFile[]="elemlist1.e";//单元文件
24
25 fp=fopen(NodeFile,"r");
26
27 //计算节点数目
28 while((ch=fgetc(fp))!=EOF)
29 {
30 if(ch=='\n')linenum++;
31 }
32 rewind(fp);
33 printf("节点数是%d\n",linenum);
34 //读入节点的坐标
35 while(fgets(buffer,200,fp)!=NULL)
36 {
37 xx[0]='\0';
38 yy[0]='\0';
39 sscanf(buffer,NodeFormat,xx,yy);
40 printf("%f,%f\n",atof(xx),atof(yy));
41 }
42 fclose(fp);
43
44 fp=fopen(ElemFile,"r");
45 int linenum=0;//计算单元的数目
46 while((ch=fgetc(fp))!=EOF)
47 {
48 if(ch=='\n')linenum++;
49 }
50 rewind(fp);
51 printf("单元数是%d\n",linenum);
52 //读入单元
53 while(fgets(buffer,200,fp)!=NULL)
54 {
55 sscanf(buffer,EshortFormat,&a[0],&a[1],&a[2],&mat);
56 printf("%8d %8d %8d %8d\n",a[0],a[1],a[2],mat);
57 }
58 fclose(fp);
59
60
61 return 0;
62 }
63
