开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 界面编程 > 查看源码
MNQM.C
资源名称：C.rar　[点击查看]
上传用户：qq5388545
上传日期：2022-07-04
资源大小：29849k
文件大小：7k
源码类别：

界面编程

开发平台：

C/C++

MNQM.C：源码内容
  1. #include "string.h"
  2. #include "alloc.h"
  3. #include "stdio.h"
  4. #include "stdlib.h"
  5. #include "io.h"
  6. int nelement,nnode,nsum_node;
  7. struct elems
  8. {
  9. int nedge;
  10. float x[4],y[4];
  11. float s11[4],s22[4],s12[4];
  12. }*elem_stress;
  13. void main(argc,argv)
  14.   int argc; char *argv[];
  15.  {
  16.  FILE *fs,*fc,*fj;
  17.  int i,j,nelement,nnode,l,have=0,num,numk,numl,fast=0,ielement,nmat,gap=0;
  18.  char str[100],coname[80],cname[80],sname[80],str1[20];
  19.  float y1,z1;
  21.  float *coora[2];
  23.  int *element[4];
  24.  struct coors{
  25.     int numj;
  26.     float y,z;
  27. /*    float s11,s22,s12;*/
  28.     }*coor_stress;
  31.  if(argc<4)
  32. {
  33. for(i=0;i<3;i++)
  34.  {
  35.  switch(i)
  36. {
  37. case 0:strcpy(str1,"control");break;
  38. case 1:strcpy(str1,"element");break;
  39. case 2:strcpy(str1,"stress");break;
  40. }
  41.  for(;;)
  42. {
  43. /*/gotoxy(wherex()-15,wherey()-1);*/
  44. printf("input %s data filename (<CR> TO QUIT):",str1);
  45. gets(str);
  46. if(str[0]=='')
  47. /*/exit(0);*/
  48. return;
  49. else if(access(str,0)==0)
  50.   {
  51. if(i!=2) break;
  52.   else
  53. {
  54. /*/gotoxy(wherex()-15,wherey());*/
  55. printf("               ");
  56. break;
  57. }
  58.   }
  59. else printf("afile not found!");
  60. }
  61.  switch(i)
  62. {
  63. case 0:strcpy(coname,str);break;
  64. case 1:strcpy(cname,str);break;
  65. case 2:strcpy(sname,str);break;
  66. }
  67.  }
  68. }
  69.  else
  70. {
  71. strcpy(coname,argv[1]);
  72. strcpy(cname,argv[2]);
  73. strcpy(sname,argv[3]);
  74. }
  76. /*/ clrscr();*/
  77.  /*read nnodem,nelement*/
  79.  if((fj=fopen(coname,"r"))==NULL)
  80. {
  81. printf("aopen failed");
  82. exit(0);
  83. }
  85.  fscanf(fj,"%d%d",&nnode,&nelement);
  86.  fclose(fj);
  88.  if((fj=fopen(cname,"r"))==NULL) /*cooresponding node and element file*/
  89. {
  90. printf("aopen failed");
  91. exit(0);
  92. }
  93.   if((fc=fopen(cname,"r"))==NULL) /*corradate file*/
  94. {
  95. printf("aopen failed");
  96. exit(0);
  97. }
  98.  if((fs=fopen(sname,"r"))==NULL) /*stress file*/
  99. {
  100. printf("aopen failed");
  101. exit(0);
  102. }
  105.  /*allocate memory*/
  106.  fgets(str,100,fj);
  107.  fscanf(fj,"%d",&nsum_node);
  108.  fgets(str,100,fj);
  110.  for(i=0;i<2;i++)
  111. {
  112. if((coora[i]=(float *)malloc(sizeof(float)*nsum_node))!=NULL)
  113. fast=1;
  114. else
  115. {
  116. printf("antoo many node; Do you want to continue(y/n)");
  117. if((getchar())=='y')
  118. fast=0;
  119. else
  120. exit(0);
  121. }
  122. }
  123.  if(fast)
  124. {
  125. printf("nreading cooradation information......");
  126. for(i=0;i<nsum_node;i++)
  127. {
  128. float y1,z1;
  129. if(fscanf(fj,"%d",&num)==EOF)
  130. {
  131. break;
  132. }
  133. fgets(str,30,fj);
  134. fscanf(fj,"%f%f",&y1,&z1);
  135. coora[0][i]=y1;coora[1][i]=z1;
  136. /*              printf(" %d, %.5f, %.5fn",i,coora[0][i],coora[1][i]);*/
  137. fgets(str,100,fj);
  138. }
  139. }
  140.  if(!fast)
  141. {
  142. if((coor_stress=(struct coors *)malloc(sizeof(struct coors)*nnode))==NULL)
  143. {printf("anot enough memory (coor_stress)");
  144.  exit(0);
  145. }
  146. }
  148.  for(i=0;i<4;i++)
  149.  if((element[i]=(int *)malloc(sizeof(int)*nelement))==NULL)
  150. {printf("anot enough memory (element %d)",i);
  151.  exit(0);
  152.  }
  156.   nnode=0;
  158. if(fast)
  159. {
  160. fscanf(fj,"%d",&ielement);
  161. fscanf(fj,"%d",&ielement);
  162. fscanf(fj,"%d",&nmat);
  163. fgets(str,100,fj);
  164. for(i=0;i<nmat;i++)     /*nmat number of material class*/
  165. {
  166. printf("nreading material %d information of group 1",i+1);
  167. for(j=0;j<3;j++)
  168. fgets(str,100,fj);
  169. }
  170. printf("nreading element information of group 1");
  171. for(i=0;i<ielement+4;i++)
  172. {
  173. fgets(str,100,fj);
  174. }
  176. fscanf(fj,"%d",&ielement);
  177. if(ielement==14)        /*exist gap element*/
  178. {
  179. gap=1;
  180. printf("nreading gap data......");
  181. fscanf(fj,"%d",&ielement);
  182. fgets(str,100,fj);
  183. for(i=0;i<ielement+6;i++)
  184. fgets(str,100,fj);
  185. }
  187. if(gap)
  188. fscanf(fj,"%d",&ielement);
  189. fscanf(fj,"%d",&ielement);
  190. for(i=0;i<8;i++)
  191. fgets(str,100,fj);
  192. }
  193. else
  194. {
  195. rewind(fj);
  196. fgets(str,100,fj);
  198. fscanf(fj,"%d",&ielement);
  199. fgets(str,100,fj);
  200. for(i=0;i<ielement+1;i++)
  201. fgets(str,100,fj);
  203. fscanf(fj,"%d",&ielement);
  204. fscanf(fj,"%d",&ielement);
  205. fscanf(fj,"%d",&nmat);
  206. fgets(str,100,fj);
  207. for(i=0;i<nmat;i++)     /*nmat number of material class*/
  208. for(j=0;j<3;j++)
  209. fgets(str,100,fj);
  210. for(i=0;i<ielement+4;i++)
  211. fgets(str,100,fj);
  213. fscanf(fj,"%d",&ielement);
  214. fscanf(fj,"%d",&ielement);
  215. fgets(str,100,fj);
  216. for(i=0;i<ielement+6;i++)
  217. fgets(str,100,fj);
  219. fscanf(fj,"%d",&ielement);
  220. fscanf(fj,"%d",&ielement);
  221. for(i=0;i<8;i++)
  222. fgets(str,100,fj);
  223. }
  224.   if(ielement!=nelement)
  225. {
  226. printf("anlining element number error!");
  227. exit(0);
  228. }
  230.   printf("nreading element information of group 2");
  232.   for(ielement=0;ielement<nelement;ielement++)
  233.   {
  234.    if(fscanf(fj,"%d",&i)==EOF) break;
  235.    for(i=0;i<4;i++)
  236. {
  237. fscanf(fj,"%d",&element[i][ielement]);
  238. if(!fast)
  239.   {
  240.   coor_stress[nnode].numj=element[i][ielement];
  241.   have=0;
  242.   /*decide the existance  this nodes cooradates:YES:have=1;NO: have=0*/
  243.   for(l=0;l<nnode;l++)
  244. {
  245.  numk=(coor_stress[nnode].numj);
  246.  numl=(coor_stress[l].numj);
  247.  if(numk==numl)
  248.    {
  249.     have=1;
  250.     break;
  251.     }
  252.  }
  253.   /*if no this node's  cooradates,then search and record it*/
  254.   if(!have)
  255.       {
  256.       for(;;)
  257.   {
  258.   /*read node */
  259. fgets(str,100,fc);
  260. fgets(str,100,fc);
  261.   if(fscanf(fc,"%d",&num)==EOF)
  262. {
  263. rewind(fc);
  264. fgets(str,100,fc);
  265. fgets(str,100,fc);
  266. fscanf(fc,"%d",&num);
  267. }
  268.  /*decide if needed*/
  269.  if(num==coor_stress[nnode].numj)
  270.    {
  271.    fgets(str,30,fc);
  272.    fscanf(fc,"%f%f",&y1,&z1);
  273. /*   printf(" %d, %.5f, %.5fn",coor_stress[nnode].numj,y1,z1);*/
  274.    (coor_stress[nnode].y)=y1;(coor_stress[nnode].z)=z1;
  275.    fgets(str,100,fc);
  276.    nnode++;
  277.    break;
  278.    }
  279. /*skip the rest string*/
  280.  else fgets(str,100,fc);
  281.  }/*end search*/
  282.        }/*end record*/
  283.       }
  284.   }/*end of the current element*/
  285.    fgets(str,100,fj);
  286.    }/*end of all element*/
  287.  /*printf("%d %d",nnode,ielement);
  288.  getch();                         */
  289.  fclose(fj);
  290.  fclose(fc);
  291.  if((elem_stress=(struct elems *)malloc(sizeof(struct elems)*nelement))==NULL)
  292. {printf("annot enough memory");
  293.  exit(0);
  294.  }
  295.  rewind(fs);
  296.  numk=0;
  297.  printf("nsearching stress data.....");
  298.  for(;;)
  299. {
  300. fgets(str,100,fs);
  301. if(stricmp(str,"               Sigma-Intn")==0)
  302. numk++;
  303. if(numk==2) break;
  304. }
  305.   fgets(str,100,fs);
  307.   printf("nreading stress data......");
  309.  for(numk=0;numk<nelement;numk++)
  310. {
  311. if(element[2][numk]==element[3][numk]) have=3;
  312. else have=4;
  314. elem_stress[numk].nedge=have;
  316. for(i=0;i<have;i++)
  317. {
  318. int inode,in;
  319. inode=element[i][numk];
  320. if(!fast)
  321.   {
  322.   for(in=0;in<nnode;in++)
  323. {
  324. if(inode==coor_stress[in].numj)
  325. {
  326. elem_stress[numk].x[i]=coor_stress[in].y;
  327. elem_stress[numk].y[i]=coor_stress[in].z;
  328. break;
  329. }
  330. }
  331.     }
  332. else
  333.    {
  334.     elem_stress[numk].x[i]=coora[0][inode-1];
  335.     elem_stress[numk].y[i]=coora[1][inode-1];
  336.    }
  337. {
  338. printf(" %d %d",numk,elem_stress[numk].nedge);
  339. printf(" %.2f %.2f",elem_stress[numk].x[i],elem_stress[numk].y[i]);
  340.  }
  341. fscanf(fs,"%d",&numl);
  342. if(numl!=numk+1)
  343. {
  344. printf("anstress data file read error!");
  345. exit(0);
  346. }
  347. fgets(str,9,fs);
  348. fscanf(fs,"%f%f",&y1,&z1);
  350. elem_stress[numk].s11[i]=y1;
  351. elem_stress[numk].s22[i]=z1;
  352. printf(" %f %f",elem_stress[numk].s11[i],elem_stress[numk].s22[i]);
  353. fscanf(fs,"%f%f",&y1,&z1);
  355. elem_stress[numk].s12[i]=z1;
  356. printf(" %f",elem_stress[numk].s12[i]);
  357. fgets(str,100,fs);
  358. fgets(str,100,fs);
  359. printf("n");
  360. }
  361. }
  362.  fclose(fs);
  364.  if((fj=fopen("test.dat","wb"))==NULL)
  365. {
  366. printf("aopen failed");
  367. exit(0);
  368. }
  369.  fwrite(&nelement,sizeof(int),1,fj);
  370.  fwrite(elem_stress,sizeof(struct elems)*nelement,1,fj);
  371.  fclose(fj);
  373.  if(!fast)
  374. free((void *)elem_stress);
  375.  else
  376. {
  377. for(i=0;i<2;i++)
  378. free((void *)coora[i]);
  379. }
  380.  for(i=0;i<4;i++)
  381. free((void *)element[i]);
  383.  free((void *)coor_stress);
  384.  }