24.12.2013 Views

Lab 3: Stress, strain, stiffness, and time-dependent properties of ...

Lab 3: Stress, strain, stiffness, and time-dependent properties of ...

Lab 3: Stress, strain, stiffness, and time-dependent properties of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Lab</strong> 3: <strong>Stress</strong>, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong><br />

<strong>properties</strong> <strong>of</strong> biomaterials<br />

Bio427 — Biomechanics<br />

This lab explores experimental methods for quantifying the mechanical<br />

<strong>properties</strong> <strong>of</strong> biological materials (recall that material <strong>properties</strong><br />

are in<strong>dependent</strong> <strong>of</strong> structural size <strong>and</strong> shape). In addition we<br />

will take our first foray into measurements <strong>of</strong> the dynamic mechanical<br />

characteristics <strong>of</strong> natural systems. We will be using both plant <strong>and</strong><br />

animal material in this lab.<br />

Goals<br />

• Learn how to measure force <strong>and</strong> displacement so that you can<br />

construct a stress-<strong>strain</strong> curve<br />

• Underst<strong>and</strong> data acquisition using a computer <strong>and</strong> an Arduino<br />

• Examine the <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biological materials<br />

• Explore the concept <strong>of</strong> anisotropy<br />

• Explore resonance <strong>and</strong> damping in natural systems<br />

Conceptual Basis<br />

Last week we explored the concepts <strong>of</strong> stress (σ), <strong>strain</strong> (ɛ ), <strong>stiffness</strong><br />

(E, Young’s modulus) <strong>and</strong> strength (σ max ). These are all material<br />

<strong>properties</strong> <strong>and</strong> represent a biomechanical ways to characterize materials<br />

living systems. These physical quantities are summarized here:<br />

Figure 1: <strong>Stress</strong>-<strong>strain</strong> relationship<br />

for ligament: from Lee <strong>and</strong> Hyman,<br />

Modeling <strong>of</strong> failure mode in knee<br />

ligaments depending on the <strong>strain</strong> rate.<br />

BMC Musculoskeletal Disord. 2002 doi:<br />

10.1186/1471-2474-3-3<br />

σ = F A<br />

(1)<br />

ɛ = ∆L<br />

L o<br />

(2)<br />

E = dσ<br />

(3)<br />

dɛ<br />

where F is the applied force, A is the cross-sectional area over which<br />

that force operates, ∆L is the change in length in response to that<br />

force, <strong>and</strong> L o is the initial length <strong>of</strong> the specimen. Young’s modulus<br />

(E) is the slope <strong>of</strong> the stress-<strong>strain</strong> relationship which is constant for<br />

Hookean materials (that is, they obey Hooke’s Law). However, as we<br />

may learn in this lab, many biological materials are non-Hookean<br />

<strong>and</strong> the relationship between stress <strong>and</strong> <strong>strain</strong> is nonlinear. Moreover,<br />

biomaterials may also show anisotropic behavior – that is, the relationship<br />

between the stress <strong>and</strong> <strong>strain</strong> depends on the direction in


lab 3: stress, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biomaterials 2<br />

which the load is applied. For example, any biomaterial containing<br />

fibers with some average axial orientation will show some level <strong>of</strong><br />

anisotropy. Examples <strong>of</strong> biomaterials with containing fibers include<br />

skin, intestine, muscle epimysium, wood, plant stems <strong>and</strong> many others.<br />

In addition to nonlinear <strong>and</strong> anisotropic behaviors, many biomaterials<br />

are viscoelastic (as we have discussed in lecture), showing a<br />

<strong>time</strong>-dependence in their their response to loads. There are a host <strong>of</strong><br />

methods available for quantifying <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> including<br />

creep tests (constant force: isotonic loads), stress relaxation tests<br />

(constant length: isometric loads) as well as cyclic dynamic loads.<br />

One fascinating application <strong>of</strong> dynamic loads is called a free vibration<br />

