!こっちを使う
!２種材料のせん断補正係数を調べる。
!片持ち梁のティモの理論値を出すとか
implicit real*8(a-h, o-z)
!入力するファイルの名前
open(8,file='tawami.dat')
!出力するファイルの名前
open (7,file='2syukekka.dat')
tyusin=50.0d-2
!b1=50.0d-2!角材の幅
b2=0.0!鋼板の幅
!h1=0.0!角材の高さ
!h2=0.500d0!鋼板の高さ
do i2=100,180,10
b2=b2+1.0d-2
h1=(2.0d0*tyusin*b2)/(tyusin-3.0d0*b2)
b1=tyusin-b2
h2=tyusin+h1
h3=h2-h1*2.0d0!真ん中の合板の長さ
ell=10.0d0!スパン
!木材
e=206.0d9!ヤング率
poa=0.3
g=e/(4.0d0*(1.0d0+poa))!せん断弾性係数
!鋼材
poa2=0.30d0
e2=206.0d9
g2=e2/((1.0d0+poa2)*2.0d0)
pkata=1.000d3!荷重(片持ち)
kakun=1.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!換算断面積
ea1=b1*h1*kakun*2.0d0*e
ea2=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*
hakokea=hakok*ea!bankのk
syotoukata=(pkata*ell**3)/(3.0d0*ei)!初等梁理論値(片持ち)
timokata=syotoukata+&
&(pkata*ell)/hakokea!ティモシェンコ理論値(片持ち)
timowarusyotoukata=timokata/syotoukata
read(8,*) atusa,ccxtawami
ccxk=(pkata*ell)/((ccxtawami-syotoukata)*ea)
sigumamoku=((e*pkata*ell/2.0d0)*(h2/2.0d0+yz))/ei
sotaigosa=((timokata-ccxtawami)/timokata)*100
!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*,"換算せん断補正係数bank"
!print*,hakok
!print*,"せん断補正係数bank"
!print*,kea
!print*,"初等梁理論値(片持ち)"
!print*,syotoukata
!print*,"ティモシェンコ理論値(片持ち)"
!print*,timokata
!print*,"ティモシェンコ/初等梁(片持ち)"
!print*,timowarusyotoukata
!print*,"相対誤差ティモシェンコ(%)"
!print*,sotaigosa

write(7,'(1pe13.5,a,1pe13.5,a,1pe13.5,a,1pe13.5)')atusa&
&,' ',sotaigosa,' ',ccxk,' ',hakok
end do
close(8)
end
