Biaxial Mechanical Behavior of Swine Pelvic Floor Ligaments: Experiments and Modeling

Size: px
Start display at page:

Download "Biaxial Mechanical Behavior of Swine Pelvic Floor Ligaments: Experiments and Modeling"

Transcription

1 Biaxial Mechanical Behavior of Swine Pelvic Floor Ligaments: Experiments and Modeling Winston Reynolds Becker Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the degree of Master of Science in Engineering Mechanics Raffaella De Vita, Chair Sunghwan Jung John J. Socha May 6, 2014 Blacksburg, Virginia Keywords: Cardinal Ligament, Uterosacral Ligament, Pelvic Floor Disorders, Viscoelasticity, Stress Relaxation

2 Biaxial Mechanical Behavior of Swine Pelvic Floor Ligaments: Experiments and Modeling Winston Reynolds Becker ABSTRACT Although mechanical alterations to pelvic floor ligaments, such as the cardinal and uterosacral ligaments, are one contributing factor to the development and progression of pelvic floor disorders, very little research has examined their mechanical properties. In this study, the first biaxial elastic and viscoelastic tests were performed on uterosacral and cardinal ligament complexes harvested from adult female swine. Biaxial elastic testing revealed that the ligaments undergo large strains and are anisotropic. The direction normal to the upper vagina was typically stiffer than the transverse direction. Stress relaxation tests showed that the relaxation was the same in both directions, and that more relaxation occurred when the tissue was stretched to lower initial strains. In order to describe the experimental findings, a threedimensional constitutive model based on the Pipkin-Rogers integral series was formulated and the parameters of such model were determined by fitting the model to the experimental data. In formulating the model, it was assumed that the tissues consist of a ground substance with two embedded families of fibers oriented in two directions and that the ligaments are incompressible. The model accounts for finite strains, anisotropy, and strain-dependent stress relaxation behavior. This study provides information about the mechanical behav-

3 ior of female pelvic floor ligaments, which should be considered in the development of new treatment methods for pelvic floor disorders. iii

4 Acknowledgments I would like to thank my advisor, Dr. Raffaella De Vita, for her guidance and mentorship during my undergraduate and graduate studies. I am grateful to Dr. Jung and Dr. Socha for serving on my committee and providing feedback on my work. I would also like to thank the other members of the Mechanics of Soft Biological Systems Laboratory, Matt Webster and Ting Tan, for all of their help over the past few years. This material is based upon work supported by the National Science Foundation under CAREER Grant No and the National Science Foundation Graduate Research Fellowship Program. iv

5 Contents 1 Introduction 1 2 Methods Sample Preparation Biaxial Testing Biaxial Test Protocol Theoretical Formulation Constitutive Model Biaxial elastic deformation Biaxial stress relaxation response Identification of Model Parameters v

6 4 Results Biaxial Elastic Properties Biaxial Viscoelastic Properties Modeling Discussion Experimental Results Modeling Conclusions Summary and Future Directions Bibliography 36 vi

7 List of Figures 2.1 a) Bath of Instron planar biaxial tensile testing system containing load cells, actuator arms, and specimen. b) Orientation at which samples were cut from swine uterosacral cardinal ligament complexes Protocol used during biaxial mechanical testing of pelvic floor ligaments. The displacement was controlled and was always the same in both directions Stress-stretch curves and normalized relaxation curves of 9 samples taken from a single sow. The locations the specimens were taken from are indicated by the schematic in the top left of the top two graphs. The stress relaxation curves represent the same samples presented in the stress-stretch curves. In both the relaxation curves and the stress-stretch curves, the filled symbols represent data collected from stretching in the direction normal to the upper vagina/cervix while the unfilled symbols denote data collected from stretching in the direction transverse to the upper vagina/cervix vii

8 4.2 Stress relaxation response for two specimens. Three tests were performed on each specimen with a 1 hour resting period between each test. a) The specimen was stretch to λ 1 = and λ 2 = during the first relaxation test, λ 1 = and λ 2 = during the second relaxation test, and λ 1 = and λ 2 = during the third relaxation test. b) The specimen was stretch to λ 1 = and λ 2 = during the first relaxation test, λ 1 = and λ 2 = during the second relaxation test, and λ 1 = and λ 2 = during the third relaxation test Continuum model fit to experimental data of 13 specimens taken from 3 sows. Filled symbols represent data collected from stretching in the direction normal to the upper vagina/cervix, the unfilled symbols denote data collected from stretching in the direction transverse to the upper vagina/cervix, and the solid lines indicate the constitutive model Dependency of α 1 (Left) and α 2 (Right), which are measures of the percent relaxation of the specimens, on the first invariant. The lines represent the functions α 1 (I 1 (τ)) = a 1 e a 2(I 1 3) and α 2 (I 1 (τ)) = a 3 e a 4(I 1 3) with a 1 = 1.078, a 2 = 2.303, a 3 = 2.175, and a 4 = Predicted stress relaxation behavior of ligaments stretched to different levels based on fitting the model to all 22 data sets simultaneously. The predictions were computed by fitting the constitutive model to all 22 sets of relaxation data simultaneously viii

9 List of Tables 4.1 Mean and standard deviation of the parameters of all 22 specimens considered in this study ix

10 Chapter 1 Introduction Recent research predicts that the number of American women with pelvic floor disorders (PFDs), including fecal incontinence, urinary incontinence, and pelvic organ prolapse (POP), will increase from 28.1 million in 2010 to 43.8 million in 2050 (38). Studies have suggested that approximately 11% of women will undergo surgery for POP in their lifetime, and that around 30% of these women will require additional procedures due to recurrence (22). Alhough the causes of PFDs are not completely understood, mechanical changes to pelvic floor ligaments appear to contribute to their development and progression (21). In particular, PFDs may result from damage to the cardinal (CL) and uterosacral ligaments (USL), which provide apical support to the uterus and upper vagina (6; 18). The CL and USL are visceral ligaments that connect the upper vagina or cervix to the pelvic sidewall and provide support to the internal pelvic organs. These visceral ligaments, which 1

11 Winston R. Becker Chapter 1. Introduction 2 connect an organ to the body wall, are different from skeletal ligaments, which connect bone to bone, in both structure and composition. The CL and USL form a thin, membranelike complex that is composed of smooth muscle, blood vessels, nerve fibers, collagen, and elastin. This structure allows the ligaments to provide support to the pelvic organs while accommodating for highly mobile structures, such as the uterus (29). The loading conditions experienced by these ligaments in vivo are likely biaxial, with larger forces occurring along the axis normal to the upper vagina/cervix. Although mechanical damage to pelvic floor connective tissue is closely linked to the development of PFDs, no experiments have attempted to investigate the mechanical response of the CL, and only a few have described the mechanical properties of the USL. To this point, the stress relaxation behavior of USLs in the cynomolgus monkey (Macaca fascicularis) (35; 32) and the tensile properties of USLs taken from cadavers (17) have been examined. While these tests provided the first description of the mechanical properties of these tissues, many things remain unknown about their mechanical behavior. To the authors knowledge, these studies were the first to show that the USL and CL undergo large strains. This finding could be expanded by employing more accurate strain measurement techniques. This includes the use of optical methods to measure strain as well as not assuming that the strain is linear. Optical measurements are preferred to grip-to-grip measurement because they do not include boundary effects near the grips and they measure the strain of the actual specimen. Furthermore, assuming linearity when computing the strain can lead to errors in strain calculations, particularly at the high strains that the CL and

12 Winston R. Becker Chapter 1. Introduction 3 USL experienced in the previous studies. The studies on the viscoelastic properties of the USL were the first to examine the relaxation behavior of these tissues (35; 32). The time-depedent mechanical properties of the CL and USL are highly relevant to the physiology because they experience loading for long periods of time in vivo. In the previous experiments, stress relaxation was performed multiple times on each specimen without any time for the specimen to recover between tests. This revealed that significant relaxation occurs in these tissues, but many details about the relaxation behavior of the CL and USL have not been studied. The effect of the mechanical history on the relaxation behavior, which has been shown to alter the viscoelastic properties of other ligaments (34), has not been studied. It is also unknown if the stress relaxation of the CL and USL is nonlinear and dependent on strain, which is the case in many articular ligaments (24; 14). Another important contribution of the previous mechanical work on the USL was the quantification of the uniaxial mechanical properties. However, because these ligaments support complex, biaxial loads in vivo, additional information about their physiologic mechanical behavior could be collected with biaxial mechanical testing (13; 36; 11; 20; 39). This would reveal if the CL and USL are anisotropic, and show how they respond to biaxial loading conditions. Although a few studies have examined the mechanical properties of the CL and USL, no attempt has been made to model the elastic or viscoelastic behavior of these ligaments. Mechanical models for these tissues should account for the biaxial loads that they can ex-