test (colloquially called a "kaboing-wonga-wonga" test) in<br />

which you impulsively load a structure <strong>and</strong> observe the oscillations<br />

that follow that load. Those oscillations have a particular frequency<br />

for that structure (called the natural frequency) <strong>and</strong>, in the presence<br />

<strong>of</strong> any visco-elasticity, they damp out in <strong>time</strong>. The natural frequency<br />

for a given structure will depend on the <strong>stiffness</strong> <strong>and</strong> density <strong>of</strong> the<br />

material from which it is made: stiffer materials increase the natural<br />

frequency, but denser materials decrease the natural frequency.<br />

For a cylindrical beam <strong>of</strong> uniform section, there is a simple relationship<br />

that relates the natural frequency to the geometry <strong>of</strong> that beam<br />

<strong>and</strong> the <strong>stiffness</strong> <strong>and</strong> density <strong>of</strong> the material <strong>of</strong> which that beam is<br />

composed:<br />

ω n = 1.7<br />

√<br />

ER 2<br />

ρL 4 (4)<br />

where E is Young’s modulus for the material from which the cylinder<br />

is made, R is the beam radius, L is the beam length <strong>and</strong> ρ is the<br />

density <strong>of</strong> the beam. For nearly all biological materials, we take the<br />

density to be 1000 kg m −3 , which is the approximate density <strong>of</strong> water.<br />

Thus, if you knew the length <strong>and</strong> radius <strong>of</strong> the beam <strong>and</strong> its density,<br />

you could observe its natural frequency <strong>and</strong> compute the <strong>stiffness</strong> <strong>of</strong><br />

the material that constitutes the beam!<br />

Figure 2: On November 7, 1940 at<br />

11:00 am the Tacoma Narrows bridge<br />

was excited by wind vortices <strong>and</strong><br />

establish large amplitude vibrations.<br />

This vibrating beam failed <strong>and</strong> fell into<br />

the sound. Though no lives were lost, it<br />

was a narrow escape<br />

Measuring forces<br />

In the lab exercises you have done so far, you have used measurements<br />

<strong>of</strong> length <strong>and</strong> <strong>time</strong> together with mathematical models to<br />

investigate biomechanical performance <strong>of</strong> humans in jumping <strong>and</strong><br />

walking. In many cases, however, it is not possible to directly obtain<br />

the data we want, <strong>and</strong> in these cases we rely on transducers to<br />

convert signals from one type <strong>of</strong> energy to another. 1 If we wish to<br />

measure changes in force, which is a signal stored as mechanical<br />

1<br />

A common example <strong>of</strong> a transducer is<br />

the microphone in your cell phone: this<br />

is an electroacoustic transducer which<br />

converts sound into an electrical signal.<br />

Your voice produces propagating<br />

changes in air pressure, which cause the<br />

motion <strong>of</strong> a coil within a magnetic field<br />

inside the microphone. This results in<br />

a potential difference across the conductor,<br />

which constitutes an electrical<br />

signal which may be transmitted.


lab 3: stress, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biomaterials 3<br />

energy, we must convert this mechanical energy to electrical energy<br />

using a force transducer. These devices output a changing voltage<br />

proportional to the force that is applied to them. One <strong>of</strong> the most<br />

common elements used in force transducers is the foil <strong>strain</strong> gauge,<br />

which is what you will be using in this lab: these devices, as you will<br />

see if you look closely at one, are simply a very thin etched metal foil<br />

that essentially forms a very long, very thin wire that doubles back<br />

on itself many <strong>time</strong>s – see Figure 3.<br />

A force transducer may be constructed using <strong>strain</strong> gauges by<br />

bonding the gauges to the surface <strong>of</strong> a linearly elastic (Hookean)<br />

beam. As the beam is deformed, the <strong>strain</strong> gauge is stretched very<br />

slightly. You may recall from physics that the resistance <strong>of</strong> a wire is<br />

