!矩形断面のたわみを読み込んで相対誤差を計算するプログラム
!Cowperによる長方形断面のk
implicit real*8(a-h,o-z)
!入力するファイルの名前
open(8,file='hakotawami.dat')
!出力するファイルの名前
open (7,file='soutaigosa1.dat')
!数値(断面に応じて変えること！)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ell=1.0d-2
poa=0.3
p=1.0d3
young=206.d9
!全体の幅
b=2.0d-2
!全体の高さ
h=2.0d-2
!左右の板厚
t=2.0d-3
!上下の板厚
t1=2.0d-3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!中
b1=b-2*t
h1=h-2*t1
sendandansei=young/(2*(1+poa))
danmenniji=((b*h**3)-(b1*h1**3))/12.0d0
danmenseki=(b*h)-(b1*h1)
emu=(b1*t1)/(h*t)
enu=b1/h
hosei=((10*(1+poa)*(1+3*emu)**2)/&
&((12+72*emu+150*(emu**2)+90*(emu**3))+&
&poa*(11+66*emu+135*(emu**2)+90*(emu**3))+&
&10*(enu**2)*((3+poa)*emu+3*(emu**2))))
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do i2=100,280,10
ell=ell+1.0d-2
syotou=((p*ell**3.0d0)/(3.0d0*young*danmenniji))
timo=((p*ell**3.0d0)/(3.0d0*young*danmenniji))+&
&((p*ell)/(sendandansei*hosei*danmenseki))
read(8,*) supan,suti
syosou=((suti-syotou)/syotou)*100
timosou=((suti-timo)/timo)*100
erueiti=ell/2.0d-2
write(7,'(1pe13.5,a,1pe13.5,a,1pe13.5)')erueiti,' ',syosou,' ',timosou
end do
close(8)
end
