附錄一:受試者同意書
研究題目:國小學童不同軟硬表面著地下肢勁度調節之機轉及其影響 指導教授:黃長福博士
研究生:許太彥
單位:國立台灣師範大學體育研究所 聯絡電話:02-29312901 轉 13
依實驗研究的規定與保護受試者權益,研究者應將研究過程可能發生 的危險和法律上的責任等事宜,向受試者說明清楚。而研究者應盡其所能 保護受試者的健康與權益,並隨時回答受試者的問題。受試者如有改變意 願時應通知實驗者,可隨時退出實驗而不受任何限制。參與實驗研究的受 試者必須明瞭並注意下列事項:
一、 實驗時間自民國九十二年二月十日至二月十五日。
二、在指定時間至台灣師範大學體育研究所力學實驗室,參與三種不同緩 衝時間、三種不同材質著地動作的實驗,每種材質測試 6-9 次。
三、所有的測驗過程,使用 KISTLER 9287 型的測力板、PEAK 高速攝影 機及 Biovision 肌電儀器來收集動力學、運動學及肌電的資料。
四、參與本實驗的受試者可深入了解到自己在不同著地條件下著地動作的 特性。
本實驗須要您的參與合作,才能圓滿順利。現在請您填寫一份您的個人資 料表。請在下表姓名欄內簽名,並請監護人簽名,表示同意並遵守同意書 內的各項規定。
自願者簽名: _____________ 監護人簽名:______________
日 期:92 年 月 日
附錄二:變異數分析摘要表
表 4-1-1 緩衝期的時間
變異數分析摘要表變異來源 SS df MS F
受試者間 146547.38 9.00
受試者內 124241.78 80.00
材質(A) 22671.29 2.00 11335.64 10.89 *
殘差 (A°S) 18733.16 18.00 1040.73
緩衝時間(B) 59469.16 2.00 29734.58 55.36 *
殘差 (B°S) 9668.62 18.00 537.15
材質°緩衝時間(A°B) 2512.71 4.00 628.18 2.02
殘差 (A°B°S) 11186.84 36.00 310.75
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-1-2 緩衝期下肢勁度
變異數分析摘要表變異來源 SS df MS F
受試者間 1717270497.11 9.00
受試者內 737668787.78 80.00
材質(A) 134779262.82 2.00 67389631.41 5.15 *
殘差 (A°S) 235464750.96 18.00 13081375.05
緩衝時間(B) 209996690.02 2.00 104998345.01 47.75 *
殘差 (B°S) 39578811.76 18.00 2198822.88
材質°緩衝時間(A°B) 7912140.04 4.00 1978035.01 0.65
殘差 (A°B°S) 109937132.18 36.00 3053809.23
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-1-3 緩衝期髖關節角勁度
變異數分析摘要表變異來源 SS df MS F
受試者間 367.44 9.00
受試者內 402.41 80.00
材質(A) 39.52 2.00 19.76 6.17 *
殘差 (A°S) 57.66 18.00 3.20
緩衝時間(B) 93.31 2.00 46.66 8.37 *
殘差 (B°S) 100.32 18.00 5.57
材質°緩衝時間(A°B) 12.71 4.00 3.18 1.16
殘差 (A°B°S) 98.89 36.00 2.75
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-1-4 緩衝期膝關節角勁度
變異數分析摘要表變異來源 SS df MS F
受試者間 21.84 9.00
受試者內 12.53 80.00
材質(A) 1.83 2.00 0.92 8.74 *
殘差 (A°S) 1.89 18.00 0.10
緩衝時間(B) 3.89 2.00 1.94 39.37 *
殘差 (B°S) 0.89 18.00 0.05
材質°緩衝時間(A°B) 0.26 4.00 0.06 0.62
殘差 (A°B°S) 3.77 36.00 0.10
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-1-5 緩衝期踝關節角勁度
變異數分析摘要表變異來源 SS df MS F
受試者間 33.70 9.00
受試者內 22.96 80.00
材質(A) 5.89 2.00 2.94 7.72 *
殘差 (A°S) 6.87 18.00 0.38
緩衝時間(B) 0.24 2.00 0.12 1.02
殘差 (B°S) 2.16 18.00 0.12
材質°緩衝時間(A°B) 0.44 4.00 0.11 0.53
殘差 (A°B°S) 7.37 36.00 0.20
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-1 之事後比較達顯著差異的組別
緩衝時間*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;b3a2-b3a3;
a1b1-a1b2;a1b1-a1b3;a1b2-a1b3;a2b1-a2b2;a2b1-a2b3;
a2b2-a2b3;a3b1-a3b2;a3b1-a3b3;a3b2-a3b3
下肢勁度*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;
a1b1-a1b2;a1b1-a1b3; a1b2-a1b3;a2b1-a2b3;a3b1-a3b2;a3b1-a3b3 a3b2-a3b3
髖*+ b1a1-b1a2;b1a2-b1a3;
a1b1-a1b3;a3b1-a3b3;a3b2-a3b3 膝*+ b1a1-b1a3;b2a1-b2a3;
a1b1-a1b3;a2b1-a2b3;a3b1-a3b3;a3b2-a3b3
踝+ b1a1-b1a2;b1a1-b1a3;b2a1-b2a3;b3a1-b3a2;b3a1-b3a3
表 4-2-1 緩衝期下肢長度變化量
變異數分析摘要表變異來源 SS df MS F
受試者間 0.1543 9.00
受試者內 0.0764 80.00
材質(A) 0.0169 2.00 0.0085 12.97 *
殘差 (A°S) 0.0118 18.00 0.0007
緩衝時間(B) 0.0274 2.00 0.0137 37.75 *
殘差 (B°S) 0.0065 18.00 0.0004
材質°緩衝時間(A°B) 0.0015 4.00 0.0004 1.12
殘差 (A°B°S) 0.0123 36.00 0.0003
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-2-2 緩衝期髖關節角度變化量
變異數分析摘要表變異來源 SS df MS F
受試者間 12683.63 9.00
受試者內 8072.41 80.00
材質(A) 1550.79 2.00 775.39 7.65 *
殘差 (A°S) 1824.39 18.00 101.36
緩衝時間(B) 1808.28 2.00 904.14 13.30 *
殘差 (B°S) 1223.95 18.00 68.00
材質°緩衝時間(A°B) 167.45 4.00 41.86 1.01
殘差 (A°B°S) 1497.55 36.00 41.60
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-2-3 緩衝期膝關節角度變化量
變異數分析摘要表變異來源 SS df MS F
受試者間 6326.85 9.00
受試者內 4051.97 80.00
材質(A) 655.17 2.00 327.59 13.12 *
殘差 (A°S) 449.44 18.00 24.97
緩衝時間(B) 1895.19 2.00 947.60 43.46 *
殘差 (B°S) 392.47 18.00 21.80
材質°緩衝時間(A°B) 47.69 4.00 11.92 0.70
殘差 (A°B°S) 612.01 36.00 17.00
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-2-4 緩衝期踝關節角度變化量
變異數分析摘要表變異來源 SS df MS F
受試者間 1630.63 9.00
受試者內 1961.77 80.00
材質(A) 119.71 2.00 59.85 1.45
殘差 (A°S) 740.99 18.00 41.17
緩衝時間(B) 58.75 2.00 29.37 1.46
殘差 (B°S) 362.68 18.00 20.15
材質°緩衝時間(A°B) 16.15 4.00 4.04 0.22
殘差 (A°B°S) 663.49 36.00 18.43
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-2 之事後比較達顯著差異的組別
下肢長度*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;;b3a1-b3a3
a1b1-a1b2;a1b1-a1b3;a1b2-a1b3;a2b1-a2b3;a3b1-a3b3;a3b1-a3b2;
a3b2-a3b3
髖*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;
a1b1-a1b3;a1b2-a1b3;a3b1-a3b3;a3b2-a3b3
膝*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;b3a2-b3a3;
a1b1-a1b2;a1b1-a1b3; a1b2-a1b3;a2b1-a2b3;a2b2-a2b3;a3b1-a3b2;
a3b1-a3b3; a3b2-a3b3
踝
表 4-4-1 著地瞬間髖關節角度
變異數分析摘要表變異來源 SS df MS F
受試者間 3329.03 9.00
受試者內 1673.67 80.00
材質(A) 22.44 2.00 11.22 0.55
殘差 (A°S) 368.16 18.00 20.45
緩衝時間(B) 259.04 2.00 129.52 5.61 *
殘差 (B°S) 415.32 18.00 23.07
材質°緩衝時間(A°B) 43.55 4.00 10.89 0.69
殘差 (A°B°S) 565.16 36.00 15.70
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-4-2 著地瞬間膝關節角度
變異數分析摘要表變異來源 SS df MS F
受試者間 3644.61 9.00
受試者內 1012.89 80.00
材質(A) 83.04 2.00 41.52 1.94
殘差 (A°S) 384.45 18.00 21.36
緩衝時間(B) 62.80 2.00 31.40 5.36 *
殘差 (B°S) 105.37 18.00 5.85
材質°緩衝時間(A°B) 7.87 4.00 1.97 0.19
殘差 (A°B°S) 369.36 36.00 10.26
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-4-3 著地瞬間踝關節角度
變異數分析摘要表變異來源 SS df MS F
受試者間 1907.09 9.00
受試者內 713.70 80.00
材質(A) 109.97 2.00 54.98 4.36 *
殘差 (A°S) 227.10 18.00 12.62
緩衝時間(B) 6.54 2.00 3.27 0.57
殘差 (B°S) 103.92 18.00 5.77
材質°緩衝時間(A°B) 16.85 4.00 4.21 0.61
殘差 (A°B°S) 249.31 36.00 6.93
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-4 之事後比較達顯著差異的組別
髖* a1b1-a1b3;a1b2-a1b3;a3b2-a3b3 膝* a1b1-a1b2
踝+ b2a1-b2a2; b2a1-b2a3;b3a1-b3a3
表 4-5-1 身體重心最低時髖關節的角度
變異數分析摘要表變異來源 SS df MS F
受試者間 26268.48 9.00
受試者內 13321.69 80.00
材質(A) 1619.97 2.00 809.99 4.58 *
殘差 (A°S) 3181.71 18.00 176.76
緩衝時間(B) 3124.75 2.00 1562.38 11.43 *
殘差 (B°S) 2459.82 18.00 136.66
材質°緩衝時間(A°B) 355.77 4.00 88.94 1.24
殘差 (A°B°S) 2579.67 36.00 71.66
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-5-2 身體重心最低時膝關節的角度
變異數分析摘要表變異來源 SS df MS F
受試者間 15389.97 9.00
受試者內 4849.62 80.00
材質(A) 718.67 2.00 359.34 8.80 *
殘差 (A°S) 734.78 18.00 40.82
緩衝時間(B) 2068.99 2.00 1034.50 49.42 *
殘差 (B°S) 376.77 18.00 20.93
材質°緩衝時間(A°B) 74.15 4.00 18.54 0.76
殘差 (A°B°S) 876.25 36.00 24.34
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-5-3 身體重心最低時踝關節的角度
變異數分析摘要表變異來源 SS df MS F
受試者間 4101.48 9.00
受試者內 2206.04 80.00
材質(A) 407.82 2.00 203.91 5.74 *
殘差 (A°S) 639.69 18.00 35.54
緩衝時間(B) 58.25 2.00 29.13 1.52
殘差 (B°S) 343.93 18.00 19.11
材質°緩衝時間(A°B) 9.20 4.00 2.30 0.11
殘差 (A°B°S) 747.15 36.00 20.75
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-5 之事後比較達顯著差異的組別
髖*+ b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;
a1b1-a1b3;a1b2-a1b3;a3b1-a3b3;a3b2-a3b3
膝*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b3a1-b3a3;a1b1-a1b2;a1b1-a1b3;
a1b2-a1b3; a2b1-a2b3;a2b2-a2b3;a3b1-a3b3;a3b2-a3b3 踝+ b1a1-b1a3;b2a1-b2a3;b3a1-b3a3
表 4-6-1 緩衝期中髖關節的角位移
變異數分析摘要表變異來源 SS df MS F
受試者間 12879.70 9.00
受試者內 8208.25 80.00
材質(A) 1504.91 2.00 752.46 7.13 *
殘差 (A°S) 1898.54 18.00 105.47
緩衝時間(B) 1781.06 2.00 890.53 12.56 *
殘差 (B°S) 1276.64 18.00 70.92
材質°緩衝時間(A°B) 158.93 4.00 39.73 0.90
殘差 (A°B°S) 1588.16 36.00 44.12
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-6-2 緩衝期中膝關節的角位移
變異數分析摘要表變異來源 SS df MS F
受試者間 5997.61 9.00
受試者內 4127.50 80.00
材質(A) 550.51 2.00 550.51 9.79 *
殘差 (A°S) 505.94 18.00 505.94
緩衝時間(B) 1868.29 2.00 1868.29 32.82 *
殘差 (B°S) 512.31 18.00 512.31
材質°緩衝時間(A°B) 61.32 4.00 61.32 0.88
殘差 (A°B°S) 629.13 36.00 629.13
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-6-3 緩衝期中踝關節的角位移
變異數分析摘要表變異來源 SS df MS F
受試者間 2063.74 9.00
受試者內 2051.51 80.00
材質(A) 97.49 2.00 48.74 1.28
殘差 (A°S) 684.97 18.00 38.05
緩衝時間(B) 71.61 2.00 35.80 1.70
殘差 (B°S) 379.15 18.00 21.06
材質°緩衝時間(A°B) 2.39 4.00 0.60 0.03
殘差 (A°B°S) 815.91 36.00 22.66
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-6 之事後比較達顯著差異的組別
髖*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;
a3b1-a3b3;a1b1-a1b3
膝*+ b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;a1b1-a1b2;a1b1-a1b3;
a1b2-a1b3;a2b1-a2b3;a2b2-a2b3;a3b1-a3b2;a3b1-a3b3;a3b2-a3b3;
踝
表 4-7-1 緩衝期中髖關節的最大角速度
變異數分析摘要表變異來源 SS df MS F
受試者間 601574.41 9.00
受試者內 280321.06 80.00
材質(A) 18032.51 2.00 9016.25 1.36
殘差 (A°S) 118985.40 18.00 6610.30
緩衝時間(B) 21280.31 2.00 10640.15 5.75 *
殘差 (B°S) 33292.39 18.00 1849.58
材質°緩衝時間(A°B) 16856.63 4.00 4214.16 2.11
殘差 (A°B°S) 71873.84 36.00 1996.50
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-7-2 緩衝期中膝關節的最大角速度
變異數分析摘要表變異來源 SS df MS F
受試者間 288628.22 9.00
受試者內 151783.24 80.00
材質(A) 15347.73 2.00 7673.86 2.68
殘差 (A°S) 51458.47 18.00 2858.80
緩衝時間(B) 18992.07 2.00 9496.03 7.58 *
殘差 (B°S) 22539.81 18.00 1252.21
材質°緩衝時間(A°B) 7420.52 4.00 1855.13 1.85
殘差 (A°B°S) 36024.66 36.00 1000.68
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-7-3 緩衝期中踝關節的最大角速度
變異數分析摘要表變異來源 SS df MS F
受試者間 137344.05 9.00
受試者內 834975.06 80.00
材質(A) 22495.45 2.00 11247.73 0.05
殘差 (A°S) 407626.48 18.00 22645.92
緩衝時間(B) 50990.85 2.00 25495.42 4.53 *
殘差 (B°S) 101314.61 18.00 5628.59
材質°緩衝時間(A°B) 13466.15 4.00 3366.54 0.51
殘差 (A°B°S) 239081.52 36.00 6641.15
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-7 之事後比較達顯著差異的組別 髖* a1b1-a1b3;a2b1-a2b2 膝* a1b1-a1b2;a1b1-a1b3 踝* a1b1-a1b3;a1b2-a1b3
表 4-8-1 緩衝期中髖關節最大角速度產生的時間
變異數分析摘要表變異來源 SS df MS F
受試者間 4424.71 9.00
受試者內 9720.89 80.00
材質(A) 212.27 2.00 106.13 0.65
殘差 (A°S) 2959.29 18.00 164.40
緩衝時間(B) 653.87 2.00 326.93 2.92
殘差 (B°S) 2016.36 18.00 112.02
材質°緩衝時間(A°B) 378.67 4.00 94.67 0.97
殘差 (A°B°S) 3500.44 36.00 97.23
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-8-2 緩衝期中膝關節最大角速度產生的時間
變異數分析摘要表變異來源 SS df MS F
受試者間 4163.73 9.00
受試者內 12739.56 80.00
材質(A) 3317.69 2.00 1658.84 8.79 *
殘差 (A°S) 3395.20 18.00 188.62
緩衝時間(B) 525.16 2.00 262.58 3.55 *
殘差 (B°S) 1334.40 18.00 74.13
材質°緩衝時間(A°B) 658.84 4.00 164.71 1.69
殘差 (A°B°S) 3508.27 36.00 97.45
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-8-3 緩衝期中踝關節最大角速度產生的時間
變異數分析摘要表變異來源 SS df MS F
受試者間 2153.07 9.00
受試者內 1912.89 80.00
材質(A) 57.96 2.00 28.98 0.88
殘差 (A°S) 596.27 18.00 33.13
緩衝時間(B) 15.29 2.00 7.64 0.37
殘差 (B°S) 372.27 18.00 20.68
材質°緩衝時間(A°B) 223.64 4.00 55.91 3.11 *
殘差 (A°B°S) 647.47 36.00 17.99
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-8 之事後比較達顯著差異的組別 髖
膝*+ b2a1-b2a2;b2a1-b2a3;b3a1-b3a2;
a2b1-a2b2;a2b1-a2b3 踝 交互作用
表 4-9-1 地面垂直反作用力的第一峰值
變異數分析摘要表變異來源 SS df MS F
受試者間 17.71 9.00
受試者內 25.90 80.00
材質(A) 8.83 2.00 4.42 14.55 *
殘差 (A°S) 5.47 18.00 0.30
緩衝時間(B) 1.92 2.00 0.96 4.21 *
殘差 (B°S) 4.10 18.00 0.23
材質°緩衝時間(A°B) 0.88 4.00 0.22 1.68
殘差 (A°B°S) 4.70 36.00 0.13
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-9-2 地面垂直反作用力第一峰值的產生時間
變異數分析摘要表變異來源 SS df MS F
受試者間 283.79 9.00
受試者內 1715.33 80.00
材質(A) 559.62 2.00 279.81 23.06 *
殘差 (A°S) 218.38 18.00 12.13
緩衝時間(B) 15.29 2.00 7.64 0.70
殘差 (B°S) 196.04 18.00 10.89
材質°緩衝時間(A°B) 84.31 4.00 21.08 1.18
殘差 (A°B°S) 641.69 36.00 17.82
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-9-3 著地到第一峰值的衝量
變異數分析摘要表變異來源 SS df MS F
受試者間 294.93 9.00
受試者內 610.78 80.00
材質(A) 111.18 2.00 55.59 11.20 *
殘差 (A°S) 89.30 18.00 4.96
緩衝時間(B) 10.87 2.00 5.44 1.32
殘差 (B°S) 73.88 18.00 4.10
材質°緩衝時間(A°B) 69.01 4.00 17.25 2.42
殘差 (A°B°S) 256.53 36.00 7.13
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-9 之事後比較達顯著差異的組別
第一峰值*+ b1a1-b1a2;b1a1-b1a3;b2a1-b2a2;b2a1-b2a3;b3a1-b3a2;b3a1-b3a3;
a1b1-a1b3;a3b1-a3b3
產生時間+ b1a1-b1a2;b3a1-b3a2;b3a1-b3a3
產生衝量+ b3a1-b3a2;b3a1-b3a3
表 4-10-1 地面垂直反作用力的第二峰值
變異數分析摘要表變異來源 SS df MS F
受試者間 64.99 9.00
受試者內 55.04 80.00
材質(A) 3.37 2.00 1.69 3.98 *
殘差 (A°S) 7.64 18.00 0.42
緩衝時間(B) 9.51 2.00 4.76 9.83 *
殘差 (B°S) 8.71 18.00 0.48
材質°緩衝時間(A°B) 2.38 4.00 0.60 0.92
殘差 (A°B°S) 23.43 36.00 0.65
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-10-2 地面垂直反作用力第二峰值的產生時間
變異數分析摘要表變異來源 SS df MS F
受試者間 2707.07 9.00
受試者內 6578.22 80.00
材質(A) 8.42 2.00 4.21 0.03
殘差 (A°S) 2757.13 18.00 153.17
緩衝時間(B) 188.02 2.00 94.01 2.14
殘差 (B°S) 791.53 18.00 43.97
材質°緩衝時間(A°B) 264.44 4.00 66.11 0.93
殘差 (A°B°S) 2568.67 36.00 71.35
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-10-3 著地到第二峰值的衝量
變異數分析摘要表變異來源 SS df MS F
受試者間 18313.34 9.00
受試者內 13515.90 80.00
材質(A) 105.03 2.00 52.51 0.28
殘差 (A°S) 3317.10 18.00 184.28
緩衝時間(B) 194.98 2.00 97.49 0.71
殘差 (B°S) 2463.10 18.00 136.84
材質°緩衝時間(A°B) 831.90 4.00 207.98 1.13
殘差 (A°B°S) 6603.78 36.00 183.44
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-10 之事後比較達顯著差異的組別
第二峰值*+ b3a1-b3a3;a3b1-a3b3;a3b2-a3b3
表 4-11-1 50 毫秒內的被動衝量
變異數分析摘要表變異來源 SS df MS F
受試者間 7293.59 9.00
受試者內 4352.24 80.00
材質(A) 187.66 2.00 93.83 1.14
殘差 (A°S) 1486.66 18.00 82.59
緩衝時間(B) 199.06 2.00 99.53 2.41
殘差 (B°S) 742.09 18.00 41.23
材質°緩衝時間(A°B) 43.47 4.00 10.87 0.23
殘差 (A°B°S) 1693.30 36.00 47.04
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-11-2 緩衝期內的全部衝量
變異數分析摘要表變異來源 SS df MS F
受試者間 185634.77 9.00
受試者內 37035.77 80.00
材質(A) 8234.65 2.00 4117.32 10.80 *
殘差 (A°S) 6863.38 18.00 381.30
緩衝時間(B) 12942.50 2.00 6471.25 25.73 *
殘差 (B°S) 4527.98 18.00 251.55
材質°緩衝時間(A°B) 538.93 4.00 134.73 1.23
殘差 (A°B°S) 3928.33 36.00 109.12
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-11 之事後比較達顯著差異的組別
緩衝期內*+ b1a1-b1a3;b1a2-b1a3;b2a1-b2a3;b2a2-b2a3;b3a1-b3a3;b3a2-b3a3;
全部衝量 a1b1-a1b2;a1b1-a1b3; a1b2-a1b3;a2b1-a2b3;a2b2-a2b3;
a3b1-a3b2;a3b1-a3b3;a3b2-a3b3;
表 4-12-1 緩衝期髖關節最大水平受力
變異數分析摘要表變異來源 SS df MS F
受試者間 26.00 9.00
受試者內 89.27 80.00
材質(A) 37.51 2.00 18.76 55.17 *
殘差 (A°S) 6.12 18.00 0.34
緩衝時間(B) 4.59 2.00 2.30 2.77
殘差 (B°S) 14.92 18.00 0.83
材質°緩衝時間(A°B) 1.04 4.00 0.26 0.37
殘差 (A°B°S) 25.08 36.00 0.70
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-12-2 緩衝期膝關節最大水平受力
變異數分析摘要表變異來源 SS df MS F
受試者間 21.85 9.00
受試者內 72.31 80.00
材質(A) 25.87 2.00 12.93 62.97 *
殘差 (A°S) 3.70 18.00 0.21
緩衝時間(B) 4.64 2.00 2.32 3.04
殘差 (B°S) 13.75 18.00 0.76
材質°緩衝時間(A°B) 2.40 4.00 0.60 0.98
殘差 (A°B°S) 21.96 36.00 0.61
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-12-3 緩衝期踝關節最大水平受力
變異數分析摘要表變異來源 SS df MS F
受試者間 18.50 9.00
受試者內 72.95 80.00
材質(A) 24.44 2.00 12.22 75.25 *
殘差 (A°S) 2.92 18.00 0.16
緩衝時間(B) 4.58 2.00 2.29 2.99
殘差 (B°S) 13.79 18.00 0.77
材質°緩衝時間(A°B) 3.91 4.00 0.98 1.51
殘差 (A°B°S) 23.30 36.00 0.65
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-12 之事後比較達顯著差異的組別
髖+ b1a1-b1a2;b1a1-b1a3;b2a1-b2a2;b2a1-b2a3;b3a1-b3a2;b3a1-b3a3 膝+ b1a1-b1a2;b2a1-b2a2;b2a1-b2a3;b3a1-b3a2;b3a1-b3a3
踝+ b1a1-b1a2;b2a1-b2a2;b2a1-b2a3;b3a1-b3a2;b3a1-b3a3
表 4-13-1 緩衝期髖關節最大垂直受力
變異數分析摘要表變異來源 SS df MS F
受試者間 486.90 9.00
受試者內 807.61 80.00
材質(A) 27.41 2.00 13.71 1.67
殘差 (A°S) 147.97 18.00 8.22
緩衝時間(B) 136.37 2.00 68.18 12.79 *
殘差 (B°S) 95.97 18.00 5.33
材質°緩衝時間(A°B) 69.96 4.00 17.49 1.91
殘差 (A°B°S) 329.93 36.00 9.16
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-13-2 緩衝期膝關節最大垂直受力
變異數分析摘要表變異來源 SS df MS F
受試者間 510.23 9.00
受試者內 765.54 80.00
材質(A) 30.05 2.00 15.02 2.62
殘差 (A°S) 103.24 18.00 5.74
緩衝時間(B) 214.23 2.00 107.12 20.57 *
殘差 (B°S) 93.75 18.00 5.21
材質°緩衝時間(A°B) 59.78 4.00 14.94 2.03
殘差 (A°B°S) 264.49 36.00 7.35
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-13-3 緩衝期踝關節最大垂直受力
變異數分析摘要表變異來源 SS df MS F
受試者間 533.61 9.00
受試者內 818.17 80.00
材質(A) 25.22 2.00 12.61 1.99
殘差 (A°S) 114.26 18.00 6.35
緩衝時間(B) 217.89 2.00 108.94 18.38 *
殘差 (B°S) 106.69 18.00 5.93
材質°緩衝時間(A°B) 62.25 4.00 15.56 1.92
殘差 (A°B°S) 291.87 36.00 8.11
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-13 之事後比較達顯著差異的組別 髖* a3b1-a3b3;a3b2-a3b3
膝* a2b2-a2b3;a3b1-a3b3;a3b2-a3b3 踝* a2b2-a2b3;a3b1-a3b3;a3b2-a3b3
表 4-14-1 緩衝期髖關節伸肌的最大力矩
變異數分析摘要表變異來源 SS df MS F
受試者間 12.15 9.00
受試者內 86.76 80.00
材質(A) 24.67 2.00 12.33 18.42 *
殘差 (A°S) 12.05 18.00 0.67
緩衝時間(B) 4.20 2.00 2.10 2.75
殘差 (B°S) 13.77 18.00 0.76
材質°緩衝時間(A°B) 6.38 4.00 1.59 2.23
殘差 (A°B°S) 25.70 36.00 0.71
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-14-2 緩衝期膝關節伸肌的最大力矩
變異數分析摘要表變異來源 SS df MS F
受試者間 7.20 9.00
受試者內 9.08 80.00
材質(A) 0.47 2.00 0.24 2.73
殘差 (A°S) 1.57 18.00 0.09
緩衝時間(B) 0.72 2.00 0.36 5.83 *
殘差 (B°S) 1.11 18.00 0.06
材質°緩衝時間(A°B) 0.49 4.00 0.12 0.93
殘差 (A°B°S) 4.72 36.00 0.13
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-14-3 緩衝期踝關節伸肌的最大力矩
變異數分析摘要表變異來源 SS df MS F
受試者間 3.22 9.00
受試者內 8.70 80.00
材質(A) 2.69 2.00 1.35 29.50 *
殘差 (A°S) 0.82 18.00 0.05
緩衝時間(B) 0.15 2.00 0.07 0.85
殘差 (B°S) 1.54 18.00 0.09
材質°緩衝時間(A°B) 0.17 4.00 0.04 0.47
殘差 (A°B°S) 3.33 36.00 0.09
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-14 之事後比較達顯著差異的組別
髖+ b1a1-b1a2;b2a1-b2a2;b3a1-b3a2;b3a1-b3a3
膝* a3b1-a3b3
踝+ b1a1-b1a2;b1a1-b1a3;b3a1-b3a2;b3a1-b3a3
表 4-15-1 緩衝期髖關節伸肌的作功量
變異數分析摘要表變異來源 SS df MS F
受試者間 5.85 9.00
受試者內 5.77 80.00
材質(A) 0.92 2.00 0.46 7.98 *
殘差 (A°S) 1.03 18.00 0.06
緩衝時間(B) 0.17 2.00 0.08 1.20
殘差 (B°S) 1.26 18.00 0.07
材質°緩衝時間(A°B) 0.07 4.00 0.02 0.27
殘差 (A°B°S) 2.32 36.00 0.06
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-15-2 緩衝期膝關節伸肌的作功量
變異數分析摘要表變異來源 SS df MS F
受試者間 2.90 9.00
受試者內 10.18 80.00
材質(A) 0.17 2.00 0.09 0.66
殘差 (A°S) 2.32 18.00 0.13
緩衝時間(B) 1.06 2.00 0.53 4.60 *
殘差 (B°S) 2.08 18.00 0.12
材質°緩衝時間(A°B) 0.50 4.00 0.13 1.12
殘差 (A°B°S) 4.04 36.00 0.11
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-15-3 緩衝期踝關節伸肌的作功量
變異數分析摘要表變異來源 SS df MS F
受試者間 0.92 9.00
受試者內 1.85 80.00
材質(A) 0.23 2.00 0.1170 4.47 *
殘差 (A°S) 0.44 18.00 0.0247
緩衝時間(B) 0.07 2.00 0.0345 1.58
殘差 (B°S) 0.39 18.00 0.0218
材質°緩衝時間(A°B) 0.01 4.00 0.0024 0.12
殘差 (A°B°S) 0.70 36.00 0.0194
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-15 之事後比較達顯著差異的組別 髖+ b2a1-b2a3
膝* a2b1-a2b3;a2b2-a2b3
踝+ b1a1-b1a2
表 4-16-1 緩衝期股直肌的 IEMG
變異數分析摘要表變異來源 SS df MS F
受試者間 0.00208 9.00
受試者內 0.00152 80.00
材質(A) 0.00016 2.00 0.00008 2.15
殘差 (A°S) 0.00067 18.00 0.00004
緩衝時間(B) 0.00021 2.00 0.00011 11.79 *
殘差 (B°S) 0.00016 18.00 0.00001
材質°緩衝時間(A°B) 0.00002 4.00 0.00001 0.65
殘差 (A°B°S) 0.00030 36.00 0.00001
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-16-2 緩衝期股二頭肌的 IEMG
變異數分析摘要表變異來源 SS df MS F
受試者間 0.00020 9.00
受試者內 0.00034 80.00
材質(A) 0.00001 2.00 0.0000027 0.47
殘差 (A°S) 0.00010 18.00 0.0000058
緩衝時間(B) 0.00004 2.00 0.0000199 7.15 *
殘差 (B°S) 0.00005 18.00 0.0000028
材質°緩衝時間(A°B) 0.00000 4.00 0.0000011 0.30
殘差 (A°B°S) 0.00013 36.00 0.0000037
全 體 89.00 F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-16-3 緩衝期脛骨前肌的 IEMG
變異數分析摘要表變異來源 SS df MS F
受試者間 0.000372 9.00
受試者內 0.000347 80.00
材質(A) 0.000032 2.00 0.000016 2.50
殘差 (A°S) 0.000117 18.00 0.000006
緩衝時間(B) 0.000013 2.00 0.000006 1.70
殘差 (B°S) 0.000069 18.00 0.000004
材質°緩衝時間(A°B) 0.000014 4.00 0.000004 1.24
殘差 (A°B°S) 0.000102 36.00 0.000003
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-16-4 緩衝期腓腸肌的 IEMG
變異數分析摘要表
變異來源 SS df MS F
受試者間 0.000684 9.00
受試者內 0.000638 80.00
材質(A) 0.000010 2.00 0.000005 0.99
殘差 (A°S) 0.000088 18.00 0.000005
緩衝時間(B) 0.000134 2.00 0.000067 4.72 *
殘差 (B°S) 0.000254 18.00 0.000014
材質°緩衝時間(A°B) 0.000005 4.00 0.000001 0.28
殘差 (A°B°S) 0.000148 36.00 0.000004
全 體 89.00
F.95(2,18)=3.55 F.95(4,36)=2.64 *P<.05
表 4-16 之事後比較達顯著差異的組別
股直肌* a1b1-a1b2;a1b1-a1b3;a3b1-a3b2;a3b1-a3b3 股二頭肌* a2b1-a2b3;a3b1-a3b2
脛骨前肌
腓腸肌* a1b1-a1b2;a1b1-a1b3;a2b1-a2b2;a2b1-a2b3;a3b1-a3b2;a3b1-a3b3
附錄三:計算關節水平受力、垂直受力和關節肌肉力矩的 C 語言程式
#define max 100
#include <stdio.h>
main() {
int m,n,o,p,q,s,t,i,j,x,y,z;
float a,b,c,d,g,dtr;
float A[max][3],K[max][3],H[max][3];
float sc[max][7][2],sla[max][3][2],aa[max][6], sd[max][3],spd[max][3],r[max][9];
char ch,ch1,ch2,ch3;
static float BW[10]={{55.1},{51},{38.9},{48.5},{46},{38},{33.4},{42.9},{68.6},{57.1}};
/*受試者體重*/
float IOT=0.53764; /*大腿近端旋轉半徑與大腿長比值*/
float IOL=0.50382; /*小腿近端旋轉半徑與小腿長比值*/
float IOF=0.47784; /*腳掌近端旋轉半徑與腳掌長比值*/
static char name[10][6][6][30]/*分別讀入 10 位受試者的比例尺轉換後的座標
={ 檔、肢段線加速度檔、角加速度檔、肢段長度檔、處理後的測力板檔*/
/*輸出結果檔(包含各關節的 X、Y 方向受力及淨肌肉力矩)*/
{{{"1121.sc"},{"1121sla"},{"1121.aa"},{"1121.sd"},{"1121spd"},{"r1121"}}, {{"1122.sc"},{"1122sla"},{"1122.aa"},{"1122.sd"},{"1122spd"},{"r1122"}}, {{"1123.sc"},{"1123sla"},{"1123.aa"},{"1123.sd"},{"1123spd"},{"r1123"}}, {{"1124.sc"},{"1124sla"},{"1124.aa"},{"1124.sd"},{"1124spd"},{"r1124"}}, {{"1125.sc"},{"1125sla"},{"1125.aa"},{"1125.sd"},{"1125spd"},{"r1125"}}, {{"1126.sc"},{"1126sla"},{"1126.aa"},{"1126.sd"},{"1126spd"},{"r1126"}}}, {{{"1221.sc"},{"1221sla"},{"1221.aa"},{"1221.sd"},{"1221spd"},{"r1221"}}, {{"1222.sc"},{"1222sla"},{"1222.aa"},{"1222.sd"},{"1222spd"},{"r1222"}}, {{"1223.sc"},{"1223sla"},{"1223.aa"},{"1223.sd"},{"1223spd"},{"r1223"}}, {{"1224.sc"},{"1224sla"},{"1224.aa"},{"1224.sd"},{"1224spd"},{"r1224"}}, {{"1225.sc"},{"1225sla"},{"1225.aa"},{"1225.sd"},{"1225spd"},{"r1225"}}, {{"1226.sc"},{"1226sla"},{"1226.aa"},{"1226.sd"},{"1226spd"},{"r1226"}}}, {{{"1321.sc"},{"1321sla"},{"1321.aa"},{"1321.sd"},{"1321spd"},{"r1321"}},
{{"1322.sc"},{"1322sla"},{"1322.aa"},{"1322.sd"},{"1322spd"},{"r1322"}}, {{"1323.sc"},{"1323sla"},{"1323.aa"},{"1323.sd"},{"1323spd"},{"r1323"}}, {{"1324.sc"},{"1324sla"},{"1324.aa"},{"1324.sd"},{"1324spd"},{"r1324"}}, {{"1325.sc"},{"1325sla"},{"1325.aa"},{"1325.sd"},{"1325spd"},{"r1325"}}, {{"1326.sc"},{"1326sla"},{"1326.aa"},{"1326.sd"},{"1326spd"},{"r1326"}}}, {{{"4121.sc"},{"4121sla"},{"4121.aa"},{"4121.sd"},{"4121spd"},{"r4121"}}, {{"4122.sc"},{"4122sla"},{"4122.aa"},{"4122.sd"},{"4122spd"},{"r4122"}}, {{"4123.sc"},{"4123sla"},{"4123.aa"},{"4123.sd"},{"4123spd"},{"r4123"}}, {{"4124.sc"},{"4124sla"},{"4124.aa"},{"4124.sd"},{"4124spd"},{"r4124"}}, {{"4125.sc"},{"4125sla"},{"4125.aa"},{"4125.sd"},{"4125spd"},{"r4125"}}, {{"4126.sc"},{"4126sla"},{"4126.aa"},{"4126.sd"},{"4126spd"},{"r4126"}}}, {{{"5121.sc"},{"5121sla"},{"5121.aa"},{"5121.sd"},{"5121spd"},{"r5121"}}, {{"5122.sc"},{"5122sla"},{"5122.aa"},{"5122.sd"},{"5122spd"},{"r5122"}}, {{"5123.sc"},{"5123sla"},{"5123.aa"},{"5123.sd"},{"5123spd"},{"r5123"}}, {{"5124.sc"},{"5124sla"},{"5124.aa"},{"5124.sd"},{"5124spd"},{"r5124"}}, {{"5125.sc"},{"5125sla"},{"5125.aa"},{"5125.sd"},{"5125spd"},{"r5125"}}, {{"5126.sc"},{"5126sla"},{"5126.aa"},{"5126.sd"},{"5126spd"},{"r5126"}}}, {{{"6121.sc"},{"6121sla"},{"6121.aa"},{"6121.sd"},{"6121spd"},{"r6121"}}, {{"6122.sc"},{"6122sla"},{"6122.aa"},{"6122.sd"},{"6122spd"},{"r6122"}},
{{"6123.sc"},{"6123sla"},{"6123.aa"},{"6123.sd"},{"6123spd"},{"r6123"}}, {{"6124.sc"},{"6124sla"},{"6124.aa"},{"6124.sd"},{"6124spd"},{"r6124"}}, {{"6125.sc"},{"6125sla"},{"6125.aa"},{"6125.sd"},{"6125spd"},{"r6125"}}, {{"6126.sc"},{"6126sla"},{"6126.aa"},{"6126.sd"},{"6126spd"},{"r6126"}}}, {{{"7121.sc"},{"7121sla"},{"7121.aa"},{"7121.sd"},{"7121spd"},{"r7121"}}, {{"7122.sc"},{"7122sla"},{"7122.aa"},{"7122.sd"},{"7122spd"},{"r7122"}}, {{"7123.sc"},{"7123sla"},{"7123.aa"},{"7123.sd"},{"7123spd"},{"r7123"}}, {{"7124.sc"},{"7124sla"},{"7124.aa"},{"7124.sd"},{"7124spd"},{"r7124"}}, {{"7125.sc"},{"7125sla"},{"7125.aa"},{"7125.sd"},{"7125spd"},{"r7125"}}, {{"7126.sc"},{"7126sla"},{"7126.aa"},{"7126.sd"},{"7126spd"},{"r7126"}}}, {{{"8121.sc"},{"8121sla"},{"8121.aa"},{"8121.sd"},{"8121spd"},{"r8121"}}, {{"8122.sc"},{"8122sla"},{"8122.aa"},{"8122.sd"},{"8122spd"},{"r8122"}}, {{"8123.sc"},{"8123sla"},{"8123.aa"},{"8123.sd"},{"8123spd"},{"r8123"}}, {{"8124.sc"},{"8124sla"},{"8124.aa"},{"8124.sd"},{"8124spd"},{"r8124"}}, {{"8125.sc"},{"8125sla"},{"8125.aa"},{"8125.sd"},{"8125spd"},{"r8125"}}, {{"8126.sc"},{"8126sla"},{"8126.aa"},{"8126.sd"},{"8126spd"},{"r8126"}}}, {{{"9121.sc"},{"9121sla"},{"9121.aa"},{"9121.sd"},{"9121spd"},{"r9121"}}, {{"9122.sc"},{"9122sla"},{"9122.aa"},{"9122.sd"},{"9122spd"},{"r9122"}}, {{"9123.sc"},{"9123sla"},{"9123.aa"},{"9123.sd"},{"9123spd"},{"r9123"}},
{{"9124.sc"},{"9124sla"},{"9124.aa"},{"9124.sd"},{"9124spd"},{"r9124"}}, {{"9125.sc"},{"9125sla"},{"9125.aa"},{"9125.sd"},{"9125spd"},{"r9125"}}, {{"9126.sc"},{"9126sla"},{"9126.aa"},{"9126.sd"},{"9126spd"},{"r9126"}}}, {{{"1021.sc"},{"1021sla"},{"1021.aa"},{"1021.sd"},{"1021spd"},{"r1021"}}, {{"1022.sc"},{"1022sla"},{"1022.aa"},{"1022.sd"},{"1022spd"},{"r1022"}}, {{"1023.sc"},{"1023sla"},{"1023.aa"},{"1023.sd"},{"1023spd"},{"r1023"}}, {{"1024.sc"},{"1024sla"},{"1024.aa"},{"1024.sd"},{"1024spd"},{"r1024"}}, {{"1025.sc"},{"1025sla"},{"1025.aa"},{"1025.sd"},{"1025spd"},{"r1025"}}, {{"1026.sc"},{"1026sla"},{"1026.aa"},{"1026.sd"},{"1026spd"},{"r1026"}}}, };
FILE *fp0,*fp1,*fp2,*fp3,*fp4,*fp5;
clrscr();
printf("\n\n\n\n\n\n");
g=9.81;
dtr=57.29578; /* change degree to radius */
/*************************************************************************/
for(o=0;o<10;o++) {
for(p=0;p<6;p++)
{
if ((fp0 = fopen(name[o][p][0],"r")) == NULL) {
printf("Can't open the file named %s.",name[o][p][0]);
delay(600);
while ( getchar() !='\n' );
return;
} m=0;
while(!feof(fp0)) {
for(i=0;i<7;i++) {
fscanf(fp0,"%f%c%f%c%",&a,&ch1,&b,&ch2);
sc[m][i][0]=a;sc[m][i][1]=b;
}
fscanf(fp0,"\n");
m+=1;
}
fclose(fp0);
/***********************************************************/
if ((fp1 = fopen(name[o][p][1],"r")) == NULL) {
printf("Can't open the file named %s.",name[o][p][1]);
delay(600);
while ( getchar() !='\n' );
return;
} m=0;
while(!feof(fp1)) {
for(i=0;i<3;i++) {
fscanf(fp1,"%f%c%f%c",&a,&ch1,&b,&ch2);
sla[m][i][0]=a;sla[m][i][1]=b;
}
fscanf(fp1,"\n");
m+=1;
}
fclose(fp1);
/***************************************************/
if ((fp2 = fopen(name[o][p][2],"r")) == NULL) {
printf("Can't open the file named %s.",name[o][p][2]);
delay(600);
while ( getchar() !='\n' );
return;
} m=0;
while(!feof(fp2)) {
for(j=0;j<6;j++) {
fscanf(fp2,"%f%c",&d,&ch1);
aa[m][j]=d;
}
fscanf(fp2,"\n");
m+=1;
}
fclose(fp2);
/****************************************************/
if ((fp3 = fopen(name[o][p][3],"r")) == NULL) {
printf("Can't open the file named %s.",name[o][p][3]);
delay(600);
while ( getchar() !='\n' );
return;
} m=0;
while(!feof(fp3)) {
for(j=0;j<3;j++) {
fscanf(fp3,"%f%c",&d,&ch1);
sd[m][j]=d;
}
fscanf(fp3,"\n");
m+=1;
}
fclose(fp3);
/********************************************************/
if ((fp4 = fopen(name[o][p][4],"r")) == NULL) {
printf("Can't open the file named %s.",name[o][p][4]);
delay(600);
while ( getchar() !='\n' );
return;
} m=0;
while(!feof(fp4)) {
for(i=0;i<3;i++) {
fscanf(fp4,"%f%c",&a,&ch1);
spd[m][i]=a;
}
fscanf(fp4,"\n");
m+=1;
}
fclose(fp4);
/********************************************************************/
fp5=fopen(name[o][p][5],"w");
for(s=0;s<m;s++) {
A[s][0]=BW[o]*0.020*sla[s][2][0]-spd[s][0]/2;
/*force of X direction at the ankle*/
A[s][1]=BW[o]*0.020*(sla[s][2][1]+g)-spd[s][1]/2;
/*force of Y direction at the ankle*/
A[s][2]=BW[o]*0.020*sd[s][2]*IOF*sd[s][2]*IOF *aa[s][2]/dtr+spd[s][0]/2*(sc[s][3][1]
-sc[s][4][1]+0.182)
+spd[s][1]/2*(sc[s][5][0]+spd[s][2]-0.45-sc[s][3][0]) -BW[o]*0.020*g*(sc[s][4][0]-sc[s][3][0])*0.41278;
/*moment of muscle at the ankle*/
K[s][0]=BW[o]*0.053*sla[s][1][0]+A[s][0];
/*force of X direction at the knee*/
K[s][1]=BW[o]*0.053*(sla[s][1][1]+g)+A[s][1];
/*force of Y direction at the knee*/
K[s][2]=BW[o]*0.053*sd[s][1]*IOL*sd[s][1]*IOL
*aa[s][1]/dtr-A[s][2]
+A[s][0]*(sc[s][2][1]-sc[s][3][1]) -A[s][1]*(sc[s][2][0]-sc[s][3][0])
-BW[o]*0.053*g*(sc[s][2][0]-sc[s][3][0])*0.4166;
/*moment of muscle at the knee*/
H[s][0]=BW[o]*0.110*sla[s][0][0]+K[s][0];
/*force of X direction at the hip*/
H[s][1]=BW[o]*0.110*(sla[s][0][1]+g)+K[s][1];
/*force of Y direction at the hip*/
H[s][2]=BW[o]*0.110*sd[s][0]*IOT*sd[s][0]*IOT
*aa[s][0]/dtr-K[s][2]
-K[s][0]*(sc[s][1][1]-sc[s][2][1]) -K[s][1]*(sc[s][2][0]-sc[s][1][0])
-BW[o]*0.110*g*(sc[s][2][0]-sc[s][1][0])*0.462;
/*moment of muscle at the hip*/
/**********************************************************************/
fprintf(fp5,"%.2f,%.2f,%.2f,%.2f,%.2f,%.2f,%.2f,%.2f,%.2f",
H[s][0]/BW[o],H[s][1]/BW[o],H[s][2]/BW[o], K[s][0]/BW[o],K[s][1]/BW[o],K[s][2]/BW[o], A[s][0]/BW[o],A[s][1]/BW[o],A[s][2]/BW[o]);
fprintf(fp5,"\n");
}
fclose(fp5);
}
} }