13 Winston R. Becker Chapter 1. Introduction 4 perience in vivo, as well as the experimentally observed stress relaxation behavior. Many constitutive theories have been applied to the relaxation behavior of other biological tissues. For cases where the stress relaxation behavior is not strain dependent, quasi linear viscoelasticity (QLV) (10) is applicable. Models similar in form to QLV have been employed to model the biaxial relaxation of epicardium (1), the urinary bladder wall (20), and aortic elastin (39). For cases where relaxation is strain dependent, modifications to QLV have been applied to model the one dimensional relaxation behavior of collagen gels (26). Other constitutive theories, such as nonlinear superposition (9) and the Pipkin-Rogers integral series (23), also capture the strain-dependent viscoelastic behavior commonly demonstrated by biological specimens. These formulations have been used to successfully model the uniaxial nonlinear relaxation behavior of ligaments and tendons (24; 14; 7; 8; 4). A few studies have considered biaxial loading conditions with these models. In particular, Rajagopal and Wineman (28) discuss the application of the Pipkin-Rogers integral series to a variety of loading conditions without specifying a particular material, and Davis and De Vita (5) formulated a three-dimensional model based on the Pipkin-Rogers integral series for articular ligaments. However, these models have not been applied to materials for which biaxial mechanical data is available. In this study, biaxial elastic and stress relaxation tests were performed on CLs and USLs taken from a swine model. The strain dependency of stress relaxation was studied by stretching specimens to different levels of strain. A three-dimensional constitutive model accounting for anisotropy and nonlinear stress relaxation was formulated and the parameters were de-

14 Winston R. Becker Chapter 1. Introduction 5 termined by fitting the model to the experimental data.

15 Chapter 2 Methods 2.1 Sample Preparation All animal protocols were approved by the Institutional Animal Care and Use Committee (IACUC) at Virginia Polytechnic Institute and State University. Four adult female domestic swine (Sus scrofa domesticus) from a slaughterhouse were dissected and the pelvic floor ligaments were extracted. The swine model was selected because its pelvic floor tissues are histologically similar to humans and because prolapse can naturally occur in swine (12). The CL and USL were carefully cleaned of extraneous fat and muscle tissue, cut into approximately 3 cm squares, and frozen (orientation of specimens shown in Fig. 2.1). Previous research has shown that the proper freezing of ligaments has little or no effect on their mechanical properties (37). Prior to testing, the specimens were thawed at room temperature. 6

16 Winston R. Becker Chapter 2. Methods 7 The length and width of each specimen was determined optically by analyzing pictures taken with a CMOS camera (DCC1645C, Thor Labs), and the thickness of each specimen was determined by averaging four thickness measurements made with calipers containing a force gauge (accuracy ± 0.05 mm, Mitutoyo ABSOLUTE Low Force Calipers Series 573, Japan) by applying a 50 gf compressive load. 2.2 Biaxial Testing Four bent safety pins connected to fishing line (4 lb Extra Tough Trilene, Berkley Fishing Company, USA) were inserted into each side of each specimen. The specimens were then mounted into an Instron Planar Biaxial Tensile Testing Device (Partially shown in Fig. 2.1, Instron, Norwood, MA, USA) by wrapping the fishing line around custom grips. These grips consisted of two pulleys and a bearing of negligible friction that could rotate to ensure that the tensions in each of the four lines attached to each side of the specimen were equivalent. During testing, the specimens were submerged in a bath containing phosphate buffered saline (PBS, 0.5 M, ph 7.4) at 37 o C. Load was recorded with four 20 N load cells (accuracy 0.25% or better) simultaneously during mechanical testing. The force along each axis was taken as the average of the loads recorded by the two load cells located along each axis (which were approximately equal). Stress was calculated by dividing the average force measured along each axis by the cross-sectional area of the sample along the given axis. Prior to testing, four poppy seeds were glued to the surface of the ligaments to produce

17 Winston R. Becker Chapter 2. Methods 8 suitable contrast for non-contact strain measurements. Images were captured with a CMOS camera (DCC1645C, Thor Labs) and a 25 mm fixed focal length lens (25 mm compact fixed focal length lens, TECHSPEC). The displacements of the poppy seeds were tracked with ProAnalyst Software (Xcitex Inc. Woburn, MA, USA). The displacement gradient was estimated using an interpolation method originally introduced by Humphrey (16) and later implemented by other researchers (30). From the displacement gradient, the deformation gradient tensor, F, and the right Cauchy-Green deformation tensor, C, were calculated as follows: F = u X + I (2.1) C = F T F (2.2) where u is the displacement vector, X is the position vector of the undeformed configuration, and I is the identity tensor. 2.3 Biaxial Test Protocol A schematic of the testing protocol is presented in Fig All specimens (n=22) were preconditioned with the same protocol: First the specimens were preloaded to 0.04 N, which was taken to be the point of zero load for the test. Then the specimens were loaded from 0.1N to 0.6N for 10 cycles at a displacement rate of 0.1 mm/s in both directions. Following preconditioning, the specimens were displaced at a constant rate (0.1 mm/s) in both directions until the force along one axis reached a force between 2 and 12 N that was selected

18 Winston R. Becker Chapter 2. Methods 9 a) b) RECTUM E2 E1 V A G I N A Figure 2.1: a) Bath of Instron planar biaxial tensile testing system containing load cells, actuator arms, and specimen. b) Orientation at which samples were cut from swine uterosacral cardinal ligament complexes. prior to the test. More specifically, 9 specimens were displaced to a maximum load between 2 and 4.9N, 7 specimens were displaced to a maximum load between 5 and 7.9 N, and 6 specimens were displaced to a maximum load between 8 and 12 N. The data from these ramps were recorded and used to determine the elastic properties of the specimens. Once the desired force was reached, the ligaments were allowed to relax at a constant displacement for at least 50 minutes while the force and time data were recorded. The testing device controlled the displacement during these tests, and stopped the displacement once a preset load was reached. By displacing the specimens to a range of forces from 2 to 12 N, we were able to ensure that a range of maximum displacements would be considered so that we could examine the strain dependent relaxation behavior of the ligaments. The displacement rate was selected to ensure that the elastic response would be quasi static,

19 Winston R. Becker Chapter 2. Methods 10 Displacement Ramp (Elastic Response) Preconditioning Relaxation Additional Relaxation (7 Samples) Time Figure 2.2: Protocol used during biaxial mechanical testing of pelvic floor ligaments. The displacement was controlled and was always the same in both directions. because similar displacement rates have been employed by others for biaxial stress relaxation tests on collagenous tissues (1), and because it is unlikely that high strain rates would be experienced by these tissues in vivo. Multiple stress relaxation tests at different displacements were performed on 7 of the samples in this study. For each sample, a relaxation test was performed as described above, followed by two additional relaxation tests immediately (n=2 samples), after a 1 hour (n=3 samples) resting period between each test, or after a 12 hour (n=2 samples) resting period between each test. For half of the specimens, three initial displacement levels were applied in increasing order (e.g. displacement at 3 N, displacement at 6 N, displacement at 9 N), and for the other half of the specimens three initial displacement levels were applied in decreasing order (e.g. displacement at 9 N, displacement at 6 N, displacement at 3 N).

20 Chapter 3 Theoretical Formulation It was assumed that the cardinal and uterosacral ligaments consist of an isotropic matrix of elastin and proteoglycans as well as a fibrous component and smooth muscle cells oriented both transversely and orthogonally to the loading direction. For this reason, the ligaments are assumed to be orthotropic. Due to the high water content of ligaments, they are assumed to be incompressible. 3.1 Constitutive Model The integral series formulation, first introduced by Pipkin and Rogers (23) was employed to model the three-dimensional stress relaxation behavior of these ligaments. The first term of the integral series is used, resulting in a Piola-Kirschoff Stress tensor, P(t), of the form (28) 11