given by<br />

R = ρ l (5)<br />

A<br />

where R is the resistance <strong>of</strong> the wire in ohms (Ω), ρ is the electrical<br />

resistivity <strong>of</strong> the wire material in ohm-meters (Ω m) 2 , l is the length<br />

<strong>of</strong> the wire, <strong>and</strong> A is the cross sectional area <strong>of</strong> the wire. So, as the<br />

<strong>strain</strong> gauge is stretched, its l increases, thereby increasing its resistance.<br />

This change in resistance can be measured using a Wheatstone<br />

bridge, which is a st<strong>and</strong>ard device used to measure unknown resistance<br />

which may be familiar from physics.<br />

By hanging known masses from the Hookean beam to which<br />

the gauge is bonded, we may calibrate the force transducer – by<br />

measuring the voltage output <strong>of</strong> the transducer when various known<br />

forces are applied to it, we can obtain a calibration curve which will<br />

then allow us to measure unknown forces by reading the voltage<br />

output from the transducer.<br />

Data acquisition<br />

Once we have a tr<strong>and</strong>ucer capable <strong>of</strong> converting a force signal to an<br />

electrical signal, we next need to take that electrical signal <strong>and</strong> record<br />

it for later analysis. In this lab course, we will use small, low-cost<br />

microcontroller boards called Arduinos to measure data from our<br />

sensors <strong>and</strong> relay this information to our laptops.<br />

Arduinos are extraordinarily versatile devices with myriad applications,<br />

as you will see as this lab course progresses. In this particular<br />

lab, we have configured the Arduino boards to read in a single<br />

analog voltage on one <strong>of</strong> its analog input channels, convert this analog<br />

signal to a digital representation <strong>of</strong> that signal 3 , <strong>and</strong> then relay a<br />

message containing this digital representation to the laptop via USB<br />

serial. We can then configure the laptop to listen for these messages<br />

<strong>and</strong> store these as data.<br />

Figure 3: A typical foil <strong>strain</strong> gauge<br />

(Source: Wikimedia Commons)<br />

2<br />

Note that ρ here is not the density <strong>of</strong><br />

the material, it’s the resistivity. They<br />

somewhat annoyingly share the same<br />

symbol.<br />

Figure 4: Schematic for a wheatstone<br />

bridge. R 1 , R 2 <strong>and</strong> R 3 are all known<br />

resistors, but R X is unknown. In a<br />

single-gauge <strong>strain</strong> gauge application,<br />

R 1−3 are precision resistors, <strong>and</strong> R X is a<br />

<strong>strain</strong> gauge.<br />

3<br />

This is roughly analogous to recording<br />

sound on a computer through a<br />

microphone <strong>and</strong> sound card.


lab 3: stress, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biomaterials 4<br />

Methods<br />

There are two parts to this lab: stress-<strong>strain</strong> relationships <strong>and</strong> dynamic<br />

vibration tests <strong>of</strong> whole structures. In the former you will use<br />

force measurement <strong>and</strong> data acquisition to plot the stress <strong>strain</strong> relationship<br />

for some biomaterials <strong>and</strong> use that relationship to examine<br />

the <strong>stiffness</strong> <strong>and</strong> the extent to which that material has an anisotropic<br />

behavior. In the latter part, you will use free vibration tests to extract<br />

the <strong>stiffness</strong> <strong>of</strong> the material that constitutes a plant stem.<br />

It’s a stretch<br />

For this part <strong>of</strong> the lab you will use cyanoacrylate glue (superglue)<br />

to mount a piece <strong>of</strong> animal tissue on a force transducer. The animal<br />

tissue we will be using is pig intestine, which is a wonderfully thin<br />

<strong>and</strong> stretchy material typically used by humans as a sausage casing.<br />

We will measure the stress-<strong>strain</strong> relationship <strong>of</strong> pig intestine by<br />

