非線性大作業(yè)(編輯修改稿)
2024-09-01 17:59 本頁面
【文章內(nèi)容簡介】 參考文獻[1][D].武漢：華中科技大學，2012[2][D].長沙：國防科學技術大學，2006附錄[]function y=conc(s) %混凝土應力應變關系（開裂點前用）fc=。s0=。su=。ft=。Ec=*10^4。 if s=0 if s=s0 y=fc*(1(1s/s0)^2)。 elseif s=su y=fc。 else y=0。 endelseif Ec*s=ft y=Ec*s。else y=0。endend[]function y=reinf(s) %普通鋼筋應力應變關系Es=2*10^5。 %MPafy=335。if abs(s)=fy/Es y=Es*s。else y=s/abs(s)*fy。end[]function y=pres(s) %預應力鋼筋應力應變關系Es=*10^5。 %MPafy=。if abs(s)=fy/Es y=Es*s。else y=s/abs(s)*fy。end[]function y=find_ini(x) %初始點求解器str_t=x(1)。phi=x(2)。global h b as h0 As1 As2 As3 h1 ps sigma Es Np str_t0 phi0n=1000。 %混凝土截面劃分為n層hs=h/n。 %混凝土分層厚度c=0。mcc=0。for k=1:n strainc(k)=str_t(k1/2)*hs*phi。 %混凝土第 k 層應變（壓為正）曲率以下緣受拉為正 stressc(k)=conc(strainc(k))。 c=c+stressc(k)*hs*b。 %混凝土合力計算（以受壓為正） mcc=mcc+stressc(k)*hs*b*(k1/2)*hs。 %對頂部取矩endstrain1=str_tphi*as。 %上部鋼筋應變（壓為正）stress1=reinf(strain1)。 %上部鋼筋應力strain2=str_tphi*h0。 %下部鋼筋應變（壓為正）stress2=reinf(strain2)。 %下部鋼筋應力c=c+stress1*As1+stress2*As2Np。 %平衡條件1：合力（以受壓為正）m1=stress1*As1*as。 %對頂部取矩m2=stress2*As2*h0。mp=Np*h1。m=mcc+m1+m2mp。 %平衡條件2：合力矩y(1)=c。y(2)=m。[]function y=find_m(s,phi)global h b as h0 As1 As2 As3 h1 ps sigma Es Np str_t0 phi0n=1000。hs=h/n。c=0。mc=0。fo
化學相關推薦
鄂ICP備17016276號-1