21 Winston R. Becker Chapter 3. Theoretical Formulation 12 ( t P(t) = p(t)f T (t) + F(t) R[C(t), 0] + R[C(τ), t τ] ) (t τ) (3.1) where F(t) is the deformation gradient, C(t) is the right Cauchy-Green deformation tensor, R[C(τ), t τ] is the tensorial relaxation function, p(t) is the Lagrange multiplier that accounts for incompressibility, and t is time. Because the ligaments are taken to be orthotropic, the tensorial relaxation function can be written as R[C(τ), t τ] = k k 2 C(τ) + k 3 M M + k 4 [M (C(τ)M) + (C(τ)M) M] + k 5 N N + k 6 [N (C(τ)N) + (C(τ)N) N] (3.2) where M is the unit vector that defines the direction normal to the upper vagina/cervix, N is the unit vector that defines the direction transverse to the upper vagina/cervix, and k 1, k 2, k 3, k 4, k 5, and k 6 are constitutive functions that depend on the invariants of C and t τ. The invariants of C are defined as I 1 (τ) = tr(c(τ)), I 2 (τ) = 1 2 (I 1(τ) 2 tr(c(τ) 2 ), I 4 (τ) = M C(τ)M I 5 (τ) = M C(τ) 2 M, I 6 (τ) = N C(τ)N. (3.3) The relaxation function is assumed to consist of components accounting for the contribution of the isotropic matrix, the family of fibers normal to the upper vagina/cervix, and the family of fibers transverse to the upper vagina/cervix. This can be written in the following form:

22 Winston R. Becker Chapter 3. Theoretical Formulation 13 R[C(τ), t τ] = R gs [C(τ), t τ] + R cfn [C(τ), t τ] + R cft [C(τ), t τ] (3.4) where R gs is the relaxation function that accounts for the contribution of the ground substance, R cfn is the relaxation function that accounts for the contribution of the family of fibers normal to the upper vagina/cervix, and R cft is the relaxation function that accounts for the contribution of the family of fibers transverse to the upper vagina/cervix. Now, it is assumed that k 1 is a function of I 1 (τ) and t τ, k 3 is a function of I 4 (τ) and t τ, k 5 is a function of I 6 (τ) and t τ, and k 2, k 4, and k 6 are zero. These assumptions allow k 1 to account for the behavior of the isotropic matrix, k 3 to account for the behavior of the fibers normal to the upper vagina/cervix, and k 5 to account for the behavior of the fibers transverse to the upper vagina/cervix. This leads to relaxation functions, R gs, R cfn, and R cft of the form: R gs = k 1 1, (3.5) R cfn = k 3 M M, (3.6) R cft = k 5 N N. (3.7) Finally, the constitutive functions, k 1, k 3, and k 5, were selected to capture the behavior of these ligaments: k 1 = c 1 (3.8)

23 Winston R. Becker Chapter 3. Theoretical Formulation 14 c 2 (e c 3(I 4 (τ) 1) 1)r 1 (I 1 (τ), t τ), I 4 (τ) > 1 k 3 = 0, I 4 (τ) 1 c 4 (e c 5(I 6 (τ) 1) 1)r 1 (I 1 (τ), t τ), I 6 (τ) > 1 k 5 = 0, I 6 (τ) 1 (3.9) (3.10) where c 1, c 2, c 3, c 4, and c 5 are constant parameters. It is assumed that the fibrous component does not support compressive loads in the definition of k 3 and k 5. The relaxation function, r 1, was selected to model the isotropic stress relaxation response of the ligament: r 1 (I 1 (τ), t τ) = 1 + α 1(I 1 (τ))e (t τ)β 1 + α 2 (I 1 (τ))e (t τ)β α 1 (I 1 (τ)) + α 2 (I 1 (τ)) (3.11) where α 1 (I 1 (τ)) and α 2 (I 1 (τ)) are functions accounting for the strain dependent viscoelastic behavior of the ligaments and β 1 and β 2 are constant parameters. Note that r 1 goes to 1 when t = τ. The functions α 1 (I 1 (τ)) and α 2 (I 1 (τ)) were defined as: α 1 (I 1 (τ)) = a 1 e a 2(I 1 3) (3.12) α 2 (I 1 (τ)) = a 3 e a 4(I 1 3) (3.13) where a 1, a 2, a 3, and a 4 are constant parameters.

24 Winston R. Becker Chapter 3. Theoretical Formulation Biaxial elastic deformation Let {E 1, E 2, E 3 } and {e 1, e 2, e 3 } be two sets of unit vectors that define the reference and deformed configurations, respectively. The ligament is assumed to undergo an isochoric deformation described by x 1 = λ 1 (t)x 1, x 2 = λ 2 (t)x 2, x 3 = 1 λ 1 (t)λ 2 (t) X 3. (3.14) This deformation was selected because a negligible amount of shear was observed in the biaxial experiments. It follows that the deformation gradient tensor, F(t), and the right Cauchy-Green strain tensor, C(t), are given by: F(t) = λ 1 (t)e 1 E 1 + λ 2 (t)e 2 E 2 + C(t) = λ 1 (t) 2 e 1 E 1 + λ 2 (t) 2 e 2 E λ 1 (t)λ 2 (t) e 3 E 3, 1 λ 1 (t) 2 λ 2 (t) 2 e 3 E 3. (3.15) To derive the first-piola Kirchhoff stress tensor that defines the instantaneous elastic response, Eq. 3.1 is evaluated at t = τ by substituting Eqs and Eq It is assumed that the surface of the ligament is traction-free so that P 33 = 0. Enforcing this boundary condition to solve for the Lagrange multiplier, p, P 11 and P 22 are found to be: P 11 = c 2 ( 1 + e c 3(λ 2 1 1) )λ 1 + c 1 (λ 1 1 ) λ 3 1λ 2 2 (3.16) P 22 = c 4 ( 1 + e c 5(λ 2 2 1) )λ 2 + c 1 (λ 2 1 ) λ 3 2λ 2 1 (3.17)

25 Winston R. Becker Chapter 3. Theoretical Formulation Biaxial stress relaxation response To model the stress relaxation response, constant stretches, λ 1 (t) = λ 1 and λ 2 (t) = λ 1, are applied. The stress relaxation response is then given by solving Eq. 3.1 by substituting Eqs , Eq. 3.15, and setting P 33 = 0. This results in P 11 and P 22 of the form: P 11 = c 1 (λ 1 1 )+ λ 3 1λ 2 2 c 2 λ 1 ( 1 + e c 3(λ 2 1 1) ) 1 + α 1(I 1 (τ))e (t τ)β 1 + α 2 (I 1 (τ))e (t τ)β α 1 (I 1 (τ)) + α 2 (I 1 (τ)) (3.18) P 22 = c 1 (λ 2 1 )+ λ 3 2λ 2 1 c 4 λ 2 ( 1 + e c 5(λ 2 2 1) ) 1 + α 1(I 1 (τ))e (t τ)β 1 + α 2 (I 1 (τ))e (t τ)β 2. (3.19) 1 + α 1 (I 1 (τ)) + α 2 (I 1 (τ)) 3.4 Identification of Model Parameters The built in minimization function in MATLAB, fmincon, was utilized to compute the model parameters by minimizing the sum of the squares function for the error between the model and the data. To determine the elastic parameters c 1 -c 5, the elastic experimental data in both directions was reduced and Eqs and 3.17 were simultaneously fit to the data. Although a small amount of shear was observed in all tests, the specimens were assumed to undergo no shear when calculating the model parameters. To compute the viscoelastic parameters, the experimental data from each direction of the relaxation tests were reduced and the elastic parameters were held constant. First, β 1, β 2, α 1, and α 2 were computed by

26 Winston R. Becker Chapter 3. Theoretical Formulation 17 minimizing the sum of the squares function for the error between one set of biaxial relaxation data of a sample stretched to an intermediate stretch level and Eqs and Then, the parameters β 1 and β 2 were held constant while α 1 and α 2 were computed for each individual data set by fitting Eqs and 3.19 to each data set. Next, plots of α 1 and α 2 versus I 1 were generated and Eqs and 3.13 were fit to the data to determine the parameters a 1 -a 4. All parameters were constrained to be greater than zero. An alternative approach to fitting the model was also considered for the relaxation. In this case the relaxation function was defined a priori and the constants a 1, a 2, a 3, a 4, β 1, and β 2 were computed by fitting the model to the data collected from all 22 relaxation tests simultaneously.

27 Chapter 4 Results 4.1 Biaxial Elastic Properties Nonlinear, anisotropic mechanical behavior was observed during the ramp portion of the mechanical tests. The ligaments experienced large strains as high as 75% during the tests conducted. While anisotropy was typically observed, it was much more pronounced in some cases than others. The results of the biaxial elastic tests for a single sow are presented in Fig. 4.1 with the anatomical locations that the specimens were taken from indicated in the inset in the top corner of each plot. This figure presents 9 of the 22 specimens tested in this study. The specimens are separated into two plots based on stress level to make it easier to see each individual curve. As shown in the figure, the specimens exhibited a greater stiffness in the direction normal to the upper vagina/cervix. It is also clear that a large amount 18

28 Winston R. Becker Chapter 4. Results 19 of variability was observed, even though the specimens were taken from a single sow. The results suggest that the mechanical behavior of the tissue may be location dependent. In particular, it appears that USL specimens and specimens located closer to the USL within the USL-CL complex are stiffer in the direction normal to the upper vagina/cervix than the specimens further from the USL are in the same direction. 4.2 Biaxial Viscoelastic Properties Multiple stress relaxation tests were performed on 7 of the 22 specimens considered in this study. In each case, regardless of whether the highest or lowest stretch was applied first, the first stress relaxation test on a single specimen always resulted in the greatest amount of relaxation (Fig. 4.2). The same result was observed when the resting time between relaxation tests was 0 hours, 1 hour, or 12 hours. These tests revealed that the strain history significantly effects the relaxation behavior of the CL and USL. Therefore, the second and third relaxation tests on these specimens were not considered to determine the relaxation behavior of these ligaments. After selecting only the first relaxation tests on the specimens tested multiple times and conducting 15 single relaxation tests on 15 samples, a total of 22 relaxation tests from 22 samples were obtained. Relaxation tests for 9 specimens taken from a single sow are presented in Fig The 9 relaxation tests in this figure correspond to the 9 stress-stretch curves. The end point for each of the stress-stretch curves is the starting point for each relaxation

29 Stress (MPa/MPa) Stress (MPa/MPa) Winston R. Becker Chapter 4. Results 20 a) c) b) 1 d) Time (s) Time (s) Figure 4.1: Stress-stretch curves and normalized relaxation curves of 9 samples taken from a single sow. The locations the specimens were taken from are indicated by the schematic in the top left of the top two graphs. The stress relaxation curves represent the same samples presented in the stress-stretch curves. In both the relaxation curves and the stress-stretch curves, the filled symbols represent data collected from stretching in the direction normal to the upper vagina/cervix while the unfilled symbols denote data collected from stretching in the direction transverse to the upper vagina/cervix.

