31293 01706 9901 This is to certify that the thesis entitled "A Rheological Model Of Apple Flesh Failure During The Magness-Taylor Puncture Test" presented by Sanghyup Jeong has been accepted towards fulfillment of the requirements for M.S. degree in Agricultural Engineering Major professor Date “.éi, / C?“ 47 0 76,710 MSU is un Altirmutim’ Action/Equal Opportunity Institution LIBRARY Mlchlgan State Unlvorslty PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. DATE DUE DATE DUE DATE DUE am my 1/98 c-JCIRC/DateDueipes-p.“ ~..-— ~.-....._ -4...‘ __ _ ,..._ ’h‘l A RHEOLOGICAL MODEL OF APPLE FLESH FAILURE DURING THE MAGNESS-TAYLOR PUNCTURE TEST By Sanghyup Jeong A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Agricultural Engineering 1997 ABSTRACT A RHEOLOGICAL MODEL OF APPLE FLESH FAILURE DURING THE MAGNESS-TAYLOR PUNCTURE TEST By Sanghyup J eong In an attempt to assess firmness, rheological behavior of apple flesh was studied by developing a laminated cell layer model which can be used to emulate Magness-Taylor puncture test. Each layer of the model consisted of a compressive and a shear-tensile element. Compressive component consisted of a Maxwell model with a fracture element. The shear-tensile component consisted of an elastic element. The annular pumping effect was also included in model. Layers were connected in series and each shear-tensile component was grounded. Simulation was performed on this model with element parameters and the postulated deformation sequence. The model depicted fairly well the Magness-Taylor test results both quantitatively as well as qualitatively. Nondestructive techniques using Hertz contact stress theory and acoustic impulse technique (Armstrong, 1990) were used for estimating model parameters. These methods were fairly accurate in predicting the model parameters. Also, a strong correlation with the Magness-Taylor firmness reading was found. ACKNOWLEDGMENTS I would like to thank Dr. Aj it K. Srivastava for his direction, support, encouragement, and patience. I also would like to thank the other members of my committee, Dr. Dan Guyer and Dr. Randolph M. Beaudry for their suggestions and assistance. I thank Mr. Dick Wolthius for his technical advice. I also give my thanks to my parents for their financial and spiritual support for me. All of the Paul Harris fellows of the Rotary Foundation deserve my heartfelt thanks for their financial support for student like me. I would like to thank all the colleagues of graduate office 9. iii TABLE OF CONTENTS ACKNOWLEDGMENT ii: LIST OF FIGURES . v: CHAPTER 1 INTRODUCTION 1.1 Justification 1.2 Objectives CHAPTER 2 REVIEW OF LITERATURE ‘ 2.1 Rheological Models Representing the Textural Characteristics of Food Materials 4 2.2 Physiological Considerations of Apple ( CHAPTER 3 THEORY i 3.1 Modeling Principles and General Assumptions 2 3.2 Modeling Compression Component 2 3.3 Modeling Shear-Tensile Component 15 3.4 Modeling Annular Pumping Component 7 3.5 Postulated Deformation Sequence 7] 3.6 Compression, Shear—Tensile, and Annular Pumping Components allied Model 24 3.7 Parameter Determination Using Rubber “Thumb Test” ................................. 26 3.8 Parameter Determination Using Acoustic Impulse Technique ...................... 3( CHAPTER 4 EXPERIMENTAL METHODS AND PROCEDURES 3] 4.1 Methods 3] 4.1.1 Compression Test 3] 4.1.2 Shear Test 34 4.1.3 Tensile Test 34 4.1.4 Magness-Taylor Puncture Test 34 4.1.5 Nondestructive Tests 4.1.5.1 Contact Pressure Measurement Using Rubber “Thumb Test” 38 4.1.5.2 Elastic Modulus Measurement Using the Acoustic Impulse Response 41 4.2 Materials 44 4.3 Computer Simulation of Laminated Cellular Model ...................................... 44 CHAPTER 5 RESULTS AND DISCUSSIONS 47 5.1 Results 47 5.1.1 Model Prediction and Magness-Taylor Puncture Test ..................... 47 5.1.2 Performance of Rubber “Thumb Test” 53 5.1.3 Performance of Acoustic Impulse Technique 59 5.2 Discussion 61 5.2.1 Quantitative Aspects of Model Prediction 61 5.2.2 Qualitative Aspects of Model Prediction 6] 5.2.3 Evaluation of Nondestructive Methods... .. 63 CONCLUSIONS 6% SUGGESTIONS FOR FUTURE RESEARCH ..66 APPENDICES ......68 LIST OF REFERENCES 7m LIST OF FIGURES FIGURE 1 - Force-deformation relationship of compression test for cylindrical apple flesh ............................................................ 9 FIGURE 2 - Compression model consisted of a Maxwell and a fracture element 1 1 FIGURE 3 - Modeling concept of cylindrical apple flesh sample ................................... 12 FIGURE 4 - Modeling concept of shear-tensile component 16 FIGURE 5 - Annular pumping effect during mastication 18 FIGURE 6 - Model representing annular effect 19 FIGURE 7 - Postulated deformation sequence of Magness-Taylor puncture test .......... 23 FIGURE 8 - Final model of apple flesh for Magness-Taylor puncture test ..................... 25 FIGURE 9 - Hertz contact stress of two spheres 77 FIGURE 10 - Stress distribution of Hertz contact problem 79 FIGURE 11 - Sampling direction apple flesh for compression test ................................. 32 FIGURE 12 - Cutting guide for preparing cylindrical apple flesh ................................... 33 FIGURE 13 - Apparatus for shear test 35 FIGURE 14 - Specimen for tensile test 36 FIGURE 15 - Apple holding apparatus 37 FIGURE 16 - Location of the Magness—Taylor puncture test 39 vi FIGURE 17 - Rubber thumb test 40 FIGURE 18 - Device for measuring the radius of curvature (Radius of Curvature Meter) ..................... 42 FIGURE 19 - Device for measuring acoustic response of apple .......................... FIGURE 20 - Structure of the simulation program ............ 43 45 FIGURE 21 - Predicted versus measured first peak force FIGURE 22 - Predicted versus measured second peak force .48 49 FIGURE 23 - Predicted versus measured stiffness constant FIGURE 24 - Comparison of model prediction with experiment ......................... FIGURE 25 - Predicted versus measured Magness-Taylor firmness FIGURE 26 - Predicted versus measured compressive failure stress ................... FIGURE 27 - Predicted versus measured shear failure stress FIGURE 28 - Magness-Taylor firmness versus failure force of rubber thumb test FIGURE 29 - First peak force versus failure force .................................... FIGURE 30 - Measured first peak force versus Magness-Taylor firmness FIGURE 31 - Measured versus predicted modulus of elasticity by using acoustic impulse FIGURE 32 - Precision of model prediction 50 ........... 51 .52 ............ 54 55 56 ........... 57 .58 60 6? Chapter 1 INTRODUCTION 1.1 Justification Fresh market apples are sorted by color, size, weight and fiminess. Apples classed by their physical properties and texture are processed for the various apple products. The texture of apples is an important quality attribute. Apple firmness is an indicator of textural qualities. Boume (1982) defined fimmess as a texture characteristic during mastication that displays moderate resistance to breaking. The Magness-Taylor pressure test developed in 1925 to measure firmness, has been accepted as an industry standard. The Magness-Taylor pressure test consists of forcing a cylindrical metal tip into the apple flesh without skin and measuring the maximum force. The penetration depth is about 8 mm (5/16 in) and the diameter ofthe tip is approximately 11 mm (7/16 in) or 8 mm (5/16 in). The plunger tip is round and the penetration depth is marked on the plunger. The destructive nature of the Magness-Taylor pressure test results in fruit loss which contributes to economic loss for the apple industry. The firmness assessment is based on sampling. Export markets require a better quality assessment of the product. A non- destructive technique that can make it possible to test individual fruit would be beneficial. The Magness-Taylor method relies on only one indicator to describe the textural quality of the apple, the maximum force to masticate apple flesh. Other parameters such as elastic modulus also contribute to texture. 2 An in-depth understanding of the Magness—Taylor test would be very helpful in developing a nondestructive technique assessing fruit firmness to meet the needs of high quality produce. 1.2 Objectives The objectives of the research were 1. to develop a rheological model of apple flesh and to examine the capability of the model for predicting the Magness-Taylor puncture test result. 2. to evaluate the potential of a rubber spherical indentor as a means of estimating the compression and shear failure stresses of apple flesh. 3. to investigate the capability of the spherical indentor as a non-destructive technique to measure apple firmness. Chapter 2 REVIEW OF LITERATURE 2.1 Rheological Models Representing the Textural Characteristics of Food Materials Much effort have been devoted to developing rheological models to describe complex texture qualities of food materials. New rheological elements have been suggested along with a combination of the conventional rheological elements e. g. spring and dashpot. These models have been subjected to different types of loading conditions and their responses have been compared to the experimental results. Maxwell and Kelvin models which are the classical models in rheology have already been studied and evaluated under compression, extension, creep, and relaxation tests. Maxwell and Kelvin models are composed of a spring and a dashpot in a serial and in a parallel pattern, respectively. Peleg (1976) proposed a model system which was basically a parallel array of generalized Maxwell body. He added two fracture elements which were suggested by Drake (1971) to the model. The general elements of this model had an elastic stage and a Maxwellian stage when the model worked as a Maxwell model until fracture. The number of each stage and the parameters of the elements determine the pattern of rheological behavior of a given material. This model can describe various rheological behavior qualitatively. Chen and Rosenberg (1977) suggested a nonlinear viscoelastic model which consisted 5 of Maxwell and Kelvin models containing an yield element. The model exhibited Maxwellian behavior and then it showed the behavior of a four element Burgers model which was a serial combination of a Maxwell and a Kelvin model afier the yield element was activated. He indicated that the model was able to represent the behavior of American cheese by adjusting model parameters. Dickinson and Goulding (1980) tested various types of cheese and found that the model suggested by Chen and Rosenberg (1977) could be used to fit individual compression curve but the model was unable to predict general load deformation behavior. Johnson and et a1. (1981) modified the model suggested by Peleg (1976) by incorporation of a contact and a fracture element. In this study the contact and the fracture element were functions of time and temperature to describe the rheological behavior of raw and cooked fish meat. They demonstrated that the model worked well to describe the rheological behavior qualitatively. McLaughlin (1987) had a different approach to develop a model describing apple tissue failure phenomenon. Yield stress data for 800 cylindrical samples was studied. Statistical test failed to reject (p=0.85) the hypothesis that the data followed a three parameter Weibull distribution. This statistical model was suggested as an adequate probability model for predicting the failure stress of apple flesh. During the past decade most of the efforts were concentrated on any developing property measurement techniques rather than developing a model having any physical meaning. 2.2 Rheological Considerations of Apple Much research has been conducted about the phenomenological and morphological characteristics of the failure of fruits. Mechanical and structural characteristics of fruits were observed and discussed. An understanding the nature of failure will help develop appropriate models having physical meaning. Peleg and et al. (1976) studied compressive failure patterns of several juicy fruits. Compression tests were performed to the cylindrical samples by using an Instron Universal Testing Machine. Fruits showed their characteristic failure patterns. Effect of the juiciness and dimension of the sample were studied. In this research it was found that structural characteristics, liquid content, and geometry of the sample affect the force— deformation relationship. Mohsenin (1977) studied the failure of food materials from the point view of solid mechanics. Working with apple tissues he found that the shear resistance was a result of the cementing agent (pectic substance) in the middle lamella holding together the cells in the parenchyma tissue. He also studied the potential usage of the Hertz contact stress theory to find failure parameters from the force-deformation relationship. Pitt (1982) considered two different types of failure, cell wall rupture and cell debonding. Sample was considered as the sum of a large number of cell layers normal to the direction of the applied force. In case of constant rate strain due to the rupture of the weakest cell of a layer, neighboring cells had higher stress level than before. If the stress is higher than cell strength then failure occurs and this process propagates through the layers of the sample. Due to this fractured transverse plane enzymatic browning occurred 7 in that cross section. This catastrophic failure was observed in the fresh apple tissue. Lin and Pitt (1986) reported that turgor pressure and the strain rate affect failure mode of fruit and vegetables and tissue strength. Khan and Vincent (1990) observed that cells of apple flesh are arranged in radial column and the cell size and shape changes along with the radial direction by using scanning electron microscope, indicating that apple flesh should be considered as an anisotropic material in the radial direction. Khan and Vincent (1993) observed that if a radial cylindrical apple flesh sample was subjected to a compressive force then it generally fractured by collapse of a single layer of cells at right angle to the direction of force. Tangential sample failed in shear. This phenomenon is due to the anisotropic property of apple parenchyma which was observed by Khan and Vincent (1990). Chapter 3 THEORY 3.1 Modeling Principles and General Assumptions To develop a model representing Magness-Taylor puncture test it is necessary to observe and analyze the mastication phenomenon. Basically Magness—Taylor puncture test involves three major factors which are compression, shear-tension and annular pumping effects. These factors will be combined into one system representing the apple flesh afier being studied and modeled individually. Throughout this research several simplifying assumptions were made. Even though the actual Magness-Taylor firmness instrument has round tip the radius of curvature is large enough so that the plunger tip was considered as flat. The second assumption was about the cylindrical apple flesh sample. Actually the layer of the sample is curved because of its radial arrangement of cells. However, the curvature was relatively large enough to be treated as flat because the sample diameter was considerably less than that of the whole apple. 3.2 Modeling Compression Component If several ideal conditions of the compression test are satisfied the saw-tooth shaped force deformation relationship as shown in Figure 1 could be obtained. The first ideal condition is to make perfect cylindrical apple flesh sample which means the flat cutting FORCE [N] 140 120 ~— 100 — 804 60« 40 -- 20- P1 <~ P2 0. 0 0.5 1.0 1.5 2.0 2.5 DEFORMATION [mm] FIGURE 1 - Force-deformation relationship of compression test for cylindrical apple flesh 3.0 10 surface are normal to the direction of the compressive force. The second ideal condition is lack of friction between cutting surface and the pressing die. The last ideal condition is satisfied if the sample is obtained in the radial direction from the whole apple because of the anisotropic characteristics of apple parenchyma studied by Khan and Vincent (1990). In these ideal conditions sample shows two or three fractures caused by collapse of a single layer of cells at right angles to the force (Figure 1). A simple model representing this phenomenon was developed. This model was the building block of the final model. The first compression region (b, Figure l) was represented by a Maxwell model and a fracture element having 6c as its yield stress (Figure 2). The model was used to represent a single layer of cells. Mathematical expression for the model at constant strain rate is, E o=nV(1—e rl) , OSoom - 2 mMDOE 7 m4‘ Total Energy <—-—— ‘ / - 7"-\\\ leferentiation/- Total Force 1 Force from Annular Pumping Final Force-Deformation Relationship FIGURE 20 - Structure of the simulation program 46 seconds of simulation, numerical differentiation of energy for the deformation was performed and it produced force-deformation relationship. To the force-deformation result, annular pumping resistance was added gradually after the second peak force. Chapter 5 RESULTS AND DISCUSSION 5.1 Results 5.1.1 Model Prediction And Magness-Taylor Puncture Test From the force-deformation relationship of the Magness-Taylor puncture test the first peak force was determined at a point showing significant drOp of force. Figure 21 illustrates the result of the comparison between predicted and measured first peak force of Magness-Taylor test. The correlation coefficient was 0.76. The measured second peak force of Magness-Taylor puncture test was compared with that of model prediction. Figure 22 shows the correlation coefficient of 0.79 which is a little higher than that of the first peak force shown in Figure 21. In some cases, the second peak force could not be observed and those cases were not considered. Even though the first peak forces are same different stiffness indicates different texture of apple flesh. Figure 23 illustrates that the correlation coefficient was 0.83. Stiffness constant was defined as the slope of the secant drawn from the origin to the first peak of the force-deformation curve. Figure 24 is an example of the simulation result. Magness-Taylor firmness predicted by the model was compared with that measured by experiment. Figure 25 shows that the correlation coefficient is 0.82. Even though the juice expelling effect produced constant resistance of 9.67 N which was calculated by using the eq. [1 6] it took time until this resistance was fully developed. 47 Measured First peak of Magness-Taylor test [N] 120 100 80 60 40 20 48 Y = 1.02 X + 4.00 Correlation Coefficient=0.76 Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold —— Fitted 20 40 60 80 100 Predicted Flrst peak [N] FIGURE 21 - Predicted versus measured first peak force Measured 2nd Peak [N] 100 90 80 70 60 5O 40 30 20 1O 49 =0.83X+10.67 Correlation Coefficient=0.79 1/ . Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold — Fitted 0 20 40 60 80 100 120 Predicted 2nd Peak [N] FIGURE 22 - Predicted versus measured second peak force Measured Stiffness Constant [Nlm] 60000 50000 40000 30000 20000 11 0000 50 Y=0.83X+3657 Correlation Coefficient=0.83 1 Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold —- Fitted 0 10000 20000 30000 40000 50000 Predicted Stiffness Constant [Nlm] 60000 FIGURE 23 - Predicted versus measured stiffiiess constant r.- FORCE[N] 51 WWfiimmv , wv,,flfi,, , .7 70 LL ,L,7,L, 7, , 2nd Peak Force mT~M%“AAAA\> Force , 1 30~ Stiffness 20 I " Constant — Experiment —Model Prediction 0 1 2 3 4 DEFORMATION [mm] FIGURE 24 - Comparision of model prediction with experiment "I Model Predicted Magness Taylor Firmness [N] 120 100 80 60 40 20 Y = 0.7 X + 15.78 Correlation Coefficient=0.82 ii Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold Measured Magness Taylor Firmness [N] — Fitted FIGURE 25 - Predicted versus measured Magness Taylor F irrnness "lf 53 The resistance increased linearly after the first peak until it reached a value of 9.67 N and then it was fixed at 9.67 N. Because the viscous element was connected in parallel to the model the resistance was just added to the result. 5.1.2 Performance of Rubber “Thumb Test” The maximum normal stress at the center of the contact area between rubber ball and the apple was calculated by using eq. [23]. Speed of the rubber thumb test was 0.00212 m/s (5 in/min). Figure 26 shows the correlation coefficient of 0.76. Figure 27 illustrates the correlation coefficient of 0.64 between measured shear failure stress and the predicted maximum shear stress. According to the Hertz contact stress theory the maximum shear stress occurred below the contact surface and it was approximately 1/ 3 of the maximum normal stress. Magness-Taylor firmness was compared directly with the failure force which was observed when the contact pressure dropped significantly. In this research this moment of “give” was detected and the testing machine was stopped at that moment. Figure 28 shows a correlation coefficient of 0.9. In rare cases, no significant drop of the pressure was observed due to tissue failure. The first peak force showed a much higher correlation coefficient of 0.90 with the failure force than that of the Magness-Taylor firmness (Figure 29). The high correlation coefficient implies that the failure force represents the first peak force. The relationship between the first peak force and the Magness-Taylor firmness was observed in Figure 30. More than a half of the data is almost identical in the range of error. 54 Y = 0.55 X + 35249 Correlation Coefficient=0.76 5E+05 4E+05 . w l 9:. i a 4E+05 ' E 1 ‘3 35+05 2 .2 35+05 1 1 . o . I .2 . a . 9:9 2E+05 a g . 8 ZE+05 1 B * 1 0 ‘5 1 E+05 2 m 5E+04 0E+00 Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold —— Fitted 0E+00 1E+05 2E+05 3E+05 4E+05 5E+05 6E+05 7E+05 Measured Compressive Failure Stress [Pa] FIGURE 26 - Predicted versus measured compressive failure stress Predicted Shear Failure Stress [Pa] 1 E+05 1E+05 1 E+05 8E+04 6E+04 4E+04 2E+04 0E+00 FIGURE 27 - Predicted versus measured shear failure stress 55 Y = 0.2 X + 40408 Correlation Coefficient=0.64 1 i l Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold — Fitted 1 Measured Shear Failure Stress [Pa] OE+00 5E+04 1E+05 2E+05 2E+05 3E+05 3E+05 4E+05 56 Y = 1.81 X-58.9 Correlation Coefficient=0.77 250 0 Red Delicious 1 1 I Red 1 Rome 1 A Fuji 1 x Granny Smith 150 -~ x Golden Delicious . e Braeburn: 100 «— + Ida Red Magness-Taylor Firmness [N] - Gala 50 —» - Jona Gold — Fitted Failure Force [N] FIGURE 28 - Magness-Taylor firmness versus failure force of rubber thumb test .«rr Measured Failure Force [N] 250 200 150 100 50 57 Y = 2.44 X - 78.41 Correlation Coefficient=0.90 20 40 60 80 1 00 11- Measured 1st Peak [N] FIGURE 29— First peak force versus Filure Force 0 Red Delicious I Red Rome A Fuji x Granny Smith at Golden Delicious 0 Braebum + Ida Red - Gala - Jona Gold 4— Fitted lL 58 Y = 0.67 X + 12.83 Correlation Coefficient=0.78 1 . e 100 1 O z . § 80 e u x (I R .. 60 .2 IL 3 a 401— __ fl 0 E 20 0 0 20 40 60 80 1 00 Measured Magness Taylor Firmness [N] FIGURE 30 - Measured first peak force versus Magness-Taylor Firmness qr 120 Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold —— Fitted 59 5.1.3 Performance of Acoustic Impulse Technique Modulus of elasticity of apple tissue was predicted using eq. [27]. In the compression mode, a Poisson’s ratio of 0.3 and the density of 800 kg/m3 were assumed for the resonant sphere model. Measured modulus of elasticity was obtained from the compression test of cylindrical apple flesh. From the compression test results initial tangent modulus of elasticity was found. In the compression test cylindrical apple flesh sample was bored in the radial direction which was different from the sampling direction used by Armstrong (1989). Figure 31 illustrates 0.7 as the correlation coefficient between the measured and predicted modulii of elasticity. IL Predicted Modulus of Elasticity [MPa] 60 Y = 0.0027 X + 2.43 Correlation Coefficient=0.7 ___ ____1 944% 1 ! 200 400 600 800 1000 Measured Modulus of Elasticity [MPa] FIGURE 31 - Measured versus predicted modulus of elasticity by using acoustic impulse Red Delicious Red Rome Fuji Granny Smith Golden Delicious Braebum Ida Red Gala Jona Gold — Fitted 1200 I P 61 5.2 DISCUSSION 5.2.1 Quantitative Aspects of The Model Prediction The laminated cell layer model showed fairly good correlation for the first peak force and the second peak force prediction. The correlation coefficient of the second peak (0.79) was little higher than that of the first peak force (0.76). Thus, the model demonstrated its prediction capability for the first two peak forces. In addition, Figure 23 showed a strong correlation (0.83) in predicting stiffness constant which meant that the model worked not just for the peak forces but also for the first slope of the force- deformation relationship. After the second peak it was hard to compare because unidentified disorder occurred as the mastication progressed. However, the information about the early period of the mastication period was enough to determine the texture of apple flesh. In Figure 32, it was found that the percentages of the sample having less than 20% error with its measured values were 51% for the first peak prediction, 78% for the second peak prediction, 73% for the stiffness prediction, and 84% for the Magness-Taylor firmness prediction. From the results explained above the laminated cell layer model was considered as a good means to depict the Magness-Taylor puncture test quantitatively. 5.2.2 Qualitative Aspects of The Model Prediction The characteristics of the force-deformation relationship was a saw-tooth shape which inferred the significant change of force or the catastrophic destruction of material. The model represented well this characteristics of the force-deformation relationship. Figure 23 shows the correlation coefficient of 0.83 for predicting the stiffness which indicates 62 ezerror AI 1st Peak Force 2nd Peak Force “10% 920% 19% 22% o e<10% c2336 45% 0 10<¢<20% io20% e>20% 16% 27% e MT Area 6.91E-06 [mA2] Velocity 0.004233 [m/s] Yield stress 1145150 [Pa] Predicted: Yield Strength 0.79 [N] Force [N] \ 0 1 2 3 4 5 Defamation [mm] FIGURE D1 - Shear test result for Y01_04 166 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00241 [m] MT thick 0.0002 [m] Sample area 0.000132 [m"2] —> MT Area 6.91 E-06 [mA2] Velocity 0.004233 [mls] Yield stress 299244 [Pa] Predicted: Yield Strength 2.07 [N] 60 50 //\\ 4o / \ E 3 30 E / \ 2. // \\ M MM, \NWM~ 1o 0 1 2 3 4 5 Defamation [mm] FIGURE D2 - Shear test result for Y02_04 167 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00137 [m] MT thick 0.0002 [m] Sample area 7.53E-05 [m"2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [mls] Yield stress 185405 [Pa] Predicted: Yield Strengt 1.28 [N] 20 / 15 g \ 3 10 — :«2 ’ \M W! 5 M o 0 1 2 3 4 5 Defamation [mm] FIGURE D3 - Shear test result for Y03__04 168 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0021 [m] MT thick 0.0002 [m] Sample area 0.000115 [m"2] —> MT Area 6.91E-06 [m‘2] Velocity 0.004233 [mls] Yield stress 248737 [Pa] Predicted: Yield Strength 1.72 [N] Farce [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D4 - Shear test result for Y04_04 _. 1.5 169 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00197 [m] MT thick 0.0002 [m] Sample area 0.000108 [m"2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [m/s] Yield stress 228580 [Pa] Predicted: Yield Strength 1.58 [N] 40 35 . / \ .. / \ .. / \ Force [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D5 - Shear test result for Y05_04 53E 170 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00193 [m] MT thick 0.0002 [m] Sample area 0.000106 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 191686 [Pa] Predicted: Yield Strength 1.32 [N] 30 25 n 2. / Farce [N] a: / 10 0 1 2 3 4 5 Defamation [mm] FIGURE D6 - Shear test result for YO6_O4 171 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00188 [m] MT thick 0.0002 [m] Sample area 0.000103 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 219752 [Pa] Predicted: Yield Strength 1.52 [N] Force [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D7 - Shear test result for Y07_04 172 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00194 [m] MT thick 0.0002 [m] Sample area 0.000107 [mA2] » MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 181504 [Pa] Predicted: Yield Strength 1.25 [N] 25 20 15 / WK HM. Force [N] 8 Defamation [mm] FIGURE D8 - Shear test result for Y08_04 173 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00164 [m] MT thick 0.0002 [m] Sample area 9.01E-05 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [m/s] Yield stress 176118 [Pa] Predicted: Yield Strength 1.22 [N] 30 . A .. l 1 Force [N] 6‘. in O 1 2 3 4 5 Defamation [mm] FIGURE D9 - Shear test result for Y09_04 174 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0022 [m] MT thick 0.0002 [m] Sample area 0.000121 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 248325 [Pa] Predicted: Yield Strength 1.72 [N] 35 30 25 / \ _ 20 / E 8 / ‘6 “' 15 I \NV/VR 1o / W V‘ ‘vw. m/ 5 0 1 2 3 4 5 Defamation [mm] FIGURE D10 - Shear test result for Y10__04 175 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00173 [m] MT thick 0.0002 [m] Sample area 9.51 E-05 [m"2] —> MT Area 6.91E-06 [m02] Velocity 0.004233 [mls] Yield stress 312775 [Pa] Predicted: Yield Strength 2.16 [N] 30 25 \ 20 Force [N] a: \ WM 0 1 2 3 4 5 Defamation [mm] FIGURE D11 - Shear test result for Y11_04 Fr 176 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00256 [m] MT thick 0.0002 [m] Sample area 0.000141 [m"2] —> MT Area 6.91E-06 [m‘2] Velocity 0.004233 [mls] Yield stress 224689 [Pa] Predicted: Yield Strength 1.55 [N] 35 30 A 25 20 Force [N] a: \- / Defamation [mm] FIGURE D12 - Shear test result for Y12_O4 177 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00277 [m] MT thick 0.0002 [m] Sample area 0.000152 [m"2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [mls] Yield stress 242944 [Pa] Predicted: Yield Strength 1.68 [N] 40 35 /\ 30 25 Force [N] N O / ’ \ \wf\ 0 1 2 3 4 5 Defamation [mm] FIGURE D13 - Shear test result for Y13_04 178 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0028 [m] MT thick 0.0002 [m] Sample area 0.000154 [m"2] ——> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 201034 [Pa] Predicted: Yield Strength 1.39 [N] 3., /A\ 2. /\ .. \ Force [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D14 - Shear test result for Y14_04 179 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0029 [m] MT thick 0.0002 [m] Sample area 0.000159 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 212239 [Pa] Predicted: Yield Strength 1.47 [N] Force [N] 40 35 25 1 2 3 4 5 Defamation [mm] FIGURE D15 - Shear test result for Y15_04 180 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00268 [m] MT thick 0.0002 [m] Sample area 0.000147 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 300421 [Pa] Predicted: Yield Strength 2.08 [N] Force [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D16 - Shear test result for Y16_04 181 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0026 [m] MT thick 0.0002 [m] Sample area 0.000143 [m"2] “—> MT Area 6.91E-06 [mA2] Velocity 0.004233 [m/s] Yield stress 189325 [Pa] Predicted: Yield Strength 1.31 [N] 30 . /\ .. J \ 1 10 t 5 W] \Mfmmmv Defamation [mm] FIGURE D17 - Shear test result for Y17_04 AI" 182 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00384 [m] MT thick 0.0002 [m] Sample area 0.000211 [m"2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [m/s] Yield stress 235456 [Pa] Predicted: Yield Strength 1.63 [N] 60 50 /"\ 30 Force [N] .., \ / \ 0 1 2 3 4 5 MWVv/W‘ Defamation [mm] FIGURE D18 - Shear test result for Yl 8_04 183 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.003 [m] MT thick 0.0002 [m] Sample area 0.000165 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [m/s] Yield stress 160527 [Pa] Predicted: Yield Strength 1.11 [N] 3.00E+01 2.5OE+01 A /\1 2.00E+01 1.50E+01 Farce [N] 1.00E+01 5.00E+00 \‘ 0.00E+OO l O 1 2 3 4 5 Defamation [mm] FIGURE D19 - Shear test result for Y19_04 184 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00221 [m] MT thick 0.0002 [m] Sample area 0.000121 [m42] ~> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 288385 [Pa] Predicted: Yield Strength 1.99 [N] 40 35 .. A .. /\ Force [N] / \Www 0 WWWWA 0 1 2 3 4 5 Defamation [mm] FIGURE D20 - Shear test result for Y20_04 ..J'I. 185 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00306 [m] MT thick 0.0002 [m] Sample area 0.000168 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 102475 [Pa] Predicted: Yield Strength 0.71 [N] 25 ! 20 l 1 1 Force [N] P/j WWWNMWl 0 1 2 3 4 5 Defamation [mm] FIGURE D21 - Shear test result for Y21_04 EAL 186 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00336 [m] MT thick 0.0002 [m] Sample area 0.000185 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 145849 [Pa] Predicted: Yield Strength 1.01 [N] 2. A /\ Force [N] a: Defamation [mm] FIGURE D22 - Shear test result for Y22_04 187 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00302 [m] MT thick 0.0002 [m] Sample area 0.000166 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 130190 [Pa] Predicted: Yield Strength 0.90 [N] 2.501901 2.00E+01 /\ 1.50E+01 Farce [N] / 1 DOE-1'01 I \ 5.005000 1 \ 0.00E+00 , 0 1 2 3 4 5 Defamation [mm] FIGURE D23 - Shear test result for Y23_O4 188 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00292 [m] MT thick 0.0002 [m] Sample area 0.00016 [mAZ] » MT Area 6.91E-06 [m42] Velocity 0.004233 [m/s] Yield stress 89790 [Pa] Predicted: Yield Strength 0.62 [N] 15 1o 2 WMMMWW 0 1 2 3 4 5 Defamation [mm] FIGURE D24 - Shear test result for Y24_04 189 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00376 [m] MT thick 0.0002 [m] Sample area 0.000207 [mA2] —> MT Area 6.91E—06 [m"2] Velocity 0.004233 [mls] Yield stress 112771 [Pa] Predicted: Yield Strength 0.78 [N] 25 ..i /1 Farce [N] Wme/Wwv 0 1 2 3 4 5 Defamation [mm] FIGURE D25 - Shear test result for Y25_04 ,. 190 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00248 [m] MT thick 0.0002 [m] Sample area 0.000136 [m"2] 9 MT Area 6.91E-06 [m‘2] Velocity 0.004233 [mls] Yield stress 211617 [Pa] Predicted: Yield Strength 1.46 [N] 30 . 1 25 /\ Force [N] a: \ i W 0 "‘~\, 0 1 2 3 4 5 Defamation [mm] FIGURE D26 - Shear test result for Y26_O4 at l i A 191 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00263 [m] MT thick 0.0002 [m] Sample area 0.000145 [m"2] » MT Area 6.91E-06 [mA2] Velocity 0.004233 [m/s] Yield stress 293539 [Pa] Predicted: Yield Strength 2.03 [N] Force [N] 0 1 2 3 4 5 Defamation [mm] FIGURE D27 - Shear test result for Y27_O4 192 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00263 [m] MT thick 0.0002 [m] Sample area 0.000145 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 200621 [Pa] Predicted: Yield Strength 1.39 [N] 30 \ Farce [N] a 1 ‘ l \Wl. Defamation [mm] FIGURE D28 - Shear test result for Y28_04 Fa 193 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00251 [m] MT thick 0.0002 [m] Sample area 0.000138 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [m/s] Yield stress 318823 [Pa] Predicted: Yield Strength 2.20 [N] Force [N] 1 2 3 4 5 Defamation [mm] FIGURE D29 - Shear test result for Y29_04 f. A. 194 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0024 [m] MT thick 0.0002 [m] Sample area 0.000132 [m‘2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 293341 [Pa] Predicted: Yield Strength 2.03 [N] 40 /\ . ::: 1 \ .. / l Force [N] Defamation [mm] FIGURE D30 - Shear test result for Y30_04 ?' 195 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00316 [m] MT thick 0.0002 [m] Sample area 0.000174 [m"2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [mls] Yield stress 177696 [Pa] Predicted: Yield Strength 1.23 [N] 35 30 A /\ 25 20 Farce [N] 6‘. \ / waWM—Mt 0 1 2 3 4 5 Defamation [mm] FIGURE D31 - Shear test result for Y31_04 196 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00318 [m] MT thick 0.0002 [m] Sample area 0.000175 [m"2] ——> MT Area 6.91E-06 [mAZ] Velocity 0.004233 [mls] Yield stress 120155 [Pa] Predicted: Yield Strength 0.83 [N] 25 A Farce [N] I 5 \R WVWV‘WwW" 0 0 1 2 3 4 5 Defamation [mm] FIGURE D32 - Shear test result for Y32_04 197 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00285 [m] MT thick 0.0002 [m] Sample area 0.000157 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 129266 [Pa] Predicted: Yield Strength 0.89 [N] 25 20 A Force [N] \Mtjkww- 0 1 2 3 4 5 Defamation [mm] FIGURE D33 - Shear test result for Y33_04 198 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00229 [m] MT thick 0.0002 [m] Sample area 0.000126 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 111836 [Pa] Predicted: Yield Strength 0.77 [N] 15 .o \ Z 5 / \ 0 0 1 2 3 4 5 Defamation [mm] FIGURE D34 - Shear test result for Y34_04 LL. 199 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00202 [m] MT thick 0.0002 [m] Sample area 0.000111 [mA2] —> MT Area 6.91E-06 [mA2] Velocity 0.004233 [mls] Yield stress 132484 [Pa] Predicted: Yield Strength 0.92 [N] Force [N] WWwW 0 1 2 3 4 5 Defamation [mm] FIGURE D35 - Shear test result for Y3 5_04 F. LIL. 200 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0024 [m] MT thick 0.0002 [m] Sample area 0.000132 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 155948 [Pa] Predicted: Yield Strength 1.08 [N] 25 20 '\ 15 Z / MWWWW O 1 2 3 4 5 Defamation [mm] FIGURE D36 - Shear test result for Y3 6_04 201 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0024 [m] MT thick 0.0002 [m] Sample area 0.000132 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [m/s] Yield stress 159658 [Pa] Predicted: Yield Strength 1.10 [N] 25 15 r E 0 1 2 3 4 5 Defamation [mm] FIGURE D37 - Shear test result for Y37_04 202 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00223 [m] MT thick 0.0002 [m] Sample area 0.000123 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 146990 [Pa] Predicted: Yield Strength 1.02 [N] 20 15 ‘1 E 8 1° 1 “8- 1 . \ Defamation [mm] FIGURE D38 - Shear test result for Y3 8_04 203 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.0021 [m] MT thick 0.0002 [m] Sample area 0.000115 [mA2] —> MT Area 6.91E-06 [m42] ‘ Velocity 0.004233 [mls] Yield stress 165813 [Pa] Predicted: Yield Strength 1.15 [N] 20 ' .5 1 I I _ 1 a 1 g 10 fl — — — — ‘1 S 0 \LW/Lr/MAA ,5. Id 0 1 2 3 4 5 Defamation [mm] FIGURE D39 - Shear test result for Y39_04 ‘1. ‘ 204 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00154 [m] MT thick 0.0002 [m] Sample area 8.46E-05 [m"2] —> MT Area 6.91E—06 [m"2] Velocity 0.004233 [mls] Yield stress 168136 [Pa] Predicted: Yield Strength 1.16 [N] n—y—np Force [N] Defamation [mm] FIGURE D40 - Shear test result for Y40_04 205 SHEAR TEST MT Dia. 0.011 [m] "Thickness 0.00214 [m] MT thick 0.0002 [m] Sample area 0.000118 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 146777 [Pa] Predicted: Yield Strength 1.01 [N] 20 .1 41 . ° ’ 57 \ L—m .. i l Mimi/WWW Defamation [mm] FIGURE D41 - Shear test result for Y41_04 206 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00237 [m] MT thick 0.0002 [m] Sample area 0.00013 [m42] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [m/s] Yield stress 116952 [Pa] Predicted: Yield Strength 0.81 [N] 20 15 -\ 4 E 3 10 E '\/\/‘ MAM/WW1 0 1 2 3 4 5 Defamation [mm] FIGURE D42 - Shear test result for Y42_04 207 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00262 [m] MT thick 0.0002 [m] Sample area 0.000144 [mA2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 134993 [Pa] Predicted: Yield Strength 0.93 [N] 25 .0 r1 1 / Force [N] . / Wm O 1 2 3 4 5 Defamation [mm] FIGURE D43 - Shear test result for Y43_04 fix,“ 208 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00256 [m] MT thick 0.0002 [m] Sample area 0.000141 [m"2] ~—*> MT Area 6.91E-06 [mAZ] Velocity 0.004233 [mls] Yield stress 123545 [Pa] Predicted: Yield Strength 0.85 [N] 20 15 l 1 z 1 g 10 12 l 5 / . \ 1 l \J\/ 0 0 1 2 3 4 5 Defamation [mm] FIGURE D44 - Shear test result for Y44_O4 209 SHEAR TEST MT Dia. 0.011 [m] Thickness 0.00229 [m] MT thick 0.0002 [m] Sample area 0.000126 [m"2] —> MT Area 6.91E-06 [m"2] Velocity 0.004233 [mls] Yield stress 130808 [Pa] Predicted: Yield Strength 0.90 [N] Force [N] 8 WWMMWW‘ _L ‘——. . —-"i4.4 4 2 3 4 5 Defamation [mm] FIGURE D45 - Shear test result for Y45_04 r- LL APPENDIX E 210 APPENDIX E Stress-Strain Plots of Tensile Test 211 TENSILE TEST Sample Thickness 0.00206 [m] Elasticity 2219400 [Pa] Width 0.009 [m] Area 1.9E-05 [m"2] Length 0.02 [m] 3.0E+06 2.5E+06 // 2.0E+06 . // 1.5E+06 - M "’ 1.0E+06 ' A I" 1», 5.0E+06 (j 0.0E+06 111 feg 1:1/\/\/\/\0/\:: -5.0E+06 . o 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E1 - Tensile test result for Y01_05 212 TENSILE TEST Sample Thickness 0.00208 [m] Elasticity 2893126 [Pa] Width 0.009 [m] Area 1.9E-05 [m"2] Length 0.02 [m] 4.0E+06 3.5E+06 / 3.0E+06 // 2.5E+06 / / 7.. 2.0E+06 ' A n/ a. /f/ f 5% 1.5E+06 ’ /\/nv// 1.0E+06 - f,// 5.0E+06 /\/: /\/V\/VA\/\/\/ 0.0E+06 .4 -5.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E2 - Tensile test result for Y02_05 Us. 213 TENSILE TEST Sample Thickness 0.00191 [m] Elasticity 2661102 [Pa] Width 0.009 [m] Are a 1.7E-05 [mA2] Length 0.02 [m] Stress [Pa] 3.5E+06 3.0E+06 2.5E+06 2.0E+06 1.5E+06 1 .0E+06 5.0E+06 ’ 0.0E+06 -5.0E+06 ’ -1 .0E+06 [ivy/11.. 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlm] FIGURE E3 - Tensile test result for Y03_05 214 TENSILE TEST Sample Thickness 0.00261 [:11] Elasticity 2019565 [Pa] Width 0.009 [m] Area 2.3E-05 [m"2] Length 0.02 [m] 3.0E+06 2.5E+06 _ // 2.0E+06 - / - ll 4 1.5E+06 A t . V (I) 1.0E+06 . //\. 1. 1W1 0.0E+06 i f + 1 . 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E4 - Tensile test result for Y04_05 0...“! L AI" 215 TENSILE TEST Sample Thickness 0.00178 [m] Elasticity 2603817 [Pa] Width 0.009 [111] Area 1.6E-05 [m"2] Length 0.02 [m] 3.5E+06 J. 3.0E+06 ” / 2.5E+06 2.0E+06 1.5E+06 ’ffl), 1.0E+06 D 5.0E+O6 : [\[7fl/ 0.OE+06334.+1 +1: t1+ 1 1 ctr 1%i/t/ii -5.0E+06 0 Stress [Pa] -1.0E+06 ” 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E5 - Tensile test result for Y05_05 216 TENSILE TEST Sample Thickness 0.00264 [In] Elasticity 2186119 [Pa] Width 0.009 [m] Area 2.38E-05 [mA2] Length 0.02 [m] 3.0E+06 2.5E+06 v v \ 2.05.06 / 1.55.06 : MAJ/l 10906 N/ / 5.0E+06 AV/\,/ 0.0E+06 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [/. 1 1 -5.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] Stress [Pa] FIGURE E6 - Tensile test result for Y06_05 TENSILE TEST Sample Thickness Width Area Length 217 0.00273 [m] 0.009 [m] 2.46E-05 [mAZ] 0.02 [m] Elasticity 2644603 [Pa] Stress [Pa] -5.0E+06 3.58-1-06 3.0E+06 2.55-+06 2.0E+06 1.5E+06 ' 1W 1 .OE+06 5.0E+06 0.0E+06 1 1 1 1 1 1/ .l 1 UV 0.02 0.04 0.06 strain [mlrn] 0.08 0.1 0.12 0.14 FIGURE E7 - Tensile test result for Y07_05 218 TENSILE TEST Sample Thickness 0.0027 [111] Elasticity 3603414 [Pa] Width 0.009 [m] Area 2.43E-05 [mAZ] Length 0.02 [111] 5.0506 .. 4.5506 / Stress [Pa] 4.0506 " / 3.5506 0 / 3.0506 " / 2.5506 " / 2.0506 " M 1.5506 0 , 1.0506 ” / 5.0E+06 / 0.0E+0611+111 ‘11111‘1' 111 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E8 - Tensile test result for Y08_05 TENSILE TEST Sample Thickness Width Area Length 0.00253 [m] 0.009 [m] 2.28E-05 [mA2] 0.02 [m] 219 Elasticity 3262707 [Pa] Stress [Pa] 4.5E+06 > 4.0506 3.5506 3.0506 : 2.5506 2.0506 1.5506 ' 1.0506 . 5.0E+06 0.0E+06 i/ y/, j?” AA ,\/V\ 0 0.02 0.04 0.06 0.08 Strain [mlrn] 0.1 FIGURE E9 - Tensile test result for Y09_05 TENSILE TEST Sample Thickness Width Area Length 0.00272 [m] 0.009 [m] 24513-05 [In/‘2] 0.02 [m] 220 Elasticity 2382488 [Pa] Stress [Pa] -5.0E+06 3.5E+06 3.0E+06 2.5E+06 2.0E+06 1 .5E+O6 1 .0E+06 5.0E+06 0.0E+06 1Mwl V'fi 0.04 0.06 0.08 0.1 Strain [mlm] FIGURE E10 - Tensile test result for Y10_05 221 TENSILE TEST Sample Thickness 0.00207 [m] Elasticity 3038228 [Pa] Width 0.009 [in] Area 1.86E-05 [mAZ] Length 0.02 [m] 4.0506 . // 3.0E+06 A “/V Stress [Pa] I V : / . 2...... : Jfi/W V /J .e/ / 0.0506 MFA—f . -1.0E+06 -2.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E11 - Tensile test result for Y11_05 222 TENSILE TEST Sample Thickness 0.00276 [m] Elasticity 2365463 [Pa] Width 0.009 [m] Area 2.48E-05 [mA2] Length 0.02 [111] 3.5506 3.0E+06 2.5E+06 . / 2.0E+06 .v/\/- I ffd/ 1.5E+06 . /A 1.0506 . fl / 5.0506 If, 0.OE+06111111111111111111/\./11 Stress [Pa] -5.0E+06 - 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E12 - Tensile test result for Y12_05 223 TENSILE TEST Sample Thickness 0.00314 [m] Elasticity 2669880 [Pa] Width 0.009 [m] Area 2.83E-05 [mA2] Length 0.02 [m] 3.5E+06 1. / 3.0E+06 ...... 3 .WV Strain [mlrn] 11 W 1 - 11 1 WJVV FIGURE E24 - Tensile test result for Y24_05 235 TENSILE TEST Sample Thickness 0.00292 [m] Elasticity 1423046 [Pa] Width 0.009 [m] Area 2.63E-05 [mAZ] Length 0.02 [m] 2.0E+06 Stress [Pa] -5.0E+06 1.5E+06 / 1.0E+06 V" V 5.0E+06 . /\// 1 0.0E+0611—11111111111./\1'*r/\1/\111 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlm] FIGURE E25 - Tensile test result for Y25_05 236 TENSILE TEST Sample Thickness 0.00248 [m] Elasticity 3227850 [Pa] Width 0.009 [111] Area 2.23E-05 [mA2] Length 0.02 [m] 5.0E+06 Stress [Pa] -5.0E+06 4.5E+06 4.0E+06 3.5E+06 3.0E+06 2.5E+06 2.0E+06 1 .5E+06 1.0E+06 5.0E+06 0.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlm] FIGURE E26 - Tensile test result for Y26_05 237 TENSILE TEST Sample Thickness 0.00274 [m] Elasticity 3545828 [Pa] Width 0.009 [m] Area 2.47E—05 [mA2] Length 0.02 [m] 5.0E+06 / 4.0506 . / Stress [Pa] R 3, 3.0506 / 7 2.0E+06 W 1 .0E+06 Kc 0.0E+06 1 1 1 1 1 1 if 1 1 1 ”W -1.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E27 - Tensile test result for Y27_05 238 TENSILE TEST Sample Thickness 0.00259 [m] Elasticity 3437427 [Pa] Width 0.009 [111] Area 2.33E—05 [m"2] Length 0.02 [m] 5.0506 4.0506 — // 3.0506 . / 3 2.0506 3 - / (D -1 .0E+06 1.0E+06 / 0.0E+06 1 1 1 1 1 1 5+ . 1 1 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E28 - Tensile test result for Y28_05 239 TENSILE TEST Sample Thickness 0.0024 [in] Elasticity 3034021 [Pa] Width 0.009 [111] Area 2.16E-05 [mAZ] Length 0.02 [m] Stress [Pa] 5.0E+06 4.0E+06 / 2.0E+06 .. / 1.0E+06 /V//.\/—/ °'°E*°6_"""""’T*\'/V'\/\ -1.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Straln[mlm] FIGURE E29 - Tensile test result for Y29_05 240 TENSILE TEST Sample Thickness 0.00202 [m] Elasticity 3415703 [Pa] Width 0.009 [m] Area 1.82E-05 [mAZ] Length 0.02 [m] Stress [Pa] 5.0E+06 4.0E+06 // 3.0E+06 / 2.0E+06 /l/ / 1.0E+06 . [V 0.0E+061-51111 5.555.- .5 WW ’ J\/\/‘ —1.0E+06 ' /\/J -2.0E+06 _ 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E30 - Tensile test result for Y30_05 241 TENSILE TEST Sample Thickness 0.0019 [m] Elasticity 3085208 [Pa] Width 0.009 [111] Area 1.71E-05 [m"2] Length 0.02 [m] 4.5506 Stress [Pa] 4.0E+06 3.5E+06 3.0E+06 2.5E+06 2.0E+06 1 .5E+06 1.0E+06 5.0E+06 0.0E+06 0.04 0.06 0.08 Strain [mlrn] 0.1 FIGURE E31 - Tensile test result for Y31_05 242 TENSILE TEST Sample Thickness 0.00305 [m] Elasticity 1748530 [Pa] Width 0.009 [m] Area 2.75E—05 [m"2] Length 0.02 [m] 2.5E+06 / 2.0506 / 1.5E+06 1.0E+06 ’ \V 5.0E+06 / Stress [Pa] 0.0E+06 -5.0E+06 -1.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlm] FIGURE E32 - Tensile test result for Y32_05 243 TENSILE TEST Sample Thickness 0.00228 [m] Elasticity 2268548 [Pa] Width 0.009 [m] Area 2.05E-05 [mA2] Length 0.02 [m] 3.5506 3.0506 ' / 2.5506 / 2.0506 _ / E I / 1.5E+06 A i Mf/l/ m 1.05106 N / L 5...... ff/ W 0.0E+06 1 1 1 1 1 1 1 1 5 1 . : WV VIA/\VA AVAV/\ -5.0E+06 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlrn] FIGURE E33 - Tensile test result for Y3 3_05 244 TENSILE TEST Sample Thickness 0.00267 [m] Width 0.009 [111] Area 2.4E-05 [mAZ] Length 0.02 [m] »_———--.——-u-vu-— Elasticity 3514630 [Pa] Stress [Pa] ‘1 -1 .0E+06 5.0E+06 4.0E+06 / 3.0E+06 / 2.0E+06 / 1 .0E+06 Amp—RA . , r AAVAV/ 41- ooe+oe . V_7I 1_V‘T a» l I _T I 1 V V 0 0.02 0.04 0.06 0.08 0.1 Straln [mlrn] 0.12 0.14 FIGURE E34 - Tensile test result for Y34_05 245 TENSILE TEST Sample Thickness 0.00272 [m] Elasticity 2438245 [Pa] Width 0.009 [m] Area 2.45E—05 [mA2] Length 0.02 [m] 3.5E+06 Stress [Pa] -5.0E+06 3.0E+06 / 2.5E+06 , ,/ g h / l 2.0E+06 , 1.5E+06 fl'/ 1.0E+06 5.0E+06 P/ 0.0E+06 1 1 1 1 1 1 ' 1 .1 w [\V/VJV\/\VIV\VIL\VAVM 7 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Straln [mlrn] FIGURE E35 - Tensile test result for Y3 5_05 _... ~--__ , 246 TENSILE TEST Sample Thickness 0.00209 [m] Elasticity 1687081 [Pa] Width 0.009 [111] Area 1.88E-05 [mA2] Length 0.02 [m] 2.5E+06 2.0906 /, <1 1.5E+06 / 1.0906 \ APA / 5.OE'I*()6L V/ 0.0E+06l111111411 AAAAAA . . . ..A WVUVWWV / Stress [Pa] 1: l -5.0E+06 Straln [mlrn] 0 0.02 0.04 0.06 0.08 0.1 0.12 O. 14 FIGURE E36 - Tensile test result for Y36_05 247 TENSILE TEST Sample Thickness 0.00273 [m] Elasticity 1794603 [Pa] Width 0.009 [m] Area 2.46E-05 [mA2] Length 0.02 [m] 3.0E+06 2.5906 ' / 2.0E+06 . // 'a‘ 1.5906 WA / E. 1 f / 1.0906 , n/AV /1 5.0906 1 / / \/\f\/J\/\/\/\ 0.0E+06 1 V/\Cf 19' 1 ~19 1 1 /\/1\/1\/ -5.0906 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Strain [mlm] FIGURE E37 - Tensile test result for Y37_05 248 TENSILE TEST Sample Thickness 0.00282 [m] Elasticity 2176547 [Pa] Width 0.009 [m] Area 2.54E-05 [mAZ] Length 0.02 [m] 3.5E+06 3.0E+06 4? /j 2.5E+06 0 // 2.0906 0 // g .55... fl // (I) :r / + “ A 1.05 06 0 /\/\/r, / 5.0906 ”(V // :1 0.0E+06 ”A 1 1 4 A A1" "1 1Av/MAA/Hv/\1fr\w .. . W V vw V V V 5.0906 l 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Stralnlmlm] FIGURE E38 - Tensile test result for Y3 8_05 249 TENSILE TEST Sample Thickness 0.00272 [m] Elasticity 1578980 [Pa] Width 0.009 [m] Area 2.45E-05 [m"2] Length 0.02 [m] 2.5906 2.0906 / / 1.5906 // 1.0906 - /‘ / H / Stress [Pa] -5.0E+06 l l Strain [mlrn] 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 FIGURE E39 — Tensile test result for Y39_05 250 TENSILE TEST Sample Thickness 0.00225 [m] Elasticity 2231997 [Pa] Width 0.009 [m] Area 2.03E-05 [m"2] Length 0.02 [m] 3.5906 3.0906 // 2.5906 - / Stress [Pa] -5.0E+06 2.0E+06 / 1.5E+06 / 1.0E+06 / 0 0.02 0.04 0.06 0.08 Straln [mlrn] 0.1 0.12 0.14 FIGURE E40 - Tensile test result for Y40_05 251 TENSILE TEST Sample Thickness 0.00234 [m] Elasticity 1916969 [Pa] Width 0.009 [m] Area 2.11E-05 [mAZ] Length 0.02 [m] 3.0906 2.5906 / / M... ’\/[l // VA// M/ M... MWU/hu/IVA .m U'V Stress [Pa] -5.0E+06 0 0.02 0.04 ' 0.06 0.08 0.1 0.12 0.14 Straln [mlrn] FIGURE E41 - Tensile test result for Y4l_05 252 TENSILE TEST Sample Thickness 0.00241 [m] Elasticity 2252100 [Pa] Width 0.009 [111] Area 2.17E-05 [mA2] Length 0.02 [m] 3.5906 3.0906 ’ // 2.5906 ’ / Stress [Pa] -5.0E+06 -1.0E+06 - 2.0906 F / 1.5906 ’ / 1.0906 F . .. / 5.0E+06> V , 0.0E+06b111.~11\-1111k VVNM fW/Ab 0 0.02 0.04 0.06 0.08 Straln [mlrn] 0.1 0.12 0.14 FIGURE E42 - Tensile test result for Y42_05 253 TENSILE TEST Sample Thickness 0.00289 [m] Elasticity 2116101 [Pa] Width 0.009 [111] Area 2.6E-05 [mAz] Length 0.02 [m] 3.0E+06 2.5E+06 . // 2.0E+06 / 1.5906 / Stress [Pa] > > \ \ 1 .0E+06 0.0E+06111111111'r1 5.0E+06 / ‘ “J Wm -5.0E+06 0 0.02 0.04 0.06 0.08 Straln [mlrn] 0.1 FIGURE E43 - Tensile test result for Y43_05 253 TENSILE TEST Sample Thickness 0.00289 [m] Elasticity 2116101 [Pa] Width 0.009 [111] Area 2.6E-05 [mAZ] Length 0.02 [111] 3.0906 2.5E+06 ' // 2.0E+06 i / Stress [Pa] -5.0E+06 1.5906 / 1 .0E+06 5.0E+06 JJ / ‘ ; Wm 0,0E+06111 1111111111 o 0.02 0.04 0.06 0.08 Strain [mlm] 0.1 0.14 FIGURE E43 - Tensile test result for Y43_05 254 TENSILE TEST Sample Thickness 0.00269 [111] Width 0.009 [m] Area 2.42E—05 [mA2] Length 0.02 [m] Elasticity 1938822 [Pa] 3.0E+06 2.5E+06 2.0E+06 1 .5E+06 / 1.0E+06 Stress [Pa] > ‘7 / > / Straln [mlrn] fl - 5.0906 /\// \ 0.0E+06 1 1 1 1 1 1 1 111/'\/\./1\)/\1r/1\ 1 1 W . V V V v -5.0E+06 l 0 0.02 0.04 0.06 0.08 0.1 0.12 FIGURE E44 - Tensile test result for Y44_05 255 TENSILE TEST Sample Thickness 0.00266 [m] Elasticity 2625809 [Pa] Width 0.009 [m] Area 2.39E-05 [m"2] Length 0.02 [m] 4.0E+06 3.5906 ' / 3.0E+06 * / 2.5E+06 / 2.0906 / 1.5906 ; AMW 1 .0E+06 / NV \ 0.0E+06111 .1. ..1.-\f\[\nA../\. VVVVVV Stress [Pa] -5.0E+06 ’ 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Straln [mlm] FIGURE E45 - Tensile test result for Y45_05 APPENDIX F 36 APPENDIX F Computer Simulation Program #include #include #include #define V1 1 /* 1 in/min */ #define V10 10 /* 10 in/min */ #define r 0.997 #define N 100.0 /* # of layers for grand simulation */ #define lambda 1 #define alpha 3.0 #define freq 100.0 /* sampling rate [Hz]*/ #define Th 0.0002 /* Thickness of unit layer */ #define Area 3.14*pow(0.011,2)/4 /* Area of the plunger */ #define Area2 3.14*0.011*Th double Vchl, Vchlo; float F[500], iFC[500],iFS[500],FC[500],pFC[500],FS[500],U[500],iF[500],FC omp[500],Fshear[500]; float dtot[500], Frac; float D_max; double d_max, f_max; double kf_s, kd_s,pc,mo; float Fvisco=9; int m, mt; float ip[45][6]={ {2.36, 1.59, 338438, 0 79, 2219400, 0.42}, {3.85, 1 10, 441721, 2.07, 2893126, 0.33}, {4.20, 1.45, 443736, 1.28, 2661103, 0.40} PMPMAAAAAAAAAAAAAAAAAAAAAAAA ONIbUTUTNl—‘WNNNHNNWih-ibih-LUthWWI-bibww PM 14>. AAAF’HAF‘HAAAAAAAAA NWNNNWWWNNWWWWW .32, 0.97, .99, 1.62, .35, 1.84, .35, 1.93, .91, 1 25, .41, 1 52, .49, 1.57, .85, 1.65, .11, 2.57, .14, 1 25, .87, 1.54, .46, 1 58, .74, 3.10, .89, 1.66, .74, 2.41, .83, 11.6, .90, 3.71, .50, 1.30, .90, 1.33, .09, 1.21, .98, 1.83, .02, 1.81, .17, 1.13, .54, 2.01, .49, 1.99, .72, 1.92, .69, 2.86, .94, 0.88, .29, 1.76, .56, 0.89, .92, 1.00, .55, 0.97, .95, 4.34, .50, 3 62, .05, 2.62, .04, 154 , .34, 2.81, .17, 2.65, .79, 1.43, .38, 1 61, .32, 0.86, .31, 1 25, 335967, 428455, 446245, 473311, 352807, 331709, 278642, 518319, 434385, 476922, 506877, 500491, 497526, 353833, 358053, 364249, 396447, 334104, 388122, 379075, 282443, 284230, 444801, 631107, 469624, 653307, 554547, 355620, 412565, 300310, 348930, 288335, 437958, 396105, 421764, 513378, 426326, 331139, 306164, 300994, 313349, 313121, }; N N H m H o o o H o N H H H N H H H H N H H H H H H H OOOOHHHI—‘l—‘l—‘OOOOH 72, 58, 32, 52, 25, 22, 72, 16, 55, 68, 39, 47, 08, 31, 63, ll, 00, 71, 01, 90, 62, 78, 46, O3, 39, 20, .03, .23, .83, .89, .77, 92, 08, .10, 02, .15, .16, 01, .81, 93, 85, 90, 2019565, 2603817, 2186119, 2644603, 3603414, 3262707, 2382488, 3038228, 2365463, 2669880, 3140829, 3262553, 3204485, 2495790, 2556493, 2106323, 2878381, 1488604, 1718797, 1910229, 1113635, 1423046, 3227850, 3545828, 3437427, 3034021, 3415703, 3085208, 1748530, 2268548, 3514630, 2438245, 1687081, 1794609, 2176547, 1578980, 2231997, 1916969, 2252100, 2116101, 1938822, 2625809, O OOOOOOOOOOOOOOO_ OOOOOOOOOOOOOOOOOOOOOOOOOO 7'51" 258 void main(void) { int x, i, j, k, c, h,z,q,qc,y[30],s[3],minj; float pd[150],pt,A,tt,sFC,sFS,ppt,max,min,Fv; float a[lSO], as[150], t, dl, d[150]; float ds[150],dss,dds,dx,temp,pop,force,pc; float Eu,Nu,SIGMAC,F_max,Es1,ct; char *fn[45]={"m:\\m01","m:\\m02","m:\\m03","m:\\m04","m:\\m05", "m:\\m06","m:\\m07","m:\\m08","m:\\m09","m:\\m10", "m:\\mll","m:\\m12","m:\\ml3","m:\\ml4","m:\\m15",'m:\\m16", "m:\\m17","m:\\ml8","m:\\m19","m:\\m20","m:\\m21",'m:\\m22", "m:\\m23","m:\\m24","m:\\m25","m:\\m26","m:\\m27","m:\\m28", "m:\\m29","m:\\m30","m:\\m31","m:\\m32","m:\\m33","m:\\m34", "m:\\m35","m:\\m36","m:\\m37","m:\\m38","m:\\m39","m:\\m40", "m:\\m41","m:\\m42","m:\\m43","m:\\m44","m:\\m45"}; FILE *destin, *M_T; Vchl=Vl*0.0254/60; VCth=VlO*0.0254/60; M_T=fopen("m:\\MagT","wt"); for(x=0;x<=44;x++) { Eu=ip[x][0]*pow(10,10); Nu=ip[x][1]*pow(10,10); SIGMAc=ip [x] [2]; F_max=ip[x][3]; Es1=ip[x][4]; ct=ip[x][5]; destin=fopen(fn[x],"wt"); q=li FV=O; a[l]=(l-r)*Vcth/(l—powl(r,N)); for(i=0;i<=l49;i++) { 2” dsli]=0.0; } Frac=SIGMAc*Area; dx=Vch10*(l/freq); Ppt=0; Pt=0; C=l; for(i=0;i<=499;i++) { iFCli]=0; iFSlil=0; Ulil=0; pFCli]=O; } for(i=0;i<=N;i++) { pd[i]=0; } as [1] =61 [1]; for(j=2;j<=N;j++) { a[j]=a[j-l]*r; aS[j]=a[j]; } for(i=0;i<=499;i++) { t=i*(l/freq); /* time */ dtot[i]=t*Vch10; for(j=l;j<=N;j++) { d[j]=pd[j]+a[j]*(t-pt); } for(j=c;j<=N;j++) { /* tt=(dtot[i]-A*(c-l))/Vch10;*/ FC[j]=Area*a[j]*Nu*(l—1/exp(Eu*((t—pt)+0.6*ppt)/Nu)); if(j==c) FC[j]=Area*a[j]*(Nu*(1-0.06*c))*(1- l/exp(Eu*((t-pt)+0.6*ppt)/(Nu*(l-0.06*c)))); 260 if(c== ) FC[j]=Area*a[j]*Nu*(1—1/exp(Eu*t/Nu)); if(FC[j]>=(Frac*(l+0.03*(c—1)))) y[c]=i; for(k=c+1;k<=N;k++) { pFC[k]=FC[k]; Pt=t; if(c==1) ppt=t; C+=l; FC[j]=O; d[j]=Th; A=Th; pdljl=Ai if((N—c+1)== ) goto SSS; a[c]=(1-r)*Vch10*(l-0.09*c)/(l—powl(r,(N—c+l))); for(k=(c+l);k<=N;k++) { a[k]=a[k-l]*r; for(k=c;k<=N;k++) { pd[k]=a[k]*(dtot[i]—A*(C-l))/Vch10; d [k] =pd [k] ; } goto FFF; } } FFF: ch=0; for(j=l;j<=N;j++) { sFC=sFC+FC[j]*(a[j]*(1/freq)); iFC[i]=iFC[i-l]+SFC; if(i== ) iFC[i]=SFC; a[0]=Vch10; as[0]=a[0]; SFS=O; z=l; for(j=q;j<=N;j++) { 261 dSS=0; ddS=O; for(k=j;k<=N;k++) { if(i== ) Z=0; dss+=a[k]*(l/freq)*z; dds+=a[k]*(l/freq)*z; } if(qc!=q) { mo=0.6; } else { mo=1;} ds[j]=ds[j]*mo+dds; force=(Esl*(ds[j]*0.7)/0.025)*Area2; if(force>F_max) q+=1; sFS=sFS+dss*force; } PC=C; iFS[O]=O; if(i>0) iFS[i]=iFS[i-l]+sFS; U[i] =iFC [i] +iFS [i] ; F[0]=O; min=200; for(j=y[l];j<=Y[2];j++) { if(Flj]=minj) { h=l; Fv=h*(i-minj)*(1/freq)*Fvisco/(minj*0.01); if(Fv>=Fvisco) Fv=Fvisco; } . -_.-——u I. i.) 262 F[i]=(U[i+l]—U[i-1])/(2*dx)+Fv; fprintf(destin,"%f %f %f %f %f\n",i*(1/freq),iFC[i],iFS[i],U[i],F[i]); } for(i=l;i<=2;i++) { max=0; for(j=(y[i]—10);j<=(y[i]+10);j++) { if(F[j]>max) { max=F[j]; s[i]=j; } 1 1 } {1 max=0; . for(i=l;i<=l89;i++) { /* 189: Depth of MT [8 mm]*/ if(F[i]>maX) { max=F[i]; mt=i; } fprintf(M_T,"%d %g %g %g %d %g\n",x+l,F[s[l]],F[s[2]],F[s[l]]/(s[l]*0.0l*Vcth),mt,F[mt] ); fclose(destin); printf("%2d\n",x+l); } fclose(M_T); } * This program was written by using C language. LIST OF REFERENCES 263 LIST OF REFERENCES Armstrong, RR. 1989. Measurement of apple firmness using the acoustic impulse response. Ph. D. Dissertation. Michigan State University, East Lansing, MI. Armstrong, RR. and GK. Brown. 1992. Non-destructive firmness measurement of apple. ASAE Paper No. 936023. St. Joseph, Mich.:ASAE. Bourne, MC. 1982. Considerations of a general rheological model for the mechanical behavior of viscoelastic solid food materials. J. of Texture Studies 7:243-255. Chen, Y. and J. Rosenberg. 1977. Nonlinear viscoelastic model containing a yield element for modeling a food material. J. of Texture Studies 8:477-485. Crandall, S.H., N.C. Dahl, and T.J. Lardner. 1978. An Introduction to the Mechanics of Solids. New York: McGraw-Hill Book Co. Dickinson, E. and LC. Goulding. 1980. Yield behavior of crumbly English cheeses in compression. J. of Texture Studies 11:51-63. Drake, B. 1971. A quasi-rheological model element for fracture. J. of Texture Studies 2:365-372. Johnson, E.A., M. Peleg, R.A. Segars, and J .G. Kapsalis. 1981. A generalized phenomenological rheological model for fish flesh. J. of Texture Studies 12:413-425. Khan, AA. and J .F.V. Vincent. 1990. Anisotropy of apple parenchyma. J. Sci. Food and Agric. 52:455-466. 264 Khan, AA. and J .F.V. Vincent. 1993. Compressive stiffness and fracture properties of apple and potato parenchyma. J. of Texture Studies 24:423-435. Lin, T. and RE. Pitt. 1986. Rheology of apple and potato tissue as affected by cell turgor pressure. J. of Texture Studies 17:291-313. McLaughlin, N.B. 1987. Statistical models for failure of apple tissue under constant- strain-rate loading. J. of Texture Studies 18:173-186. Mohsenin, N.N. 1977. Characterization and failure in solid foods with particular reference to fruits and vegetables. J. of Texture Studies 8:169-193. Mohsenin, N.N. and H. Gohlich. 1962. Techniques for determination of mechanical properties of fruits and vegetables as related to design and development of harvesting and processing machinery. J. of Agricultural Engineering Research 7:300-315. Mohsenin, N.N. 1986. Physical Properties of Plant and Animal Materials. New York: Gordon and Breach Science Publishers. Peleg, M. 1976. Consideration of general rheological model for the mechanical behavior of viscoelastic solid food materials. J. of Texture Studies 72243-255. Peleg, M. 1976. Compressive failure patterns of some juicy fruits. J. of Food Sci. 41:1320-1324. Pitt, RE. 1982. Models for the rheology and statistical strength of uniformly stressed vegetative tissue. Transaction of ASAE 1776-1784. Shigley, J .E. and CR. Mischke. 1984. Mechanical Engineering Design. New York: McGraw-Hill Book Co. White, RM. 1986. Fluid Mechanics. New York: McGraw-Hill Book Co. "71111111111111.1111