cutting a small strip <strong>of</strong> the material <strong>and</strong> stretching it by applying a<br />

force to one end. This force will create stress within the material <strong>and</strong><br />

cause it to elongate (<strong>strain</strong>). As you modulate the amount <strong>of</strong> stress<br />

applied to the material, you will measure its deformation <strong>and</strong> thereby<br />

derive its stress-<strong>strain</strong> curve.<br />

Preparing your sample The pig intestine is obtained from a local<br />

butcher <strong>and</strong> arrives sleeved around a strip <strong>of</strong> plastic. The material is<br />

difficult to work with once it is removed from the plastic core sheet,<br />

so you will want to do as much preparation <strong>of</strong> the sample as possible<br />

before cutting away the plastic.<br />

First, determine whether you want to test the material along the<br />

axial or circumferential direction – you will test the gut in both directions<br />

for this lab, so the order doesn’t matter. Next, find a clean<br />

section <strong>of</strong> the intestine <strong>and</strong> glue two pieces <strong>of</strong> stiff plastic to the material<br />

with a ∼ 2 cm wide section <strong>of</strong> material between them. Try not<br />

to glue your fingers to the intestine, or to each other. Once you have<br />

glued the plastic tabs to the gut, cut away material so that you are left<br />

with a thin strip <strong>of</strong> material glued at either end to the tabs.<br />

Measuring stress <strong>and</strong> <strong>strain</strong> In this section, you will obtain a stress<strong>strain</strong><br />

relationship for your sample by stretching your sample with<br />

known weights. Begin by ensuring that you have an Arduino board<br />

connected via USB cable to your laptop, <strong>and</strong> start the yodaq application.<br />

When you hit the start button on the GUI, you will see a live<br />

display <strong>of</strong> the voltage output from the force transducer. Hang your<br />

sample from the force transducer using a thin loop <strong>of</strong> cotton thread,


lab 3: stress, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biomaterials 5<br />

<strong>and</strong> adjust the <strong>of</strong>fset knob on the transducer control box so that the<br />

output voltage from the transducer is 0 V.<br />

Next, you will hang weights from the free end <strong>of</strong> the sample to<br />

apply a stress. Our weights in this case are 1/4-20 nuts, which each<br />

weigh 12.6 g, <strong>and</strong> you will suspend these from your sample using a<br />

hooked paperclip weighing 0.2 g.<br />

As you add more mass to the free end <strong>of</strong> the sample, measure<br />

the extension <strong>of</strong> the sample with calipers or a ruler. Plot your results<br />

on the worksheet for both an axially-loaded sample <strong>and</strong> a<br />

circumferentially-loaded sample.<br />

Free vibration<br />

For this part <strong>of</strong> the lab we will be measuring the natural frequency<br />

<strong>of</strong> vibration for a biomaterial (a stick), <strong>and</strong> use these data to estimate<br />

the <strong>stiffness</strong> <strong>of</strong> bamboo.<br />

Recall from Eqn 4 that the natural frequency <strong>of</strong> a vibrating cylindrical<br />

beam is determined partly by its <strong>stiffness</strong>, as well as its density<br />

<strong>and</strong> shape. Therefore, if we can measure the natural frequency <strong>of</strong><br />

vibration for a structure where we know the dimensions <strong>of</strong> the structure<br />

<strong>and</strong> the density <strong>of</strong> the material from which it is composed, we<br />

can estimate the <strong>stiffness</strong> <strong>of</strong> the material.<br />

We will measure the natural vibration frequency <strong>of</strong> these beams<br />

by applying a foil <strong>strain</strong> gauge (see Section ) directly to the surface <strong>of</strong><br />

the beam, loading that beam <strong>and</strong> then releasing that load quickly (the<br />

kaboing-wonga-wonga test). As the beam vibrates, the <strong>strain</strong> gauge<br />

will be cyclically stretched <strong>and</strong> compressed, <strong>and</strong> we will be able to<br />