30 Winston R. Becker Chapter 4. Results 21 test. For example, the endpoint of the stress-stetch curves denoted by green squares in Fig. 4.1a (λ 1 = 1.08 and λ 2 = 1.14) represents the constant stretch for the relaxation test denoted by green squares in Fig. 4.1b. Substantial stress relaxation was observed in all specimens tested. In general, the relaxation was very fast for the first few hundred seconds, and the samples continued to relax slowly through the end of the test. The specimens that relaxed a greater percentage, typically had greater rates of relaxation throughout the tests. A wide range of relaxation behavior was observed in the tests conducted. For example, some samples retained about 65% of the initial load at 3000 seconds (red squares in Fig. 4.1), while others maintained as little as 20% of the initial load at 3000 seconds (green squares in Fig. 4.1). The observed differences in relaxation rate and percent relaxation were weakly dependent on the magnitude of the stretch applied to the specimens. In general, relaxation was higher for specimens stretched to lower strain levels when compared to specimens stretched to higher strain levels. Interestingly, in all of the stress relaxation tests performed, the normalized stress relaxation was identical in both directions even if the elastic behavior was anisotropic. This can be seen in Fig. 4.1, where many of the relaxation curves in the two different orientations overlap (i.e. filled purple triangles overlap unfilled purple triangles). 4.3 Modeling Representative results of the model fitting to the biaxial data are presented for thirteen specimens taken from three sows in Fig The model was able to capture the elastic

31 Normalized Stress (MPa/MPa) Normalized Stress (MPa/MPa) Winston R. Becker Chapter 4. Results 22 a) 1 b) Time (s) Time (s) Figure 4.2: Stress relaxation response for two specimens. Three tests were performed on each specimen with a 1 hour resting period between each test. a) The specimen was stretch to λ 1 = and λ 2 = during the first relaxation test, λ 1 = and λ 2 = during the second relaxation test, and λ 1 = and λ 2 = during the third relaxation test. b) The specimen was stretch to λ 1 = and λ 2 = during the first relaxation test, λ 1 = and λ 2 = during the second relaxation test, and λ 1 = and λ 2 = during the third relaxation test.

32 Stress (MPa) Normalized Stress (MPa/MPa) Stress (MPa) Normalized Stress (MPa/MPa) Stress (MPa) Normalized Stress (MPa/MPa) Winston R. Becker Chapter 4. Results 23 a) 1.6 b) Stretch Time (s) c) 1.6 d) Stretch Time (s) e) 2.5 f) Stretch Time (s) Figure 4.3: Continuum model fit to experimental data of 13 specimens taken from 3 sows. Filled symbols represent data collected from stretching in the direction normal to the upper vagina/cervix, the unfilled symbols denote data collected from stretching in the direction transverse to the upper vagina/cervix, and the solid lines indicate the constitutive model.

33 α2 α1 Winston R. Becker Chapter 4. Results 24 Parameter Mean Standard Deviation c 1 (Pa) c 2 (Pa) c c 4 (Pa) c Table 4.1: Mean and standard deviation of the parameters of all 22 specimens considered in this study. Figure 4.4: Dependency of α 1 (Left) and α 2 (Right), which are measures of the percent relaxation of the specimens, on the first invariant. The lines represent the functions α 1 (I 1 (τ)) = a 1 e a 2(I 1 3) and α 2 (I 1 (τ)) = a 3 e a 4(I 1 3) with a 1 = 1.078, a 2 = 2.303, a 3 = 2.175, and a 4 =

34 Winston R. Becker Chapter 4. Results 25 and stress relaxation behavior of the ligaments. The mean and standard deviation of the elastic parameters computed for 22 data sets are presented in Table 4.1. In many cases, the elastic constant c 1 in Eqs , which represents the contribution of the ground substance, approached zero when fitting the data. This suggests that the ground substance had a relatively small mechanical contribution when compared to the fibrous components of the tissues. The elastic constants c 2 and c 3, which account for the contribution of the fibers normal to the upper vagina/cervix, were typically larger than the constants c 4 and c 5, which account for the contribution of the fibers transverse to the upper vagina/cervix. This is a result of the specimens being stiffer in the direction normal to the upper vagina/cervix. Differences in these parameters allowed the model to capture the anisotropic behavior of the specimens. The values of α 1 and α 2 versus I 1 are presented in Fig In general, the magnitude of α 1 and α 2 were higher for lower levels of stretch. Incorporating the dependency of α 1 and α 2 on I 1 allowed the model to capture the strain dependent relaxation behavior of the tissue. When the alternative approach of fitting the model to all 22 data sets simultaneously was considered, the model predicts the general trend in the relaxation, but cannot fit every data set because of the noise in the data. The predicted stress relaxation behavior of ligaments stretched to different initial stretches when the model parameters are determined with this method is presented in Fig. 4.5.

35 Normalized Stress (MPa/MPa) Winston R. Becker Chapter 4. Results lx=ly=1.05 lx=ly=1.1 lx=ly=1.15 lx=ly= Time (s) Figure 4.5: Predicted stress relaxation behavior of ligaments stretched to different levels based on fitting the model to all 22 data sets simultaneously. The predictions were computed by fitting the constitutive model to all 22 sets of relaxation data simultaneously.

36 Chapter 5 Discussion 5.1 Experimental Results The anisotropic mechanical response with the direction normal to the upper vagina/cervix being stiffer is not surprising considering that the majority of the physiological loading on these ligaments likely occurs in this direction. It is noted that the anisotropy was less pronounced for a few of the specimens in this study. This suggests that in some cases the collagen fibers may be less organized, or that there may be a similar quantity of collagen fibers oriented in each direction. The viscoelastic tests on these tissues revealed the amount that they relax, if they relax different amounts in each direction, and if the stress relaxation is nonlinear (i.e. strain-dependent). Significant stress relaxation was observed in this study, which is similar to the results of the 27

37 Winston R. Becker Chapter 5. Discussion 28 previous stress relaxation tests. In the final step of Vardy et al s (2005) representative USL load response, the stress relaxed to about 70% of its initial value at just over 2000 seconds. Although this was after relaxing multiple times, the percent relaxation is similar to the relaxation observed in some of the tests in this study. The relaxation tests also showed that, even though the elastic behavior was anisotropic, the normalized stress relaxation response was always the same in both directions. It is unclear why this was observed; however, it may be because all of the individual collagen fibers relaxed approximately the same amount. While the amount of collagen fibers oriented in each direction is likely different, the amount that each individual fiber relaxes may be the same. In this case, normalizing the relaxation by dividing by the initial stress would normalize each direction for the concentration of collagen fibers and result in identical normalized stress relaxation in each direction. The idea that relaxation takes place within individual collagen fibers has been speculated by previous studies that have shown that the spacing and orientation of collagen fibers do not change during stress relaxation (27). Alternatively, the observed relaxation behavior may be due to relaxation taking place within the isotropic matrix rather than within the collagen fibers. If this is the case the model would need to be reformulated to include a relaxation term in the term accounting for the ground substance. This mechanism may be less likely because collagen is stiffer than the other components of the tissue, and likely supports more force than the matrix of the tissue. Considering that the stress has decreased to as little as 20% of the initial stress, it is likely that the force supported by the collagen fibers must decrease for this much relaxation to occur.

