\‘’—ΝŠw@\‘’HŠw“Α˜_

‚a‚`‚r‚h‚bƒvƒƒOƒ‰ƒ€

•½”i‚S•Σ’PƒŽxŽA“™•ͺ•z‰Χdj

10 REM save"plate",a
20 REM plate 4-simple-supported
30 REM uniform load
40 REM 1998.6.26
50 DEF FNSINH#(X)=(EXP(X)-EXP(-X))/2
60 DEF FNCOSH#(X)=(EXP(X)+EXP(-X))/2
70 DEF FNTANH#(X)=(EXP(X)-EXP(-X))/(EXP(X)+EXP(-X))
80 K=10
90 READ E: REM young modulus
100 READ H: REM thickness of plate
110 READ NU: REM poisson' ratio
120 READ A
130 READ B
140 READ Q: REM uniform load
150 READ X
160 READ Y
170 '“ό—Νƒf[ƒ^
180 DATA 2.1e+6
190 DATA 0.2
200 DATA 0.3
210 DATA 16
220 DATA 16
230 DATA 0.05
240 DATA 8
250 DATA 0
260 D=E*H^3/12/(1-NU^2)
270 PAI=3.14159
280 W=0
290 FOR M=1 TO K STEP 2
300 MPY=M*PAI*Y/A
310 MPX=M*PAI*X/A
320 ALM=M*PAI*B/2/A
330 W=W+(1-(ALM*FNTANH#(ALM)+2)/2/FNCOSH#(ALM)*FNCOSH#(MPY)+1/2/FNCOSH#(ALM)*MPY*FNSINH#(MPY))*SIN(MPX)/M^5
340 NEXT M
350 W=W*4*Q*A^4/PAI^5/D
360 PRINT
370 PRINT "d=" ; D
380 PRINT "a=" ; A
390 PRINT "b=" ; B
400 PRINT "x=" ; X
410 PRINT "y=" ; Y
420 PRINT "q=" ; Q
430 PRINT "Ί³½³";K;"Ί³ΟΓή"
440 PRINT"w="; W
450 END