see this variation in <strong>strain</strong> as a function <strong>of</strong> <strong>time</strong>.<br />

Begin by ensuring that you have an Arduino board connected via<br />

USB to your laptop, <strong>and</strong> that the Arduino is supplying power to the<br />

Industrologic SGAU <strong>strain</strong> gauge amplifier. The <strong>strain</strong> gauge will<br />

form the fourth arm <strong>of</strong> a Wheatstone bridge, where the other three<br />

arms are precision 350 Ω resistors.<br />

Start yodaq again, <strong>and</strong> you should see an output voltage from the<br />

<strong>strain</strong> gauge amplifier that varies as <strong>strain</strong> is applied to the gauge.<br />

Adjust the VRO trimpot on the <strong>strain</strong> gauge amplifier to bring the<br />

output voltage up to the mid-point <strong>of</strong> the Arduino’s range.<br />

Apply a force to the end <strong>of</strong> the beam, <strong>and</strong> release this force to<br />

allow the beam to vibrate freely. You should see the output voltage<br />

from the <strong>strain</strong> gauge oscillate accordingly. Hit the “Stop” button<br />

after one such oscillatory period, <strong>and</strong> measure the <strong>time</strong> between<br />

successive peaks <strong>of</strong> the oscillation by clicking on the plot area. This<br />

is the period T <strong>of</strong> the natural vibration, <strong>and</strong> from this value you can<br />

Strain (%)<br />

Strain (%)<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

1<br />

2<br />

3<br />

4<br />

0 2 4 6 8 10<br />

4<br />

3<br />

2<br />

1<br />

0<br />

1<br />

2<br />

3<br />

0 2 4 6 8 10<br />

Time (s)<br />

Figure 5: Kaboing-wonga-wonga tests<br />

for two otherwise identical cylindrical<br />

beams <strong>of</strong> differing length. A single<br />

period <strong>of</strong> oscillation is highlighted for<br />

each beam.


lab 3: stress, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong> <strong>of</strong> biomaterials 6<br />

determine the natural frequency:<br />

ω n = 2π T<br />

Use this measured value along with Eqn 4 to estimate the Young’s<br />

modulus <strong>of</strong> the beam.<br />

Next, try <strong>and</strong> clamp the beam at a different point, so that its effective<br />

length is increased. Does the natural frequency change? Do<br />

you still get the same value for the Young’s modulus? If not, can you<br />

speculate as to why?


<strong>Lab</strong> 3: <strong>Stress</strong>, <strong>strain</strong>, <strong>stiffness</strong>, <strong>and</strong> <strong>time</strong>-<strong>dependent</strong> <strong>properties</strong><br />

Note: H<strong>and</strong> in one worksheet per lab group<br />

<strong>Lab</strong> Section:<br />

Name 1: Name 2:<br />

<strong>Stress</strong> (σ)<br />

Strain (ɛ)<br />

In the graph above, plot the stress <strong>strain</strong> curve for your samples. Use the same graph for the stress <strong>strain</strong><br />

curves for both axial <strong>and</strong> circumferential loading. Remember to label your axis ticks!<br />

<strong>Stress</strong>-<strong>strain</strong> data<br />

1 Sample length (m)<br />

2 Sample thickness (m)<br />

3 Sample width (m)<br />

4 Cross-sectional area (m 2 )<br />

Longitudinal<br />

Circumferential<br />

Free vibration data<br />

1 Sample length (m)<br />

2 Sample radius (m)<br />

3 Natural frequency (Hz)<br />

4 Young’s modulus<br />

Length 1 Length 2<br />

Questions<br />

1. What evidence do you have that intestinal lining is anisotropic? What do you think the function <strong>of</strong><br />

anisotropy is in the function <strong>of</strong> an intestine?<br />

2. If the beam you used in the free vibration test was anisotropic, how would its natural frequency vary<br />

with the direction <strong>of</strong> loading?

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!