38 Winston R. Becker Chapter 5. Discussion 29 In order to determine the dependeny of relaxation on the initial strain, multiple stress relaxation tests were first completed on 7 of the specimens. Performing multiple tests on one sample overcomes the problems with identifying a trend that may result from the variability between different samples. Furthermore, stress relaxation tests at multiple stretch levels have previously been performed for this ligament (35). The problem with this method is that the tissue has a mechanical history, which can affect the specimens viscoelastic properties (34). In this study, the order of the stretch levels was varied to determine if the same trend in the normalized stress relaxation would be observed regardless of which stretch level was applied first. This revealed that the mechanical history is important for determining the relaxation behavior, which aided in the design and interpretation of this study. One limiting factor that may have affected the relaxation results in these experiments is that the hydration level of the tissue cannot be controlled. The longer each test and resting period was, the longer the specimen remained submerged in fluid, which could potentially alter the mechanical properties. It is possible that the change in hydration over this time period affected the result as much or more than the mechanical history. In particular, if a resting period of 12 hours is actually needed so that the mechanical history does not effect the relaxation results, it may be necessary to find a different way to preserve the specimen rather than submerging it in a bath during the resting period. After determining that multiple stress relaxation tests could not be performed on individual samples, single stress relaxation tests on different samples were compared to determine the relationship between relaxation rate and the initial strain. The rate of stress relaxation de-

39 Winston R. Becker Chapter 5. Discussion 30 creased with increasing strains, which is consistent with previous studies on other ligaments, such as the medial collateral ligament (14; 25). In each of those studies specimens stretched to lower initial strains relaxed more than specimens stretched to higher initial strains. The underlying cause of this trend in stress relaxation behavior is unknown (14); however, it has been speculated that it may result from water loss (25). Alternatively, others have speculated that relaxation may occur within individual collagen fibers (27), suggesting that this trend may have to do with molecular behavior within individual collagen fibers. Tests aimed at understanding the molecular changes that occur during stress relaxation in these tissues could advance our understanding of this experimental result. Understanding the mechanical behavior of these tissues is especially important to disease pathology because many other researchers have noted that significant remodeling takes place in these tissues. In particular, collagen fibers become disorganized under certain physiologic conditions, especially during the development of PFDs. For example, CL samples from 8/10 women with uterine prolapse had loosely arranged connective tissue and less dense fiber bundles (31), and samples of the USL had significantly reduced organization in individuals who underwent hysterectomies (3). Other studies have determined that diminished collagen levels are associated with the development of PFDs (2). These studies suggest that the mechanical properties of the tissues likely change as the tissue becomes diseased, which should be examined in future studies. One important limitation that must be considered in the context of the experimental results is the limited information available for the history of the sows. Because they were obtained

40 Winston R. Becker Chapter 5. Discussion 31 from a slaughterhouse, we do not have detailed information on their health and history. Nevertheless, within the same sow, we are confident of our findings. 5.2 Modeling After collecting data on the elastic and relaxation behavior of the CL and USL, a threedimensional finite strain nonlinear viscoelastic constitutive relation based on the Pipkin- Rogers integral series was formulated to model their mechanical behavior. The model was capable of describing the elastic and relaxation responses quite well. A large amount of variability was observed in the elastic parameters, c 1 -c 5. For c 1, the values of the mean and standard deviations are not representative of the overall data set. In all but three cases, the value of c 1 was much less than one meaning that the term containing c 1 did not contribute to the relaxation function. In the three cases where the value of c 1 was larger, the term containing c 1 likely had a slightly more significant contribution, but it was still much less than the contribution of the other components. While the values of c 2 -c 5 were all approximately the same order of magnitude, the large standard deviations are representative of the significant variability in the shape of the stress-strain curves. The function r 1, which was selected to be identical in both fiber directions, depends only on the first invariant because the stress relaxation response was isotropic for all ligaments tested. If the parameters depending on the strain invariant, I 1, in the relaxation function are taken to be constant, our model collapses to quasi-linear viscoelasticity. Simplifying the model to one

41 Winston R. Becker Chapter 5. Discussion 32 dimension would result in a relaxation function similar in form to the nonlinear viscoelastic constitutive model proposed by Pryse el at (26). In particular, they employed a functional form, A i (ɛ)e t/τi(ɛ), that is similar to the form of r 1 (I 1 (τ), t τ) (Eq. 3.11) to account for changes in relaxation due to the magnitude of the initial strain. The model presented here is highly versatile and could be employed in the modeling of many biological tissues. In this case, it was used to capture the behavior of a tissue with two primary collagen fibers orientations. However, the model could easily be extended to model the behavior of other collagenous tissues with fibers embedded at multiple orientations. Other extensions, such as altering the relaxation function to account for the active mechanics of biological tissues, could be done following the work of previous studies (19; 33), and may allow this model to provide additional information. This extension is particularly relevant when studying pelvic floor ligaments because they include a significant amount of smooth muscle and their smooth muscle content changes during disease development. The modeling of the CL and USL remains limited because only biaxial experimental data can be collected. While biaxial mechanical tests are useful for characterizing the anisotropic behavior of soft tissues, it is not possible to fully characterize the three-dimensional elastic behavior of anisotropic materials with biaxial mechanical tests because the mechanical behavior in the third direction is not characterized (15). In the tests performed, it was also impossible to completely prevent a small amount of shear from occurring. Finally, the test could be improved if we could determine the orientation of the fibers within the specimen prior to and during testing. This would aid us in explaining why we see anisotropy in the

42 Winston R. Becker Chapter 5. Discussion 33 mechanical response.

43 Chapter 6 Conclusions 6.1 Summary and Future Directions This work presents the first biaxial elastic and viscoelastic tests on the CL and USL. The results revealed that the specimens are anisotropic, and stiffer in the direction normal to the upper vagina/cervix. The normalized relaxation behavior of each specimen was the same in both directions and the amount of relaxation was dependent on the initial strain. Specimens stretched to lower initial strains relaxed more than specimens stretched to higher strains. A novel three dimensional constitutive model for multiple fibrous groups embedded in a matrix that accounts for anisotropy and strain dependent stress relaxation was formulated to describe the experimental observations. It was assumed that the specimens were composed of two families of fibers oriented perpendicular to each other and that the 34

44 Winston R. Becker Chapter 6. Conclusions 35 tissues were incompressible. The resulting model effectively fit the elastic and viscoelastic data collected through mechanical experiemnts. The results of this study provide detailed mechanical information about pelvic floor ligaments, which should aid in the development of new treatment methods for PFDs in the future. This study provides the basis for additional work on pelvic floor connective tissues. First, the condition of the sow and the effects of parity on the mechanical properties of pelvic floor connective tissue should be carefully studied. Second, this study would benefit from collecting the structural information of the tissues. In particular, the orientations of the fibrous components should be quantified before and during mechanical testing. This would reveal the unloaded molecular configuration as well as identify how this configuration changes during loading and relaxation. Performing these studies would aid in the connection of many of the mechanical observations to their underlying structural causes. For example, this may explain why the elastic response was isotropic while the relaxation response was anisotropic. Finally, future studies should consider the active mechanical properties of these tissues. Previous histological examinations have revealed that a significant amount of smooth muscle is present in these ligaments. Computing the active mechanical properties of these tissues would reveal the physiological role of this smooth muscle.

45 Bibliography [1] Baek, S., Wells, P., Rajagopal, K., Humphrey, J.: Heat-induced changes in the finite strain viscoelastic behavior of a collaagenous tissue. Journal of Biomechanical Engineering 127(4), (2005) [2] Campeau, L., Gorbachinsky, I., Badlani, G., Andersson, K.: Pelvic floor disorders: linking genetic risk factors to biochemical changes. BJU International 108(8), (2011) [3] Cole, E., Leu, P., Gomelsky, A., Revelo, P., Shappell, H., S, H.M., Dmochowski, R.: Histopathological evaluation of the uterosacral ligament: is this a dependable structure for pelvic reconstruction? BJU International 97(2), (2006) [4] Davis, F., De Vita, R.: A nonlinear constitutive model for stress relaxation in ligaments and tendons. Annals of Biomedical Engineering 40(12), (2012) [5] Davis, F., De Vita, R.: A three-dimensional constitutive model for the stress relaxation of articular ligaments. Biomechanics and Modeling in Mechanobiology pp (2013) 36

Biaxial mechanical properties of swine uterosacral and cardinal ligaments

Biaxial mechanical properties of swine uterosacral and cardinal ligaments Biomech Model Mechanobiol DOI 10.1007/s10237-014-0621-5 ORIGINL PPER Biaxial mechanical properties of swine uterosacral and cardinal ligaments Winston R. Becker Raffaella De Vita Received: 12 June 2014

More information

Lectures on. Constitutive Modelling of Arteries. Ray Ogden

Lectures on. Constitutive Modelling of Arteries. Ray Ogden Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics

More information

Full-field measurements and identification for biological soft tissues: application to arteries in vitro

Full-field measurements and identification for biological soft tissues: application to arteries in vitro Centre for Health Engineering CNRS UMR 5146 INSERM IFR 143 Prof. Stéphane Avril Full-field measurements and identification for biological soft tissues: application to arteries in vitro using single-gage

More information

AN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS EFFECT IN ARTERIAL TISSUE

AN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS EFFECT IN ARTERIAL TISSUE XI International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XI E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) AN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS

More information

