!応急橋のせん断補正係数と3点曲げ試験のたわみの理論値、
!等曲げ試験のたわみの理論値を出すプログラム
implicit real*8(a-h, o-z)
b1=0.120d0!角材の幅
b2=0.009d0!鋼板の幅
h1=0.120d0!角材の高さ
h2=0.500d0!鋼板の高さ
h3=h2-h1*2.0d0!真ん中の合板の長さ
ell=7.00d0!スパン
ell2=ell/2.0d0!スパンの半分
ell3=1.30d0!等曲げ試験の荷重までの距離
!e=6.921d9!ヤング率
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!断面積
ea=b1*h1*kakun*2.0d0*e+b2*h2*koun*e2!換算断面積
!yz
yz=(e2*(b2*(h2-h1)*(h2/2-(h2-h1)/2)))/(e2*(b2*(h2-h1))+e*2*kakun/2*b1*h1)
!木材部分断面2次モーメント
dnijimoku=b1*kakun*(h2**3)/12.0d0-b1*kakun*(h3**3)/12.0d0
!鋼板部分断面2次モーメント
dnijikouhan=b2*koun*(h2**3)/12.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=8.65473226E-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*,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
