非线性源程序

1、#include #include #include #include void main()int i,ii,j,k,sn=30,jmax;int concrete=0,sbar=0,pbar=0; /混凝土、钢筋破坏标志int J_kailie,J_xiaoya,J_chushi,J_squfu,J_pqufu,J_pohuai;float s0,s01,s02,dc,ds; /截面中线应变，S开头表应变，f开头表应力float scy,scu; /混凝土屈服应变、极限压应变float ssy,ssu; /非预应力钢筋屈服应变、极限应变float spy,spu,fps,sp; /预应力钢

2、筋屈服应变、极限应变、应力、应变float sci,fci,fc; /混凝土条应变、应力和抗压强度float fpy,fpu,fpe; /预应力筋屈服强度、极限强度、有效预应力float ss,fss,fsy; /非预应力钢筋应变、应力、屈服强度float ep,ep1; /预应力钢筋弹性模量float es,ec; /非预应力钢筋和混凝土弹性模量float em,sf1,sf2,mm,mic,mis,mip,ffc,ffp,ffs,dsf;float c500,m500,ao2;float h=600,b=300,aas,as1=0,as2=402,ap=784,dd=180,z;FILE

3、*file1;/*赋初值（根据算例）*for(i=0;i500;i+)ci=0.0;mi=0.0; j=0;ii=0;fc=35;fsy=400;fpe=1000;ep=200000;ec=28000;fpy=1540;fpu=1800;es=200000;scy=-0.002;scu=-0.0033;ssy=0.002;ssu=0.01;spy=fpy/ep;spu=0.05;ep1=(fpu-fpy)/(spu-spy);/*M-曲线*dc=-0.0000001;s0=-0.0001;ds=0.00001;doj+;cj=cj-1+dc;doii+;/*求混凝土应变应力* ffc=0.0;

4、mic=0.0;for(int i=0;i=0)fci=0.0; else if(sciscy)fci=-fc*(2*sci/scy-(sci/scy)*(sci/scy); else if(sciscu)fci=-fc; else concrete=1; ffc=ffc+fci*b*h/sn;mic=mic+fci*b*h/sn*z;/*求非预应力钢筋应变应力*ffs=0.0;mis=0.0;for(int k=0;k2;k+)z=aok-h/2;ss=s0+z*cj;if(abs(ss)=ssy)fss=es*ss;else if(abs(ss)=ssu)fss=ss/abs(ss)*fs

5、y;else sbar=1;if(z0.0)aas=as1;else aas=as2;ffs=ffs+fss*aas;mis=mis+fss*aas*z;/*求预应力钢筋应变应力*ffp=0.0;mip=0.0;sp=s0+dd*cj+fpe/ep;if(spspy)fps=ep*sp;else if(sp1.0);mj=mm;/*各阶段破坏标志*if(mj=0) J_chushi=j; if(sp-fpe/ep)=0) J_xiaoya=j; if(ss-ssy)=0) J_squfu=j; if(sp-spu)=0) J_pqufu=j; if(pbar=1|concrete=1) J_p

6、ohuai=j;while(concrete=0&pbar=0);jmax=j;/*输出结果*file1=fopen(MC1.dat,w);for(j=0;jjmax;j+)fprintf(file1,%18.8g,%18.8g n,cj,mj);fclose(file1);/*求P-曲线*float a,x,ln,da,cc,dc,d0,d1,dd1;float p500,d500,mom100,coc 100;for(i=0;i500;i+)pi=0.0;di=0.0;for(i=0;i 100;i+)momi=0.0;coci=0.0;dc=0.00000002;cc=c0;j=0;a=

7、4000;x=2000;ln=50;da=a/ln;docc=cc+dc;j=j+1;i=-1;doi+;if(ci=cc&ccci+1)mm=mi+(cc-ci)*(mi+1-mi)/(ci+1-ci);while(ijmax1);pj=mm/a;d0=0.0;for(int n=1;n=ln;n+)momn=mm/a*n*da;jj=-1;dojj+;if(mjj=momn&momnmjj+1)cocn=cjj+(momn-mjj)*(cjj+1-cjj)/(mjj+1-mjj);while(jjjmax1);dd1=cocn*da*(n*da-da/2);d0=d0+dd1;d1=cc*b*(a+x);dj=d0+d1;while(cc=cjmax