home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e004 / 1.img / WINDFALL.MCD < prev   
Encoding:
Text File  |  1987-10-23  |  3.4 KB  |  153 lines
  1. .MCD 20000 0
  2. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  3. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  4. .CMD SET ORIGIN 0
  5. .CMD SET TOL 0.001000
  6. .CMD MARGIN 0
  7. .CMD LINELENGTH 78
  8. .CMD SET PRNCOLWIDTH 8
  9. .CMD SET PRNPRECISION 4
  10. .CMD FILENAME data1 data1
  11. .CMD FILENAME data2 data2
  12. .TXT 1 0 1 45 
  13. a1,45,43,44
  14. MOTION OF FALLING BODY WITH WIND RESISTANCE
  15. .TXT 2 0 3 70 
  16. a3,70,68,185
  17. A parachutist jumps out of an airplane and opens his chute 5 seconds 
  18. into his fall. Examine the characteristics of the flight, including 
  19. position, velocity, acceleration, and "jerk."
  20. .TXT 4 0 2 45 
  21. a2,45,43,64
  22. To start, set up N short time intervals dt. 
  23. Also define units:
  24. .EQN 3 1 1 9 
  25. N:20
  26. .EQN 0 16 1 12 
  27. i:0;N
  28. .EQN 0 23 1 11 
  29. sec:1T
  30. .EQN 0 11 1 10 
  31. lb:1M
  32. .EQN 1 -50 1 14 
  33. dt:.5*sec
  34. .EQN 0 16 2 12 
  35. t[i:i*dt
  36. .EQN 0 23 1 10 
  37. ft:1L
  38. .EQN 0 11 2 19 
  39. g:-32*ft*sec^-2
  40. .EQN 2 -11 1 13 
  41. lbf:lb*g
  42. .TXT 1 -40 1 32 
  43. a1,32,30,31
  44. Initial conditions as follows:
  45. .EQN 2 1 2 15 
  46. y[0:2000*ft
  47. .TXT 0 18 1 46 
  48. a1,46,44,45
  49. Parachutist jumps out of plane at 2000 feet.
  50. .EQN 2 -18 3 13 
  51. v[0:0*ft/sec
  52. .TXT 1 18 1 40 
  53. a1,40,38,39
  54. Parachutist starts with zero velocity.
  55. .EQN 3 -18 1 13 
  56. m:200*lb
  57. .TXT 0 18 1 33 
  58. a1,33,31,32
  59. Mass of parachutist plus chute.
  60. .TXT 2 -19 3 72 
  61. a3,72,70,141
  62. Force on parachutist results from gravity and air resistance. 
  63. Gravitational force is constant; air resistance is proportional to 
  64. velocity.
  65. .EQN 4 1 1 11 
  66. Fg:m*g
  67. .EQN 0 13 1 14 
  68. Fg=?lbf
  69. .EQN 2 -13 5 23 
  70. α_nochute:.5*lbf/((ft/sec))
  71. .TXT 0 29 2 35 
  72. a2,35,33,54
  73. The two α constants are a measure 
  74. of air resistance.
  75. .EQN 6 -29 5 20 
  76. α_chute:5*lbf/((ft/sec))
  77. .TXT 2 29 3 47 
  78. a3,47,45,132
  79. The Ff function represents the frictional 
  80. force for air resistance over time, which 
  81. changes after 5 seconds when the chute opens.
  82. .EQN 4 -29 1 64 
  83. Ff(t):α_nochute*(Φ(5*sec-t))+α_chute*(1-Φ(5*sec-t))
  84. .TXT 2 -1 1 69 
  85. a1,69,67,68
  86. Now compute velocity, acceleration, position, and "jerk" over time.
  87. .EQN 2 1 4 32 
  88. v[(i+1):v[i+(Fg+Ff(t[i)*v[i)/m*dt
  89. .TXT 2 39 2 34 
  90. a2,34,32,66
  91. Change in velocity computed from 
  92. acceleration (force over mass).
  93. .EQN 4 -39 2 20 
  94. y[(i+1):y[i+v[i*dt
  95. .TXT 0 39 1 34 
  96. a1,34,32,33
  97. Position computed from velocity.
  98. .EQN 3 -39 4 19 
  99. a[(i+1):(v[(i+1)-v[i)/dt
  100. .TXT 2 39 1 38 
  101. a1,38,36,37
  102. Acceleration computed from velocity.
  103. .EQN 3 -39 2 9 
  104. a[0:g
  105. .EQN 2 -1 4 17 
  106. j[i:(a[(i+1)-a[i)/dt
  107. .TXT 1 40 3 37 
  108. a3,37,35,81
  109. "Jerk" computed from acceleration. 
  110. (Jerk is the third derivative of 
  111. position.)
  112. .EQN 4 -40 3 12 
  113. j[0:0*g/dt
  114. .TXT 4 5 1 20 
  115. a1,20,18,19
  116. POSITION OVER TIME
  117. .TXT 0 40 1 20 
  118. a1,20,18,19
  119. VELOCITY OVER TIME
  120. .EQN 1 -45 17 30 
  121. y[0&y[N&y[i{1,1,15,25,l}@t[N&t[0&t[i
  122. .EQN 0 38 17 33 
  123. 0*ft/sec&v[11&v[i{1,1,15,25,l}@t[N&t[0&t[i
  124. .TXT 19 -31 1 24 
  125. a1,24,22,23
  126. ACCELERATION OVER TIME
  127. .TXT 0 44 1 18 
  128. a1,18,16,17
  129. "JERK" OVER TIME
  130. .EQN 1 -51 17 34 
  131. a[12&a[0&a[i,0*g{1,1,15,25,l}@t[N&t[0&t[i
  132. .EQN 0 39 17 36 
  133. j[11*1.1&j[12*1.1&j[i,0*g/sec{1,1,15,25,l}@t[N&t[0&t[i
  134. .TXT 10 -6 1 3 
  135. a1,3,1,2
  136. 0
  137. .TXT 0 41 1 3 
  138. a1,3,1,2
  139. 0
  140. .TXT 8 -74 1 77 
  141. a1,77,75,76
  142. ---------------------------------------------------------------------------
  143. .EQN 1 44 28 11 
  144. a[i/(ft/sec^2)=
  145. .EQN 0 16 28 11 
  146. j[i/(ft/sec^3)=
  147. .EQN 1 -33 27 12 
  148. v[i/(ft/sec)=
  149. .EQN 2 -26 25 7 
  150. t[i/sec=
  151. .EQN 0 10 25 12 
  152. y[i/ft{19075}=
  153.