数値組込み関数の桁あふれ

 投稿者:idx  投稿日:2014年 1月16日(木)14時35分14秒
  REM 部分荷重下の単純支持板の中央のたわみ
REM 「板とシェルの理論」130頁143式 参照
REM  a:板のx方向の幅
REM  b:板のy方向の長さ
REM  u:荷重のx方向の幅
REM  v:荷重のy方向の長さ
REM  q:単位面積当たりの荷重(kN/cm^2)
REM  t:板厚
REM  mu:ポアソン比(0.3)
REM  E:縦弾性係数 20500(kN/cm^2)

INPUT PROMPT "a(cm) =":a
INPUT PROMPT "b(cm) =":b
INPUT PROMPT "u(cm) =":u
INPUT PROMPT "v(cm) =":v
INPUT PROMPT "t(cm) =":t
INPUT PROMPT "q(kN/cm^2) =":q

LET mu=0.3
LET E=20500
LET D=E*t^3/12/(1-mu^3)
LET GS=0
FOR M=1 TO 1000 '100回程度ならOK
   IF MOD(M ,2) = 0 THEN GOTO 10

   LET AM=M*PI*b/2/a
   LET GM=M*PI*v/4/a
   LET s=1/M^5*SIN(M*PI*u/2/a)*(1-1/COSH(AM)*(COSH(AM-2*GM)+GM*SINH(AM-2*GM)+AM*SINH(2*GM)/2/COSH(AM)))
   LET GS=GS+S
10 NEXT M
   LET W=GS*4*q*a^4/d/PI^5
   PRINT "Sigma=";GS
   PRINT "w(cm)=";W



END
 

戻る