at. ‘ .- .:.;...,.3._ ‘ . ...ar 13.], 53:: : 1.. c i he i I s . . ax “Hmtki. $59.... I. :1 3.! .5 ”3 fr .w....... 5:. x .23.: .3 .I xtfifiafl. 3:. 8321‘. ink {A .. a E“: .l... a: flu . 55.1... Earn, ‘3! E. x: . a x 3.9:"...lssfl5 :l r} . :8: .I. 7 q . 1; 30¢ I 1...! 3:: 4‘ :..n...i..1 x I. :4 .f s :5... 1:! TREE:- LIBRAR‘s’ i iiyfiacl'aigan shah . University ( r3 "W5 1 h ‘ This is to certify that the thesis entitled KINETIC MODELING AND SEAT FACTORS RELATING TO MANIKlN-SEAT INTERACTION IN AUTOMOBILES presented by Nathan J. Radcliffe has been accepted towards fulfillment of the requirements for MS Engineering Mechanics degree in Major professor 8-23-01 Date 0-7639 MS U is an Affirmative Action/Equal Opportunity Institution PLACE IN RETURN Box to remove this checkout from your record. TO AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 8/01 cJClFlC/DatoDuepes-sz KINETIC MODELING AND SEAT FACTORS RELATING TO MANIKIN—SEAT INTERACTION IN AUTOMOBILES By Nathan J. Radcliffe AN ABSTRACT OF A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Materials Science and Mechanics 2001 Dr. Robert Hubbard ABSTRACT KINETIC MODELING AND SEAT FACTORS RELATING TO MANIKIN-SEAT INTERACTION IN AUTOMOBILES By Nathan J. Radcliffe Automotive seats are integral components in the interactions between occupants and vehicles, and ultimately the road surface traveled. Seats are important in determining the comfort of the ride, and therefore the user’s impression of the total vehicle. The occupant’s static deflected position, supported posture and comfort may be related to a set of quantitative values known as Seat Factors. Seat Factors are quantitative measures and calculated values of structural properties for complete seats related to biomechanical interactions. Using specific procedures, measurable Seat Factors can be attained which relate to the physical properties of each seat pan. These measured quantities may then be used in a kinetic model, which calculates additional Seat Factors that can be used to predict the static deflection and rotation response of the buttocks and thigh region in contact with the seat pan, due to imposed loads at the hip joint center and knee. This thesis provides information regarding the procedures used to determine Seat Factors, assumptions used in modeling, kinetic models and components, model validation techniques, applications and future work. Copyright by Nathan J. Radcliffe © 2001 ACKNOWLEDGEMENTS I would like to thank Dr. Hubbard for positively influencing my engineering career over the past six years. Your positive guidance, mentoring, fi’iendship, and support have guided me in the right direction through the often tough but always encouraging years of engineering school at Michigan State University. Thank you to my master’s committee members, Dr. Reid-Bush, and Dr. Rosenberg. Your help and guidance has helped this project and my education progress successfully through both my undergraduate and graduate programs. I would like to thank my family for pushing me to succeed from the very beginning. Mom, Dad, Isa, and Eliott, I thank you for all of your loving support. And finally, I would especially like to thank my fiancé Terri for her loving patience and support without which I would be further challenged to succeed in this project and beyond. TABLE OF CONTENTS LIST OF TABLES VII LIST OF FIGURES Ix NOMENCLATURE XI INTRODUCTION 1 EXPERIMENTAL DATA ACOUIsrrION METHODS 2 SAE 1826 W’ Seating Manikin 2 ASPECT Seating Manikin 5 KINETIC MODELS a Model I. Two-spring kinetic model with support forces under the H-pt. and knee. 6 Model I]. Two-spring kinetic model with single spring constant, support forces under the H-pt. and thigh at a fixed distance ETSL from the H-pt. 7 Static Mathematical Analysis for Physical Simulation of Model 11 using a Dynamic Solver Software 8 Results from Model I] 11 Model III. Two-spring kinetic model with variable spring constant support forces under the H-pt. and thigh at a variable distance ETSL from the H-pt. 13 Results from Model III 14 Model IV. Coupled force and moment kinetic model 16 Static Mathematical Analysis for Mathematical Simulation of Model IV using MS Excel Solver to Optimize Calculated Seat Factors 17 Mathematical Modeling Methods 19 Step 1. Experimental Data Input 19 Step 2. Static Reaction Force and Moment Calculation 20 Step 3. Optirnizatim of Modified Taylor Series-Based Seat Factor Matrices 21 Step 4. Solving Equations to Determine Simulation Results Based on Calculated Seat Factor Matrices 23 Step 5. Graphical Analysis of Optimization and Simulatim Results 24 Step 6. Occupant Positim Approxrmat' ions 29 Results and Conclusions from Model IV 31 CONCLUSION 39 Future Work 40 APPENDICES APPENDIX A 45 APPENDIX B 168 APPENDIX C 170 REFERENCES 183 vi LIST OF TABLES Table 1. Input data acquired fiom experimental incremental loading of the rmnikin is shown in white font within dark cells, along with calculated values based on experimental data shown in grey. 20 Table 2. Individual Deflection-Rotation Solver for Occupant Position Approximations. 29 Table 3. Average H-pt. deflection, knee deflection and relative thigh rotation difference across experimental loading range for linear, second, and third order Seat Factor Matrices. 32 Table 4. RMS and Percent Error vs. Order of Taylor Series Expansion. 37 Table L-l . Seat Key for seats analyzed within the Seat Factor Solver Spreadsheet.____47 Table A-1. Seat A SFS Analysis Worksheet. 48 Table A-2. Seat B SFS Analysis Worksheet. 52 Table A-3. Seat C SF S Analysis Worksheet. 56 Table A-4. Seat D SFS Analysis Worksheet. 60 Table A-5. Seat E SFS Analysis Worksheet. 64 Table A-6. Seat F SF S Amlysis Worksheet. 68 Table A-7. Seat G SFS Analysis Worksheet. 72 Table A-8. Seat H SF S Analysis Worksheet. 76 Table A-9. Seat I SFS Analysis Worksheet. 80 Table A-lO. Seat J SFS Analysis Worksheet. 84 Table A-1 1. Seat K SFS Analysis Worksheet. 88 Table A-12. Seat L SFS Analysis Worksheet. 92 Table A-13. Seat M SFS Analysis Worksheet. 96 Table A-14. Seat N SFS Analysis Worksheet. 100 Table A-15. Seat 0 SFS Analysis Worksheet. 104 Table A-l6. Seat P SF 8 Analysis Worksheet. 108 vii Table A-17 Table A-1 8. Table A-1 9. Table A-20. Table A—21. Table A-22. Table A-23. Table A-24. Table A-25. Table A-26. Table A-27. Table A-28. Table A-29. Table A-30. . Seat Q SF S Analysis Worksheet. Seat R SFS Analysis Worksheet. Seat S SFS Analysis Worksheet. Seat 4 SFS Analysis Worksheet. Seat 8 SF S Analysis Worksheet. Seat 9 SF S Analysis Worksheet. Seat Neon/ASPECT SFS Analysis Worksheet. Seat TandC/ASPECT SFS Analysis Worksheet. Seat Saturn/ASPECT SFS Analysis Worksheet. Seat TandC/J826 SF S Analysis Worksheet. Seat LHtan/ASPECT SF S Analysis Worksheet. Seat Saturn/1826 SFS Analysis Worksheet. Seat Neon/J826 SFS Analysis Worksheet. Seat LHtan/J826 SFS Analysis Worksheet. viii 112 116 120 124 128 132 136 140 144 148 152 156 160 164 LIST OF FIGURES Figure 1. SAE 1826/0scar manikin (Unloaded - Loading Step 1). 4 Figure 2. ASPECT manikin (Unloaded). 5 Figure 3. Two-spring kinetic model with support forces under the H-pt. and knee. 6 Figure 4. Two-spring kinetic model with single spring constant, support forces under the H- pt and thigh at distance ETSL fi'om the H-pt. 7 Figure 5. Descriptive two-dimensional figure of H-pt.-knee rotation about the ETSL for similar triangle-based analysis of the ETSL. 9 Figure 6. Descriptive two-dimensional plot of H-pt.-knee spring constant linear magnitude relationship for calculation of ETSS at the ETSL. ' 10 Figure 7 a-d. Examples of Working Model results fi'om physical simulations of Model H. Plots indicate good experimental-simulation results for a linear seat (Seat G) in Figures 7a and 7b and poor experimental—simulation results fi'om a non-linear seat (Seat A) in Figures 7c and 7d. 11 Figure 8. Two-spring kinetic model with variable spring constant support forces under the H-pt. and thigh at a variable distance ETSL from the H-pt. 13 Figure 9 a-d. Examples of Working Model generated results fi'om physical simulations of Model III with results of Model 11 shown for comparison. Plots indicate experimental-simulation results for a linear seat (Seat G) in Figures 9a and 9b and experimental-simulation results from a non-linear seat (Seat A) in Figures 9c and 9d. ' 15 Figure 10. Coupled force and moment kinetic model. 16 Figure 11 a,b. Plots for comparison of experimental results (diamonds) to simulation results based upon third order Taylor series expansion matrices (line) (Seat A). 25 Figure 12. Thigh position as a fimction of incremental H-pt. and knee loading. The graphical point of thigh rotation (zero deflection) during H-pt. loading is approximately 225mm fi'om the H-pt. along the thigh for this seat (Seat A). 26 Figure 13 a,b. Third order surface models of reaction force and moment as a function of H- pt. deflection and relative thigh rotation (Seat A). 28 Figure 14. Interpolated and extrapolated H-pt. force surface Plot, shown with extrapolation and interpolation beyond the path of experimental data (Seat 0). 33 Figure 15. Two-dimensional contour of the SAE 1826 Manikin [4]. A ‘/4 inch, 2-D plastic template with rods extending from the H-pt. was used as a reference for data acquisition by Hubbard and Gedratis [1] and Hubbard and Radclifl‘e [6] within this project [4]. 35 Figure 16. Proposed two-dimensional linked system of reaction forces and moments (dark arrows) with input loading (white arrows) at crucial joints 42 Figure C-l . Descriptive two—dimensional figure of H-pt.-knee deflection and rotation for analysis of the ETSS. 182 NOMENCLATURE Arabic Symbok and Acronyms ETSL = efl‘ective thigh support location ETSS = effective thigh support stimress F = generic reaction force under the H-pt (could be RH or F") Fn = simulation reaction force under the H-pt. FH = input force at the H-pt. FK = input force at the knee FT = input force at the thigh n K fl, = nm order reaction force-deflection stifiness coeficient n K ,9 = nth order reaction force-rotation stimtess coeficient n K m 6 = 11‘” order reaction moment-deflection stiffness coeflicient " K m 9 = 11‘11 order reaction moment-rotation stifliress coefficient Kn = spring constant of the seat under the H-pt. KK = spring constant of the seat under the knee (physically non-existent) LK = distance from the H-pt. to the knee along the thigh xi Mn = SF = Greek Symbob 5" = .9" = 9 THIGH = large male (90th percentile for height and weight, 180 lbs.) generic reaction moment around the H-pt (could be Mu or Mn). simulation reaction moment around the H-pt. experimental reaction moment around the H-pt. mid male (50* percentile for height and weight, 170 lbs.) experimental reaction force under the H-pt. small female (10til percentile for height and weight, 115 lbs.) relative deflection of the H-pt. (same as AH) relative rotation of the thigh (referenced to the 2-D template) relative deflection of the H-pt. (referenced to 2-D template) relative deflection of the knee (referenced to 2-D template) thigh angle referenced to horizontal xii INTRODUCTION The pln'pose of this investigation was to develop an experimental data interpretation technique coupled with a kinetic modeling method to quantitatively describe and evaluate the mechanical properties of automotive seats. These quantitative measures and calculations are known as Seat Factors [lfit was theorized that ifa modeling method could be developed which produced similarresults to data acquired experimentally when identical loading procedures were followed, the model and a calculated set of Seat Factors were validated and could be used to define the mechanical characteristics of a particular seatét‘The evolution of this project included several different kinetic models, and their validatioh using experimental results fiom 30 production and prototype automotive seats, which is discussed in upcoming sections. For the initial kinetic modeling, static equations were derived to solve for calculated Seat Factors, and physical modeling was attempted for the purpose of validation using Working ModelTM (WM) [2]. {I Working Model“ is a two-dimensional, dynamic simulation software, which has the abrlrtjito predict motion of rigid bodies with imposed forces and constraints. Though this physical seat factor validation technique provided much insight, the most successful kinetic modeling method involved experimental data acquisition coupled with mathematical analysis using Microsoft (MS) Excel Solver [3] to optimize modified Taylor series expanded, two by two descriptive seat factor matrices based upon the experimental response of the manikin loaded incrementally into the seat] Traditionally, seating design is based upon analysis of comfort, occupant posture, contact surface pressure distribution, loaded manikin behavior, and software model positioning timing the computer-aided design process. Though these techniques provide designers with valuable insights into the interactions between people and seats, these methods do not lead to a complete, quantitative set of Seat Factors that defines the structural characteristics ofthe seat}; was theorized that a simplified, condensed set of Seat Factors coupled with a kinetic modeling approach could be used as a predictive tool for occupant positioning, and may also be used to correlate quantitative Seat Factors to subjective responses in terms of comfort. 6he latter of the two proposed applications will _ ”may (‘ be considered for future work, and will not be investigated within this thesis. ‘1‘ 1 Previous work in seating biomechanics and kinetic modeling by team members of the Biomchanical Design Research Laboratory at Michigan State University has lead to the definition of Seat Factor data acquisition techniques [1] and the development of the ASPECT manikin [5] through support of the ASPECT program. The ASPECT Manikin is a replacement for SAE’s J826 "Oscar" maniklnEThis manikin has an articulated back section that can be used to take various manikin-seat interaction measurements including hip joint center location, thorax, pelvis and thigh angle, seat back and seat pan angle, and lumbar prominencei'i EXPERIMENTAL DATA ACQUISITION METHODS SAE J826 “Oscar” Seating Manikin Gil a previous work by Hubbard and Gedraitis [l ], an experimental technique for measured seat factor data acquisition with the SAE J 826 [4] manikin was used. This technique is presented in detail within SAE Paper Number 1999-01-09'5’31‘ but can be ._..——J Ar" summarized in the following steps: 1. Using the SAE J826 “Oscar” seating manikin (shown in Figure 1), incremental loads were measured at the H-point (H-pt.) and knee axes. Note: H-point (H-pt.) is a term commonly used within the automotive industry, which relates to the location of the H-pt. axis on seating manikins when a specific loading procedure is followed [4]. The H-pt. axis is analogous to the hip joint center (HJC) of a human. References to the H-pt. within this report refer to the location of the physical H-pt. axis on the mnikin, and not the position of the H-pt. axis according to the established H-pt. procedure [4]. This was accomplished by suspending the seating manikin along the H-pt and knee axes by tensile lines with forces measured using voltage-output load cells. Incrementally loading the H-pt axis while holding a constant, minimal knee load, and subsequently, incrementally loading the knee axis with no additional applied H-pt load investigated a range of loads similar to a person sitting down into a seat. The measured loads could then be known for that particular incremental loading sequence. I” u“. ,doJ *'.-_ Figure l. SAE J826/Oscar manikin (Unleaded - Loading Step 1) f”. I, ; (o; II(A reference position in each seat was measured using the 2-D plastic SAE J826 Iemplate [1,4])The manikin without additional weight was then placed on top of the seat pan surface with the rear of the manikin seated firmly against the lower portion of the seatbackCPrior to measuring the position of the loaded SAE 1826 manikin in the seat, the torso/back section of the manikin was adjusted to vertical for each , wa incremental load. This minimized interaction between the back of the manikin and the seat back, and ensured consistency between trials. Known loads were incrementally added, and the positions of the H-pt and knee axes were measured afier each applied load. H-pt. axis measurements were made on the left and right side by measuring the position at the tip of a rod extending out fiom the H-pt. axis center. Left and right H-pt and knee axis measurements were measured at the same distance from the vertical plane of symmetry of the manikin. By averaging the left and right measures, the variation in level on each side of the manikin was averaged. Average measures, which corresponded to deflections at the plane of symmetry of the manikin, were used for simulation data input purposes. Measurements of H-pt and knee positions could be acquired using a motion analysis system or any method of vertical position measurement with a fixed reference i.e. floor [1]. ASPECT Seating Manikin Data acquisition using the prototype ASPECT manikin [5], butt-thigh segment (ABT), was accomplished in a similar manner as the procedure discussed for the SAE 1826 manikin. Only the butt-thigh region of the ASPECT manikin was used for data acquisition? order to remove interaction between the back section of the manikin and , , ‘ the scat [6T Data acquisition techniques and results ofASPECT manikin testing are presented in detail within Appendix B: Seat Pan Load-Deflection Measurements For Determination of Seat Factors with the ASPECT Butt-Thigh (AB7) Segment of the ASPECT Manikin Prototype at MSU and the SAE J826 Manikin by Robert Hubbard and Nathan Radcliffe. Results of ASPECT manikin data acquisition and simulation are discussed with comparisons to the SAE 1826 within Appendix B. ABT frOI'lt _ f 3 w, p: .> _' d Figure 2. ASPECT manikin (Unleaded) KINETIC MODELS During the course of work on this project, four separate kinetic models were successively developed. Each of the models are discussed in this section with varying levels of description based upon their relevance to the final outcome of the project. Model I. Two-spring kinetic model with support forces under the H-pt. and knee. Slider prohibits horizontal translation of the H-pt. Flt H-pt Figure 3. Two-spring kinetic model with support forces under the H-pt. and knee Model I, a two-spring kinetic model with support forces under the H-pt. and knee, was the first model developed. It separated the complex, distributed load interaction between an automotive seat and a manikin into two input forces and two support forces at the H-point and knee. This model was a simple way to represent experimental behavior, as it was trivial to calculate the spring constants under the H-pt. and knee with a measured force-deflection response. It was realized prior to the onset of analysis and modeling that a support force directly under the knee did not allow the thigh to rotate about a point between the H-pt. and knee as observed in measurement by Hubbard and Gedraitis [1]. In order to allow for realistic rotation behavior of the thigh in a kinetic model, the support must be positioned at some distance known as the effective thigh support location (ETSL) away fi'om the H-pt. along the thigh, which lead to the development of Model II. Model II. Two-spring kinetic model with single spring constant, support forces under the H-pt. and thigh at a fixed distance ETSL from the H-pt. 1 Slider prohibits horizontal translation of the H-pt.F K FH L K . | Knee H-Pt ’ 0; ETSL ETSS ;_ Kn — Figure 4. Two-spring kinetic model with single spring constant, support forces under the H-pt and thigh at distance ETSL from the H—pt. In Model H, the first support force was located at the same point under the H-pt. as seen in Model I. This model differed fiom Model I, in that the second support force was located at some distance known as the eflective thigh support location (ETSL) away fi'om the H-pt., along the thigh In designing this kinetic model, it was hypothesized that the second support would exist at the point of zero deflection (ETSL; thigh pivot point referenced to original thigh angle) during H-pt. loading: Based upon experimental measurements [1] and calculations, the distance fi'om the H-pt. to the point of rotation along the thigh (ETSL), rennined relatively constant during H-pt. loading and typically stabilized at larger imposed loads. Therefore, a single ETSL was chosen based upon a t middle point in the incremental H-pt. loading range to simulate the hill-range, 9 ' incre -loaded response of the manikin. The single ETSL used for all imposed loads in simulation was calculated based upon the point of rotation for loads of494N on the H-pt., and 29N on the mi? Based upon experimental data acquired by Hubbard and Gedraitis [1], one spring constant for each of the two support springs was calculated to represent the behavior of the seat throughout all imposed loads in the software simulation. These spring constants were referred to as K" (under H-pt.) and ETSS (under thigh) or effective thigh support stiflhess. The H-pt. spring constant (Kn) was calculated based upon the relative force- deflection response at 494N on the H-pt., and 29N on the knee. The ETSS used for simulation was an average of the calculated ETSS during the full range of incremental knee loadings. Static Mathematical Analysis for Physical Simulation of Model H using a Dynamic Solver Software Using static analysis, the following equations were derived to calculate individual Seat Factors relevant to the kinetic simulation of Model I] (and Model 111) below. The spring constants KH and KK were based upon relative, incremental characteristics of the force-deflection response of a particular seat. Each of these spring constants were the incremental input force (in this case referenced to the SAE 2-D template [4], which was assumed to have negligible mass) divided by the relative deflection referenced to the initial position of the H-pt. and knee as established by the SAE 2-D template. The effective thigh support location (ETSL) was at the point of zero deflection (instantamous pivot) during H-pt. loading with constant nominal knee loads, referenced to the initial position of the SAE J826 template [4]. When the H-pt. deflected into the seat and the knee elevated during incremental H-pt. loading, there was always a point of zero deflection. Calculation of the ETSL was done for H—pt. loading only. The ETSL used dining H-pt. loading was assumed to be the constant ETSL for incremental knee loading. This point was calculated by similar triangle-based analysis as derived below. — Initial Position — Position after H-pt. Loading Figure 5. Descriptive two-dimensional figure of H—pt-knee rotation about the ETSL for similar triangle-based analysis of the ETSL LK _AH-AK ETSL AH AH*LK ETSL:— AH—AK (1) In order to calculate the efi‘ective thigh support stifi‘ness at a point of support between the H-pt. and knee, the equation derived below using proportional analysis of similar triangles was used. This ETSS equation was only valid during knee loading, but the average ETSS value calculated from experimental knee loading was used as a constant for incremental, simulated H-pt. loading. Assumption: Linear change in spring constant magnitude between H-pt. and knee Magnitude of K“ H-pt. + Lx AH- AK ETSL Figure 6. Descriptive two-dimensional plot of H-pL-knee spring constant linear magnitude relationship for calculation of ETSS at the ETSL. K” —Erss _ K” -KK ETSL LK ETSS: K” — ETSL“? 'K") (2) K Equation 2 for ETSS was based upon the assumption that the spring constant at the location of the thigh support was linearly related to the spring constants under the H- pt. and knee (virtual, as the knee of the manikin extends beyond the support of the seat) according to the effective location of the thigh support between the H-pt. and knee along the thigh. Typically, automotive seats decrease in stiffness as you move from the support under the H-pt. toward the fiont of the seat underthethigh The constant value used for KH was the last calculated spring constant during H-pt. loading with constant nominal knee loads. An alternative equation for determining ETSS proposed after Model II analysis, which was not fully investigated due to the success of Model IV, is given in Appendix C. 10 Results from Model H Whenasinglespringconstantwasusedforeachspringineth[2]trial, along with a constant ETSL, the results of all physical simulations for Model H yielded linearly varying deflections at the H-pt. and knee as a function of applied load. Two sample sets of typical results for Model H are shown below in Figure 7 a-d. Simulation results for an approximately linear seat (Seat G) and a stiffening, non-linear seat (Seat A) are shown. These two seats represent good (Seat G) and poor (Seat A) experimental- sirnulation correlations and they were indicative of the range of results of Model H. mefi) KnoeDclIectlonIG) ms 235 335 435 585 see 0 50 100 150 200 250 300 o . A 40 E .5- +ModaIIl v; 20‘ +Modalll : -1: \ +Exp : o j—\.—I- , +59 E ' E \ é -20 5 -20 \ -25 g ,0 x t g 4° \. -35 l -80 i .40 ‘ a0 1 H-pt.Loadino Fema(N) KneeLoad‘ng Force (N) Figure 7a Figure7b WWNW KnceDcsectlontAI 135 235 ass 485 585 ass 0 so 100 150 200 250 300 A 0 ' A 30 -5 +M°°°'" ‘ 20. +Modalll 3 _ \4‘\ g 10.. I 10 \ \ +5“) I O \ +531) E45 ‘Q . £540 \ \ 1 cm \ A? .20 \ \ l -25 -30 \ \_ so \ \‘ '35 H-ptl-owimFomem) - KnceLoadanoroerN) Figure 7c Figure 7d Figure 7 a-d. Examples of Working Model results from physical simulations of Model H. Plots indicate good experimental-simulation results for a linear seat (Seat G) in Figures 7a and 7b and poor experimental-simulation results from a non-linear seat (Seat A) in Figures 7c and 7d. For seats that exhibited a linear force-deflection response in experimental testing (as with Seat G), simulation results were similar to the measured responses. However, seats that exhibited non-linear (typically stiffening) forceodeflection responses during experimental manikin loading (as with Seat A) did not yield results that were acceptably similar between the experimental response and kinetic model physical simulation. This was due to the choice of a single spring constant calculated fiom the force-deflection results of a middle loading step (494 N of H-pt. force and 29 N of knee force). A single spring constant seat factor could not effectively simulate the soft behavior of the stiffening seats when small, initial loads were applied. In Figure 7c, the large load stifliless under the H-pt. of Seat A is similar for both the experimental data and physical simulation, as indicated by approximately parallel lines between 385 N and 585 N OfH- pt. loading force. The offset difference rs due to the imbility of a single (linear) spring constant to effectively simulate a stiffening (non-linear) seat with soft properties at small, . impqsed loads- .. Further downward deflection of the H-pt. during H-pt. loading in experimental results caused the knee to deflect upward by means of rotation prior to the onset of knee loading. This phenomenon accounted for the offset in knee deflection during knee loading. In other words, differences between Model H simulation and experimental results of H-pt. deflection influenced Model H physical simulation of knee deflection. 12 Model HI. Two-spring kinetic model with variable spring constant support forces under the H-pt. and thigh at a variable distance ETSL from the H-pt. Slider prohibits horizontal translation of the H-pt. H-pt ETSS(FH) Slider allows for variation of ETSL Figure 8. Two-spring kinetic model with variable spring constant support forces under the H-pt. and thigh at a variable distance ETSL from the H-pt. Analysis of Model H simulations suggested that a variable spring constant model, with variable ETSL would likely produce better results since Seat Factors varied with incremental loads. Model HI was investigated extensively with physical simulations in Working Model”. Using WM’s data import abilities, simulations were run with variable spring constants KH(FH), ETSS (FH, Fk), and ETSL (FH ) as a function of the applied loads on the H-pt. and knee. The choice of variable spring constants and a variable ETSL required data exchange between WM and a spreadsheet, along with the construction of a sliding thigh support under the thigh Each of these factored dramatically into the complexity of this physical model. Mathematically, the same equations presented in the description of Model H were applied directly to the physical simulation of Model IH. The difference existed in the selection of variable K", ETSL and ETSS as fimctions of incremental loads. Since KH ,.- 13 l “ F“ and ETSL were calculated based only upon experimental deflection results of H-pt. -m-..’_—_~_.r-. .~_ loading with constant knee loads, the final coefficient of each of these terms was used as a constant coefficient during knee loading. Comparably, for the ETSS, the first - coefficient calculated fi‘om incremental knee loading with constant H-pt. loading was i used to represent the ETSS durinng-pt. loading with. constant knee loading. fl’-‘—-—‘——-v~- —_ I ""‘ Results from Model III When variable spring constants were used for each spring in each trial, along with a variable ETSL, the results of all physical simulations for Model 1H yielded non-linear deflections at the H-pt. as a fimction of applied load. Two sample sets of typical results for Model ID is shown below in Figure 9 a-d. Simulation results for an approximately linear seat (Seat K) and a stiffening, non-linear seat (Seat A) are shown. The Model HI response under the H-pt. was acceptably similar to experimental deflections for each of the seats. Problems occurred in simulating the deflection of the knee or rotation of the thigh. l4 H-potrlDeIledimoo WWMW 135 235 335 435 535 335 0 50 100 150 200 250 300 A 0 T +Modalll A 2 g -5 +Exp g 104% I‘m V 04 - \ . E +Modalm g \ \o g -15 g -10 \ J '2° 3: +Modclll \ l .25 40‘»— +Em x 40 .50.L_ H-DLLoadimForcam +Modellll KmaLoadingFomem) Figure9a H-poIruDeIIectIonw 185 285 385 485 585 685 g'K ' t ' +W" E -10 X\ +Exp S 4 \ v 15 A +ModaIlll L. ‘ g .20 \I E - . 5 a s -20 -30 w— 8 .30. -35 g 40« ~40 4] .50.L_ WWWFUOBM +m'" KneeLoad'noForce(N) Figure9c Figure9d Figure 9 a-d. Examples of Working Model generated results from physical simulations of Model 111 with results of Model 11 shown for comparison. Plots indicate experimental-simulation results for a linear seat (Seat G) in Figures 9a and 9b and experimental-simulation results from a non-linear seat (Seat A) in Figures 9c and 9d. Though Model IH could effectively represent H-pt. deflection during physical simulation, significant differences occurred between simulation results and experimental data for knee deflection. The differences were attributed to the dual-dependency of thigh rotation and knee deflection upon the ETSL and ETSS. If either of these parameters contained errors, the response was typically unacceptable. Since the equations presented for ETSL and ETSS could not be used directly for all types of loading as previously discussed, error due to assumptions of constant ETSS for H-pt. loading, and constant 15 ETSL for knee loading caused the rotational response of the thigh to be less than acceptable. Model IV. Coupled force and moment kinetic model Slider prohibits horizontal translation of die H-pt.F H-pt Figure 10. Coupled force and moment kinetic model Model IV was a simplified development from Models II and HI for modeling the support of an automotive seat. Mode] IV used only two Seat Factor terms to describe the structural properties of each seat. This model represented the support structure of the seat using a reaction force-under the H-pt. and a coupled moment about the H-pt. axis of the seating ' The generic support force Fjfilied on the same basic assumptions as spring constant-based support force (Kn) in the previous models (I-III). Model IV differed from Models I through 111 since the behavior of the seat between the H-pt. and the knee was described by a reaction moment M, and a rotation of the thigh between the H-pt. and knee. The generic reaction moment about the H-pt. axis M, was coupled to the support force F. This approach exchanged the two previous Seat Factors, ETSL and ETSS for a Seat Factor M based on the moment angle response. By reducing the number 16 of calculated Seat Factors, the model and corresponding Seat Factor calculation procedures were simplified. The coupled force and moment system of Model IV was extremely effective for modeling the effects of the seat on manikin posture. However, the separation ofthe moment into a force and distance increased the level of understanding of thigh interaction with a seat. Initially, Model IV was investigated using a similar physical kinetic modeling approach within WM as Models I through III. During initial attempts to simulate the response of the seat using Model IV, significant dynamic instabilities were present within the static physical modeling in WM. [It was then realized that the use of a dynamic solver to determine static results was not an ideal approach to solving the proposed problem. An important point to be made is that this is currently a static investigation into the force- deflection characteristics of automotive seats. A dynamic investigation could be pursued N, .fl in future work based upon a successful static analysis of mechanical seat characteristics. Due to these dynamic instabilities, a mathematical modeling approach using static analysis techniques was investigated. As will be shown from the results of Model IV, the mathematical modeling approach provided the best results during simulations. Mathematical modeling is also software independent, which increases applicability. For these reasons, this mathematical approach will be the basis of material discussed in the remainder of this thesis. Static Mathematical Analysis for Mathematical Simulation of Model IV using MS Excel Solver to Optimize Calculated Seat Factors Model IV represents two input forces at the H-pt. and knee, and two reactions: a coupled reaction force at the H-pt. axis and reaction moment around the H-pt. axis. To 17 efl‘ectively simulate the response of an automotive seat to two variable input forces, the reaction force and moment were calculated using the force-deflection results and physical parameters of the experimental data acquisition process. Calculation of Seat Factors was accomplished using the MS Excel Solver optimization tool and a MS Excel [3] spreadsheet designed for this analysis that will be referred to as the Seat Factor Solver (SFS). The processes for determining calculated Seat Factors related to Model IV, along with mathematical simulations using the calculated seat factor parameters, are described f), in the following sections( This procedure used modified Taylor series expansion maize?) to approximate the behavior ‘of the seat during experimental data acquisition. Modified i ‘ Taylor series equations were used to successively decrease the difieremes between I experimental data and simulation results with the addition of error correcting, higher order terms. These Taylor series expansions are modified in that the 0‘” order, linearization term is neglected, along with terms which are coefficients multiplied by both variable terms (K*[6 " * .9" ])(see equations 5 and 6). First, second, and third order modified Taylor series expansions were investigated within the Seat Factor Solver to {r/educe the RMS difl‘erenceaéween the experimental and simulation results. The ms of each modified Taylor series term increased the complexity of the kinetic model, but potentially reduced the differences between the experimental and model responses. Following the computation of the modified Taylor series seat factor matrices, several tables and plots were generated to graphically interpret the matrix-based simulated response of the kinetic model to imposed loads. Each of these steps are documented within this section. 18 Mathematical Modeling Methods The objective of this mathematical modeling procedure is to determine a set of two simulation equations (F .. and MD) with optimized parameters ("K13 ,nKfi, ’nKnuS ”KM ), which represent the behaviors ofautomotive seats based upon experimental results and can be used to simulate these behaviors. 3 Step 1. Experimental Data Input Prior to optimization and simulation, experimental data must be acquired according to the experimental data acquisition procedure documented by Hubbard and Gedraitis [l] and Hubbard and Radcliffe [6]. Data acquired previously by Hubbard and Gedraitis [1] was'used for simulation with Models 11 through IV. An additional set of data acquired by Hubbard and Radcliffe [6] (Appendix B) was added to Model IV mathematical simulations using the SFS to investigate the response of the ASPECT vs. SAE J826 manikins. 3’ An example of the experimental input region within the Seat Factor Solver spreadsheet is shown below in Table l. The first four columns (denoted by dark cells with white font) are the only experimental inputs necessary for simulation. Deflections of the H-pt. and knee relative to the position of the 2-D template [4] (dlh, dzk) were chosen to have a positive down (into seat) coordinate system for Model IV simulations. Positive relative angular deflections (Def Ang) are knee-elevating rotations. l9 Table I. Input data acquired from experimental incremental loading of the manikin is shown in white font within dark cells, along with calculated values based on experimental data shown in grey. H-PT KNEE Step 2. Static Reaction Force and Moment Calculation The static reaction forces and moments, which were based on experimental loadings and measurements, were calculated for each incremental loading step according to the equations of statics shown below. Figure 10 is repeated below as a free-body diagram for statics-based analysis and derivation of equations 3 and 4. 20 '7 . . . Slider prohibits horizontal translation of the l-I-pt.F H-pt Figure 10 repeated. Coupled force and moment kinetic model R” = F” + FK (3) M” =FK *LK *(cos(@m,c,,)) (4) The set of reaction forces and moments varied slightly with the variation in angle (9mm, ) between experimental results for each seat. These experimental reaction forces and moments (equations 3 and 4 based upon experimental, incremental input forces F H and FK) were then used as a basis of comparison to optimize the stifiiiess parameters (K’s) within the modified Taylor series expansion matrices (equations 5 and 6) using the Seat Factor Solver spreadsheet within MS Excel [3]. Step 3. Optimization of Modified Taylor Sena-Based Seat Factor Matrica The modified Taylor series-based equations used to calculate the simulation reaction force and moment based upon optimized seat factor parameters (K parameters) were as follows: F, =Z(,,Kfi, *5"+,,K,, *3") (5) 0 21 .9 *3") (6) M, =Z(,,KM. *6"+,,K 0 The K parameters were Seat Factor stifiiress coefiicients used to create the coupbd 2x2 matrices optimized in the Seat Factor Solver spreadsheet. These K parameters are the calculated Seat Factors that define the mechanical interaction between themanikinandtheseat. A f 7', ‘ . , . I ._ _) <5 MS Excel Solver was used to reduce tthMS difference een the experimental reaction forces and moments (equations 3 and 4) and simuhtion reaction forces and moments (equations 5 and 6) by optimizing the K parameters within the modified Taylor series-based, Seat Factor matrices. The linear model was optimized first, and the K coeflicient results following optimization were used as a guess for the first matrix in the second order expansion. The two matrices optimized for the second order expansion were used as a guess for the first two terms of the third order expansion. Using previous results as guesses for higher order expansions ensured proper convergence (for all optimization procedures discussed within this thesis of models higher than first order, guesses were used to ensure proper convergence to a solution or set of optimized parameters). The Seat Factor matrices optimization process could be done regardless of the order of modified Taylor series expansion (number of 2x2 seat factor matrices). For this project, models were generated using linear, second order, and third order modified Taylor series expansions. As will be discussed in upcoming sections, increasing the order of modified Taylor series expansion added terms to the simulation equations 22 (equations 5 and 6), which typically reduced the RMS difference between the calculated, and simulation static reaction force and moment. The number of sets of 2x2 matrices was dependent upon the order (n) of modified Taylor series expansion i.e. Linear (first order): One, two by two matrix with a total of form K Seat Factor parameters, Second order: Two, two by two matrices with four K parameters per matrix, totaling eight K Seat Factor parameters, etc. RMS error, or difference, was chosen for optimization procedures to minimize the difi'erences throughout the data points as a mean, and to remove the effect of signs of the differences. Optimization of the Seat Factor matrices was accomplished by calculating theRMS difl’erence betweentheexperimentalreactionforcesandmomentsandthe simulation reaction forces and moments. This number was then minimized, allowing the K coeflicients of the Seat Factor matrices to change (optimize) using MS Excel Solver. Step 4. Solving Equations to Determine Simulation Raids Based on Calculated Seat Factor Matrica \ ‘. \ l 4 i ._.—l [1:0 solve for simulated deflections (5 ") and rotations (.9") using the optimized set of Seat Factor matrices based upon equations 5 and 6, a MS Excel Solver optimization process similar to that discussed in Step 3 was used. By minimizing the RMS difference between the experimental and simulation static reaction forces and moments with a fixed set of Seat Factor matrices and a set of H-pt. deflection and thigh rotation values to be ' optimized, the simulated H-pt. deflection and thigh rotation values could be solved within a specific, RMS error range. In other words, this is the opposite procedure of Step 3 as the Seat Factor Matrices were held constant to optimize the H-pt. deflection and thigh rotation values, which solved equations 5 and 6 within a specific RMS error. From the 23 simulated H-pt. deflection and thigh rotation values, H-pt and m position could be calculated to compare to the original set of experimental data. Based on numerical analysis and graphical comparisons the models could then be validated within a known error range. Step 5. Graphical Analysis of Optimization and Simulation Raults Within the SFS analysis worksheets, several plots were generated to display the results of optimization and simulation of Seat Factor matrices. Graphical results for each seat are shown in Appendix A. In the section below, select figures shown in each table within Appendix A are explained. Plots of H-pt. Force vs. H-pt. Deflection, and H-pt. Moment vs. Relative Thigh Rotation, with both experimental and simulation paths present for comparison provided a graphical method of error analysis for each order of modified Taylor series expansion. Examples of these plots are shown for a third order simulation below in Figure 11a and 11b. Graphically, the differences between the experimental (points) and simulation (solid line) responses could be analyzed. These plots also provided a graphical representation of the non-linear (typically stifi‘ening) response of each seat. For most seats as seen in Figure 11a, the seat stifi‘ened under the H-pt with increasing applied load. When knee loads were applied following incremental H-pt. loading, which increased the total H-pt. reaction force without adding additional load at the H-pt, the H-pt. position remained stable or elevated (deflection decreased) slightly. As seen in Figure 1 lb, incremental H-pt. loading with constant nominal knee loading lead to an increase in relative rotation with constant reaction moment. Knee 24 loading with constant H-pt. loading increased the reaction moment and caused the manikin to deflect into the seat and reduce the relative angle between the H-pt. and the knee along the thigh. 300 Order Nbdel 3rd Order Nbdel Hpt Reaction Force vs. Absolute Reaction Moment vs. Relative Rotation Deflection 900 150000 J 800 ‘ 130000 c Mi(stat) :33 . ’2‘ 110000 2 ,E 90000 “—M‘a'd) V 500 z 8 400 : 70000 3 8 o Rh(stat) o 200 2 30000 100 — F(3rd) 10000 - 0 a , . _1 t 2 1 . 0 20 40 60 0000 0 2 4 6 8 10 Deflection (mm) Relative Rotation (deg) Figure 11a Figure llb Figure 11 a,b. Plots for comparison of experimental results (diamonds) to simulation results based upon third order Taylor series expansion matrices (line) (Seat A). Figure 12 is an example plot of thigh position as a fimction of H-pt and knee loading for the third order model of Seat A. This plot was created for each seat to graphically represent the deflection path of the linked H-pt and knee (analogous to a human femur) during simulated incremental loading, and to determine the point of rotation throughout H-pt. loading (similar to ETSL fi'om Models H and 111; see conclusion 7). For most seats (excluding Seat C, a suspended mesh prototype), a finite region of rotation during H-pt. loading was observable. 25 THIGH POSITION AS A FUNCTION OF H-PT AND KNEE LOADING 3RD ORDER MODEL -200 Will: (N): H-pL-Knee +20 +185-29 +26429 Position Referenced to SAE Jszs Tern plate (nun) O 100 200 300 400 500 Figure 12. Thigh position as a function of incremental H-pt. and knee loading. The graphical point of thigh rotation (zero deflection) during H—pt. loading B approximately 225m from the H-pt. along the thigh for this seat (Seat A). Three-dimensional surface plots were generated to graphically display the H-pt. reaction force and moment as a function of the H-pt. deflection and relative thigh rotation. These plots are interpolations and extrapolations of the reaction forces and moments beyond the path of experimental deflection and rotation based upon the modified Taylor series expanded seat factor matrices for each seat (see Figure 13a,b). Graphically, these plots provided individualized insight for each seat into the quality of modified Taylor series expansion model-based interpolations and extrapolations within and outside the range of experimental loading and deflection. Models that generated regions of large, variant local curvature outside the range of experimental loading and deflection had likely lost correlation to experimental results. Models with a continuous trend are better for predicting manikin position outside the experimental path. 26 In most H-pt. reaction force surface plots, a slight to dramatic stifl‘ening trend was observable within and around the range of experimental loading as the reaction force typically increased with an increasing slope as H-pt. deflection increased. This is expected due to th{stifl‘ening properties of seating fog Reaction moment response was less consistent than the reaction force surface plots, and required individual analysis for each seat. 27 RelativeThigh v r . , Rotation (699.) q. 53 FORCE SURFACE MODEL-3RD ORDER H-pt Deflection (mm) WSURFACEMGRDGIDER Nag' Relaivenigh Rm (deg) Figure 13 a,b. Third order surface models of reaction force and moment as a function of H-pt. deflection and relative thigh rotation (Seat A). 28 Step 6. Occupant Position Appron'mations Using average loads based upon a large set of human subject data gathered in dissertation work by Bush [7], sets of typical reaction forces and moments were calculated for small female, medium male, and large male body types. Using these loads as inputs into the Individual Deflection-Rotation Solver module in the Seat Factor Solver, third order Seat Factor Matrices interpolated and extrapolated mathematically to solve for the approximate H-pt. deflection and thigh rotation due to these particular imposed loads. Since the thigh surface contours of the SAE J826 and ASPECT manikins are intended to represent 50% males (body dimensions and deflected contours), results for a small female or large male based upon SAE 1826 or ASPECT manikin experimental data acquisition are only estimated. This module of the Seat Factor Solver is shown below in Table 2. Table 2. Individual Deflection-Rotation Solver for Occupant Position Approximations. HmflolvaeLawF __ . .. flfidOOCLPANTPosrnmAPI-‘RCNMAHOBS “MIOW'I—CW .V - . 2550mm FH FT F M |sr= 271 112 am 25378 [MM 43 129 537 23230 ILM 4a 135 as 30817 _ worms :2le 051mg INDIVIDUAL DEFLECTION-RO'I'ATION SOLVER SF 358 4.9 F M dzh DerArg Enu’? RSE w 456 68 25373 3581 4.9, 1. 000m * um 47.2 7.0 The terms F n and F1 within the table above were the average measured forces in Newtons under the H-pt. and thigh (235 mm forward of the H-pt.) respectably, for each of the three body types observed in experimental data acquisition by Bush [7]. The terms F and M were the true experimental reaction force and moment components in Newtons 29 based upon F H and F1, which were used as inputs into third order Model IV Seat Factor matrices to calculate the H-pt. deflection (dzh) and the relative thigh rotation (Def Ang). In the lower right comer of the table, results fiom average loads for three difl‘erent body typesarelisted. This module proved the ability of the Seat Factor matrices to interpolate and extrapolate results outside the path of experimental loading. Though only approximate, this module provided comparative insight into the deflected positions of the manikin when loading criteria indicative of transmitted loads outside the experimental loading path for three different body types were used as inputs. As expected, loads indicative of a small female deflected the manikin into the seat less than loads indicative of a medium male, which deflected the manikin into the seat less than loads indicative of a large male. Results for each set of deflections and rotations due to loads indicative of these three body types are shown in Appendix A for each of the 30 seats analyzed within this project. 30 Results and Conclusions from Model IV Model IV and corresponding mathematical methods are an effective way to determine quantitative Seat Factors that define the static interactions between seating manikins and automotive seat pans. The ability to generate equations that define the structural properties of seats as a system of multiple components based upon experimental data acquisition provides insight into the individualized behaviors of seats when subjected to varying loads. The following statements are conclusions based upon mathematical modeling analysis during the course of this project: 1. Simulation error for Model IV is significantly less than the current H-pt. positioning tolerances and less than the certainty of data acquisition. Current automotive H-pt. positioning tolerances allow error in H-pt. positioning of up to 20mm during occupant positioning and design procedures. As observed in Table 3 below, Model IV defined the structural properties of seats and predicted deflections and rotations within a mean error much less than 20mm for each of three orders of modified Taylor series expansion models tested. As shown below in Table 3, Model IV had a mean error less than that of experimental data acquisition certainty at the second and third order modified Taylor series expansion level (data acquisition has an uncertainty of approximately i2mm). The ability to model and predict the static deflected position of a manikin’s buttocks and thigh region based upon experimental data acquired fiom industry standard manikins is a significant advancement in the study of biomehanics and occupant positioning. 3] Table 3. Average H-pt. deflection, knee deflection and relative thigh rotation difference across experimental loading range for linear, second, and third order Seat Factor Matrices. VG. DIFFERENCE BETWEEN EXPDATA AND SIMULATION ALONG EXP. LOADING PATH 1. 1 .75 0.67 2.37 1 .82 0.74 2.40 2.1 1 0.88 1.1 2.35 4.07 1.74 1.75 1 .68 1.04 1 d d d ‘ I0 d ooopppopooo b «at 0b are pooopoo .3 to ..n 1 1. 1. O. 6 1 1 O O 1 1 1 1 1 O. 1 0 O. 1 0 0. 7 0.63 ' F is a statistical outlier due to convergence error in SFS. and therfore is not included in statistical calculations 2. Calculated Seat Factor equations can be used to interpolate or extrapolate results outside the path of experimental data acquisition. Using a set of H-pt. deflections and thigh rotations within a realistic range of deflection as inputs, the Seat Factor matrices can be used to interpolate or extrapolate reaction force and moment 32 results outside the path of experimental data acquisition. These data sets were displayed as three-dimensional surface plots in Figure 13a,b, and for each seat within Appendix A. An example of this interpolated and extrapolated surface with experimental path overlaid on the surface is shown below in Figure 14. The experimental path appears to be segmented since it weaves above and below the suriace generated by the Seat Factor matrices. This is due to the slightly inexact surfiwe fit of the modified Taylor series expanded, Seat Factor Matrices to the line of experimental data acquisition H-pt. Deflection (mm) D -5 Relative Thigh Rotation (deg) Figure 14. Interpolated and extrapolated H-pt. force surface Plot, shown with extrapolation and interpolation beyond the path of experimental data (Seat 0). The coupled set of equations (equations 5 and 6) could also be solved for H-pt deflection and thigh rotation results based upon any input force and moment system. 33 This was described within Step 6 of the Mathematical Modeling Methods section of this thesis. Simple static equations (equations 3 and 4) are necessary to determine input reaction force and moment based upon a system of forces at the H-pt. and another point along the thigh. For example, Bush measured forces under the buttocks region corresponding to the H-point, and at a region centered 235 mm forward of the H-pt. along the thigh [7]. These measurements can easily be transformed into an equivalent reaction force and moment set at and around the H-point. . Interpolation and Extrapolation accuracy may decrease for loading inputs that are different from the load/deflection path of experimental data acquisition. For this reason, it is important that the path of experimental data acquisition extends beyond the range of plausible simulation loading. This means that if simulation results and calculated Seat Factors rmtrices are required for a particular body type, the experimental loading protocol should be similar or preferably span below and above those loads necessary for simulation. Therefore, interpolation of small fermle or medium male results using seat factor matrices calculated from a large male experimental loading procedure should produce more accurate deflection results than extrapolation of small female experimental loading to a large male simulation. Another fundamental reason for acquiring data points beyond the range necessary for simulation is that the Taylor series expansion results can drift by loss of correlation at the extreme experimental data points at each end of the experimental loading and deflection path. 34 4. H-pt. deflection and thigh rotation can be referenced to any reference position. Data acquired for this project used the SAE J826 [4] two dimensionaL plastic, drafling template as the initial location to which all other manikin deflections were referenced i.e. relative thigh rotation. Within Model IV, H-pt. deflection and thigh rotation were relative measures that were referenced to a set of initial points. This procedure was chosen to attempt to locate the position of the manikin with respect to an approximately “weightless” buttocks-thigh contour in contact with the seating surface. Though this was an effective procedure, H-pt. deflection and thigh rotation could be referenced to any known reference position that would be convenient for analysis with modifications to the Seat Factor Solver spreadsheet. Alternative reference positions include the loaded SAE J826 manikin location as prescribed by the SAE 1826 H-pt. procedure [4], or the loaded ASPECT manikin position as establismd by the ASPECT program [5]. The two-dimensional contour of the SAE J826 2-D template [2] is shown below in Figure 15. NNPOINT TRAVEL PATH Figure 15. Two-dimensional contour of the SAE‘J826 Manikin [4]. A V4 inch, 2-D plastic template with rods extending from the H-pt. was used as a reference for data acquisition by Hubbard and Gedratis [l] and Hubbard and Radcliffe [6] within this project [4]. 35 5. Adding additional orders of modified Taylor series expansion to the set of Seat Factor matrices reduced the RMS error, but increased the set of possible solutions. Each additional term added to a Taylor series expansion firrther reduced the error in approximating the original function, but further increased the set of possible solutions. This was observed in the difference in mean RMS error and mean percent error as shown at the bottom of Table 4 below. A diminishing decrease in difference occurred with the addition of each order of Taylor series expansion to the set of Seat Factor Matrices i.e. the reduction in difference was greater between the linear and second order model than between the second and third order model. The MS Excel solver tool is iterative and uses an initial set of numbers to guess at a solution based upon user prescribed criteria. To ensure a coherent result, a reasonable solution must be used as a starting point for iteration. The larger the order of Taylor series expansion, the better the guess must be to allow for iteration to a logical solution set. In order to initially solve for the seat factor matrices, the linear matrix ~should be calculated first, and then used as a guess for the first order matrix in the second order Taylor series expansion. For each additional order of expansion (71), the n-l order expansion should be used for the first n-l matrices. In order to solve for the results of the simulation, a good data set to use as a guess for the solution is the experimental result. From this solution, the solver will iterate to a solution based upon the chosen calculated seat factor matrices. 36 Table 4. RMS and Percent Error vs. Order of Taylor Series Expansion. VI. 1 *FisastfisficflafliadatomgamminSFS, atdfltafueisruitdtxbdinsuisfidm . Order of Taylor series expansion should be minimized based upon error tolerances. As previously discussed, the addition of Taylor series terms to the model decreased the error, but added to the overall complexity of the seat factor set. Based upon the error tolerance of the application, the lowest possible order of modified Taylor series expansion should be used. This will reduce the number of possible solutions, and reduce the overall complexity of the model. As user application 37 tolerances are not curremly known for each application, suggesting a specific order of expansion as the best model would be an arbitrary statement. The best choice for expansion may vary with each seat modeled based upon the type of response and degree of non-linearity during experimental data acquisition. . The thigh lines of the manikin in all seats pivoted about a specific region of the seat during H-pt. loading with constant applied knee load. As a fundamental part of the meoretical basis of Models H and III, the link between the H-pt and the knee can rotate about a small region or point between the H—pt. and knee during H-pt. loading with constant kme loads. This is described for these models by the ETSL term (Equation 1). For Model IV, this assumption was verified as the phenomenon was visible using graphical analysis (see Figure 12) of the simulation results based upon a third order expansion. For each seat, it was apparent that there was an approximate point (region) of zero deflection or thigh rotation. This point is the center of thigh support during H-pt. loading with constant, nominal knee loads. This point likely changes position during knee loading of the manikin in the seat. Future work could determine whether there is a comfort correlation between subjective testing and the position of the thigh support. 38 CONCLUSION The objective of this research project was to develop models that described and simulated the force-deflection response of the buttocks-thigh region of industry standard seating manikins in interaction with automotive seats. This objective was met by successively developing and evaluating several kinetic models using validation techniques with 30 automotive seats/conditions. Models I-III proved to be inadequate in modeling the force-deflection/rotation response using physical modeling techniques. At the conclusion of this project, Model IV in conjunction with the Seat Factor Solver spreadsheet in MS Excel [3] proved to be the best model and method combination for describing and simulating the response of seating manikins in automotive seats. Model IV used modified Taylor series expansions to describe and simulate the force-deflection/rotation behavior of the buttocks-thigh region of manikins as observed in experimental testing. The ability of Model IV to model the force-deflection/rotation response of the manikin in an automotive seat was proven to have an error at the second and third order level of modified Taylor series expansion less than the uncertainty of experimental data acquisition. This constitutes extremely successful validation of this kinetic model. Model IV and corresponding mathematical modeling techniques could be very useful for software models of occupant positioning in automotive interiors. Model IV my also be of interest to occupant positioning software models for aeronautical, office, or heavy machinery applications. Currently, software models assume the surface of the seat to be rigid, and do not account for the properties of the foam and suspension within individual seats. Model IV and corresponding experimental data acquisition and 39 mathematical modeling methods could dramatically improve the ability of a software model to position an occupant correctly within a seat and interiorjml Future Work Though this work is a major step in developing models and methods for defining the static characteristics of automotive seat pans with respect to a progressively loaded manikin, thisproject isonlyafirst stepwithinthe fiillscopeofkineticmodeling of manikins and occupants in automotive seats. As with most initial investigations, analysis begins for the simplest scenarioEthis case, static deformation of an automotive seat due to imposed loads by a seating manikin was investigated, and subsequent seat factor matrices were developed and calculated to define the mohanical characteristics of each seat. The next step is to develop dynamic models and experimental data acquisition techniques that investigate and define dynamic characteristics of seats when dynamic A loads are applied indicative of an uneven road surface. This is an important topic, as K seats are an integral component of interaction between the human occupant and the vehicle, and ultimately the road surface traveled. Suspension characteristics of automotive seats may be of similar importance as the vehicle’s main suspension system ”"1 . . . . l in terms of the overall user’s impressron of the vehicle. _,__J One of the important advantages of the mathematical modeling technique used in conjunction with Model IV is the ability of the Seat Factor Matrices to predict force- deflection/rotation response outside the range of experimental loading. To validate this ‘ “v -. theory, testing could be conducted which used an initial data set to model and predict the if i: I I response of the seat to a manikin loaded outside the range of experimental loading. Subsequently, an experimental loading procedure that investigated the same loading as 40 checked in simulation outside of the experimental loading path could be pursued to validate the ability of the model to predict position of the buttocks-thigh region of the manikin. Another topic of great importance, is determining possible correlations between Seat Factors and subjective comfort ratings of individual automotive seats. This requires _.4 -w-r—y- comparative analysis of Seat Factor sets as compared to subjective responses to questions regarding individual seat comfort in various positions and regions relating to each seat.£lf L. correlations between subjective comfort ratings and quantitative Seat Factor sets could be established, the seating design process would be aided dramaticallijhis work my be of interest to seating manufacturers in automotive and other industries. This project investigated the interaction between the H-pt. to knee link and the automotive seat pan. Though this is an important interaction interms of occupant positioning and possibly comfort, many other interactions between the occupant/manikin I ' and automotive interior could be modeled. One proposition is to extend the current Model IV analysis techniques to additional segments of the occupant/manikin in a linked system. This model may appear as shown below in Figure 16. 4] fi« I f . l/ '7 FN O Fr F L w i 0“} MK 3 .. FP -333% MP i Mu FK . O F“ :3 .w MF Fr Figure 16. Proposed two-dimensional linked system of reaction forces and moments (dark arrows) with input loading (white arrows) at crucial joints. The ability to model the complete interaction of an occupant/manikin and an automotive seat/interior is contingent upon the ability of researchers to experimentally measure transmitted forces and specific deflected locations at crucial joints using a seating manikin.EheASPECT manikin [5] has an articulating back section with jointed regions similar to the kinetic model proposed in Figure l6; If incremental transmitted loads and _..t manikin joint positions could be measured experimentally, this type of kinetic model may be feasible for modeling, simulating, and gaining predictive postural positioning 42 knowledge for a two-dimensional occupant in an automotive interior. The ability to describe and predict the posture of a complete occupant in an automotive interior due to variable inputs is the ultimate objective for kinetic modeling of manikin-automotive interior interaction 43 APPENDICES APPENDIX A 45 APPENDIX A The manufacturers and types of seats corresponding to letters used within this document are provided as a legend in Table L-l. Tables A—l through A-30 represent the analysis of 30 seats/conditions using Model IV and the Seat Factor Solver spreadsheet. Tables A-l through A-22 represent data acquired by Hubbard and Gedraitis [1]. Tables A-23 through A-3O represent data acquired by Hubbard and Radcliffe [6]. Tables A-23 through A-26 were tested using the Prototype ASPECT manikin. Tables A-27 through A-30 were tested with the SAE J826 manikin for comparison to Tables A-23 through A-26. Within each table, the seats are analyzed mathematically and then graphically to describe their response to manikin loading. Colunms of analysis for 3rd order, 2"‘1 order, and Linear models show calculated Seat Factor matrices and the subsequent RMS difference associated with each modified Taylor series expansion. The Seat Factor matrices based upon equation 5 and 6 are labeled as Stifliless Matrices, and are transposed within the spreadsheet to ease display. Next to the transposed Seat Factor matrices, the variable ‘x’ refers to 6 " and ‘ang’ refers to .9" in equations 5 and 6. The ‘Significance Magnification on M Calculation’ was always one, or equal weighting on optimizing the F and M components according to the calculation of the RMS difference between experiment and simulation. For seats 0 through 9 or Table A-15 through A-22, the raw data from acquisition by Hubbard and Gedraitis [l] for the first 2-D template reference position was lost. As this was a predominantly a validation procedure, 13.64 mm was added to the H-pt., and 46 6.44 m added to the knee position of the unloaded manikin (step 2) to establish a reference position for modeling these seats. These values are the average differences betweenthe step 1 and step 2 positions measured by Hubbard and Gedraitis [1]. Table L-l. Seat Key for seats analyzed within the Seat Factor Solver Spreadsheet. gomammonozzr'x "IQTIITIUOW 0 2t is SEAT KEY Vectra Taurus BMW/Renee ruck 605 Unknown Aurora own and JCllNeon JCIILH own and JCl/Neon 47 scans? cabs—E mam < sow ._-< 2.3 unpesauao I “59:: .o n :95 48 33. Santa 25...: =55 Sac-coo 2 a o o N o 8 9 ON o . . 8&5 . o . Ear I o u see. . . . 8N 3mm m o 4. . n a. m m. M 8a m / 8+m.— W 22: I m a 8» 2.x: . 08, II I. noewN 3.3300 :2th 9.3.x .3 .85: 8.58: 2282 .s 85a can-om a: .802 Ice... .30: 89.: 3%. 8.88m its: at: 8.61.8 9 a o e N o 8 9 8 o . _ 8.3 o in: I fl 8N n 2.3m . 36o m 84 m. m. ) 9..st I I. 8.3 m a So our; . m I I IEII I 82 I I I I III I . BEN 8388 8:23. 2.51 .2, 5982 838m as; .9 8a“. 8.58m 5: Is: 820 EN 38: .85 EN 88. Site assoc Es. 8.688 o. o a e N o 8 3 ON 0 . . . Send o in: I fl 3%: . 8N coma M 8. m. / ) 8c M /. 3m... m 75: I ( a 8o 35:! o g . iIIItII . 89 I.I.1.i II . meow N __o_.n.__.n Sis: 9.13. .3 E88: can-om 2282 .9 85". 8.8.x a: too: 320 Ba .80: .620 En 62353 warez mam < new ._-< 23. 49 2:20“. 00.88.085.9095u02v. .20 88a 8.880 a... .9 .son. at I III I It IIIttII 02 938i 8. m 8381. a 580+ . B 8.98 - 82ml... 1 m 8+9 . m 8.2.. , 038-..: 8. w 848.1 i 8871 8.. m 02- m I I I I I I I I IIII I I II I go 38: .85 can 838.. 85. 05 a... a 8.68.. o 8 828.. 8.5. 8 0 8. o 0.9 l- . 0.8 0.8 0.8 0.0a adv 0.0. m , , - 0.8 , . . 7.0.0 ) , r/ 0.8 ) A” 0.0.. rm. Arm/u 0.2 M Z t 0.0M 1 4 1 0.0 It, .. tit . .. .t r I . .. II- . . .1. .. ,. I! lIitl. ..Il:.-t .r .- °.°FI .8383 £322. mam < sow ..-< 2...: 50 0.2 .55 8.880 a... -A§-i~ v All. V lu‘N $.53 A v .a\ N... .020 010.332 contam E52 0809 AN.“ .83. .83 8.33. :95 2.63m 0.8m.fl 89°- . gwgmfl . 89.9590 €5.20 .55: 123:8an .2088. 08. 88883 EEO almanac: 83.5w 86“. 8.283 £32... mam < sow ..-< 2...; 51 .8583 ones... mam 0 seem .N-< 2...; 23:8330 I 1.3%.. .04 3...: 52 53 .08. 8.8.8. 25...: .EE. 8.8800 up Or 0 D Q N o 00 0' ON 0 t t . t . oo+md . o . . . 2...... I :35..— e \ DON 3+w6 W O 8' .1 » \\.\ u . m m 8.. m 3...: I . L 000. 2...... o I.-I . 3cm... 8.81.8 c383. 9.3.01 .9 .0080! 00.803. 05.002 .0» 880 00.8.3. .0... .0001 Bee... .0002 See... .08. 8.53. 25...: .55. 883.00 N. or o o v N o on 3 on o L a L b I 8+m.° I p a O 8ch. I .55. o \\ 8N 3+mh W 00' .a 0 \\ m M So w . m J 00+w w Mr M d 80 .08.: I .I.II - It- 08. 30...: o - t .3. II: III. 00+mN 5.80.30 00.301 9.1.3. .9 .0080: 8.88m 05.0.0.1 .9 8.0.... coco-om .0: tees. 80.0 new .0003. .020 new .08. 0833. 250.3. .52. 88.0.0 N. 0. 0 o 0 N 0 00 3 0N 0 b h b Ir b 8+w.° t h h o .20... I w .55. e \t 8... 3+m6 m 8v .3 a \\ m n so I 0020 w Mt M 80 .20.: I a 20...: o . ttttt 80. 8.33.. 00+m N 00.888 3.1.3. .3 2.0602 Sauce: 05.022 .3 88". 8.8.01 .0: too: .020 0.0 .282 820 En .8583 808... mam m .80 N... 2...: (m)°m1938f3VSWMH'-’°!W 000 00. 8. 000 80 000 000 00. 8. 00 0 III I: .-.i...-I|--7- ,.-,-..--.I-. III I I .I. I III I I I Jfi one 0.0 0.00. . 0.. 000 .00 00 000 0.0.. 0.00. 8. 000 $0 0.0 0.0-000+ E 0.... 0.00. 000.0001? _000 0.... 0.0 000:? 00 00 _..00. 0.0. .0. . .- . . - 000 0.00 00 00.000 00 .00.- 0.00. 00000-1 000 3... 0.0 00.0. .- 0 0. 0 a... 00 _ .00- 0.00. 00.000 4 .. 0.. 0.0. 0.0 00.00-9- ON F.5- .v.anc 9...»? 10.. 000 0.00 0.0 00-? 00 000. 0 0.. .00 .00 0.0 _ . 00.- 00 ".00. 0.... 00. .00. :0 0 0.00. 000. III I-.-I- . - , ., -I-I-I - .- 000. 0 0.0 0.0 .08.. > x 38... .85 000 00.0.... 85. 0... 01.08.6000...80.80 00.... 0 00. 000 000 00. 000 000 00. 0 -. .000 32...?! .i-r: , .. - i. . , - -. 1.7.00. 0.00 . 9m. 890 new ....... o cm 0.0. 8050.0I 000 0.0 X 0.0. v 6.6 CNUIIOII / . . I. 000 00. .- -/ . ) 0.0.- w 000 m 0.0.- I 0.0. ( 0.00- . r 0.0“! 0 0 0 I I 0 CO 0.80 .. . -i .. -. . I I- I I. - :I..I- .. .I- I.-,I- :I II. I 0.0... 2.88“. 8.80.00 .05. .9 .20. 0.5. 8.6.0.0 .0... .9 8.0“. .0-.. .8383 0222 000 0 .80 .0-< 0.0.... 54 . 6838.36: EES_§§0 at z 0.. _ 5:: 263m .. Ag? , y ~ A!- - . V AEE.z. E0502 SEC 019.285 885m EOEOE EEO canine! gram 88“. “8333 ”are? Em m .3 .~-< 03.; 55 629.83 $315 $5 0 .sm .m-< 03.; u 3:3:qu I "583 .0... .35 56 .Qowcoiom 2,3. a o . 8+w.o 3+3. W 4. m m M m J gm; W -ll- ! .1 i xi. 8+5. con-8m 02301 .9 .5602 8:801 830.2 .9 88m San-om E: .28: .35.. .322 825 3.3 8:33. 55.x :55 838.3 v n e an 9 on o . 8§m.o . r o \ . 8n \\ 09%. 8am So IT. ..;\ro+u_~ 8.5.2 8.38m 9.50m .9 .5802 Sac-om 05.02: .9 8.0“. 5.831 fix Inc: 820 new Ito: .620 new 8%. Siam 5.2.1 :55 83230 v u a on 3 on o . 8E6 . . o . ,Y§.. «$3 8.3.1 .9 .85: 8.3!: 38: 820 En 233.1 .9 85". 8.83: E: 182 .020 En 32383 $322 Em 0 .8m .3. 23. 57 can 8.. 8' So 08 SN 8" on. 8. 8 o -I -- 8. fl - ..m 938:? a 08.08:? -_ 3780-? on 8.80 . 8.2»..- a. .33 . . m 8...: . 8.8» -.- m 938-? 1 933-..- - w gun? o9- w lull. .- itlillII-IIO II -v--.1-I---I.I nlllir ' 1|:- 11-3-1111' .1 illll-u 8N- 38: .85 can 868.. 85. us. at 3 83:3. - I 8&8“. so...» can 8a 8a 8. 8. on o 85 8o 8m 2:. 8m 8m 8. o 80:... 1.. - int-ll 1- Iii-{-11- -. 1.11-1- .. 5.!!!- 4 odm 121 .1. --.--.-s - .:!--.i-l-.--.. 1.1.1 0.2. v/ 3.4 can 820 new ....... V/ ode can 82 n can . o 2 II can . o 9 anal-OI Qmu 0.00 98 ) / ad... ) 0.3 m X-./ o? .m- . 2: / . 3. fl 0 o . ...L....-- - -- --------- .- .1: I.- -- c.2- 968.3“. 5.8230 85. .9 85m 85. 8.3.8 5+. .9 8a“. .9: .8333 £352. mum o Ham .3. 2.5 58 ..a. ’VNQC. . ‘ h\\\| Afifiozv EoEoE .‘ { . . .‘\ ...‘~\‘- NONNOMOHsM :83 8.30: :95 9.33. \\\o A‘ - .‘ \\\“ O .Q\ \\0 5 \\ 5 .. .2. 8.6". SEO omniuoz 80:...m 8.0“. 6 ‘90 .”" K ‘ v.96. ‘.....( .8353 “.32.... mum 0 saw .3. 2...; 59 3:... 3 n8.— .3 El 9.. ad 08 .24. 888 5 8.. 3 m8 8%.. n8 .2 . oi 0.8 .22 mm”. «12.-u. 9.58 fin 2 k 2.3.. o..- «.8. 38.- 2.. P932 «.8 am «950m 20.555.955qu .2505E a2 a c at. «.09 no at. 9:200 :8 >000 coda n «mm 25 docs. 2523 I u I m 00.0 I (km 2.: «EGO; kakQ Q2N .33.... .39 8. 339 to... 2. 322 ..9 9.. 880 8» 82 N9 .2. 8.... 1:sz moose-.30 Em Nana." New 2 v.03 03 x 0.3:- On x awNmN snn map 00v :2 .2 u .2 u .2 n. momma mom a: 2.“ mm 2 "8.3.3.8 :. 8 8.82.2.2: 8:886 2.833.858 7803.858 7.85.3858 .2 u. t 2.... >08 8a .583 .3 8.5 2.509: 8.3 8.9 no.0 m20F~::< ula—VI. ...~ :74]. .Vvl< 9.3:... "=32:qu E "59.3 o... :2: 32383 £322 95 m :3 .n-< 2.3 803 8.18! 9:93. :2: l . Air: . 11111 §§§§§€_ (N) .0303 8+3 iconic 8:3. (55¢ .9 E28: 5.881 23.83 .9 86“. 8.38: a: 33m. Inc... .822 50:: 3.3 8.33 2.3.x 6!. 5% 33 a o v a e on o. o 8+m.o o O . 8“ 3mm, 8: m o \\ o8 M 8m... . o8 6.5.2 I 3.}: o v..- o8 .3. 11.1. .. -. Rim.” £33 833. Bit: .2 E85: 8.681 2282 .3 86". can-om a: .822 520 new 18: van 65 88. 8.1.8. 8.3.x o v u / xv 6.3: I1 2%.: . 00+wd 02%! 02311 .9 Each! 8.8.01 .802 .020 En AEEV 5:00:00 ow on on or o §§§§§6 053 .9 8—0"— 531 EL 38: .320 PM (N) mg .8533 “iii mum m =3 .3. as: 80 001 St 000 RM 000 8M. 00.. 8F 3&8 I9. 895 I0! 57% IOI 0020!... on; 1.1 00.0; . 3.000 1... onion I0: 8.00? ..OI 0.0+ 300282005 . Eggnlaitcofau-ISESE (“I“) OWL 929“ 3V3 0! WU m gluon. 530503.386u8cx 0.0a 0.0M 0.9 0.0. 1 0.0 \ \ 0.0 0.0. 0.0... 0.0 w. v °.8I 0.3. (W) m §§808v§00u Tl'llfitli gouge“. £011.385ufii 0.0v 0.00 0.8 0.0a 0.0“ 0.0? 0.0.. 0.0 0.0 0.0- (“I“) ma .3fiio3 “.322 mum m .3 .m-< 03¢ . ..8u.co.§¢ .EE.S.§=ooa.: z Q. ... 59.: 2.5.3.. .83 8.1.3. 8.. 8.5 Site at .32.. 3.1.3. 8 o ,. . . . 8... 182.... . , . . 8.... gr u . r, . .. ' . 3 8189.0 - . 5:53.11... git... I 2.8.0 . \ . I mag/””6000 8......- ',/a..—,.—,—.w-.a.u—... b.901g 8881888.- ’Xfiflvf‘fo’.’ 1. 95.20.58: %%u 2.83.. o '55:," 0 1&qu 883.0889 . . < isx . com o . ‘ \\\“\ n. a, Zwmuwuwt. . 89 z”... §bwwwxxxawxufiu 8886881.”. _ 8., ‘ 0.888.- 888...- .:6n.mn..o8....o8t_m58§ 8605-3823382 8583 2322 mam m .sm .m-< n..efi 67 $03.83 232$ mum u 88m .e-< 2.3 «caisson 3 .583 04 :95 68 a8. 8.1.3. 2!... .5... 8.3.8 o c u o a. v. 9 8 o 8 3+3. ..car. In 8v . 7 8n . 2...... . a . . . 8+mo o a O O O A1 x - 3 / 34W». ) .\ o . 09 m can ( ..cazu // 09m.— w M 8.. .55.. . fl 8. _ 84m.“ 8.3.8 8.38”. 9.1.3. .3 :18... 8.88 33.89. .9 85h. 8.68.. a: .80: 820 2a .82.. :8... as. 8.151 8.1.3. an... 8.82.8 o v a o a- v. 9 8 o 8. 3.2.: in: l_ as. 8.3 .55. . _ . o \ o q 8N 4. . . 8.. m 3+3 ) X o m /1\ w NV 80 ( ini|_ . . 8.3 ( a .8 a!!! o 8.: 8m... 8.3.8 8:83. 9.1.01 .3 E98: 5.38... 05.0.2 .2. 8.0”. 5.8.3. 5.. .80: .020 can .80: .020 ten .8... 88.3. 3.13. .5... 8.8.3 o v u o N. v. 9 on o 8- . . _ 8. .c . O O w gym [I_ O 0 Soft o \ o 1 Ba 3m... m . 8v m. . ) 7 .\ ) I/ < 80 W I. 8+”? 9‘ 652' ( a Air: 0 . g— oom. N 8.3.8 8.30: 2.13. .9 .882 8.38". ass-2 .9 85.. 8.68m 8... 18.... .85 2a .88.. .86 en 8.9.53 “.322 Em . .sm .3. 2...; 69 O'NKB IOI 8mg IOI so p.80 I... as IOI an; .0. Pine: . 3.08 I... 3% IOI 8.03 IOI DINIOI 9“ En- .39 So a...“ a... 8.. Sn. .8. 8o «.8 a... 8. 98. an: 8.. 0.8 ..3. 8 can. 9": 8o am] cs 8 .3- s. .3 2.» Eu a... 8 0.8. 2": .9 3.. a... on .35. .313 a: 3: .2. 8 as... $3 8» 3. 2. on Hole. .23 .8 a... 2. 8 98. to? no. as a... o u 8. 3.... o oo 2. 88.. > x 8... 8.13m § l .7..- O 1 I (“I“) am; 928i" 3V8 W mu own-0d ‘\ g § 38.28.59.” piggaitsfisuuasasacg .8563 8%.... mum . sum .3. 2...... 70 EErzv €885 .620 0104082 885w 85m 3%. 8.16m 2 . _ 2959.531 829.53 WEE/x mum m :3 .o-< as: 71 "53:26 3 “583 ,3 5...: .8333 £322 mum o sow s-< oz: Agvgmim o o v N o N. v. o. - o O 4 - 8... / 8. m. / 8o w 1. 08 $5.: I— ..u o . o8. E: 36 a 5.8.3 8.181 3.3:. .9 :25: 538m .538... .9 8.3. 8.88... a: .80: i8... 38.... :8: 8.3 8.38 3.3.3. a o v u o a. v. o. . P . . . v 0 N . . 8~ / 81 m / 8c M Gail /. Bo Air: . o8 III I I! . 3m.“ 8.3.8 8.681 9.30m. .9 885: coca-om 05.8.2 .2. 8.0“. :0. .081 sf 38.... .85 EN .80: .86 En .§.S§m§.§om a o c N o a. v. v A; / / / 96.2 I 2%.: o 8+m~ 5‘81 Deg-Ix .-> E982 ERICK .000! .520 En §§§§§é 8:00:00 05.33 .9 88m 838: 1.. 1615203 (N) "‘05 .8383 ”abuse. mum o sow .3. oz: 73 nvnfioc IO. Sag I... bowglt and; . Dug IOI 1 . . . a; I01. \\ § Ex 7"— O 53 8.8» I? \ ail? . -..\.-- . . 83.... 1V DINIOI § (W)m1m3vsqmam 8.. o8. 18: ..20 can Eggulntusgu-ISZSE 8 o g 8o 8... 8w Sn .8 8. o 0.8 .85 - - -11 I. III}. I II 28 I. owe. can :20 2.8 ....... / a 2 ..20 2.8 ....... If . :8 o 3 ..2o2aIII // 98 m 8205' J/ - 28 M 5+ / o.9 ) fiulol / x , o8 ) / o 8 m Affiflh m 98 I , .. ad. I [1 < ..C V 0.0 2:. 2:. 2. 88.. 9.. 8.0.. inggguoocx 820%.388ufii .8533 v.38... mum o 2.3 .3. 2...... 74 _ 38. 8.3.3. 3...... 8:88 at m m a. .L 5...» 9.63. .E... 5.38 n... u m m. _ , a can. E . . .. . .. 83- .I... \b 33... (mm... \\\\ \\\I. .22. .55... 8.3.2 o .z. 8.... ..k‘sskRuNu» gem.— I 98w \\\w~§\\\\ 8m. 88 820 05.38: 88.8 .58: 320 85.38: 8&3 8.0.. 32.3.83 2.22:. mum o .8 .3. 0...; 75 32933858... mum : Em . w-< 03:... neg-.325 E .523 .u... :95 76 .§.co=-sm9§ox €58.38. 8 o v N o n. 8 9 8 o A; . c . 8.». g m. / 8.. .m. 2...... l /. o8 3.x: . 82 1’1I I. 3+m u no 8.18: 9.1.0! .9 £2.62 Sacco”. 05.0.2 .9 88m .5331 if .86: .08... two: .88: .§.8.!o¢ 8.8.8. o o ' u o N. AI/ / 3...! / 2...... . 5'8”. 9.11m .3 «co—:02 5.3001 Inc! .090 “EN A§.§suzo>§ox a o v N o / m 8+m..m l. 6.9: II 3%.: o .- I . - - 88.... 8888 5.16m 3.33. .9 2.8.0: 8:83. 05.022 .2. 88". 8.83: a... .805 SEC En .805 SEC En 77 .8333 .33... 95 : .sm .3 2.3 ovu.§ IOI 8a.: I.» B v.80 I... angel? . an" 3; oué; . 8.08 no: g§+ 3.02 IO: QINIQI 7--- O 53 390139.005 ggggnttgsu-flgzg § (ununmxszeravsmmum § § o 85 80 8» 8Q 08 SN 8.. o .zv 88“. 5.8.3 .91.: 83m at 339.63 m§a=< mum a sum .m-< 29¢ 78 . 3.382%: Asevcofioio at m 0.. . 5E. 25oz m , h . p . . V \ . A, . \ RE 5 V. \\\\h N x x \ \\\ «fives \\\\. \A\\\\\\\\\ ' \\\\\\s...~.\\\\\“\ {Rx§\§ §%\ ‘ a n h - Ehrzv .560! 320 05.38! 835m «:28! SEC 054082 8-2.5 080.... 80383 £322 mum : sow .w-< as: 79 829103 22.3 mum _ .sm .o-< 03.: n «coasauiu I .595 .04 3.: 80 Bomé .I. 8m: 8.13 93.1 .9 2982 8.68m :8: .85 v5 . o . 8n 9 . .1 34m m m \\ 8v m / ) 98 W 3..”— W M €32 I J ( 4 8o 33:! o 98— I... gm.“ 5.8800 8.301 9.311 .9 E28: 8.603; 830-2 .9 88“. c263: sf .80: i8... .80: .09.: 303 c9181 93.: AEEV 3:850 o o v u c or on on or o P .- > v o . 8w 8. m. / 8.. w t 8.5... M fail I m. d o8 A.‘Iz . IIIII .IIIfigF IIIIiI- . III--. . gm.“ 5:00.30 c3501 9.33. .9 E88! censom 05.8: .9 020“. 3.881 E: .80: 520 v5. .03: .020 new $035.18: 933m 0 o v w o . . 8+m.o . c . . 8a 3+m.m W 8v m. ... § w M So 98 053.988ucofilozxf 383.5205 .8333 £324.. mum _ “8m .o-< 03.; 81 8o 3v 8' 0mm Sn 0mm ova.~.8 IQ. Rang IOI 5 pg I as I on; -0. 8.0: . Qua Io. anion I0! unbow I0: (ula) om; czar avs 0: ma unused O‘NIOI 8? S“. :8: .86 am» igantucofisuuicszoafig 8 a 9: So So g 8n 8.... 8. 0 I2: 98 0.8 92 97 98 35 EN ....... / 0.0.. 85 EN ....... // 0.0N .88 RI o.» .86 SI 0.8 . / I v ed W I” 99 .BI / I fiulol 0.0. o9. XI. 0.9. odN. 288“. §o§x9gu3§ 283“. 8.62.8 at .3 85“. .9: .8363 “iii mum _ 98m .3. 03:. 82 . 3.38.8.3. .. Q. _ €32.58. l SE 3. E5 8.3.8 at a n V 1 {‘4~M\.p , Vlv‘ s ‘. ’1’. \\ b .0.§0W\\\\\\\v~ ’ ‘ .....‘l \N\\\\\§\ Q » l’ObObQ§NN\M\QvMRA\QQ\ ’ \\\\\\\Mfl . . ém ‘ ....»w.u~wu§... .9“ 2. 8.8 E22. ...-.5: . . .~\\\\ V .\§b\§ 820 05.1.0: 835w 8.0... EEO Dru-.000: coat—6 E85: aegis; £3.22. mum _ .8m .3. as: 83 "cognate E "52.3 o... ......E $39.83 ”gee mmm . :3 57¢. 2...; as. 835. 9.3.: §§§§§5 05.022 .3 090“. San-om «3.. .032 .525 (N) ”:1 §§§§§5 2:68... .9 83... Sac-om a... (N) 0010:! .80: .85 EN 38. 8.33. 23.x .5... 8.3.8 n c u o 8 8 2 o _ 8m... . o . - 8... «cm... m g m. / ) 8» w /. 8a.. m M 75.2! ( J o8 ....zz . . I ill! ..8 l . 8m. u 5.3.8 5.5... 3.1.: .9 295: 8.68: .5322 .9 8.8. 8.3.x a: 18: .85 En 38: .85 En .8383 “.32... mum 28m .o.-< as: 85 80 8v 8' can 80 0mm 8.4. cm— Sp on o 3 (ununmuza'avsoamamood ,8_ ' 0 é 300582005 igniuibcsgm-Scozng 86 gonzo“. 2.8.0“. 833.383.; 8.6.3:...9853... .8383 22:5 Em . sum .o.-< 22; _ ..oou.8.§om a, . .3...» 9.33m , L . v . \ CanasJ 0.8“»: \ \v \A‘I“\M§o. \\\§Q~0-\\\\ \ H 08.4. w- 96- gm 3 A n. \ K 88685 o8o " -~\\\\.\V~\s\s\\\ .5358: 2.88“. 88 .9? x..\ g. gr " \\\\\\\‘N~s~\ 8838.... 88. H: --:\ .030 010.382 835w 8.0m .020 $3.38: 885m E25: 87 2.2.9.53 22.22 mum H :3 ... .-< 0...: 3.2.53 .3322 95 V. .3 .. .-< 2...; "52.... ...... .3... 88 88. 8.33. 2...... ....E. 8.3.2... 0 0 c a o 9 8 on o. o . . 8+0... . o o 0 0 0 A can I a ......E . . . 8... . .1 vO+m m M o 0 8w m . m 8.. m maria ......5' . m «. 8» SIZE 0 lwl . gr III! . mo+wd 5.80000 02301 02301 .9 E0562 5:000! 05.002 .9 00.0... 8.600”. E: .000: .09.... .000: .02.... 38. 8.13 9.0.3. ...-5 8.3.6 o o v N 0 av on on or 0 . . 8+md I. c . . . can... a ...-Em . , 8m 3%.... m 8v m ) 9 \I 80 IN\ 8 mo+m.— W ( J 80 {lsli Ilr; ,. . 8.: 09m.“ 5.30.00 8.1.2.. 028.01 .3 E080: 8.6001 05.005 .3 020". 00.6001 5: .03! .020 “EN .082 .020 ucn 68. 8.08". 2.8.0... ...»... 8.60.3 o a v n o 9 on 8 e o . 8&6 o 0 o o 0 6.8“. II a .55. o 1 8“ 3+m.m W 8v W / , a. w A. 3.3m; m 80 6.9: I ( 4. 30!: o . VII! I- in: ‘1 08.. ill,[l . 39m u 5.32.00 5:83. 02.0.01 .9 E0822 8.90: 05.89. .9 020“. 8.30m 51 .0005 .020 En .0022 .020 Pm 3.2.83 22.2.... Em 0. .sm .. .-< 2...; § .3. O 8w- (“I“) ml czar 3V8 03 me: mad 03 t: .. 7.03. R .33. £3 8» 8v 8' 8» So 8m 8w 8. 8. 8 o 80 can .3 B0 3.8- 39 80 ....m ....o 928:? 8 .08.- 9.09 838.? g 0.50 .....0 hag IOI 8 32. $9 20 0.8 o... 8.80 .- . 8 0.8.. «.09 820+ 30 0.8 20.. _ 8.37? 3 .0.” .o. :4 8.0:. .. 0: Eu 2... a 0.8. an: 8.80.... 08 a...“ .20 348:? 8 o. E. .0. .3 3.8. no: (ON N... . 0.0 IO! . . a.“ a“ 08. 33 p 0.9 o... «.8. 0.03. o o.o ad 88.. > x .08: .85 can 51—. ......flmH 80803903352008.2305...» 8a 8a 8a 8. 8. 8 a 3. .005..- n '00:-..11 n 2: 320.5 ....... o.» 8208“ ....... / 153' . \qu w 0203'! 5+ on. ) fiulol 0.07 m . 0.9. 0.0m. A20 020» 8.60000 005. .9 020". 005. g 8.0“. 8.60.30 at! 83“. a... 08:95.3 «.322 mum 0. .8m .. .-< 29¢ 90 ...-:2. E83: .85 22.8: 835m .58: . 38. 8.30m .65. 8.62.8 a... m. 0.. . 59.: 2.0.3. ' . . .. 0.0.03 ‘ '§ §\\\\\\. . - \\\~ .Q“. \\\ .020 019.0002 80:30 020... .8383 “.322 mum v. 3% .. .-< 2...... 91 «coca—.53 nah—05‘ wmm .— aaom .N . -< 030,—. "8....330 : 2.3%. .o... :95 92 308 00.081 9.131 o 0 v N 0 $0.5 00:00:00 8 00 0' on o 860 00...... I 0 3...... o N 8“ 3&0. m o p 80 w. ) 80 m ... m .03.: I 3+w ( a 80 ...-...: o 08. r 30." 00.80000 0000.01 02.0.01 .9 0.00.0! 00.6001 05.00.? .9 00.0.. 00.6001 .0... .0002 .000... .000! '00... as. 8.0.3. 9.000 .85. 8.00.00 0 0 v u o u. 00 00 00 00 o L . . 8...... is. I. . . o :85". 0 \ 8N 3m... m \ 8v m . ...... 0 Gail v ( a 80 30!: 0 25.. F 8.0”“ 838600 00.301 9.0.01 .9 .0900: 00.6001 05.00.? .9 020.... 0000001 3.. .0003. .005 00a .0002 .005 00m 68. 8.06.. 9.000 00... 8.088 0 0 v a o u. 00 00 9 on 0 09m 0 0.0.1 I + . 0 30.01 o 000 3+0... m \\ 8.. m. / > \w 8. w / 3.0.. W H a 0.02:: ( a 8.. 30.0: o 89 09w.“ 00.80000 0000.01 9.1.01 .9 .0800: 00.8001 05.009. .9 00.0... 0000001 5.. .0002 .005 0.0 .000: .005 0.0 022.03 090.2 m.m .. .sm .2-.. 03.; 93 80 000 8' can Sn 00m SN 09 80 or. 0 (“I“) m1 W 3V3 01MB mend 8?. 00F 08. 3002 .020 010 23492058188210.0003? 0.0m 0.00 0.00 00' 1 0.8 XXX 0.8 ) . WM 0.2 m 4 0.0 0.0? .2. 00.01 8.03.8 at .9 851 1+. .8000; 0.0205 Em .. .sm .2... 2...; 94 .020 010.3002 003.5 .0050: a... .0855. -\\.&&r. . .090 019.0002 000.50 09.01 .08. 8.881 . 59.: 02.0.01 .8383 0322 $5 .. sow .~_-< as: 95 «000030.00 E “5%... .0... .82. 82.9.83 ......22 Em .2 8% .m . -< 2...; so... 8.881 9.1.0.. o 0 0 a 00.151 9:1 .9 E0032 00.8001 05.29. .9 88. 8.80.0 0.... .000: .000... .0002 .000... .80. 8.881 9.8.01 ...-5 0280000 0 o v n o 00 on on or 0 . 00.0.0 I o . . . .231 fi w ...-.01 0 8n 00.5.0 m 8? m. / ) 8° ., 7 m 9 202' 10 mo+m . m J 000 53!: 0 gr 11!}-1‘-.. . . mafia 00.80000 000081 900.01 .9 0.00.0: 00.8001 05.002 .9 00.01 00.8001 E: .0005. .005 00a .0005. .005 00a .08. 8.0.01 9.0.01 .85. 8.80000 0 0 v w 0 00 on on or 0 . . 8&0 o . . . 00w 00.000 W 80 m M a. w J 8.9m... m M .08.: l ( J 80 ...-...: o 08. 8m: 8.52.8 00.881 900.01 .9 0.00.03. 00.8001 05.009. .9 00.01 00.8001 .0... .000: .005 En .0005. .005 En .8383 032... $5 2 8% .m .-< 2...: A2. 00.0.... 02800000001900.0135. 28.01 8.8080 a... .3 88. 8+. 80 000 80 08 80 08 80 8. 8. 8 0 8. ... m 8 8 _ . . -.- o .. W . I. 8.- m 00. 0 lg: 8”: .08: 80.0 000 Eggninttfigtlfifioafifi 0 8. 0 0.8 0.0? 0.0. 0.00 0.0. 0.8 0.0 w ”MM m \\ v0.0 5+0 °.mF 0.0. ) /. . ) m I 00. m 0.0.. ( . 0.0 . 0.0.. 0.0 0.8- 0.0. 820.83 m.008... 000 2 .30 .2-.. as: 98 AEE.zV €060: .095 own-.03: 80:5 20:5: “Evcofiio at 9. Q. . 3%. 8.30m .55 02.0.0: I. x mm» . ‘ -.p‘r ‘ / 9“\\\\\\V\- . aw,........-...“%§. 7, vi 0. kfia \H‘fi ’ .86 9515: 835w 85”. .8333 “3%.:- mmm .2 sum .m_-< 03.: .8583 m§==< mum z sow .:-< 03: fl «co-aghu‘cu z 25:3 .0... $5 100 3.3883281... o. o o v u o 8+3 , . . o O . . 8a 3%.... m g m. / mom. M N1 8» w ..cail 1. m- a 8a 3.2: . 08. .||. 8&0 8823 5.8.3.. 021.01 .3 «clue! coaucom 05.00.? .9 090“. 8.8001 sf 08.; inc... .08: IS... .83 8.13 9.10: o. o a V a . . . . . O . 8m 8' m. / So .m- eézl / 8o 3...: . -I I; - 1- .82 V‘r.‘-}- \YI yo, .8omd c9303 8:06: 25301.9 E0502 5.800: 05.00.? .9 020.... c9600! .01 18s. .85 .5 08: .85 EN 303 8.13m 98.3.. o. o a v u o . 8m... . o . . 8... 3&3 W 8. m. 8.3% So w 1. gill ( 98 3.2: o w ‘11 l 11 3.08. i I 8a.” 8%. 5 8.1.3. 9.1.0: .2. E15! 5600! 05.8.2 .9 020m c3801 sf .82.. .85 2n :8: .85 E» .8583 £52.. Em z 3% ..._-< 03.; 101 8..- Sv 8.. 8n _ 8n 8" 8a 8. 8. 8 o 8. 8. ova-.8... m 838.? a 8.80:1 .. 8 8.80 . . . u 8.5.... . o m 3.30. o: M 8.03 - 8. 8.8.” a? m ail... I. 8.8. no. 8.. m Oanol O 8.. m, 8.... ..8: .86 oz." Sgégntufigu-Iigg 8 o 8. 8.. 8.... 8.. non 8a 8. o 2... QB .02.: - I 98 98 g EN ....... 9.9 o 8 m .86 Ball :8 W o o. ) ill 8“ ) . o... .m- . ed. M - 2..- o... 98. ad.- 2. 88“. 2. 88.. Sic-......ggusé Sfioioutésauai .0298? ...-.5... mum 2 .sm ... . -< 2...: 102 .55. .8an .690 Btu-loo! 803m E832 3%. 8..-8m 9.58.628 .9: m .- Q. .1 .2... 9.3.: m ... w. w . . 1 2.8.0 . 8......- 82.83 8988.0 2.8.3 231%.“ .020 05.382 Slim 00.0... .8333 mass... Em 2 .sm .3... 2.... 103 vv 0 v0.2 "sacs-.30 3 "59.3 o... .....E Boga—.53 max—0.2. mum 0 .3m .m T< 030L- 88. 8..-8m 3.3.... .5... 8.0.80 0 0 v N 8 9 on o b 0 1r 8.”.0 aim I F p o ...-...... . ¢ 8.. . .V\ a 09m 0 80 m 2.351.. .v 3.5 . ( 80 30...: 0 coo. new.“ go 00.33. 023.01 .3 .0080! 8300a 05.00.? .3 00.0“. 8:000”. a... 1002 80:3 .0002 I2... .8... 8..-am 2...... .se. .23.... o o v N 00 00 cu o . _ . 8m... 6...... I L . o ......E o \ 8.. 09m.» m 80 m. 09w. W .2521 / ( 8.. ...-...: o r- 2!. 8... r its-L. 8..-Nd c930 8.38m 9.13m .3 .0080: 8.303. 05.022 .9 00.0.. 3.8001 .0: .000: .020 new .000: .020 05 .08. 5.33. 3.1.8. 2...... 882.00 0 0 v N o 8 3 on o . r . 8+0... . a .05.... ..l .55. . L\\ 8.. 3+0 0 m 80 m. .. ...- 8+”. ...-...... . H 8+0.“ .8. 8.1.3. 8.60000 023.0". .9 E050! 5.8001 05.022 .9 00.0". 8:000: E... .000: .020 0.0 .032 .020 0.0 3.9.63 ”.32... mum 0 .3 .m .-< 0...... 105 8» on. 8.. 8n 8.. 8a 88 8. 8. 8 o 8. 08.80 no. 8. m 88.80.? 538.? _ 8 8.98 . . -.\\\- 8.8.. no. h o m 8.3.. r? m 8...: . 8. w 8.08 Io- -. 8+8 .... 1 8.8. I? 8.. m DIN I0! LO! 3... .m- 88. 1.8... 805 can 888.. 85. .8- .a... .o 8.8%.. 0 8 838.. so...» 88 .8 8. 8. S o 8. So 8... 8c So 88 8. o 8.8 8...... . .8 - - - - o 8 c... A - . adv 069 .320 EN ....... k 0.00 0.0. - - 0.00 . a.» m 88.0 ual... . AIM/1 8.8 m 0.0 cut I - I z / odw 06: W) _ I I 05hr MI . ......- ( I . a... ( . o...- - o...- OdN. 4 0.0 8.8. 8.0- .z. 88.. 283.. 8:00:00 005. .3 88.. 005. 8.833 8+. .3 88.. ...... .8383 8.8.8... mum 0 8% .m .-< 0...... 106 €252. EoEo! 820 Burrito: scram EoEoE .820 Dani—.62 008:5 8.6.... .8... 8..-8m .3...» 2.0.8. .8383 WEE... mum 0 wow .m .-< 2...; 107 .8383 .322... mum a .8m .2-.. as: 10,? "583 0.. :25 108 so... 8..-.3. 9.13m 2...... 8.85.. o o v a o u. 8 9 8 o b n b 8+m-o A.=au l . ° 30ch 0 \ 8N 3.2m.» m \ St W. 2.3.... I. 8+m . ( 4 8o 23...... o 08. . 8&8 5.8.3 8..-...... 9.11m .2. .228... 8.303. 228.... .... 8..... 888.. ...... .000: .02.: .03: .02.: 2...... 8..-um 9.13m .52.. 8.62.8 o o v u o u- 8 3 8 o . . . 8m... 2.8... I _ o ...-...”. . \ 8... 3+3. m \ 8.. m. :55: o L .I 08' - I -- 5L 8&8 8.8800 8.33. 028.01 .... .2028: 3.8001 05.83. .0.. 02...". 5.800”. a... .08: .020 ....N .032 .020 new . .8... 8..-.3 25-3. 2...... 8.6800 o o .. a o 8. 8 9 8 o 8m. 0 in... I. . o w 23...... o \\ 88 . o 38.. m 8.. m. w \\ m. 8.. w m / 8+w w W M 2..."... ..I a 8.. 28...... . 8... 18.1.31 888 8.688 933. .... .88.... 3.8.x .532... .... 8..... 8.68m ...: .082 .....0 .5 18: .85 ...» 8.9.33 ......5. 8m .. s8 .2-.. 0...... 109 00m. 09. 00* 000000 on“ con 03 00. On A20 020". 00300000299098.08v. 8.880 ...... .... 8..... ...... 8. .. m 83.8 .... a 838 I? 8.88 no. .. . 8 8.08 . .»\-\\\\ 8.2.7.. . o .o. 8.3. ..- m 883 . 8. 8.98 no. .. 3.3" :0: \x I. 8.8. I? 8.. m 9N ..OI 0 8.. w r 8“! 38.. .85 can 8.83 85. ...... ...... ... 88...... - ... 8.8.... 8..: 8o 88 8. 8. 8 o 8. .8 .8 8. 8m 88 8. 98 8.. - - - - 98 I25 - I I .. l 0.3 ... I I I 0.3 into 0:“ ....... 0 On .320 flcn ....... 1|? 0.00 . OS .85 ...» III / . .85 ...... .II 10,, Q8 .8 f 0.8 80+ ll... 06. ) 50+ 14V 9m. w; . . v 0.0 M X 0.02 ( o... v 0.09. T 1 0.0 0.8. ...... .2. 8..... 8.9.33 ......22 W8 .. .8 8 . -< 0...... 110 :25. 3. :55 8.850 5.: AEE.ZV E0502 .520 casino: 885$ «SE05 80.0- 80700.20 82ng 000N609 I MC (50' _ 683888... 8508808991 .. a... . 8..: 2.00m . x. ’\> I \E? x \\\\\\\\s. .620 0104282 855m 090“. 62383 2324‘ mum m =8 .2-< as: 111 38!. oi: oi 088 no. . z m #0855 08550 ...n. 8595. 8.33 «10.5 ...—...! 7.3... MIMI: 38 n. «.8- 3 8 N... s. m 088.. ...v. 90 9:8 0... x 2 u b.3350! 000550 Sow... 5.8. 80 Bow: 880 m.» «comm. 8.80 «n. 8&3 8.8 8» «9w»... 0.8. 08 8808 «.8. 8m. 8&3 3.. ~00 8.3.“ $8 o5 8&9. gnu. own 8.50.0 «Bo. 8. 8+wod o o «10:- 83.. .05... ENDED DZN 0.: No 9.9. 0.9 a... at. n. n. 2 u ..n Q; :2 VB 8.“. «<9.- 3 0.8 ... .8 «.0. as. to fin >000 .2 ... Sosa. .8. 8. 5.." 8... 8. «9 .5 9 . ~88 8.8 x 808 n8 8. 8w 2.2 2 u. 88a 39 a: In um . $000.50! 000055 .E Ik >00m man 9.» mzo_..<:_x0¢&< 20.501 .255qu 8v 35 0:... 8&8. m «9mm. N50: 03 93: «3 «guns .3 S ta 0:8 08 gun.“ at” on» n KKK men 5+wns .80" 05 «80“ who gown: was. '8 :5“. 08 5+wod 5m: «3 RR. our. .3...» 88. o: 88. 5: «8m... 88 .8 .08. .8 .93.... 03.. B8 .08. 08 «938 atom n.« 28. In 8.3... o o o o «10.:- :..cl :5: ...-...: 30...: ggiigi Eéélmu!l M RR“ 2.: .500: fifiko ka. 3.0 .0 a. u u 2 .808 02.. , .8883 2 8 83:28: 8:883 .583 .3 82.. 2 E31 2.83 85. 2 3.8.1 .8. 2: 25 E038 ». 3.388 288 ...-t! E36 film. 300: .8353 mi? mum 0 .sm 8 _ -< oz... 112 38. 8..-8.. 2!... .25. 8823 o o c u o 8 9 8 o L . p 8&6 2.3:. .II . L o 3...... . 8n 3+m.o m \\ 8' W. .55.. . o8. . 8%“ 88.8“. 8..-8m .2331 .2. .5an coco-om 33.022 .2. 88.. 8.8-3. .1: .80: I83 .03: .35 so... 8.5.... 3...... .52. 832.3 o o v u o 8 9 cu o . . . 8+w.o in... I . . o 35...: o \1 8n 3...... w 8. m / m \ ) m. [S m + . w H1 65.3.! /. 8 m . m a 8.. 3.2: o .ltrislinilitlt o8. T .. .. 3+...“ 8.8.50 8.33.. 25-31 .2. .880: 8.8.3. 83.82 .2. 8.0“. 8.881 a... .03: .620 ucm .032 820 ucu 38. 8.39.. 21.8. ....E. 8.885 o o v N o 8 ow ON 0 p \F b 8*“.0 h b earl a 3:5. . \\ 8... sou o m 8' W 0 \\ w 3 is! 2.. m /~ ( « So 3.; 0 fl _ 3.3 {Ii 1:. 82 8.888 3:53. 2.301 .2. £23: 8.8.3. 05.92 .2. 88m coco-om a: .802 820 En .03: .85 En 829%? mi? mum 0 .sm .. .-< 2...: 113 :5 8.8“. 5.6180 85. .3 88m 85. 0.8. Az. 88“. 8.883 3.... .2. 85”. .9: 8... 8c 8' 93 8a SN 8a 8. 8. 8 o 8. d 3.39.? 8. m 838.9. a 5.80.... . ‘8 g o \\ m 2.2.»... v m 363 -o: M 233 . . 889.? m «39...? I. 232.9. . m o..~IOI M 8.. w m '4' 8“- ..8: .85 oz.” £33§x9u$13835m.882822» 8 o 8. a...» 0.8 3. Q8 0.8 O 0.9 0 98 W 98 m . . m / oor W: OON WI 4 v 0.0 M 0.9 .W 1 0.2. + 9o 2:. 3.9%? ”2.52 mum 0 ”am s _-< 23. 114 i I |‘ t I ll...! u :55 Satin 3.: 10.00.. 880m 8003 80.5 93.38: 085m E25! 2.5.5 «5.5! 0.. $90 012082 885w 88... :83 8..-am :95 .253. .8583 252< mum o 38 .:-< 2%; 115 .8583 “iii mum a sum .2-< as: 3.0 3,0— n§§§8 : «583 .04 :95 116 3.3 8..-8m .21.: or o o v N o Sowd // EH 3+wd m 8...”. W 8§WN :55 5:03.00 EDI 1.. 3.x: . \\ N d (N) ”10:1 §§§§§5 F // 8.8.!” 83.x 9.3.x .9 .85.: 88...: .5322 AS .2... 8.88m 1: .802 '8: Enos. .Ioc: 8.3 5!... .21.: EE. fiuoio 2 . o v u o 8 3 8 o p . . - 8+”.O ABNVK 1 O 3.5. . 8m 8m. (N) 3310:! \\ m 6:8: 1 a Air: . . o8. 1.1111 . Bow a issue 5!... 93.: ..> 8.5.: 8.6.... .5322 .9 85. 5:8... u... .822 ..Eo 2m 18: .85 new 3.3 8.5.: 25...: :55 8.6230 2 o o v N o 8 3 8 o 8.3 earl A. m 3.5. . . 8... 3+... m 8. m. .u. 1N“ 8. ) 3+3 w m ( a 8. . 1 82 8+“ N 83.3 8.3.x 93.x 9 .85.: 8.....m .5322 .9 8a”. 5.8... a: ...8: .85 E. ....s. .86 En .8353 £322 mum m ...m ....< .3; 117 118 8. 8. 8. 08 08 80 SN 02 8. 8 0 8. 8. m 93.81 80.8.1 a $381 , V 0 8 x \‘ a 1 _ m 82...... v 0 m 848 .7 m 8..: . 0.. 8.08 I? w 3+8 1 .... 8.8. I... 8.. m OaNIOI O 0... m 80- .80: 8.20 can Eggnlaiusgualszfiz» 8 0 8. 08 08 8c 08 08 8. 0 0.8 8...... - - - 0.8 0.8 71?. 0.8 0.9 .85 05 ....... 4/ 0.9 0 8 W .85 20 1 0.8 W 0 8 ) .81 0.0. > 0.0. w m . 1!. m- . 0.0.. . . 4 q . . 0 0 0.00. 0.0.. 8838x9833: 580.081.98.003... 0809.83 8222 mum a 080 .0 _-< 20.: gpgl gum—gr. ggwl §§u Earzv 59:02 820 019182 885m Elan! Smél 28782“. Smwgp D 33p . :83 8030: 02.: 900.0: .620 9.3.382 oootam 83m .8383 £322 mum a .80 .0 _-< 0.000 119 .8583 $322. mum m .30 .0_-< as: 3.0 10.9 n 0203030 3 .5003 .0... 00.5 120 as. 8.30m 3.33m :55 8.3.0.0 0 v a 0 8 8 8 0. 0 _ . 8&0 a.” . . 0 03.5. . q 80 3&0 .\ m. o \\ 8. w. . m ... 8.. m J 84” w ( d g :50: .l . 08. :0}: o 8¢wd c2802 8.8.9.. 3.1.8. .3 .5502 Sac-om 03.82 .9 8.0“. 8.63m xi 1.0: .80... 18: .85 .80: 80.8.1.8 2.03m .55 8.3.8 0 v a 0 9 8 8 0. 0 . . 8m... earl . _ 0 Sara . 8w .2 w / 8+w.r M M g ....-. . I 89 A 3r: 0 8+WN 830m 3.13m .3 .0200: 8.68.... 05.8.2 .3 8.0“. 538m .0... .80: .86 05 .80: .85 000 38. 8..-um 90-3. 05.. 80088 0 v u 0 9 8 8 0. 0 . . 8+m.0 6.0... .l . . 0 ...-Em . 08 w 3m... \ ... 80 ) /. 8+3 m N m ( g sail a 20...: o . 80. 8m. 0 0383. 0233.. .9 .0882 coco-om 05.82 .3 88“. 3:83. sf .80: 320 En .08: .85 En 329.33 23.05. 0% m .80 .0 .-< 2.3 121 08 09 8' 03 08 08 88 8. 8. 8 0 00. 9a .080. ..03 .80 Fm 0.0 m 08 ...mm. “.03 8. 80 «.... 0.0 08.80.? m B. ..v... 3.3 Bug lo: a .... .... ...... ......+ .. 80 mm 0.... 8.80 . \\-x....\\\\\W\.\\\\fi a —°.80 68' 8g + 1 u d u 0 1. .1 \H1 . . - l\ .1... 11‘ o 0.0 2.0 ....mm 0.0 8+9 -.. \ .. . o 8 :8. w. .3 8.0.0 . . m 3' .0 F” .... .U .\l\|\ul| 8| 8 _ ..8. 0.03 8.08 no. \ . w 0.0 1.8 0.0 88.... . . I. ON 0.00. .33 0N8? 1.! 8w. m 80 ..8 0.0 8 0.3. 0.03 0.81: o vow oww 0.0 or..- ) 8 .98. 2.3. m 8. .0 n. 0.0 0 fig. 0 03 08. 0 0.0 0.0 .08 _ > x flu... 8.5 05 0:88.. 85. u:- a I :9... so... 8..-am. .0 502...... 02...... 08 08 08 8. 8. 8 0 8. 80 80 8.. 8.. 8w 8. .08... 0.3 88: adv / 0.8 .IIIPV/ 0.00 08.0 2N ....... / Qua 8.0 new ....... I 0.8 I I 0.8 I 0.8 365' // 0.09 m 3020' /> 0.0N m 3+ [MW 0.2. ) .8191 in 0.9 ) W/ . 0.0 m .-X 0.0. m. . . I /, . 0.0 . 0... be 0.0. . 0 0 0.09. 0“. 2. 8.0“. .2. 8.0.. 8080885..388“.3.x 56.381.985.31 .8553 .322 mum m .80 .03. as... 122 ET. 3700 ggwn g3“. 37%- 393?. 59530 eiFZV Elan! .520 9%! 885m E252 3.0.5.16”. REESE—.00 5.: m a? . no.5. 9:201 ~55.» Io\\ \\\\\\\ to. \1\\.\\\\\“\‘ .090 054802 serum 851 62.9.83 ”.322 mum m .30 .2-.. 2...... 123 .8583 “32" Em ¢ .3 53 oz: 3.2 .593.— .U .4 :95. 124 125 803 Sig 9.100. is. 8.80:8 9 0 0 v a o N. Y on 8 or on o h p P b b 8+w.° A.::vul - h D ...-Em o ‘ 8w 3&0 m \ a // 8v m ) 08 M 10/1 8+”. W ...fiill . ( 80 3.1: o . 08 00.33.028.01 .Bsgzi 05.82.9020“. 58001.}... .000: I003 .000: I00: 303 8.1.3. 021.01 £05 8.30000 or 0 0 v u 0 ~- 0. 3 8 or on o h n p b P gw.° AENVKI h F o 43:“ O \ 8N / am \ i I m 80 w 09m; $52!. /. ( 8° ....za . , 1:111:13 - r 8.. i:|--1.‘ v1.1: 3+m~ 3.62.00 8.1.01 9.1.01 .3 E0502 3.8001 05.00.? .0> 090“. 8:000: I: 1002 .020 new .0002 .020 EN 303 8..-8: 2.0.2. is 8.62.00 9 0 o v a o w- v. 8 8 or on o . _ . _ . 8+wd an . k o 385. o \1 8N X 3&0. m \\ 80 m. f , a. M ... / 09m W W €22 . L . iii}! sf: - , 8° ilrll :r-‘ it 3+mn 3:00:00 8.18m 050.0: .9 E0502 00.6001 05.002 .0) 020“. 03000: a! .0002 .020 Pm .0002 .020 Pa 829.33 0222 mam w :3 .o~-< as: 126 80 09 80 08 08 08 80 on. 8. 8 0 8. 9300.? 8. m 08.80.? a 5.80.? ._ 8 00.80 . 7 n 0008:? v 0 .. 00.8. f o 8.0: 1 8. m §+ w 00.30.? 1 00.02 .... 8.. m O'NIOI O 09. w; L: 8“. .80: .86 can Sgguaaitcgguaisazoasz» 8 0 g 80 08 09 80 000 8. 0 0.8. 7 11.1.: . 02...- - - - 0 2. . 0.00 08 If! g g ....... I; 0.8 0 00 m 80.0 .20 ll. /./ 0.3 m I 09 M 5 n /I 08 ) II . w VII DON / 0 00 I; w / ( I 1 0.0. ( . . x- , v 0.0 . 1 1 1 . . 0.0 0.00. 1 . -1 0.0.- 20 88h. 2 88“. 8080.085. .9 88085. 83800.01 98.0001 0809.83 £322 000 0 .30 .0~-< as: gwgn gwgr I gugw I .020 Dania: coltsw E95! 0...sz .58: _ 303800.00. .. 0.. _ 5...: 0250a 2V 020". .690 019.08: gum 020“. 0809.53 0352 mum v .80 .0~-< 0:: 127 «8..—8.003 £9205 mum u “now ._N-< 030,—; I :6 906— " 55-3030 3 «£000.. .04 :95 128 or o 303 8.30m 9.00m 0 v a 5-1 p o $.32 I 30!: o // 1. 1.03.0.5 8+md .3... m 2w gm.“ 00.30: it .9 .0080: 00.8001 §§§§§° (N) was or o .08. 80300. 930m 0 V N 6.5: I 30!: o / / 1002.020!“ 8+wd gm 0231K .0) is! c2800". :55 882.0 8 9 8 0 Earl: . 305 o K1 M 808000 03.00.? .9 00.0.... 00.00001 3.. .80: .08: £05 00:00:00 00 9 8 0 600: ..I . 1 0th o \\ \ N N so; §§§§§° 05.82 .9 8.0“. 8.08". .0: 082820000 (N) ””5 or o .003 800.3. 800.0: o v N o // 6.8: I 08:: . / 1002002005 8+md c2331 93¢ .0> Elbe! coaulom 055 808000 0' ON 0 6.8“. 1. 30!: o \\ \1 a §§§§§° 053 .9 020”. coim sf 332.0203 (M) can: .8353 0322 000 0 .80 ._~-< as: 129 8... 8' 8' can Sn 8" 8n on. 8. on o I , 8' 9a r. a. To! 80 SM 2. SN :8. .9.va 8F 80 am ed 3.80 I... B. :3- 5.8.. 838.? a 80 an... 5.6 8 8 :8. .38.. 8.80:? 80 . _. 2. as 8.80 . 8 :1”... .9. 5 $08.... a 2.0 3' 6.0 on; aO‘ W on :69. 9.0.9 . 3' s. 2' ed 8.03 a ...8- 39. 3.08... m o: :3 3 3.37? 1 Sn 3 .5 o a 9?? a 0.8- :3: o .8 w. a .3 8.. ) on 0.8. ““03 m 08 3: 06 o «.3. «.03 . o8- ad ad 88.. > as: .85 can 20.3.3.8- 13833-853852» .th 8..-Wm o8 8... Ba g 80 8... 8. Sn 8a 8. 80:: u u u f 80:... - u - : o 8 0.8 g g ....... 5850 EN ....... 0.9 .86 .5 II: . m .86 .5 .ll / 0.8 m 3+ . \II E . 0.8 \I . m 2. m . - o.o . r: - 2:. sz 080“. 8.6030005...)qu2 2v 8.0... 8:00:00 5.... .9 090m .91 .8533 mgaé mum a sum ._~-< as: 130 €22. Ens: b0v5 oxalate} 005:5 EoEoE Ass. 8:850 at m. 9 m m . a . ~\\ .020 054082 Ritam 090“. 303 8..-um _ 2%: 250m .8583 £322 95 a sum ._~-< 0:: 131 .8583 £335 mum a :3 .--< 03.; 3.0 00.9 £30328 I no 5:02:50: 2.33 .0... :95. 132 5‘81 £50.01 .9 E9100“: c231 .082 30:... 05.022 .0) 005m 8.88m it. .082 Doc: (N) 0010:! 38. 8..-8x 9:23. 2 0 c v a o u. r . 1 , 8m; infiiflfifiiliflfi'ot ASE“ 8:00:00 co 00 cc ON 0 3m; m. w 8+3 $5.21 § curs. . 5111111..- .L 08, [111.11 rL 8+3 8.888 8..-6m 933. .9 .55: 8:80: 05.8.2 .9 8.0... 8.681 a: 38: .85 EN 08: i=0 3.. $038.13 930m 9 o o v a o u p > . I o H . 8a 4. 8.. m 80 m =5: l 8» 2:52. . 1 111 1 .68. 1 l ....... y 1 8m N 5.335 Siam 9.3.x .9 .58: 8.381 as; .9 8.0... 8.63m 5: 08: .85 En .80: .85 En 329.53 2935 Em a “am .--< 2.3 133 can 09 8.. 8n 8a 8a 8a 8. 8. 8 o o9 a. m 3.8. no. m 838:? u 838.? x \ - 8 agl‘l \. . ‘ v 0 m $31-..- m 3.0; a 00. 2.89 I? m 3.3“ I? ... 3.8. no. (\V 8.. DINIOI 8.. M can. 38: .85 am» 9.98.. 85. 3- E1 3 8.65... u 8 cos-on. $1.. 8 o R: 8o 8... § o8 8a 8. o '8: - - ad. 7 is i If: s- . 92 f 2: QB 5020 g ....... 0.8 0.0.5. o 8 . 820 .5 II 99 o 9 Vr 98 0.8 5+ . ) / o 8 w. 28 w o 2 ( ; o9 ( W 4 v 06 .o . 99. o 0.0a- 0.0—3 28.0“. can-coagggugx 2.88“. 8.888 .919 88“. at .8??? $93.2 mum a How .NN-< 2...: 134 3.88.5.0 gmwgwl ggw I ...-52. E98: 820 019.80: gm Elbe! 89.0 010.80: 80:5 090". 38. 5.33. :95. 9:33. “85.53 agape mum a sum .NN-< 03.; 135 Ed 8o 9?- .3 3.? at . : a 3.3a. 2.9 «$5 .93”. «in. «85 8.2 8... «a. 8..... 8.. as. s u a u u Em 25. #302 $th ozm 8.83. 9... 8- 3.88. 2.2. 8- 3.88. 8.3. oi n 5G. :3. .562 $95 at... ...So .3. . 3.8. 8.9 x 3.8. 86. x : ... 3 u a a Eggs”. 2 .3 5885.50: .8833. 53-83.! 82:5» ...-85.! 32:5» 7.85.! 82.53 .59.: .3 SE 2.58.: 9.: 86 no.» .623 85. 2 .58.: 55 max 35 Em mix ..ESE? 2:3,: . 8&9» «Qua o 3mg» :8 as as 29528 a ii... :38 N c n 3082 8.3.. «ula 8... 5+me tan .8 .93.. 83a 8m. .89 .B 3.. m 8 m 8 Ba. «..3 cc; . S .3 «com: :8“ .3 6+3. 8v; 8... 6+3. 3:.” .3 «88 «an and n E m. 3 3.0. «63 as; 8 3 8&2 38. 3 5+3.“ o8.“ .9 .o+m~.. .83. 8: 93m 8. cm.» s o. o 9. 8.: :3 nouv 8 :3 «bout. to: .9 5.va 8.8 9w 5.3.. 0.3 v3 2.9 3v 8... o o o 3 «3. SB m «av on a: 3.3.“ 302 «3 .ofiv.‘ 98 v9 5+2. N «‘8 8w «Sn 2' mod n m m an 8.2 p.58 0 Ev 3 3n «coma... 93. v8 Sou: 8... 5» 5+3 . «3.. 2n 9.8 a EN 0.0. m on 8.2 v.3. 9.3 8 .3 5.75.0 3:. m1... 5+3.— m. . 2 3n 5+3 N «8.. 03 50:. .3 «ad in. m «n 9.0. :3 0.19 an 23 .9de n3! .8 Sound :5. .8 5+3... «moo. mam :3. 3n 8.. ..o. c an 3.9 .68 as? on 2.x «com: .8». 8a 8&2. 33 m." 6&3 83 new 38. mmu a . mac. n on 3.... 2.8 0.03 on can .gwms awn! o.“ 5+m~d 8... 8m 5&3 90.. 8a 88— XN 3.0. m.n~ o a. and. fins men: on =— uo+mnd 89. no. 5&3. 8w: .0. 5+9: EN... 3. «30. mm. 3... 03 m 9 mod. :3 0.3 as our «comm... anom— 3 8.60.1 38— 3 Emma «mac. 3 28. 8 3.7 o? o v 56 3.3 mag an 3 8+3... a o 8.8... o 0 8+8... o o o o 8.0 o.o o o 8.". 95 :3 o o «.85 .13.: .13... «.8...- :E«.: :5... «<8:- 33... :5... ...-...... ...-.5. I! I. is is .5 ...... it... it... 9.. .... .... Efizxd EMQQO DEN EWQQO Qua: «3330:. u~§.8~1 51.1 L 82.8.83 £3.22~ mum HUmmm<\Ub:a... “new .mm-< 03¢... 136 NE 22 .93 8..-um 8..-.... .52. 8388 o v w o u. v. o. 8 9. 8 o .Fo mu... . . 8&6 ..cerl 0 o w 2...... o 8.. o arm.» \ ... m \« 8. m . I w 8.. ) .55 09m. m. m 68...: o M So . 3+3 8... 8.32.8 c2181 9.38. .9 E28: San-om 3232 .3 88... coco-ox a: 38: 820 .5 38: .85.. .93 8..-8m 2!... .25 8.326 o v N o N- V. o. 8 3 ON a L . r F 8%... F + . o .ucuull ”H w 8.. 3a.... w .55. o |\ 4. 8.. w ...W m .2qu! 35. m 8.. m 38:! o g M con 7! ...--- r: lull}; . 3+m a r ; st-.- #1:? -- a g o8. con-8m 5.6800 9.1.5. .9 .350: 8.8.8. oSBB< .9 88“. est-om E: 18: 820 an 38: 8.5 new 303 8.18m 9.33. .EE. 8.62.00 o q a o a- v. o. 8 3 on o . L . . 8&3 . . . o A11 in: I 8“ arm... 3...... o an\\l 5.9) 8v m. . 8o m .w - Gail: mom. ( 8.. ....z: o L . 8+3 ...:i--. . -. «sill! I... o8. 8..-8m 8.60ch in. .9 .5802 8.681 0.23.? .9 020m cocoon”. a... .82.. .85 .5 38: .85 .5 .3553 £58.. mum SEQDEP sum .3... oz... 137 .: 3.. .3 .8 5 o 8... 84 8c So 8... SN com on. 8. 8 8. a9 8. 3 8 o 8 3. 8. .a 8' 3 a 8. ON No... ‘0? «prop: IO: 2' t a 8.3.1.... 8 a 8w. 3 8.87? so 9 o 8 :18. 5 8.27? 03 on .u 5.30 IQ. an 3. S. 3.89 . m 8» 8 o «38 .0. w. ow 3. w: 3.0km 2.: 3 2.“ R a 8.8m . m 8 a. n! a 3.8. 1? 8a a o 1 on 8. m: auburn? . M 8. 8 o 8.3 no. a a 8. a: 33 o ...: m 8. a o 8.. m on 7s. «9 w B v o ( o 0:. .3 8a. a o 38.. > x _ _ .3: i1 ... 382 .020 05 l 2.98.. 85. 2- a.: .o 8:23. a I 828a 5:: on. 8F 8 8 9. cu o 8m 8.. con 8a 8. ado ado .85---- / 92 82... --- 98 0.8 new ....... 93 820 95 ll - 98 , J/ o.8 ) 5+ 41/ . o2 ) X , , a 2. w J/ m V , I . 2: m 0.2 m J + . 4 od . . I [It 4 < o o v. 2:- 0.0“: I!!! nxlll 0.0—.0 sz 88“. 8.3.30 85. .9 88m 85. AZV 090“. 8:00:00 3.: .9 09.0“. .91 Esau; £3.92 mum HommmSonfie sum .m~-< 2.3 138 3.3 8.1.8. :95. 9.3:. 889.00st 0800—600sz EEO oxniuo! 32w E050! REE.ZV .5602 a... 3.3 8..-Ba .... :95 25.2... EEO Graig! Satan 8.6.... 63383 $522 mum 553525 :3 .mm-< 03¢ 139 2. 3 Lo: :4 as 7.3 .2: «hi on 33 mm 2. 92.8. fiu >03 . 8n. 9| macaw 00m wow «0* 24 3.: x no won sun on — 8v 2: vmvvn nun N w r :N um . 52%. .0055! 2 u t tk >OOm on“ 2.93596. < 8281 . :80 09 lg .BHB agilafiimfiflla‘i 28m E'Eiimflfiaiiigiii .II kkm Ex: QMQD: kmfiko QZN fl QRW is: HWQOE RMQQO ka. £8.38 a .a .523 .3 3.: e .58.: gleesgn‘efigi, 5.83:2: 8:528 , .523 :5. a 2.3.: :2 Sun? 2 E3: ...-82 r: §E§3§§§3 o.— a: ho °§fifififi38§888$ “8.383 “3222 mum hommmSEaam sow .vmé 93¢ EEQMLVH! 140 I41 ...-six 93:. 26.8.33 «popooc~o«.v.o. 8 3 8 o 8- LT r b b P. h|b|» 8+“.o b b o m \\ 8n 0 3&3 m K 8. m. «Sill m w -\ 8c M 3+ .. . .l. ( 3%: o m 23“. can a ( 355. o 8.: 3”.“ 5g 8.1.3. 33¢ .9 E25: 8.38m 05.32 .9 88“. coco-om a... 33:35.. 33: 30:3 3.3 8.181233. €6.88»; «popoocuovvo. 8 3 8 o 8- - rL r b p L! 8+m.° Aucpwvullllb - o fl \0\. W 633.”. 0 ix 8N Venn 8' J m an m m 8» m 6:52! warm... W. 000 ( ...-I: o ( _ i It 82 3&4. 8.688 8..-3x 23:. .9 :85: 888a .53.? .9 88”. 888: a... 38: Ito 2a 38: .85 Ba 3.3 8..-am 233. =55 8.33 «.c.»ov~o~-1o. 8 3 on o 8- hhhhh p . 8+m.° P b p o ‘m . furl: \\ m A55. . R 8a arm» u ~N o9 m. m 8.. ) Gail. 3+3 m m 23:: o M 8a 7 3m: 82 8888 8.33. (£101 .9 2.0.52 8.68m oSBB< .9 88m coco-om a... «195 92 182820 .5 382820 .5 .8533 £332 mum hummm<>=aam sum .¢~-< 033 142 «9 "8 o8 Sc 8.. 8n 8» on... 8a 8. 8. 8 a o 8. mac d o a; 373' I? cop W ”a. 8.87? 3 Ale 8.87? 3 8.2.. .9. o 3.87? 5 8.8.. 9 . .o. o «38 IT M7 “9 «9:... ,+ m 3.8... . w m. “9 aé.+ a _~m o 3.3 no. a s g.— OIOI m q o . m 8.. O" W o o ( .... i. 89 o o > . x _ F 82.... 8.38m : 5 .28: 320 93 £33 85. vi .9: B 88...... a I 828... 5.: 8. 8. 8 8 9 on o as 8. 8n 8“ 8. o 2... : 1 98 0.8 99 can . . 2.. m QS w . o... , 2.“ 2... ) call 5.? ) u , , - w 0.8- w , m 1 0.8. q ‘ 0.9. 98. AZv 83n— AZV 090“. 8823 85. .3 8.3. 85. , 5.3.8 .9... .9 83”. .9: .8333 £58... mum Swamfieaam .sm 3.... 2.3 3.3 8..-8m 5.: 9.5:. 808.0. \\ \. 88383 g? 808—.ng Sooougoomwl \ AEEV 8:00:00 5.... b2:0 mini—85 885w Eon—02 0.807- 08de ooouécowfl coco/SONG oooTooonI 08méoovfi 250698- .83 8..-Ba Es. 8.62.8 .9: m. . ca...» 2.28. E. 8.9". .520 Dani—82 835w 88m .8583 £52.. mum 533.5% 3% .3... 29¢ I43 32.9.83 £32... mum SEQEBZ sum .2-.. 03.; .58... ,3 .32.. 144 1. ...-38.1.3. 9.1.8. :5... 8888 fireouvuovvo. 8 8 9 8 o hhhhhh p r gw.° b b h o blur/Pain} o Sm o o M o o 4. vewd M 00¢ w W coo \m . II . ......Ell ( .5... 8+m . m 8» ...-...: o ( 33.5. o — gr 8+3 8888 8.1.3. it .... ...-E2. 8.603. 05.8.3 .... 83... 8.83: a: 18: lg: 13: Inc... .8... 8..-8m 9...}. .8... 8888 c.u.o.ooc~o~.vo. 8 8 9 8 o pppppp h p 8+w.° r b \hl“|k o n q 8m 3.3m \ 8.. m. I 8.. m gill 3.3% .....«II 8» ( ....z: . m ......E . F . 82 3. ~ m 882.8 8..-.3. 2.3:. .... .88.... 888.. 3.83... .3 83... 888a ...: .82. 820 new 18: 820 EN .8... 8..-.3. 2...... :5... 8.880 :Nwopoovwcw.v.o. on 8 9 ON a pppp b p gm.° P p b o m 8... 368 m -\\9\4 8.. W m, RN 8.. M I. . earl: .82... ...?m . ( 8o 3...... o 3...... o h 89 86..“ 8.8.3 8..-.8. 3.13m .9 .88: 888.. 3282 .... 8...“. 888.. ...... «195 38.... .85 .5 38... .85 ...» .8383 828.... mum Smmmscooz sum .3-.. 2.3 I45 ... ... _o9 .8 8 o 8.. 8c 8.. 8n . 8n 8“ 8N 8. 8. 8 o 8 a... .9 + 8. s 5 o «u . : .9. ...v d 3 ....- . .. . 8. "P ... ...-.8 no. m .... r... o 8.03:... a .... 9.- a: . . .8 3 o 8.87? \W 8 8 3.. .... 8.2.1? . . .. .. 8» m] o .38.? . 2“ ov o «0.08 10.. W a x... 8.2... .9. ... E iv 0 6N8" 0 8: W . C... a W 8“ 8 El! 8.8.1.- .. a 8.. 8' 8.8. I? .8.- u 8. 5 o 8.3.? a ..u 3 83 0.9. m a. o 8.. ) m a 8. 8. m - n o ( 0 _8p. 03 oom- .. . MW .83 > x F _ 8... 88.3.. .. .... 38.... .85 can 8.83 8.... u..- .9: .o 88...... a I 8.8.. ......» 8. 8. 8 8 3 8 a 8m 8' can 8.. 8. a a8 a...» '25 - I n .- 2.» I 2 ..8 - 0.0“ ‘20 EN ....... I 0.8 cm... 8203' ..-/... 2... o o 8.. . 2... 2... m 0.8 w 98- m . , . .m- . ( ‘ I .- °.°P o 8. - - - . . . 0.9.. o a OS. 0 o.- fv 020“. 5.30.30 85. .2. 8.0“. 85. .2. 88.. 3:00:00 3.1 .9 8.0m ...-I .8883 £38.. mum Sway-<82 sum .8-.. 0......- 146 So... 588: Q. _ 8..: 8.5.3. .8... 8..-8m 5.... 2...... ... m. .55. 8.3.8 .9: :55 5.52.8 a-.. , D69“. D 890 I oovéou U W\\- \\\ I ..wquwmww.-H H.385 \ 35.2. .88: 8%. 2.83". So - .\ \\\\\\\\\ \ ‘\\\\\\\I \xx- -~\\V\~\‘ o8..8on_ \.\.\\\\~. - 8.68.- 89 8.88.0 08.. .580 0193.82 scram 88m .620 051082 885w .5802 32.9.83 £9.35. mum 583282 3% .8-< as... 147 .8583 £32.. mum 58385:. .an . 0N-< 033. .58... n... :95 I48 «is 93 .8... 8..-.3. 2.3:. =5... 882.8 ~.o.uov~o~.v.9 8 9 8 o 8. .... p.» .o h . + 88... ......E-l w o o. w 2...... o o .1 8m 38.... w > 09 ... m N m M 8» w 8...... l 3m. . .m- 8» 3...... o 9 . I::! ,1--. i g 8... .ll. - l 8#mw €0.20 8..-.3. 981m .9 .88: 5.88 2233 .2. 88“. 888m ...: .802 52.3 .802 .32.... .8... 88.2... 2.8.. 8.... 8.88 «wopoov~o~-v.o. on 0‘ ON 0 ON- 2.. o .25.! .\«\ o w 2...... o - 8m 3m. m m. W 8v m. M So m . w ( ABNVEII mo+w w (M can 8...: o g g -- --- 1!! 1+!!!- 4. 8o. Fliilzrllvi: Vii}--. . . 88 a 882.8 8:83. 3.8.3.. .2. E23: 8.803. 05.82 .9 8.0.. 5881...... .822 .020 new .032 .020 new .8... 88.3. 3.18 .....5 8.82.8 ~.o.oov~o~.1o. 8 .... 8 o 8. 88 o in... I .\q\ o .68.... b\. a m - 8.. m . 8c .M I 38.. w ( 5.82 (w.- 80 33...... o .34 L 88.“ -.--.iln-tixf-{-L 8. 8.32.8 3:53. 3.301 .9 E232 5.3.x 05.3.2 .9 88... 3.831 a... 18.2 ......o .5 .28: .85 En .8533 88...... mum 59535.5 .8 .8-< 2...... 149 ... ....0. ..3 ...“. m. ”3 8 8.. 8.. 08 80 08 8.. 8. 8. 8 0 0.... 8 0 :-!-I I... 8. 8 B. 09 . d 8.. 8 0 . m. um. ”.1 “9 ...-.00 I? 8 m I 8.03.? ... .~ W .3 08 .8 0 8.08.? - 8 m. ON 8. cl: onéwvlol . . --.-312+“... 03 W p 8.08... a . - ..-..- x m «0 w.- 33 8.03 . . . . . . . , . 0 m 8.. «v F «0.08 no. N 8 ...... .8111 8.08. 9 a 08 a 0 8. I] 8.08 . w 8 w 03 2 SN N” a 3.009 :0! u an m?! .mw 3.9.2.19: 8.- u 3. ma . 3.5.5. .6.- d 8 W 8' 83 0.6. m 8 .0. .P w B P o ( 0 W 8. 1: - - - 0 0 0 8m .83 > 8.... 88.3. i... 38.2 .020 8.0 8.83 8.... 0.... ...-... .0 88...... 0 I 8.8.. 8..... 8. 8. 00 8 0.. 8 0 88 8.. 80 8... 8. 0 0.8 1.?! ! .L !I.l!-! :.! - I .. .- i. i:. n!- -z-..I v! ,-! . .i 5.!!! iv)... . .85---- . i.e.-.--! 00. 0 8 0.8 ‘20 EN ....... 0.00 .020 EN ....... 0.0m seagull 08 3202...]! 03 - 0.0. / 0.8 .8le.. 0.0 ) ..nulbl / 0 8 ) 0.0.. m - w 1 °.°Nl Odo. .2. 8.0.. c0888 35. .9 020“. 85. 0.91 .2. 85m 8.62.8 ...-... .9 8.0.. ...... .0298... ...»...5. mum 588.83.... 38 .8-< 0...... 150 A603 5'61 8..: 3.3.x .... 8 8.5 8.3.8 ...... oooomlooooo... 0.88m.- 880.00 cocoordooomn. oooomwéoooowl ooaoowéoocmw I \\. \\.K s . Q ‘\\\ \\\\\x\ I .020 Optic—BE 885w E0502 2.6.2. .88: E E C .9 '5 o E." I v s I :- sz 020". ...... .8... 8..-.3. .. 8.... 2.8.3. .090 omnivoz 855m 8.0m .8583 “.22... mum 583.52.... .8 .8-.. 28.. 151 . 2.. 2 u 0.6 5 av 22 ... 8' a... Elfil as 3 .3. .m n ttn... SE H.303. txuz... n «km 25. 4.30:. 59.0 E x o 38: .85 9..» c 00 E c: I n. . 8..-s... .. .33 5. ...... ...... 8.... n. 828.. 8..: 8a 8. 8. 8 o 8o 8m. 8. Sn RN 8. . .I-lil.- 0.05 .- I 1 J ; :l-.. .l..-----. .. f .. .-.. .- 0.8 50:3--- 32.3-:- . ode + 0.3 '30 EN ....... 502° UCN ....... 0.8 9.? , 2:. ......I ...... ......” mr, ...-W . .. / . in 98 .m- fiu / // ad... .m- , ed“ M ..... . / . a? .m- . ., , o8 . - a... II 0.0 - c - I. I - I ii ..I- .- iI- . Ode- : 8.0“. 8:00:00 85. .2. 8.0m 005. 8.888 5-: .2. 8.9. a... 82.3.83 .32.... mum 8858.; sum 8-... 2...; 154 800900 cocoopéooomu Scoop-0008*. ,8ooou-oooomwl EEO Dan-.082 8&5m .5502 ::sz .85: 3.3 Siam AEEV coauofio 3.1 c. Q. . 20F: O>=GUK ... - ... 2 8.8 . ...-9‘ .020 GEM-.082 coatzm 8.0... .8333 £322 mum 83025 3% .R-< 23 155 i E 2 ‘5 IE? fl ktm 2‘: Hwfio: khnko QZN u kt.“ Z§< HWQOE EMQQO 5km. .52... .u 4 .35 2 E39: .523 85. 2 ER... , 22 8mm? 8 2.3.: :8 .s 252:8 a :83... 2.58, i .%i% iggliguafio igfii Elgar»! El HIE—Mali Elm: iEimmi—Nliaiig ”'32 QWQQQ QZN 33.... T35. Quake GEM. AAAOJNQO: .5393}... I .6 tea-use?! 8:85:92 , 32.9.83 2322 mum 333:5 sow .w~-< 2.3 156 157 so... 8..-...... 8..-3m .55. 8.9.3... o o v N 0 av on on or o b p - §N.O F; p p h :7 0 E3... -I -\fi com d 3m: :35. . h\o\- a / ) 8? m .5?!- mo+m.. w .1 8o UN. 3...... o ( 8.. r a 8m." .. 8.: 5:08 cog 0.50.0”. .-> .5805 8.300.... 05.03... .9 8.0... 5.600". a... .82.. .85.. .98.... .22.... 303 5.33. 0.50.0”. :55 3.30.30 c o 0 u o 00 on on or o . . Sm... .- . . . o 6:3“. I. T 8N . :35. o 4. - / 13.... M J .lIlil! 8.. m A832]! memé fl Ox! 8Q \lle 3...... o g M So ...-E- ; .1- ..--um i- . mom...“ .- 35:- .- -82 cc:0 8.8300 9.3.0”. .9 2.0.5.2 8.3003 82.32 .9 00.0“. 8.800”. a... .28: .....o Em 18: .220 .....u .9... SEE 9.0.2.. 2...... 8.82.8 o c v N o 9. on ca 9 o . . 8&6 T . . . o .20... II 8... . O - 3+3 ......m \\- 4. / \l \ 00‘ m 8.9:!- . w .1 o8 m new w m ( 3...... o ( 8» will!- - “Lg-$.1- mO+wN r- L. 009 . 9.0.8050 9.3.0”. .3 5080.2 8.800.... 0.3.00.2 Q. 020.. Sim E... 38.2 ......o ...o % .28: .....o ....” .8232... .328... mum esteeam ...om an... 2...: 8... Sc 03 8m 8m 8w 8.. on. 8. o... o :1! 8. d 838:? 8. m U ”2.8... . 8. 8. P3 . 8 8. T8 9 2 .8... \ m 2 8.. 8. . -.. 1 8. 9. o 25 . m 2 2.. 4m. «38 . n .8 5 o a 3 a? 8. 8-39? . .... '51 _oo o ounovulol R an 8.. a 1. ..8 m o 3.2.... . m 8 5. 3.. a an a o 0.7... m 8 o... 8.. - w 9.... ... o m R .37 P: ( 2. : a {II}: I 1&1--- ..-L 8a. a ...... .3 o o o 88.. > [x— 182 .85 Q... 9.68.. a... a... 5.8.... u 8 .. fi 85. .o .. cogs... .3...» 8a 8. 8. 8 a 98 .85.. - - - - A can '20 EN ....... 0dr .I Vflv g. .35 Ba 0.» + v . 5+ ,1 w... 0... ) [WV 3 m . 2... m I . 98. 3%.. z 85m 8.30.30 85. .9 88.. 85. gonzo“. cocoa—.00 «a... .9 8.0m 3-... 32.9.83 .39.... mum egéaam ...om ...N-< 2...; 158 ....o... 8..-...: 6...: 9.1.... .020 010.382 Satan E055! .EE.z. 525: :8... 8..-am 8268...... :5... 8.62.8 ...-I .. my :25 93... 908..- L r c. I. m o8..ou , . 8868. .... . . 1‘oooooo/ . vcoououuuug 8838- >>>. pfisssuusswuuhhn- o8v88 a 494.49% QMMQHNAANMVMQ. 2.8.0... 88 >D.pr§.~\\MN““v , ~ 8889- 88 ’4...‘ .-..s~ .q- c8382”. v D § \§- 000 7 -§~§§NN““\ 88 ’ .020 019.085 oootaw 020“. .8553 .32.... mum 33.5% .....m .3... ......fi 159 n ERN 2.: HMQO} khfiko QZN H kkld i: dunno: khfiko ka: 2258.8 u. .3 88.2.8: 8:883 "58... .u .. 8...: 2 .58.: .58... as. e .58... , ...... . 8%.. ... .83. , =8 .... ...... 23.5.... n. .83.... 2.88 , .6032 5... :.o. 8... 8... a ._ 2.: 3.1 as... 3: 8... 8." 2... ...... 8,. 39 z . x u 8.88. 2.3. E- 928: 8.3 a... «$8.. «.8 8.88 9.2 . 8.82 8~ .. 3.82 .2 ... . ... . x . ...-85.! 82.58 ...-85.8 81.55 7393-: 32.58 a... 8.. 8.. .952... 559... .2me 8&8.» 8.8... 8.8: n a «9.39“ as: who 5+m§s Sp «8 Sowoos Sons m8 930 v8.3a Bows: ~99 08 «com»: 3°: 05 «9qu... 893 In 989 «Rd—o 8.9.0.. 23w 0?. Sowom.‘ “BR 08 5+3... :2“ '8 E 88.2... «9va most 8v «9on 85— 9m 8+9: 5o: 2... ovum. 03.3... 5&5. 8E. onv 5&8... m8: 8v 8&3. 8...: o; 38. 8.9 «0.83 §— 3n 8&9.“ 8am. «on 5&qu atnw 9n otfi unomn 3&9.“ 0:! «No «9&8... «SN. «an 8&2; 039 SN KNN. 8mm mun 8&3; nun! g «oomoud two? at rigs 83.. 8— 03: 338. 8&8. o 0 858.0 0 o 8+§d o o o o 9.3:- ...5... .5... «to... .85! :25: «135 .13! .38: ..lgl .435. Efizfiu QMQEO DIN EUQQO DEW «Ago: .8383 .1522 am 8.5.82 sum .85... 2.3 160 803 5.39.. 02301 EB. 8.8050 or o o v u a on ov ow . o h b b b 8+w.° P + 0 25K II » 8n + 3+3 235. o a / ) 8c m 25.: I. 8+m.. w \ 8c w :.z: o M So r . m 8+m~ 89 5'8. 5.8230 3.13m. .9 .5602 ego-om 82022 .9 88“. Siam E: :8: 32.: ..Bz :25 303 8:83. 02.33. :55 8:00:00 9 o o c N c on at on o r Ly b r gw.° w b P o ecu“. ll 4 8N + . _ 3:5. o ... / 19m w. W oov m 3:82 I new... W N 000 W 3%: + t 8.. _ ii. 8+3 L o8. 5.18m 5.82.00 in. .9 .5602 5:003. 05.3.2 .9 8.0m 5.38m E: .302 820 vcw two: 820 new 8.3 8..-am 938. :55 c2688 9 a o v a o 3 av ca 0 . . 09w... .F b . o 8.8“. II can . o 1 3+w.m A 37.! \4 ... . . I . w ~ 8» m an : 8+w. m ( est: 0 ( fi 08 Faiiillr€ir -. 1i} mo+wd L 89 5:20”. 8.8800 9.11m .9 .850: 8.68m 2282 .9 88“. 8.881 a: .03: 520 Pa .803 820 En $3983 82% mum £3282 sum .83. 03a. 161 288“. 8.6800 35. .9 88m 85.. A20 88“. 8.32.6 .9: .9 8.0“. 3+. 80 00.. 000 08 03 000 00. 00. 8 0 - slla amp .0 m 02.08:? 8. w 8 2 0.80 . m. w H. ”W 058 '9. m n: 8.. 09 ..T 08 8 0 3...: m 2 0:. on. «0.8... . m .8 or. o 3 3 5. g 8.000 -T w :0 8 01 8.9.0.? w a “PI ...N' I. 80 0.. _ 3.0: .0. . m ON o3. 6v .0 «a Tm E1 90.? m 8 +12. 39 00.. w 02+. 00 4.0: w R 0 0.. 3|? ( 0: a $ 4 I I. l .. 11-..-.. 000- 0 E. 03 0 0 _0 _ 008.. >_ x .08: .020 9.0 05003 85. 05 gr: .0 8.300... a 8 80.8.". 00.5 _ so: 8..-am 0 p 000 8. 8. 8 0 000 8... 00+ 80 000 00. 0 900 I l I I; .I It!!!” en I» I all! lit! 1: .1... .3. .--- ..A. 0.8 .85---- 0.8 0000..---- . M 0.9+. 0 8 ‘ ....... to 000 ”W 0.9 0.0.. III 0.2 805 En . 0.8 If 0 ... fiulol . lull: 0.0 ) 0. ) // 0.0+ W ON W 0.0.. ( 0.0. ( . 0.2. T 0.0 -.. 0.0a. O.DN. .l inl r ill? 0? II I: 0.07 329.003 050.2 mum 303.002 .80 .0~-< 0.000 I62 303 8..-um 00...... 3.301 0. .320 Graig: 885m 50:52 :55; «co—to: .900 0330: 00:: 900.0: EEO 019.285 88:6 8.5". 0809.003 0305 mum 082082 08m .0~-< 2000 163 i .8336 : =0 00.08.0000: 808.000.“ .5003 85.2.53: Ina 288% 2.50%: E =8 2. as 23.58 2 000808 20.38 .5802 all? . _ . 0.08 0.8. .2 8. Eli—ugh. i 662 ... lil 9 0 n0. 0 9 :.! 0.08 0.09. a: .8 EEE 5.33 Eii 8L8. lil 2 u we. 0 N. 8,0. 0.00... 0.8. 2 § Efiigfiigfilfimfl. 8 0 ... 8. 0 8 «0.0. m a... 0.09 3 :0 ggigligfiifimfl. 0. h 0.00. 0 8 3...: 0.30. 00: «n 80 Efiiafiaiiiimfi Sn 0. m. 0.3. m 2 3.2 0.0% 0.0: a «an EfiigliE—fifilii 20 50 0.0. 0.: no: 00.... 0.5 a $0 gfiigiiialil “2 t u 0.0. m. z 08, 0.05 0.09 a 0: g . . 8 0 0,0 0 0 00.0 0.5 0.03 0 0 0:55. .11 00a :3 ii :0. .III. I: .lli. 9c 0.5 :5 ENQEO nth. «Ago: hulk E ts. I E SW .l. RR.“ 2‘: QNQQ: RNQQO GEN fl kg” E: QWQO: RMQRO GEN .500... .3 00:: 2 2.8.2 House—.53 max—g mum 03253:.— 30m .om-< 03.0... 303 c2381 9.33... :55 5.6050 o. o o v w o oo ov on o . . . oo+w.o N p F . o m 25?. II cow 1/ 3+m.m m 3521 o \ 8* W 5515: I newé w K 80 w 330: o M 000 ( . 8.0.0 .... ..- 000. 3:88. 8:088 9.10". .9 E282 coon-om 05.3.2 .9 88m coon-om a: loos. 80:: 332 825 303 8..-8m 9.38 :55 8.888 o_. o o v u o oo ov on o b P b b gwuo W + b O m €05... .l. 80 . ....Em o ... c [I ) A0 «5 006.. m 000 w 23:: o m 8» wt izllti E 3i: -i . no+w.~ .- - ooo. 5.33. 5.62.00 9:201 .9 .5802 563m 8239‘ .9 88». 8.331 a... .0002 820 won .282 820 new 3.3 8..-um 9.13 :55 08808 o_. o o v u o oo 9 om o L oo+m.o A {v L . o En u I w . 000 3&0 :05 o \o\\. ... :88: ll 3&6 m « ooo UN, 30:: o M ooo r1: ............... new.“ 1115: Era: Ewli 1 8o— coism 3.6200 9.30m .9 .580: 5:081 8239‘ .9 8.9. 8.33. RI .322 820 En .322 .820 En 0809.33 2305. mum 022.55 000m .0m-< 2000. 165 08 09. 00.. 000 08 000 000 0.... 8. S 0 I .- III.-- ---III-iII - I III II II ..III -- .00. d ... m 00.08:? 8 0:08 . B. .3 . 8 m. «0 0 0500:... - I. I_ ... .3 .35 .9. o \I........IIII.IIIH\\\ m 3 1.: u! v o m. ml". “3 «0.000 .. w 3. 00. 050.7? 8- m 0' c .....\\I Z 8. H: 00.03:? - u on o 362+ oo? .0. 03 m 8 0 0.0.? m ... .... 3.. W. ... 09 ( n. “a III- I - - - 12.. 0 Am! >0 Ix.— 0 1002.020 0% 0.08.. 85. 00- .0.: .0 8.300.. 0 8 000.80 00.... 30.... 8..-Ba 0 80 0 000 08 8.. 000 000 8. 0 . 0.0. - I I I. I- - ,. . .. . -. -. .09 .80... I I 0.0 5000...... Off 0.9 . 0.00 2N ....... O 0 .‘EO 2N ....... ..8 3. V as 305 Hall. 0.0.. :05 0.0 % 0.0.0 , 0.8 5+ 0.0.. m fiulol AMI!» 0.0. w 0.0N. m {1 0.0.‘ w . 0.00- 0.0 0.8! 1 4 4 0.0 0.00: I: I I I I. ---- ---.III.I-ItI.|I v0.n .20 8.5.... 8.3800 005. .9 090m 005. g 88m 00.62.00 a-.. .9 88¢ 3+. 30.9.3... 0.322 000 000.20.10.80 .03. 0.000 166 .000. 833. 00.0.0 9.30.. 0. 0. .50.. 00.60.30 a... ‘ 1‘ I Q NWQH§\$§4\ I .200 02.8000 ...... z 0.. .1... .— a: ... ‘ =5 . ......0............... .000. 8..-am . 00.00 900.00. 1V ‘ i . a 005.2. .850: 0000:8006 .2. 8.00. 000. ..Fvuvpfladwgg . «HMHMHszI 80.10000 I 000. ‘11 Y .80.... \ - A. . .-....lv. . 0 . 0‘. ‘ O 04 v... 0000070000... I \\\M\ 08780.6 000 ’4‘ IHOWHWHNHWWII 710. ~\\\\\ 0000300000. I \ 80108.6 000 5 80000080.... I \ 008.I . 00080000000... 0000a Eon-800mm I 807805. I .005 9.0.3002 80.50 .820... .005 00.0.08: 80:00 85“. 000000.003 0.0.005 mum 00.5000: .000 .00-.. 0.000 167 APPENDIX B 168 APPENDIX B Appendix B contains an umnodified comparative report of results fi'om testing four seats with both the SAE 1826 [4] and ASPECT [5] manikins, entitled: ‘Seat pan load-deflection measurements for determination of seat factors with the ASPECT butt- thigh (ABT) segment of the ASPECT manikin prototype at MSU and the SAE J826 manikin’ by Robert Hubbard and Nathan Radcliffe. This original report is presented as supplemental material with description of experimental methods and results not fully discussed within the body of this thesis. References made within this report refer to the reference section of Appendix B, and not the body of this thesis. 169 Seat pan load-deflection measurements for determination of seat factors with the ASPECT butt-thigh (ABT) segment of the ASPECT manikin prototype at MSU and the SAE J826 manikin Robert Hubbard and Nathan Radcliffe A. Introduction: Seat factors were first measured and documented by Hubbard and Gedraitis [ l] in the development program for the ASPECT manikin [2]. This initial measurement of seat factors was done using the SAE J826 manikin [3]. In a recent study for her BS. in engineering at Michigan State University (MSU), Bogard measured the position of the H-point and the lumbar prominence using the new ASPECT manikin [4]. Bogard’s measurements were done in seats that were used by Gutowski [5] in posture measurements with human subjects. The purpose of the present study is to measure the vertical deflection response of the ASPECT manikin butt-thigh segment (ABT) and the SAE J826 as they are progressively loaded into selected automotive seats. The vertical load-deflection measurements are the basis for interpretations of seat pan responses, which are important seat factors B. Methods: 3.1. Loafing of the ASPECT butt-thigh segment (ABT): As in our previous work [1], weights were progressively added to the ABT of the ASEPCT manikin prototype that we have at MSU. As the weights were added, the heights of the H-point axis and ham of the ABT shell were measured above the floor. The H—point axis was measured at the same distance from the center on both the right and lefi sides of the ABT. The height of the front of the ABT was measured at the upper edge of the ABT shell on the mid-line. Weights were added first to the rear part of the ABT on or near the H-point axis, and then weights were added toward the front of the ABT. So that the loading of the ABT could be known In terms of equivalent loads at the H-point axis and a lateral axis through the front of the ABT the ABT was suspended from lead cells that were attached to these two axes. We recorded as weigWsuelyaddeQ The lateral axis through the front of the ABT was 384 mm in front of the H-pomt ax15. \ cords to front of ABT 0 Figure l. Weighting of the ASPECT butt-thigh segment (ABT) Figure 1 shows the ABT hanging by cords attached to H-point axis and a lateral axis through the fi-ont of the ABT. Threaded steel rods (5/16 inch diameter, 5 inches long) were used to extend the HJC axis of the ABT. Figure 1 also shows the ABT with all of the weights added in the loading steps shown and listed in 170 Table l. The torso weights with the ASPECT manikin prototype at Michigan State University (MSU) are built up with eitha‘ six or five steel plates. The weights of the MSU ASPECT manikin are partitioned differently than the final and standard ASPECT manikins, but the weighting of these manikins lS equivalent. (Ind step 0 in Table l was accomplished by placingthe 1826 drafting template in the fest; D, ‘ The H-point axis for the .1826 template was extended to both sides with a rod so that the height of the H- point could be measured while the J826 template was in a seat. as shown in Figure l and Table l. ues of moment about the l-l-point axis in Table 2 were computed as the load measured at the front 0 e ABT multiplied by the distance from the Hopoint axis to the front of the ABT. In actual seats, this distance is at a small angle (typically less than ten degrees) and the moment is computed with the cosine of that angle. Table 2 gives the loading of the Hflmtaxis and moment about the H-point axis for the ABT and weights Th The H-point of the J826 template is 100 mm above the bottom of the template, which is the same distance M the H-point axis to bottom of the ABT directly under the H-point axis. However, the contour of the ABT under the H-point is different fi'om the J826. The ABT contour was based on measured human contours and the region around the bottom of the pelvis is more prominent and less flat than the 1826. This is probably the reason that the H-point axis with the ABT typically is approximately 10 mm lower in seat than the with J826 once both manikins are loaded into seats, which is apparent in the results belo‘al Table 1. Weight placements for loadingstemfor the MSU ASEPCT manikin butt-thigggment. Load step q 0 J826 template for zero load reference (a 1 No weights added to ASPECT butt-thi ent 2 Two 6-plate H-point weights on H-point axis 3 Step 2 plus two S-plate H-point weights on H-point axis 4 Step 3 plus two S-plate Hyoint weights on H-point axis 5 Step 4 plus two torso weights on H-point axis just inside shell 6 Step 5 plus two torso weights on H-point axis just inside shell with slot down and back rear edge restinflgainst shell 7 Step 6 plus two torso weights behind center structure with slots forward 8 Step7 plus two torso weights behind center structure with slots forward just outside weights added in 7 9 Step 8 plus two torso weights on center structure with slots up and back lower edge against nylon busthgOAm two torso weights were not used) 10 Step 9 plus two thigh wei 11 Step 10 plus two thigh weights 12 Step 1] plus two thigh weights I71 T the Step Load on Front Moment H-poin about Figure 2. ASPECT manikin butt-thigh segment in a Saturn seat with all of the weights positioned as the ABT for the loading steps of Table l. described in Table land as shown in Figure l. 13.2. Deflections into seats of the ABT when loaded Figure 2 shows the ABT in a Saturn seat with all of the weights positioned as described in Table 1 and also shown in Figure 1. Tape was added to the weights on the thighs to keep them from sliding rearward. The positions of the H-point axis and hum of the ABT were measured up from the floor using a meter scale with millimeter graduations. 3.3. Loading of the J826 manikin: Seat pan stiffnesses were initially measured in our previous work using the SAE 1826 manikin [1,3]. Measurements using the 1826 manikin in selected seats were repeated in the present study for comparison with similar measurements using the ASPECT manikin as described above. The back angle of the 1826 manikin was adjusted to be vertical for each loading set. Figure 3 shows the 1826 manikin is a Chrysler LH seat loaded with all of the weights in it as described in Tables 3 and 4. The loading used in the present study was the same as in our previous study [3] except for the loading of the knee axis. The loading of the 1826 manikin was similar to, but not the same as, the loading with the ASPECT manikin as described in the previous sections, 8.1 and 8.2. As with the ABT the cosine of the thigh angle was not used to compute the moment in Table 4 but was used in data analysis. In the results section below, the responses Earn the ASPECT manikin and 1826 manikin will be compared on the basis of load- deflection of the H-point axis and moment-angle about the H- point axis. Figure 3. SAE .1826 manikin in a Chrysler LH seat with all of the weights positioned as described in Table 3. Table 3. Weight placements for loading steps with the SAE 1826 manikin. Load step 1826 template for zero load reference No weights added to SAE 1826 manikin with torso segment vertical Two disk torso weights on H-point axis Step 2 plus two disk torso weights on H-point axis Step 3 plus two disk torso weights on H-point axis Steg4plus two disk torso weights on H-point axis (Two 1826 torso weights were not used) Step 5 plus two cylinder thigh weights Step 6 plus 10 pormds (44.5 Newtons) hung from knee axis ”\IGMbWN-‘O Step 7 plus 10 pounds (44.5 Newtons) hung from knee axis Table 4. Loads at the H-point and knee axes of the SAE 1826 manikin for the loading steps of Table 3. Load Step Load on Load on Moment H-point Knee about H- ax'n axis point Newtons Newtons N-mm 0 0 0 0 l 170 27 12300 2 245 28 12800 3 322 28 12800 4 395 32 14600 5 474 31 14200 6 501 70 32000 7 505 113 51600 8 506 158 72200 173 C. Results and Discussion: C1. Measurements of deflections due to seat pan loading with the ASPECT butt-thigh segment (ABT)- Table 5 is the raw data for the deflection measuranents of ABT into four automotive seats: Chrysler L11 (L11), Neon, Saturn, and Chrysler Town and Country (TC). Height measurements were made relative to the flat, level floor. Similarly, Table 6 is the raw data for the deflection measurements of the SAE 1826 manikin in these same four automotive seats. Table 5. Raw data for height measurements of ABT H-point axis and front of ABT hip hip front hip Hip from hip hip front hip hip front Load Step LH LH LH Neon Neon Neon Satrn Satrn Satrn TC TC TC right left ri t lefi Right left right lefi mm Mm mm mm Mm mm mm mm mm mm mm mm 0 445 445 51 l 437 459 557 410 410 521 465 468 577 l 448 446 476 448 445 515 416 412 481 463 462 525 2 434 429 482 436 431 519 406 403 491 454 452 533 3 422 417 488 426 417 530 397 394 494 448 447 536 4 412 411 492 412 406 536 392 393 502 443 437 541 5 405 410 499 407 401 542 384 381 510 441 434 541 6 400 407 501 402 396 545 379 378 515 439 429 541 7 399 399 503 400 394 551 378 373 525 435 427 542 8 394 393 506 396 388 560 375 371 525 432 422 545 9 392 391 510 394 387 560 373 371 525 427 418 547 10 390 390 502 393 388 545 373 371 514 424 417 536 l 1 391 391 498 391 382 538 370 367 501 423 413 528 12 388 391 494 388 383 529 370 367 499 420 412 521 Table 6. Raw data for height measurements of SAE 1826 H-point axis and knee axis hip hip knee knee hip hip knee knee hip hip knee knee hip hip knee knee Load L11 L1! L11 L11 Neon Neon Neon Neon Satan Satrn Satra Satrn TC TC TC TC right left right left right nen right left right no right left right 1‘“ right k“ mm mm mm mm mm mm mm mm mm mm mm mm mm mm mm mm 0 445 445 511 511 456 441 556 556 407 411 521 521 460 462 573 573 1 430 431 517 515 429 427 559 556 400 396 522 519 456 450 561 555 2 425 429 520 519 421 415 566 560 395 390 532 527 455 442 562 554 3 419 420 526 526 413 412 568 566 388 386 535 530 441 438 564 562 4 410 414 536 532 408 403 574 570 383 381 536 532 430 426 567 563 5 403 407 538 537 402 397 573 570 381 377 537 533 419 418 570 569 6 404 402 530 529 400 396 566 565 379 375 532 527 415 413 558 554 7 403 403 521 518 400 397 554 551 380 375 523 517 414 415 538 534 8 403 403 503 503 399 395 534 530 377 374 504 499 415 416 518 516 174 For the measurements using the ASPECT manikin, the seat pan loading of Table 2 can be combined with seat pan deflection measurements of Table 5 to provide load0deflection responses. Similarly, with measurements using the SAE 1826 manikin, the seat pan loading of Table 4 can be combined with seat pan deflection measurements of Table 6 to provide load-deflection responses. Radcliffe [6] is working to interpret these measurements in terms of: 1. vertical translational deflections of the H—point axis in response to downward loading tlnrough the H- point axis and 2. rotational deflections of the line fi'om the H—point axis through the knee axis. By fitting these translational and rotational respornses with a mathematical model, the response of a seat pan can be characterized so that loadings other that those used for the measurements can be predicted. Once these models are available, they will serve as Seat Factors. Figures 4 though 7 show the measured responses of fan automotive seats (Chrysler LH, Chrysler Neon, GM Satin-n, and Chrysler Town and Country) to loading of the 1826 and ASPECT manikins referernced to the 1826 template. For each seat, the first plot (a) shows the force-deflection response of the H-point axis; the secornd plot (b) is the moment-angle plot relative to the H-point axis. Notice that the vertical scales are difl‘erent for the difl’ercnt seats. H-pt Force vs. H—pt Deflection AapectvaABJMSeatingManifina Fifi-mm -10 6% E O z 10- § 2. . +LHIBI' 30 +4.11.” 40 s x 50 60 Figure 4a. Force—deflection responses of a Chrysler LH seat with the ASPECT (ABT) and 1826 manikins 175 H-pt. WtuRelativeThighRotation Aapect vs. SAE 1826 Seating Manikin H-pt. Ital-mm 4.00 4.00 -2.00 0 20000 30000 40000 50000 60000 70000 80000 4 2.00 . 4.00 0.00 ) 8.00 (fa/K . +01.“ 10.00 12.00 RalaIva Tlnlfln Won (0 Race lav“ Weaning.) Figure 4b. Thigh moment-angle respornses of a Chrysler LH seat with the ASPECT (ABT) and 1826 manikins 176 H-ptFoscevaH-ptbeflection massajmmm Marco-m o 100 zoo 300 400 500 000 O K . - _s O :20 .2,30 -e—n~£n~nsar 340 ‘50 +n~m~um ‘ao N O Figure 5a. Force-deflection responses of a Chrysler Neorn seat with the ASPECT (ABT) and 1826 manikins H-peMamnnuRelativeThighRotaflon AspeavaSABJMSeatingManikina HUMM 10.0.».3!’ 88888 as». 88 10.00 12.00 14.00 Relative Thlgh Rotation (0 Knee lavatlna Rotation) (a...) a. 8 Figure 5b. Thigh moment-angle responses of a Chrysler Neorn seat with the ASPECT (ABT) and 1826 manikins . 177 3.9112”tuan AapectvaSAEJMSeatingManildna MM“) -10 100 200 300 400 500 0 \ +smasr iffy: -....- O _s O 8 ma M «mum B 8 Fignne 6a. Force-deflection respornses of a GM Saturn with the ASPECT (ABT) and 1826 manikins H—ptMomentvaRelativeThighRontion AapectveSAEJWSeatingManiltine III-n. worm) Memmnommm Iatatloanag.) A N 8 8 10.00 12.00 I Figure 6b. Thigh moment-angle responses of a GM Saturn seat with the ASPECT (ABT) and 1826 manikins 178 H-pt PorcevnH-ptDeflcction AspectvaSAEJMSeating Manflcina "ii-ml") 0 100 200 300 41!) 500 600 x L O / \ -..... \ -..—.... “Mlomlm 8 8 8 8 8 Figure 7a. Force-deflection responses of a Chrysler Town and Country seat with the ASPECT (ABT) and 1826 manikins H-fihionentnkelative'l'highkotation AapectuSAEJMSeatingManildna ”m0?“ .0 .N p .I'o .1. .6» 88 8 8 8 8 H-e Relative midi Rotation (0 Knee Mating Rotation) (deg) SIX) Figure 7b. Thigh moment-angle responses of a Chrysler Town and Country seat with the ASPECT (ABT) and 1826 manikins The force-deflection plots (a) are all similar. With the 1826 manikin, the respornses are nearly linear. With the ASPECT manikin, there typically is an initial offset, tlnen a region of softer respornse from approximately 100 to 200 Newtons, then the stifliness with the ASPECT and 1826 manikins are nearly the same from 300 to 400 Newtons. The initial offset and softer response with the ASPECT manikin is probably due to the diflerence in contour shape under the H-point axis. As mentioned above, the ASPECT manikin is more prominernt due to realistic pelvis shape of tlne human data. If a single seat pan stifliness were to be selected as a Seat Factor, we suggest using the stifiness measured with the ASPECT manikin between 300 and 400 Newtons on the H-point axis. However, Radcliffe [6] is 179 developing a mathematical modeling process that represents the complex and nonlinear responses of automotive seats to loading with both the ASPECT and 1826 manikins. The moment-angle plots (b) are also similar. Starting fiom the origin, which is the orientation of the 1826 template, the knee axis of the 1826 manikin (without added weights) was typically 2 degrees higher than the 1826 tennplate. The exception was with the Chrysler Town and Conmtry seat where the 1826 manikin knees wee slightly [owe than the tenplate. The Town and Conmtry seat had a crease in the midlinc of tlne seat pan, which was a problem for placement of the of tlne J826 template. As the H-point axis of the 1826 manikin was loaded without sigrnificant change in moment about the H- point axis, the lance axis of the 1826 manikin rotated upward a few degrees. After the initial loading of tlne H-point axis and with an increase in moment about the H-point axis, the lance axis rotated downward into the seat. With the ASPECT butt and thigh segment (ABT) in each seat, the float of the ABT was approximately 4 degrees lower than the knees of the 1826 manikin. As the H-point axis of the ABT was loaded without increased moment about the H-point axis, the front of the ABT rose and this upward rotation of the ABT was consistently greater than the corresponding rotation of the 1826 manikin. Again, tlnis difference is probably due to the more prominent contour of the ABT under the H-point axis; with H-point loading, the ABT moved down in the seats with H-point and the thighs rotated up more than with the 1826 manikin, which has a flatter contour than the ABT. With increasing moment about the H-point axes, the moment-angle responses of the seats were equivalent with both the ABT and 1826 manikin. For a single parameter to represent support of tlne thighs by a seat, we suggest the slopes of the moment-angle responses between 20,000 and 40,000 Newton-millimeters. However, as stated above, Radclifle [6] is developing a mathematical modeling process that represents the complex and nonlinear responses of automotive seats to loading with both the ASPECT and 1826 manikins. We believe that this development by Radclifl‘e is preferable to using simpler seat factors, such as force- deflection and moment-angle stimnesses. References: l. Hubbard, R., and C. Gedraitis, “Initial Measurenents and Interpretatiorns of Seat Factors for the ASPECT Program”, Soc. of Auto. Engin. (SAE) Paper Number 1999-01-0958, presented at the 1999 SAE - International Corngress, March, 1999. 2. Schneider, L., M. Reed, R. Roe, M. Manary, C. Flannagan, R. Hubbard, and G. Rnnpp, “ASPECT: The Next Generation H-point Machine and Related Seat Desigrn and Measurement Tools”, Soc. of Auto. Engin. (SAE) Paper Number 1999-01-0962, presented at the 1999 SAE lntenational Congress, March, 1999. 3. “Devices for Use in Defining and Measuring Vehicle Seating Accommodation”, Society of Automotive Engineers Recommended Practice 1826. 4. Bogard, H., “Additiornal Methods to Characterize Lumbar Supports Using the ASPECT Manikin”, thesis for the Bachelor of Science Degree, Department of Materials Science and Mechanies, Michigan State University, 2000. 5. Gutowski, P., “Influence of Automotive Scat Factors on Posture and Applicability to Designn Models”, thesis for the Master of Science Degree, Department of Materials Science and Mechanies, Michigan State University, 2000. 6. Radcliffe, N., thesis in preparatiorn for the Master of Science Degree, Department of Materials Science and Mechanics, Michigan State University, 2001. 180 APPENDIX C 181 APPENDIX C This section describes and derives a secondary method for determining the ETSS in Models II and 111, which was not investigated fully due to the efi‘ectiveness of Model IV in modeling the interaction between seating manikins and automotive seats. _ Initial Position — Positiorn alter H-pt. Loading Figure C-l. Descriptive two-dimensional figure of H—pt-knee deflection and rotation for analysis of the ETSS. ZF,‘=0=F,, +FK =R,, +R, ZM” =0=R, *ETSL-FK *LK E ISL K geomeny=AT=AH+(AK—AH)* ETSS=KT=£= FH+FK_RH 3 FK*LK2 AH+(AK_AH).§ZS_£ I...2 *(AK—AH)+LK *LT‘AH K 182 REFERENCES . Hubbard, R., and C. Gedraitis, “Initial Measurements and Interpretations of Seat Factors for the ASPECT Program,” Society of Automotive Engineers (SAE) Paper Nmnber 1999-01-0958, presented at the 1999 SAE International Congress, Marcln, 1999. . Working Model® 2D Software. Knowledge Revolution. Version 4.0 ©1992-1996 . Excel Software. Microsoft Corporation. Ofice 2000 Version . “Devices for Use in Defining and Measuring Vehicle Seating Accommodation”, Society of Automotive Engineers Recommended Practice 1826. . Schneider, L., M. Reed, R. Roe, M. Manary, C. Flannagan, R. Hubbard, and G. Rupp, “ASPECT: The Next Generation H-point Machine and Related Seat Desigrn and Measurement Tools”, Soc. of Auto. Engin. (SAE) Paper Number 1999-01-0962, presented at the 1999 SAE International Congress, March, 1999. . Hubbard, R., and N. Radclifi‘c, “Seat Pan Load-Deflection Measurements For Determination of Seat Factors with the ASPECT Butt-Thigh (ABT) Segment of the ASPECT Manflcin Prototype at MSU and the SAE 1826 Manikin,” submitted to TecMath Inc., reprinted in Appendix B. . Reid-Bush, T. “Posture and Force Measures of Mid-Sized Men in Seated Postures,” dissertation for the degree of PhD. Dept. of Materials Science and Mechanics Michigan State University, 2000. 183 IIImamwiljfljfluijfllfljnglrI