Simple Shear Testing of Parallel-Fibered Planar Soft Tissues

Simple Shear Testing of Parallel-Fibered Planar Soft Tissues John C. Gardiner Jeffrey A. Weiss e-mail: jeff.weiss@utah.edu Department of Bioengineering, The University of Utah, 50 South Central Campus Drive #2480, Salt Lake City, UT 84112 Simple Shear Testing of

More information

Adaptation of a Rabbit Myocardium Material Model for Use in a Canine Left Ventricle Simulation Study

Adaptation of a Rabbit Myocardium Material Model for Use in a Canine Left Ventricle Simulation Study Adaptation of a Rabbit Myocardium Material Model for Use in a Canine Left Ventricle Simulation Study Matthew G. Doyle e-mail: mdoyl009@uottawa.ca Stavros Tavoularis 1 e-mail: stavros.tavoularis@uottawa.ca

More information

Linear viscoelastic behavior

Linear viscoelastic behavior Harvard-MIT Division of Health Sciences and Technology HST.523J: Cell-Matrix Mechanics Prof. Ioannis Yannas Linear viscoelastic behavior 1. The constitutive equation depends on load history. 2. Diagnostic

More information

NOTE To improve the TA to class ratio, two groups will perform this experiment at one time. Group times are listed on the class website.

NOTE To improve the TA to class ratio, two groups will perform this experiment at one time. Group times are listed on the class website. Laboratory 3: Viscoelastic Characterization of Tendon using Harmonic Oscillations November 8/9, 2006 BIOEN 5201 Introduction to Biomechanics Instructor: Jeff Weiss TA: Heath Henninger Lab Quiz: A 10 point

More information

in this web service Cambridge University Press

in this web service Cambridge University Press CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.

More information

Effects of Aging on the Mechanical Behavior of Human Arteries in Large Deformations

Effects of Aging on the Mechanical Behavior of Human Arteries in Large Deformations International Academic Institute for Science and Technology International Academic Journal of Science and Engineering Vol. 3, No. 5, 2016, pp. 57-67. ISSN 2454-3896 International Academic Journal of Science

More information

Modelling Anisotropic, Hyperelastic Materials in ABAQUS

Modelling Anisotropic, Hyperelastic Materials in ABAQUS Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,

More information

Development of a reduced relaxation function and comparison of stress relaxation for anatomically paired tendons

Development of a reduced relaxation function and comparison of stress relaxation for anatomically paired tendons Journal of Mechanical Science and Technology 23 (2009) 1893~1898 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-009-0504-3 Development of a reduced

More information

Lecture 8 Viscoelasticity and Deformation

Lecture 8 Viscoelasticity and Deformation Read: pg 130 168 (rest of Chpt. 4) 1 Poisson s Ratio, µ (pg. 115) Ratio of the strain in the direction perpendicular to the applied force to the strain in the direction of the applied force. For uniaxial

More information

Discontinuous Galerkin methods for nonlinear elasticity

Discontinuous Galerkin methods for nonlinear elasticity Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity

More information

Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis. Supplementary Document

Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis. Supplementary Document Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis Supplementary Document Shan Yang, Vladimir Jojic, Jun Lian, Ronald Chen, Hongtu Zhu, Ming C.

More information

A new strain energy function for the hyperelastic modelling of ligaments and tendons

A new strain energy function for the hyperelastic modelling of ligaments and tendons A new strain energy function for the hyperelastic modelling of ligaments and tendons University of Manchester BMC-BAMC 2015 Anterior cruciate ligament reconstruction surgery Ligament and tendon hierarchical

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,

More information

UNIVERSITY OF CALGARY. Viscoelasticity of Articular Cartilage and Ligament: Constitutive Modeling and Experiments. Sahand Ahsanizadeh A THESIS

UNIVERSITY OF CALGARY. Viscoelasticity of Articular Cartilage and Ligament: Constitutive Modeling and Experiments. Sahand Ahsanizadeh A THESIS UNIVERSITY OF CALGARY Viscoelasticity of Articular Cartilage and Ligament: Constitutive Modeling and Experiments by Sahand Ahsanizadeh A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT

More information

Testing and Analysis

Testing and Analysis Testing and Analysis Testing Elastomers for Hyperelastic Material Models in Finite Element Analysis 2.6 2.4 2.2 2.0 1.8 1.6 1.4 Biaxial Extension Simple Tension Figure 1, A Typical Final Data Set for Input

More information

Constitutive Modeling of Biological Soft Tissues

Constitutive Modeling of Biological Soft Tissues Constitutive Modeling of Biological Soft Tissues Attila P. Nagy 1, David J. Benson 1, Vikas Kaul 2, Mark Palmer 2 1 Livermore Software Technology Corporation, Livermore, CA 94551, USA 2 Medtronic plc,

More information

Plane Strain Test for Metal Sheet Characterization

Plane Strain Test for Metal Sheet Characterization Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue

More information

Predeformation and frequency-dependence : Experiment and FE analysis

Predeformation and frequency-dependence : Experiment and FE analysis Predeformation and frequency-dependence : Experiment and FE analysis Nidhal Jridi 1,2,*, Michelle Salvia 2, Adel Hamdi 1, Olivier Bareille 2, Makrem Arfaoui 1, Mohammed Ichchou 2, Jalel Ben Abdallah 1

More information

Comparative Study of Variation of Mooney- Rivlin Hyperelastic Material Models under Uniaxial Tensile Loading

Comparative Study of Variation of Mooney- Rivlin Hyperelastic Material Models under Uniaxial Tensile Loading Comparative Study of Variation of Mooney- Rivlin Hyperelastic Material Models under Uniaxial Tensile Loading A. N. Jadhav 1, Dr. S.R. Bahulikar, N.H. Sapate 3 1 M Tech Design Engg, Mechanical Engineering,

More information

EXPERIMENTAL EVALUATION OF SHEAR STRENGTH OF WOVEN WEBBINGS

EXPERIMENTAL EVALUATION OF SHEAR STRENGTH OF WOVEN WEBBINGS EXPERIMENTAL EVALUATION OF SHEAR STRENGTH OF WOVEN WEBBINGS Kevin L. Peil +, Ever J. Barbero +, Eduardo M. Sosa* + Department of Mechanical and Aerospace Engineering, West Virginia University (WVU), Morgantown,

More information

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS 1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal

More information

Modeling of Fiber-Reinforced Membrane Materials Daniel Balzani. (Acknowledgement: Anna Zahn) Tasks Week 2 Winter term 2014

