!応急橋のせん断補正係数と3点曲げ試験のたわみの理論値、 !等曲げ試験のたわみの理論値を出すプログラム !b1=120にするとyzが変わってしまうのでb1=420にしてkakun=7から !kakun=2に変更 implicit real*8(a-h, o-z) b1=0.120d0!角材の幅 b2=0.009d0!鋼板の幅 b3=0.007d0!継手部分の幅 b4=0.068d0! h1=0.120d0!角材の高さ h2=0.500d0!鋼板の高さ h3=h2-h1*2.0d0!真ん中の合板の長さ h4=0.007d0! h5=0.125d0!元のh3 ell=7.00d0!スパン ell2=ell/2.0d0!スパンの半分 ell3=1.30d0!等曲げ試験の荷重までの距離 e=7.694d9!ヤング率 g=e/15d0!せん断弾性係数 poa2=0.30d0!ポアソン比 e2=206.0d9 g2=e2/((1.0d0+poa2)*2.0d0) p3ten=1.000d5!荷重(3点曲げ) p3ten2=p3ten/2.0d0!荷重(3点曲げ)の半分 ptoumage=1.4000d5!荷重(等曲げ) ptoumage2=ptoumage/2.0d0!荷重(等曲げ)の半分 kakun=7.0d0!角材の数(計算で使うから一応実数) koun=2.0!d0合板の数(計算で使うから一応実数) !dmenseki=b1*h1*kakun*2.0d0+b2*h2*koun!断面積 dmenseki=b1*h1*kakun*2.0d0+b2*h2*koun+b3*h5*koun+b4*h4*koun!継手がある場合の断面積 !ea=b1*h1*kakun*2.0d0*e+b2*h2*koun*e2!換算断面積 ea=b1*h1*kakun*2.0d0*e+b2*h2*koun*e2+b3*h5*koun*e2+b4*h4*koun*e2!継手がある場合の換算断面積 !足した部分-yz yz=(e2*(b3*h5*(-h2/2.0+h5/2.0)& &+b4*h4*(-h2/2.0+h4/2.0)))/(e2*(b2*h2+b3*h5+b4*h4)+kakun*e*b1*h1) !木材部分断面2次モーメント !dnijimoku=b1*kakun*(h2**3)/12.0d0-b1*kakun*(h3**3)/12.0d0!元のdnijimoku dnijimoku=b1*kakun*(((h2/2+yz)**3)/3.0d0-((h2/2+yz-h1)**3)/3.0d0)& &+(b1*kakun*((h2/2-yz)**3)/3.0d0)-(b1*kakun*((h2/2-yz-h1)**3)/3.0d0)!追加 !鋼板部分断面2次モーメント !dnijikouhan=b2*koun*(h2**3)/12.0d0!元のdnijiokuhan dnijikouhan=b2*koun*((h2/2+yz)**3)/3.0d0+b2*koun*((h2/2-yz)**3)/3.0d0& &+b3*koun*((h2/2+yz)**3)/3.0d0-b3*koun*((h2/2+yz-h5)**3)/3.0d0& &+b4*koun*((h2/2+yz)**3)/3.0d0-b4*koun*((h2/2+yz-h4)**3)/3.0d0 !dnijikouhan=(b2*koun*((h2/2-yz)**3))/3.0d0& !&+b2*koun*((h2/2-h1-yz)**3)/3.0d0& !&+b3*koun*((h2/2-yz)**3)/3.0d0-b3*koun*((h2/2-h5-yz)**3)/3.0d0& !&+b4*koun*((h2/2-yz)**3)/3.0d0-b4*koun*((h2/2-yz-h4)**3)/3.0d0!追加 !曲げ剛性 ei=dnijimoku*e+dnijikouhan*e2 hakob=b1*kakun+b2!bankのb hakoh=h2-h1!bankのh hakom=hakob*h1/(hakoh*b2)!cowperのm hakom2=hakom**2!bankのm**2 hakom3=hakom**3!bankのm**3 hakon=hakob/hakoh!bankのn hakon2=hakon**2!bankのn**2 a=e2/e!bankのα a2=a**2!bankaのα** a3=a**3!bankaのα*3 hakok=(20.0d0*(a+3.0d0*hakom)**2)& &/((e/g)*(60.0d0*hakom2*hakon2+60.0d0*a*hakom*hakon2)& &+(e/g2)*(180.0d0*hakom3+300.0d0*a*hakom2+144.0d0*a2*hakom+24.0d0*a3)& &+poa2*(30.0d0*hakom2+6.0d0*a*hakom-4.0d0*a2))!bankのk* kea=hakok*ea!bankのk syotou3ten=(p3ten2*ell2**3)/(3.0d0*ei)!初等梁理論値(3点曲げ) timo3ten=syotou3ten+& &(p3ten2*ell2)/kea!ティモシェンコ理論値(3点曲げ) syotoutoumage=(ptoumage2*ell3*(3.0d0*ell**2-4.0d0*ell3**2))& &/(24.0d0*ei)!初等梁理論値(等曲げ) timotoumage=syotoutoumage& &+(ptoumage2*ell3)/kea!ティモシェンコ理論値(等曲げ) timowarusyotou3ten=timo3ten/syotou3ten timowarusyotoutoumage=timotoumage/syotoutoumage ccxtawami=6.82379538E-03 ccxk=(p3ten*ell)/((ccxtawami-syotou3ten)*ea) sigumamoku=((e*p3ten*ell/2.0d0)*(h2/2.0d0+yz))/ei print*,"ccxk" print*,ccxk print*,"sigumamoku" print*,sigumamoku print*,"yz" print*,yz print*,"hakom" print*,hakom print*,"hakon" print*,hakon print*,"断面積" print*,dmenseki print*,"曲げ剛性木" print*,dnijimoku*e print*,"曲げ剛性鋼" print*,dnijikouhan*e2 print*,"曲げ剛性" print*,ei print*,"換算せん断補正係数" print*,hakok print*,"せん断補正係数" print*,kea print*,"初等梁理論値(3点曲げ)" print*,syotou3ten print*,"ティモシェンコ理論値(3点曲げ)" print*,timo3ten print*,"ティモシェンコ/初等梁(3点曲げ)" print*,timowarusyotou3ten print*,"初等梁理論値(等曲げ)" print*,syotoutoumage print*,"ティモシェンコ理論値(等曲げ)" print*,timotoumage print*,"ティモシェンコ/初等梁(等曲げ)" print*,timowarusyotoutoumage end