3.1.(0.01).c

混凝土结构杆件的非线性计算

#include <stdio.h>
#include <stdlib.h>
#include<math.h>
main(){
	int   e0=355;
	double  xc, s,s1,c,mc,elc,e0c,m,n,mx,el0=0.0022;
	FILE  *fp;
	fp=fopen("out(0.01).txt","w");
	for(elc=0.0001;elc<=0.0035;elc+=0.0001){
		xc=600;
		do{
			xc=xc-0.01;
			c=xc*(32665.948-1855831.659*elc-32665.948*log(1+445.324*elc)/(445.32*elc));
			mc=(32665.948-1855831.659*elc-32665.948*log(1+445.324*elc)/(445.324*elc))*xc
				*(565-xc)+(-73.3544+16333.2374*elc-1237211.63484*elc*elc+73.3544*
				log(1+445.324*elc)/(445.324*elc))*xc*xc/elc;
			if(elc/xc*(565-xc)<=0.002275)
				s=200000*elc/xc*(565-xc);
			else 
				s=455;
			if(elc/xc*(xc-35)<=0.002275)
				s1=200000*elc/xc*(xc-35);
			else 
				s1=455;
			
			n=c+s1*942-s*942;
            m=mc+s1*942*530;
			e0c=m/n-265;
			mx=n*355;

		}
		while((e0-e0c)>0.0001);
        fprintf(fp,"elc        xc        e0c         n            mx\n");
		n=n/1000;
	    mx=mx/1000000;
        fprintf(fp,"%.4lf  %.6lf  %.6lf  %.6lf  %.6lf\n",elc,xc,e0c,n,mx);
	}
	fclose(fp);
}