Modeling of Fiber-Reinforced Membrane Materials Daniel Balzani. (Acknowledgement: Anna Zahn) Tasks Week 2 Winter term 2014 Institute of echanics and Shell Structures Faculty Civil Engineering Chair of echanics odeling of Fiber-Reinforced embrane aterials OOC@TU9 Daniel Balani (Acknowledgement: Anna Zahn Tasks Week 2 Winter

More information

SIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS

SIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS Engineering MECHANICS, Vol. 18, 2011, No. 1, p. 23 32 23 SIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS Tomáš Lasota*, JiříBurša* This paper deals

More information

International Journal of Pure and Applied Mathematics Volume 58 No ,

International Journal of Pure and Applied Mathematics Volume 58 No , International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin

More information

2.1 Strain energy functions for incompressible materials

2.1 Strain energy functions for incompressible materials Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to

More information

FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Ever J. Barbero Department of Mechanical and Aerospace Engineering West Virginia University USA CRC Press Taylor &.Francis Group Boca Raton London New York

More information

Biomechanics. Soft Tissue Biomechanics

Biomechanics. Soft Tissue Biomechanics Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic

More information

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources

More information

Effect of off-axis cell orientation on mechanical properties in smooth muscle tissue

Effect of off-axis cell orientation on mechanical properties in smooth muscle tissue J. Biomedical Science and Engineering, 2011, 4, 10-17 doi:10.4236/jbise.2011.41002 Published Online January 2011 (http://www.scirp.org/journal/jbise/). Effect of off-axis cell orientation on mechanical

More information

CYTOLOGY & HISTOLOGY THE STUDY OF CELLS AND TISSUES

CYTOLOGY & HISTOLOGY THE STUDY OF CELLS AND TISSUES NAME: DATE: PARTNER: CYTOLOGY & HISTOLOGY THE STUDY OF CELLS AND TISSUES For ease of study, multicellular animals are often examined at various levels of structural organization. Starting from the most

More information

Constitutive model of brain tissue suitable for finite element analysis of surgical procedures

Constitutive model of brain tissue suitable for finite element analysis of surgical procedures Journal of Biomechanics 32 (1999 531 537 Technical Note Constitutive model of brain tissue suitable for finite element analysis of surgical procedures Karol Miller* Department of Mechanical and Materials

More information

MATERIAL MODELS FOR CRUMB RUBBER AND TDA. California State University, Chico

MATERIAL MODELS FOR CRUMB RUBBER AND TDA. California State University, Chico MATERIAL MODELS FOR CRUMB RUBBER AND TDA California State University, Chico Waste Tire Products for CE Applications Whole Tires Tire Shreds (TDA) Crumb Rubber/Tire Buffings Whole Tires TIRE DERIVED AGGREGATE

More information

MSC Elastomers Seminar Some Things About Elastomers

MSC Elastomers Seminar Some Things About Elastomers MSC Elastomers Seminar Some Things About Elastomers Kurt Miller, Axel Products, Inc. www.axelproducts.com Visit us at: axelproducts.com 2 Your Presenter Kurt Miller Founded Axel Products 1994 Instron Corporation,

More information

Numerical models of cerebral aneurysm biomechanics

Numerical models of cerebral aneurysm biomechanics ANEURISK PROJECT: LABORATORY OF BIOLOGICAL STRUCTURE MECHANICS www.labsmech.polimi.it Numerical models of cerebral aneurysm biomechanics Laura Socci, Dario Gastaldi, Giancarlo Pennati, Pasquale Vena Gabriele

More information

What is a Strain Gauge? Strain Gauge. Schematic View Of Strain Gauge

What is a Strain Gauge? Strain Gauge. Schematic View Of Strain Gauge ( ) : 1391-92 92 What is Strain? Strain is the amount of deformation of a body due to an applied force. More specifically, strain (ε) is defined as the fractional change in length. Strain can be positive

More information

Pointwise identification of elastic properties in nonlinear heterogeneous membranes, and application to soft tissues

Pointwise identification of elastic properties in nonlinear heterogeneous membranes, and application to soft tissues University of Iowa Iowa Research Online Theses and Dissertations Spring 2009 Pointwise identification of elastic properties in nonlinear heterogeneous membranes, and application to soft tissues Xuefeng

More information

Mathematical Modelling of Biological Soft Tissues

Mathematical Modelling of Biological Soft Tissues Mathematical Modelling of Biological Soft Tissues 1., William J. Parnell 1., Barbara Lynch 2., Hazel R.C. Screen 3. and I. David Abrahams 4. 1. University of Manchester 2. Ecole Polytechnique 3. Queen

More information

STRUCTURAL OPTIMIZATION OF A MATERIAL EXHIBITING NEGATIVE STIFFNESS

STRUCTURAL OPTIMIZATION OF A MATERIAL EXHIBITING NEGATIVE STIFFNESS International Conference on Engineering Vibration Ljubljana, Slovenia, 7-0 September 05 STRUCTURAL OPTIMIZATION OF A MATERIAL EXHIBITING NEGATIVE STIFFNESS Zuzana Dimitrovová*, Jan Heczo, Helder C. Rodrigues

More information

EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD

EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University

More information

Lecture 8 Viscoelasticity and Deformation

Lecture 8 Viscoelasticity and Deformation HW#5 Due 2/13 (Friday) Lab #1 Due 2/18 (Next Wednesday) For Friday Read: pg 130 168 (rest of Chpt. 4) 1 Poisson s Ratio, μ (pg. 115) Ratio of the strain in the direction perpendicular to the applied force

More information

Wave Propagation Through Soft Tissue Matter

Wave Propagation Through Soft Tissue Matter Wave Propagation Through Soft Tissue Matter Marcelo Valdez and Bala Balachandran Center for Energetic Concepts Development Department of Mechanical Engineering University of Maryland College Park, MD 20742-3035

More information

AE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF HOOP STRESSES FOR A THIN-WALL PRESSURE VESSEL

AE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF HOOP STRESSES FOR A THIN-WALL PRESSURE VESSEL Objective AE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF OOP STRESSES FOR A TIN-WA PRESSURE VESSE This experiment will allow you to investigate hoop and axial stress/strain relations

More information

On the Constitutive Models for Heart Valve Leaflet Mechanics

On the Constitutive Models for Heart Valve Leaflet Mechanics Cardiovascular Engineering: An International Journal, Vol. 5, No. 1, March 005 ( C 005) DOI: 10.1007/s10558-005-307-x On the Constitutive Models for Heart Valve Leaflet Mechanics ELI J. WEINBERG and MOHAMMAD

More information

Photoelastic Experiments with a Transmitted-Light Polariscope

Photoelastic Experiments with a Transmitted-Light Polariscope Photoelastic Experiments with a Transmitted-Light Polariscope Objectives: To become familiar with a Transmitted-Light Polariscope To physically see the effects of loading various specimens and their unique

More information

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQF LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQF LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS EDEXCEL NATIONAL CERTIICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQ LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS 1. Be able to determine the effects of loading in static engineering

More information

Lab Exercise #5: Tension and Bending with Strain Gages

Lab Exercise #5: Tension and Bending with Strain Gages Lab Exercise #5: Tension and Bending with Strain Gages Pre-lab assignment: Yes No Goals: 1. To evaluate tension and bending stress models and Hooke s Law. a. σ = Mc/I and σ = P/A 2. To determine material

More information

Size Effects In the Crushing of Honeycomb Structures

Size Effects In the Crushing of Honeycomb Structures 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 19-22 April 2004, Palm Springs, California AIAA 2004-1640 Size Effects In the Crushing of Honeycomb Structures Erik C.

More information

Numerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35

Numerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35 Numerical Model of the Influence of Shear Stress on the Adaptation of a Blood Vessel BMT 03-35 Mirjam Yvonne van Leeuwen Supervisor: Dr. Ir. M.C.M. Rutten Ir. N.J.B. Driessen TUE Eindhoven, The Netherlands

More information

Deflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis

Deflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis Rotating Machinery, 10(4): 283 291, 2004 Copyright c Taylor & Francis Inc. ISSN: 1023-621X print / 1542-3034 online DOI: 10.1080/10236210490447728 Deflections and Strains in Cracked Shafts due to Rotating

More information

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity 63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain

More information

MECHANICAL CHARACTERIZATION OF BRAIN TISSUE

MECHANICAL CHARACTERIZATION OF BRAIN TISSUE ROLE OF MOISTURE CONTENT IN MECHANICAL CHARACTERIZATION OF BRAIN TISSUE HENRY W. HASLACH, JR. DEPARTMENT OF MECHANICAL ENGINEERING CENTER for ENERGETICS CONCEPTS DEVELOPMENT UNIVERSITY OF MARYLAND COLLEGE

More information

A FAILURE CRITERION FOR POLYMERS AND SOFT BIOLOGICAL MATERIALS

A FAILURE CRITERION FOR POLYMERS AND SOFT BIOLOGICAL MATERIALS Material Technology A FALURE CRTERON FOR POLYMERS AND SOFT BOLOGCAL MATERALS Authors: William W. Feng John O. Hallquist Livermore Software Technology Corp. 7374 Las Positas Road Livermore, CA 94550 USA

More information

Understanding Frequency Domain Viscoelasticity in Abaqus

Understanding Frequency Domain Viscoelasticity in Abaqus Paper # 12 Understanding Frequency Domain Viscoelasticity in Abaqus By Saurabh Bahuguna, Randy Marlow*, and Tod Dalrymple Dassault Systèmes Simulia Corp., Great Lakes Region Presented at the Fall 172 nd

More information

DEFORMATION THEORY OF PLASTICITY

DEFORMATION THEORY OF PLASTICITY DEFORMATION THEORY OF PLASTICITY ROBERT M. JONES Professor Emeritus of Engineering Science and Mechanics Virginia Polytechnic Institute and State University Blacksburg, Virginia 240610219 Bull Ridge Publishing

More information

FEM model of pneumatic spring assembly

FEM model of pneumatic spring assembly FEM model of pneumatic spring assembly Tien Tran Xuan 1, David Cirkl 2 Department of Applied Mechanics, Faculty of Mechanical Engineering, Technical University of Liberec, Liberec, Czech Republic 1 Corresponding

More information

Characterization of Viscoelastic Behaviors in Bovine Pulmonary Arterial Tissue

Characterization of Viscoelastic Behaviors in Bovine Pulmonary Arterial Tissue University of Colorado, Boulder CU Scholar Mechanical Engineering Graduate Theses & Dissertations Mechanical Engineering Spring 1-1-2011 Characterization of Viscoelastic Behaviors in Bovine Pulmonary Arterial

More information

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation

More information

A novel elastometer for soft tissue

A novel elastometer for soft tissue Computational Methods and Experiments in Material Characterisation II 111 A novel elastometer for soft tissue S. M. Harrison 1, M. B. Bush 1 & P. Petros 1,2 1 University of Western Australia 2 Royal Perth

More information

Constitutive Relations

Constitutive Relations Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field

More information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the

More information

Fatigue in osteoporotic human trabecular bone Comparison between genders

Fatigue in osteoporotic human trabecular bone Comparison between genders loading to the load corresponding to a normalized stress level (Δσ/E 0 ). Fatigue in osteoporotic human trabecular bone Comparison between genders André do Carmo Saraiva Abstract In the present work, the

More information

06 - kinematic equations kinematic equations

06 - kinematic equations kinematic equations 06 - - 06-1 continuum mechancis continuum mechanics is a branch of physics (specifically mechanics) that deals with continuous matter. the fact that matter is made of atoms and that it commonly has some

More information

The nonlinear mechanical properties of engineered soft biological tissues determined by finite spherical indentation

The nonlinear mechanical properties of engineered soft biological tissues determined by finite spherical indentation The nonlinear mechanical properties of engineered soft biological tissues determined by finite spherical indentation BMTE 07.12 Bart van Dijk March 2007 Internship at Eindhoven University of Technology

More information

Material testing method for understanding rubber s mechanical characteristics Quasi-Static Tension Test

Material testing method for understanding rubber s mechanical characteristics Quasi-Static Tension Test Material testing method for understanding rubber s mechanical characteristics Quasi-Static ension est 1 Introduction In this working group, we conduct tension tests of rubber materials to characterize

More information

Deflections and Strains in Cracked Shafts Due to Rotating Loads: A Numerical and Experimental Analysis

Deflections and Strains in Cracked Shafts Due to Rotating Loads: A Numerical and Experimental Analysis International Journal of Rotating Machinery, 9: 303 311, 2003 Copyright c Taylor & Francis Inc. ISSN: 1023-621X DOI: 10.1080/10236210390147416 Deflections and Strains in Cracked Shafts Due to Rotating

More information

Linear Constitutive Relations in Isotropic Finite Viscoelasticity

Linear Constitutive Relations in Isotropic Finite Viscoelasticity Journal of Elasticity 55: 73 77, 1999. 1999 Kluwer Academic Publishers. Printed in the Netherlands. 73 Linear Constitutive Relations in Isotropic Finite Viscoelasticity R.C. BATRA and JANG-HORNG YU Department

More information

Chapter 7. Highlights:

Chapter 7. Highlights: Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true

More information

A Review On Methodology Of Material Characterization And Finite Element Modelling Of Rubber-Like Materials

A Review On Methodology Of Material Characterization And Finite Element Modelling Of Rubber-Like Materials IOSR Journal of Engineering (IOSRJEN) ISSN (e): 50-301, ISSN (p): 78-8719 PP 06-10 www.iosrjen.org A Review On Methodology Of Material Characterization And Finite Element Modelling Of Rubber-Like Materials

More information

Constitutive Equations

Constitutive Equations Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module

More information

SOME TOPICS FROM CONTINUUM MECHANICS RELATED TO BRAIN NEURO-MECHANICS

SOME TOPICS FROM CONTINUUM MECHANICS RELATED TO BRAIN NEURO-MECHANICS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, Series B Volume 3, Number 1, Pages 1 19 c 2012 Institute for Scientific Computing and Information SOME TOPICS FROM CONTINUUM MECHANICS RELATED

More information

NATIONAL REVIEW COURSE. Cells, Tissues, and Membranes

NATIONAL REVIEW COURSE. Cells, Tissues, and Membranes NATIONAL REVIEW COURSE Cells, Tissues, and Membranes I. Cell Types A. Prokaryote bacteria cells; a cell that does not have a nucleus in which to store its genetic material. B. Eukaryote plant or animal

More information

Lecture 7 Constitutive Behavior of Asphalt Concrete

Lecture 7 Constitutive Behavior of Asphalt Concrete Lecture 7 Constitutive Behavior of Asphalt Concrete What is a Constitutive Model? A constitutive model or constitutive equation is a relation between two physical quantities that is specific to a material

More information

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the

More information

ME 2570 MECHANICS OF MATERIALS

ME 2570 MECHANICS OF MATERIALS ME 2570 MECHANICS OF MATERIALS Chapter III. Mechanical Properties of Materials 1 Tension and Compression Test The strength of a material depends on its ability to sustain a load without undue deformation

More information

NET v1.0: A Framework to Simulate Permanent Damage in Elastomers under Quasi-static Deformations

NET v1.0: A Framework to Simulate Permanent Damage in Elastomers under Quasi-static Deformations NET v1.0: A Framework to Simulate Permanent Damage in Elastomers under Quasi-static Deformations Vahid Morovati, Roozbeh Dargazany morovati@msu.edu, roozbeh@msu.edu Michigan State University Department

More information

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3

More information

A fibre-reinforced poroviscoelastic finite element model for the Achilles tendon

A fibre-reinforced poroviscoelastic finite element model for the Achilles tendon A fibre-reinforced poroviscoelastic finite element model for the Achilles tendon Anna Gustafsson 2014 Master s Thesis Faculty of Engineering, LTH Department of Biomedical Engineering Supervisors: Hanifeh

More information

arxiv:q-bio/ v1 [q-bio.qm] 30 Nov 2003

arxiv:q-bio/ v1 [q-bio.qm] 30 Nov 2003 Material forces in the context of biotissue remodelling arxiv:q-bio/0312002v1 [q-bio.qm] 30 Nov 2003 K. Garikipati, H. Narayanan, E. M. Arruda, K. Grosh, S. Calve University of Michigan, Ann Arbor, Michigan

More information

Nonlinear Elasticity, Anisotropy, Material Stability and Residual stresses in Soft Tissue

Nonlinear Elasticity, Anisotropy, Material Stability and Residual stresses in Soft Tissue Nonlinear Elasticity, Anisotropy, Material Stability and Residual stresses in Soft Tissue R.W. Ogden Department of Mathematics, University of Glasgow Glasgow G 8QW, UK Email: rwo@maths.gla.ac.uk Home Page:

More information

ELASTOPLASTICITY THEORY by V. A. Lubarda

ELASTOPLASTICITY THEORY by V. A. Lubarda ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and

More information

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2 Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force

More information

GG303 Lab 12 11/7/18 1

GG303 Lab 12 11/7/18 1 GG303 Lab 12 11/7/18 1 DEFORMATION AROUND A HOLE This lab has two main objectives. The first is to develop insight into the displacement, stress, and strain fields around a hole in a sheet under an approximately

More information

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS

LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds

More information

Characterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam

Characterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam Characterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam Qusai Hatem Jebur 1,a, Philip Harrrison 1,b, Zaoyang Guo,c, Gerlind Schubert 1,d & Vincent Navez 3,e 1 School of Engineering, University

More information

Mechanical Properties of Polymers. Scope. MSE 383, Unit 3-1. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.

Mechanical Properties of Polymers. Scope. MSE 383, Unit 3-1. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Mechanical Properties of Polymers Scope MSE 383, Unit 3-1 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Structure - mechanical properties relations Time-dependent mechanical

More information

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 Samantha Ramirez TENSION AND COMPRESSION TESTS Tension and compression tests are used primarily to determine the relationship between σ avg and ε avg in any material.

More information

Computational Biomechanics Lecture 2: Basic Mechanics 2. Ulli Simon, Frank Niemeyer, Martin Pietsch

Computational Biomechanics Lecture 2: Basic Mechanics 2. Ulli Simon, Frank Niemeyer, Martin Pietsch Computational Biomechanics 016 Lecture : Basic Mechanics Ulli Simon, Frank Niemeyer, Martin Pietsch Scientific Computing Centre Ulm, UZWR Ulm University Contents .7 Static Equilibrium Important: Free-body

More information

Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur

Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 9 Table of Contents 1. Plasticity:... 3 1.1 Plastic Deformation,

More information

*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials

*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials 14 th International LS-DYNA Users Conference Session: Constitutive Modeling *MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials Jesper Karlsson 1, Mikael Schill 1, Johan Tryding

More information

Two problems in finite elasticity

Two problems in finite elasticity University of Wollongong Research Online University of Wollongong Thesis Collection 1954-2016 University of Wollongong Thesis Collections 2009 Two problems in finite elasticity Himanshuki Nilmini Padukka

More information

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business

More information

Computer Modeling in Bioengineering

Computer Modeling in Bioengineering Computer Modeling in Bioengineering Theoretical Background, Examples and Software Milos Kojic Harvard School of Public Health, USA University of Kragujevac, Serbia University of Texas Health Science Center

More information

3D constitutive modeling of the passive heart wall

3D constitutive modeling of the passive heart wall 3D constitutive modeling of the passive heart wall Thomas S.E. Eriksson Institute of Biomechanics - Graz University of Technology Department of Biophysics - Medical University of Graz SFB status seminar,

More information

MECHANICS LAB AM 317 EXP 3 BENDING STRESS IN A BEAM

MECHANICS LAB AM 317 EXP 3 BENDING STRESS IN A BEAM MECHANICS LAB AM 37 EXP 3 BENDING STRESS IN A BEAM I. OBJECTIVES I. To compare the experimentally determined stresses in a beam with those predicted from the simple beam theory (a.k.a. Euler-Bernoull beam

More information

Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec

Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling

More information