!ifを使って要素を設定するプログラムifyouso.f90の
!計算用プログラム
!長方形断面k=5/6
implicit real*8(a-h,o-z)
dimension x(300000),y(300000),z(300000)&
&,narax(300000),naray(300000),naraz(300000)&
&,mesh(200000,8),e66(6,6)
!材料の寸法(m)
b=1.0d0!幅
h=10.0d0!高さ
ell=3.5d1!軸長
p=(1.000d5)/2
!ヤング率の入力
ezz=7.694d9!軸方向のヤング率
!せん断弾性係数の入力
gxy=ezz/15.0d0!せん断弾性係数
!gxz=gxy
!gyz=gxya
!片持ち初等梁理論
plsanjo=p*(ell**3)
dniji=((b)*(h**3))/12
riontikatamoti=((p)*(ell**3))/(3*ezz*dniji)
!片持ちティモのせん断補正項
dmenseki=b*h!断面積
hosei=((p*ell)/(((gxy*dmenseki*5)/6)))
timokata=(((p)*(ell**3))/(3*ezz*dniji))+((p*ell)/(((gxy*dmenseki*5)/6)))
print*,"danmenniji"
print*,dniji
print*,"riontikatamoti"
print*,riontikatamoti
print*,"plsanjo"
print*,plsanjo
print*,"断面積"
print*,dmenseki
print*,"せん断弾性係数"
print*,gxy
print*,"せん断補正項"
print*,hosei
print*,"片持ち梁ティモシェンコ理論値"
print*,timokata
!print*,""
!print*,
end
