!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