非線性大作業(yè)-閱讀頁
2024-08-24 17:59本頁面
【正文】 in1)。 %下部鋼筋應(yīng)變（壓為正）stress2=reinf(strain2)。 %預(yù)應(yīng)力鋼筋應(yīng)變（壓為正）stress3=pres(strain3)。 %對頂部取矩m2=stress2*As2*h0。m=+mc+m1+m2+m3。[]function y=find_phi(s,a)global h b as h0 As1 As2 As3 h1 ps sigma Es Np str_t0 phi0n=1000。c=0。for k=1:n strainc(k)=(k1/2)*hs*a+s。 c=c+stressc(k)*hs*b。 %對頂部取矩endstrain1=as*a+s。 %上部鋼筋應(yīng)力strain2=sh0*a。 %下部鋼筋應(yīng)力strain3=sh1*a(str_t0h1*phi0)sigma/Es。 %預(yù)應(yīng)力鋼筋應(yīng)力c=c+stress1*As1+stress2*As2+stress3*As3。[]clcclearglobal h b as h0 As1 As2 As3 h1 ps sigma Es Np str_t0 phi0ntd=1000。Algorithm39。levenbergmarquardt39。ScaleProblem39。Jacobian39。Display39。off39。TolFun39。TolX39。 %迭代選項(xiàng)h=400。as=40。ps=80。L=5400。As2=339。sigma=。Np=As3*sigma。str_t0=x0(1)。 %初始曲率和初始上緣應(yīng)變（以壓為正）str_pc0=str_t0phi0*h1。 %上緣應(yīng)變列表m(1)=0。 %初始曲率for i=2:length(str_t) phi(i)=fsolve((a)find_phi(str_t(i),a),[phi(i1)],options)。endfigure。grid。{\it\phi} (m^{1})39。ylabel(39。)。M_phi關(guān)系圖39。emf39。%求構(gòu)件P_f關(guān)系，由對稱性取半跨進(jìn)行結(jié)構(gòu)計(jì)算dL=10。xx=(L/2:dL:0)。 %荷載列表mm=zeros(length(str_t),length(xx))。aa=zeros(length(str_t),length(xx))。for i=1:length(str_t) %跨中混凝土上緣應(yīng)變作為控制變量 for j=1:length(xx) %求節(jié)點(diǎn)彎矩 if xx(j)L/3 mm(i,j)=3*m(i)/L*xx(j)。 end end aa(i,1)=phi(i)。 for j=2:length(xx) aa(i,j)=interp1(m,phi,mm(i,j))。 %共軛梁法求跨中撓度 endendfigure。grid。{\itf} (mm)39。ylabel(39。)。P_f關(guān)系圖39。emf39。P_exp=[0 ]。figure。hold all。39。Marker39。*39。hold off。xlabel(39。)。{\itP} (kN)39。legend(39。,39。,39。,39。)。P_f對比圖39。emf39。15
化學(xué)相關(guān)推薦
鄂ICP備17016276號-1