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

Similar documents
INFLUENCE KINDS OF MATERIALS ON THE POISSON S RATIO OF WOVEN FABRICS

Size Effects In the Crushing of Honeycomb Structures

GG303 Lab 12 11/7/18 1

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4

BioMechanics and BioMaterials Lab (BME 541) Experiment #5 Mechanical Prosperities of Biomaterials Tensile Test

Strain Gauge Application and Measurement of Unknown Load

Photoelastic Experiments with a Transmitted-Light Polariscope

SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS

Image Acquisition and Sampling Theory

CHAPTER 6 FRICTION AND WEAR ANALYSIS FOR BUSHING

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

Numerical and Experimental Studies on Thermoforming Process. Sogang University

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

MSC Elastomers Seminar Some Things About Elastomers

December 1999 FINAL TECHNICAL REPORT 1 Mar Mar 98

Lectures on. Constitutive Modelling of Arteries. Ray Ogden

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

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

Developing Evaluation Model of Tire Pattern Impact Noise

Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation

Lecture #8: Ductile Fracture (Theory & Experiments)

HOMOGENEOUS ELECTRORHEOLOGICAL FLUIDS APPLIED TO VIBRATION CONTROL

Evaluation of in-plane orthotropic elastic constants of paper and paperboard

Mechanical Shock Testing for LIGA Materials Characterization

Experiment Two (2) Torsional testing of Circular Shafts

Lab Exercise #5: Tension and Bending with Strain Gages

STRUCTURAL OPTIMIZATION OF A MATERIAL EXHIBITING NEGATIVE STIFFNESS

Alternative speckle photography techniques for plastic deformation investigation

Special edition paper

VMS-GeoMil. Background

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

ME 243. Mechanics of Solids

Exercise: concepts from chapter 8

On Non-Ideal Simple Portal Frame Structural Model: Experimental Results Under A Non-Ideal Excitation

Schur decomposition in the scaled boundary finite element method in elastostatics

DEVELOPMENT OF DROP WEIGHT IMPACT TEST MACHINE

1 Nonlinear deformation

TE 75R RESEARCH RUBBER FRICTION TEST MACHINE

**********************************************************************

18. FAST NONLINEAR ANALYSIS. The Dynamic Analysis of a Structure with a Small Number of Nonlinear Elements is Almost as Fast as a Linear Analysis

1.105 Solid Mechanics Laboratory Fall 2003

Material Characterization of Carbon Fiber Reinforced Polymer Laminate Using Virtual Fields Method

Chapter 2. Rubber Elasticity:

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS

7. FORCE ANALYSIS. Fundamentals F C

Chapter 3. Experimentation and Data Acquisition

Prof. S.K. Saha. Sensors 1. Lecture 5 June 11, Prof. S.K. Saha. Purpose Classification Internal Sensors. External Sensors.

3-dimensional joint torque calculation of compression sportswear using 3D-CG human model

Study of Rupture Directivity in a Foam Rubber Physical Model

Fatigue-Ratcheting Study of Pressurized Piping System under Seismic Load

A rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates

A STUDY ON MULTI-AXIS FORCE MEASUREMENT OF POLYMER SKINS USING FBG SENSOR

Difference Between Fixed and Floating Reference Points AASHTO T-321

Design issues of thermal induced effects and temperature dependent material properties in Abaqus

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

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

Class XI Physics. Ch. 9: Mechanical Properties of solids. NCERT Solutions

STRAIN GAUGES YEDITEPE UNIVERSITY DEPARTMENT OF MECHANICAL ENGINEERING

Plane Strain Test for Metal Sheet Characterization

DEVELOPMENT, CONSTRUCTION, DELIVERY AND WARRANTY OF A DEVICE ABLE TO DETERMINE BOUNCE PARAMETERS OF A TABLE TENNIS BALL ON A TABLE SURFACE

The waves of damage in elastic plastic lattices with waiting links: Design and simulation

1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine Spring 2004 LABORATORY ASSIGNMENT NUMBER 6

Static and Kinetic Friction

Figure Two. Then the two vector equations of equilibrium are equivalent to three scalar equations:

FAVIMAT (AI)ROBOT2 FAVIGRAPH. Automatic Single-Fibre Testers

The strain response of silicone dielectric elastomer actuators

The Strain Gauge. James K Beard, Ph.D. Rowan Hall Auditorium November 2, 2006

A FAILURE CRITERION FOR POLYMERS AND SOFT BIOLOGICAL MATERIALS

Human Arm. 1 Purpose. 2 Theory. 2.1 Equation of Motion for a Rotating Rigid Body

Model tests and FE-modelling of dynamic soil-structure interaction

Simulation of the mechanical behavior of a rubber pumping system under large deformations.

Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures

Basic Considerations in Nonlinear Analysis

22 Which of the following correctly defines the terms stress, strain and Young modulus? stress strain Young modulus

AN LS-DYNA USER DEFINED MATERIAL MODEL FOR LOOSELY WOVEN FABRIC WITH NON-ORTHOGONAL VARYING WEFT AND WARP ANGLE

EFFECTS OF VALLEY SHAPE ON SEISMIC RESPONSES OF FILL DAMS

Because the third wire carries practically no current (due to the voltmeter's extremely high internal resistance), its resistance will not drop any

Effect of Liquid Viscosity on Sloshing in A Rectangular Tank

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.

(Refer Slide Time: 1: 19)

Exercise: concepts from chapter 5

Electronic Supplementary Information

walls, it was attempted to reduce the friction, while the friction angle mobilized at the interface in the vertical direction was about degrees under

Development of Laser Thickness Gauge in Steel Plate Shearing Line

University of Bristol - Explore Bristol Research. Publisher's PDF, also known as Version of record

Address for Correspondence

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions

ENSC387: Introduction to Electromechanical Sensors and Actuators LAB 3: USING STRAIN GAUGES TO FIND POISSON S RATIO AND YOUNG S MODULUS

Question Figure shows the strain-stress curve for a given material. What are (a) Young s modulus and (b) approximate yield strength for this material?

6 Examples of pneumatic structures and material models for membranes

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

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.

Hooke s Law. Equipment. Introduction and Theory


Influence of Back Pressure on Slab Edge Deformation Behavior during Width Reduction Pressing *1

Pullout Tests of Geogrids Embedded in Non-cohesive Soil

Work and Energy. W F s)

STRAIN ANALYSIS OF HDPE BUTT FUSION JOINT IN SIDE BEND TEST

Class XI Chapter 9 Mechanical Properties of Solids Physics

PREDICTION OF OPEN HOLE COMPRESSIVE FAILURE FOR QUASI-ISOTROPIC CFRP LAMINATES BY MMF/ATM METHOD

Transcription:

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 their history curves obtained from biaxial tension tests. hree different kinds of tests, uniaxial tension test, equi-biaxial tension test, and strip-biaxial tension test are carried out. We correct the measured displacement values by using strain fields measured from the specimen s photographic images in order to obtain the actual strain data. In this report, we firstly outline the equipment for biaxial tension tests. Secondly, we present the shapes or other characteristics of the specimens used in uniaixial tension tests, equi-biaxial tension tests and strip-biaxial tension tests. Finally, actual test procedures are explained. Experimental equipment. In this working group in JANCAE, we conduct material tests to characterize the history curves of rubber materials. A testing device for biaxial tension tests is used to elongate a specimen for three cycles. he deformation rate is controlled by an outside device. In this working group, trains calculated from measured displacements are corrected by an image processing. Images are taken by a digital CCD camera placed on the top of the biaxial tension device, and at the same time the load and displacement are measured. Fig. 1 shows the schematic of the biaxial tension device. Since very low current flew out of a PC used for controlling the devices, electric noise occurred. herefore, we blocked the noise by using an outlet with three lines, one of which is connected to the ground. Below, we show the outline of individual experimental equipments. CCD Camera ripod sustaining a camera Specimen Button for internal control Load indicator displacement indicator Fig. 1 ension test device.1 Biaxial tension test he biaxial tension test device used for this experiment moves its two arms along the center lines of a specimen to apply strains on it. he performances of the device are presented in table 1. he biaxial tension

device and the specimen are shown below. Figure is an entire view of the biaxial tensile test device. Here, the X- and Y-directions are defined as illustrated in the figure. he specimen is set so as to be pinched by chucks. Different chucks are used for specimens in tension and biaxial tests, which are shown in Figs.3 and, respectively. Y X Fig. Entire view of the biaxial tension test device Y X Fig. 3 Chuck of the uniaxial tension test device Fig. Chuck for the biaxial tension test device X-Y axis driving Speed setting Elongation speed Displacement display Distance display Device elongation range Basic elongation range Detection by load controlling load cell Speed instruction by external control Power source able 1 Performance of the biaxial tension device Usage (Control and performance) Motor: HC-KFS3G 1/0, units Independent (Internal control) vernier dial and digital display 100 100 mm/min (1.6 to 0 mm/sec) Distance detection by potentionmeter. Format: LP-00F MIDORI PRECISIONS CO., LD. Digital display: from 0 (Size within chuck: 50) to 999mm X-axis: from 50 to 00 mm, Y-axis: from 50 to 80 mm X-axis: from 60 to 360 mm, Y-axis: from 60 to 0 mm DC 0 10v with analog output 1.0 V (X-axis: scale 0 mm) 1.0 V (Y-axis scale 35.8 mm) Load cell U-BR 500N (51.0 kg), two units each Amplifier D-30A, EAC Corporation Output by an open collector with upper limit: DC 0 to 10V with analog output 10.0 V at 999.9 N, implying 100 N/1.0 V Analog input: DC0 10.0 V (0mm/sec) AC00-0V 15-0A (Current requirement less than 6A) During operation as interlocked, operation stops when top/bottom/corners are outputted for a load cell protection. Internal/external overrun sensor is arranged for the width drawn due to incorrect instruction.

. Image capturing equipment Figure 5 shows the digital CCD camera used in the image capturing to obtain displacement data. he performance of this camera is presented in able. Since the spatial resolution is set at 08 000 [pixels] in this test, the temporal resolution is expressed in the formula given below. able 1 indicates a temporal resolution in the case of spatial resolution of 08 08 [pixels]. Assume that it takes t [sec/frame] to record one image frame, 3 t (11. 0.0 y) 10 (1) where y [pixel] is the spatial resolution in the vertical direction. herefore, assuming the temporal resolution is f [frame/sec], we have 1 f () t which means the temporal resolution is f 16.0. able Performance of digital CCD camera Name Adimec-000m Spatial resolution [pixel] 08 08 emporal resolution [frame/sec] 16. Recording method Image file format Progressive scan 8bit grayscale 3 Outline of specimens Fig. 5 Digital CCD Camera 3.1 Shape and dimensions of specimen for uniaxial tension test Figure 6 present the shape and dimensions of a specimen, and Figure 7 shows the photo of an actual specimen for uniaxial tension tests. A uniaxial tension test in this context indicates a loading scheme in which a specimen is elongated only in one direction and the deformation in the direction perpendicular to the tensile direction is not constrained. Figure 8 shows an appearance of actual testing. We made a specimen for a uniaxial tension test by dividing a specimen for a biaxial tension test into two and by cutting out the chuck portions arranged in the direction perpendicular to the tensile direction. he maximum stroke of the chuck for a uiaxial tension test is 360 [mm] and is attached along the X-axis. 10 60 Cross section1 3 10 (a) Plain view R= R= 6 1 3 10 56 10 3 (b) Cross section view along the cross section 1 Fig. 6 Shape and dimensions of a specimen for a uniaxial tension test

Fig. 7 Photo of a specimen for a uniaxial tension test Fig. 8 Appearance of a uniaxial tension test 3. Shape of specimen for biaxial tensile test Figure 9 shows the shape and dimensions of a specimen, and Figure 10 is a photo of a actual specimen. he biaxial tension test device can impose arbitrary deformations on a specimen in the two direction independently (referred to as X-axis and Y-axis). his specimen is used for biaxial tension tests. In a strip-biaxial tension test, a specimen is deformed in the X-direction, and the stretch in the Y-direction is fixed. In a equi-biaxial test, a specimen is elongated in both X- and Y-directions at the same loading speed and amplitude. Figures 11 and 1 show the appearances of these tests. 10 56 10 10 Cross section 1 56 10 R= R= 3 (a) Plain view R= R= 1 3 10 56 10 3 (b) Cross section view along the cross section 1 Fig. 9 Shape of specimen for biaxial tests

Fig. 10 Photo of a specimen for biaxial tension tests Fig. 10 Appearance of a strip-biaxial tension test Fig. 11 Appearance of a equi-biaxial tension test 3.3 ypes of rubber Each rubber has its own hardness and damping characteristics. he damping values of a rubber with hardness 50 are tanδ=0.03 and tanδ=0.05, while those with hardness 65 are tanδ=0.0 and tanδ=0.7. he names of these materials are defined as in able 3. Also, tanδ is called the loss coefficient, which indicates the damping performance (internal friction) of a rubber. able 3 Names of rubbers Namel Material A Material B Material C Material D Hardness 50 65 50 65 Damping tanδ=0.03 tanδ=0.05 tanδ=0.0 tanδ=0.7 esting method.1 esting procedure As described in the beginning of this report, the picture images of the deformation process of a specimen are captured by a CCD camera placed on the top of the test device and the strain values measured are corrected with the obtained image information (details are omitted here). ime intervals for capturing picture images are changed according to tension speed and amplitude. he number of recorded images is 100 pictures per half amplitude motion for a specimen, so about 00 pictures are captured in a single test. able indicates the number of specimens used in each test. he uniaxial tension tests were conducted three times at each of the loading speed of m/sec and 0 mm/sec. he strip-biaxial tension tests were conducted once at each of the loading speed of m/sec and 0mm/sec. he equi-biaxial tests were conducted once at each of the loading speed of 0.6 mm/sec, m/sec and 0mm/sec. For image capturing, an additional test was conducted for each material, each speed and each loading speed.

. biaxial tension tests able Number of specimens (Unit: piece) Uniaxial Strip-biaxia Equi-biaxial 8 16 In this working group, we conducted three different kinds of tension tests as follows: (1) Uniaxial tension test () Strip-biaxial tension test (3) Equi-biaxial tension test Using the results of these three tests, we will characterize the energy of rubber for each test. Uniaxial tension tests were conducted with a specimen shown in Fig. 7. At the same time, in a test accompanying image capturing to correct measured strain values, we used a specimen with white dots being drawn as shown in Fig. 13 and captured its deformation by a CCD camera to measure their displacements. he strip-biaxial and equi-biaxial tension tests were conducted with a specimen shown in Fig. 10. Like in the uniaxial tests, in a test accompanying image capturing, we used a specimen with white dots being drawn as shown in Fig. 1 and captured its deformation by a CCD camera to measure their displacements. hose tests were conducted at the elongation speed of [mm/sec] and 0 [mm/sec] to study the effects of energy absorption in response to the change of the loading speed. Since the test results showed a similar trend irrespective of loading methods, we additionally conducted only for a equi-biaxial test at the loading speed of 0.6 [mm/sec]. he initial value and increment of the loading amplitude for one of the three uniaxial tension tests are 1.5 and 0.5, respectively, while those of the other two are 1.5 and 0.5, respectively. he tests were continued until a specimen is broken. On the other hand, the initial value and increment of the loading amplitude the strip-biaxial and equi-biaxial tension tests are 1.5 and 0.5, respectively, and the tests were continued until a specimen is broken. Also, an image capturing at every test was conducted for the amplitude smaller than the stretch value at fracture by one. Sampling frequencies were Hz, 5Hz and 50Hz for 0.6 [mm/sec], [mm/sec] and 0 [mm/sec], respectively. emperatures were maintained within the range of about 18 to C in all the tests. able 5 indicates a list of conducted tests and able 6 indicates a list of sampling frequencies at each loading speed. Fig. 13 Uniaxial test s specimen for image capturing Fig. 1 Biaxial tests specimen for image capturing

able.5 Loading conditions of the tests Name of Loading emp. Humidity Loading Specimen method (C) (%) Speed(mm/sec) Amplitude (mm) Equi-biaxial 0 37 90, 105, 10, 135, 150, 165 Strip-biaxial 0 37 90, 105, 10, 135, 150, 165 Uniaxial 0 3 0mm/sec 90, 10, 150, 180, 10, 0 Equi-biaxial 0 3 90, 105, 10, 135, 150, 165, 180 Strip-biaxial 0 1 90, 105, 10, 135, 150, 165, 180, 195 0 36 90, 105, 10, 135, 150, 165, 180, 195, 10, 5, 0 0 36 90, 10, 150, 180, 10, 0 Uniaxial 0 36 mm/sec 90, 10, 150, 180, 10, 0 A Equi-biaxial 0 36 0.6mm/sec 90, 105, 10, 135, 150, 165, 180 Equi-biaxial 19 1 90, 105, 10, 135, 150 Strip-biaxial 0 6 90, 105, 10, 135, 150, 165, 180, 195 Uniaxial 0 3 0mm/sec 90, 10, 150, 180 Equi-biaxial 17 90, 105, 10, 135, 150, 165 Strip-biaxial 18 90, 105, 10, 135, 150, 165, 180, 195, 10 0 33 90, 105, 10, 135, 150, 165, 180, 195, 10, 5 0 33 90, 10, 150, 180, 10 Uniaxial 0 33 mm/sec 90, 10, 150, 180 B Equi-biaxial 0 36 0.6mm/sec 90, 105, 10, 135, 150, 165, 180 Equi-biaxial 18 90, 105, 10, 135, 150, 165 Strip-biaxial 18 90, 105, 10, 135, 150, 165, 180, 195, 10, 5, 0, 55, 70, 85, 300 33 90, 105, 10, 135, 150, 165, 180, 195, 10, 5, 0, 55, 70, 85, 300 1 35 90, 10, 150, 180, 10, 0, 70, 300 Uniaxial 0 35 0mm/sec 90, 10, 150, 180, 10, 0, 70, 300 Equi-biaxial 0 1 90, 105, 10, 135, 150, 165, 180, 195 Strip-biaxial 0 35 90, 105, 10, 135, 150, 165, 180, 195, 10 19 50 90, 105, 10, 135, 150, 165, 180, 195, 10, 5, 0, 55, 70, 85, 300 19 50 90, 10, 150, 180, 10, 0, 70, 300 Uniaxial 19 50 mm/sec 90, 10, 150, 180, 10, 0, 70, 300 C Equi-biaxial 0 0.6mm/sec 90, 105, 10, 135, 150 Equi-biaxial 19 55 90, 105, 10, 135, 150, 165, 180, 195 Strip-biaxial 0 50 90, 105, 10, 135, 150, 165, 180, 195, 10, 5 19 51 90, 10, 150, 180, 10, 0, 70 19 51 90, 10, 150, 180, 10, 0, 70, 300 Uniaxial 19 51 0mm/sec 90, 10, 150, 180, 10, 0, 70, 300 Equi-biaxial 0 51 90, 105, 10, 135, 150, 165, 180 Strip-biaxial 0 51 90, 105, 10, 135, 150, 165, 180, 195, 10, 5 19 5 90, 105, 10, 135, 150, 165, 180, 195, 10, 5, 0, 55, 70 19 5 90, 10, 150, 180, 10, 0, 70, 300 Uniaxial 19 5 mm/sec 90, 10, 150, 180, 10, 0, 70, 300 D Equi-biaxial 1 6 0.6mm/sec 90, 105, 10, 135, 150, 165, 180 able 6 Sampling frequency at each loading speed (Unit: Hz) Loading speed 0.6 mm/sec mm/sec 0 mm/sec Sampling frequency 5 50 5 reating method of experimental results 5.1 How to calculate strain field By using the coordinate data at each point in a specimen obtained by an image processing (details are omitted here), we follow the displacements at the center of a white dot marked on the specimen, and calculate strains by comparing the images before and after the deformation. he center point of a square consisting of four coordinate data points in a specimen is defined as one node of a four-node element. Strains in the element were calculated by interpolation with appropriate interpolation functions.

(1) Calculation of strain field (a) Interpolation of displacement A displacement value known at kth point (hereinafter referred to as node k) by the image processing is denoted by u i, and the position of a dot in the undeformed state is denoted by X i. hen, the displacement is interpolated according to the following formula: n u1 N u1 { N} { u1} (3a) k 1 n k 1 u N u { N} { u } (3b) where i=1 indicates the X-axis and i= the Y-axis. Here, u i is a value of the displacement u i at node k, and N is an interpolation function associated with u 1 and ( k u ). Each vector in formulae (3a) and (3b) are defined as follows; (1) () ( ) { 1} n u { u1, u1,, u1 } (a) (1) () ( ) { } n u { u, u,, u } (b) (1) () ( ) { } n N { N, N,, N } (c) where n indicates the number of nodal points of elements being used. When a square consisting four points are used as one element of an interpolation function in the above formulae, we have n= and a specific interpolation function associated with each nodal point is expressed as follows; (1) 1 N (1 )(1 ) (5a) () 1 N (1 )(1 ) (5b) (3) 1 N (1 )(1 ) (5c) () 1 N (1 )(1 ) (5d) In the above formula, and are within 1, 1, which are called the natural coordinates (Fig. 15) by finite element method. In this working group, the center of a square consisting four points are used, we used 0. hese are regarded as a kind of parameters and have the following relations with spatial variables ( X1, X ) : X1 N X1 (6b) X N X (6b) where ( X1, X ) are a specific numerical values of ( X1, X ) at a node k. With these interpolation functions, we treated the experimental results. he target of the interpolation is a deformation gradient tensor. node (-1,1) node3 (1,1) (-1,-1) node1 (1,-1) node Fig. 15 Natural coordinate system (b) Calculation of the deformation gradient tensor Denoting the position of a material point in the reference configuration (undeformed state) by X i and the current position (after deformation) by x i, the deformation gradient tensor is defined as

F x i ij (7) X i where the current position x i is xi ui Xi (8) hus, when substituted to formula (10), this can be written as ( ui Xi) ui Fij ij (9) X j X j Substitution of this expression into (6) yields { N} { ui} { N} Fij ij ij { ui } (10) X j X j From this, if each component of (1) () ( n) { N} N N N, (11) X j X j X j X j can be calculated, the deformation gradient tensor component F ij can be obtained. For the details of how to calculate formula (1) is explained in (3) later. (c) Calculation of the right Cauchy-Green deformation tensor and Green_Lagrange strain According to the definition, the right Cauchy-Green deformation tensor is calculated by the formula C F F (1) and the Green_Lagrange strain is defined as 1 E ( C I ) (13) where, I is a unit matrix. When each component of the deformation gradient tensor F is calculated from the formula (13), the components of the above two tensors can be easily obtained. () Spatial derived function of interpolation function As indicated by the formula (13), to obtain the specific values of strain and etc., one have to evaluate derivatives of interpolation function N with respect to X 1 or X. However, the interpolation function is directly a function of parameter (, ) and directly associated with ( X1, X ) by the formula (9). Accordingly, when the chain rule of differentiation is applied in a general manner, we have ( k) N N N (1a) X1 X1 X1 ( k) N N N (1b) X X X When the above formulas are summarized and described in a matrix, we have N N X 1 X1 X 1 ( k) ( k) (15) N N X X X in which N N (16) can easily be calculated from formula (5). On the other hand, the right side matrix of formula (15) is calculated as explained below. Firstly, assuming that interpolation function N is a function of (X 1, X ), we apply the chain rule to obtain ( k) N N X1 N X (17a) X1 X ( k) N N X1 N X (17b) X1 X he corresponding matrix notation is given as

N X1 XN X1 ( k) ( k) (18) N X1 XN X When the above formula is compared with formula (15), we can identify the following relationship: 1 X1 X X1 X 1 (19) X1 X X X Finally, by substituting formula (6) into to the right side of the above equation, we obtain: X1 X N ( k) N ( k) X1 X (0) X1 X N ( k) N ( k) X1 X As a result, the matrix in the right side of formula (15) can be calculated from formula (19). hus, the derivative of interpolation function N with respect to X 1 or X can easily be calculated. (3) Examples of strain field (stretch) in actual tests Figures 16, 17 and 18 show the distributions of strain fields obtained from the image processing in a uniaxial tension test for material A at the loading speed of mm/sec, the one in a strip-biaxial tension test for material D at mm/sec, and the one in a equi-biaxial test for material D at 0mm/sec. Figure 16 for the uniaxial tension test and Fig. 17 for the strip-biaxial tension test show that when elongated in the horizontal direction, the material is contracted in the vertical direction as the loading increases. his is because the volume of the material remain constant even though it is deformed.

(a) Before deformation (b) At stretch 1.5 (c) At stretch (d) At stretch.5 (e) At stretch 3 (f) At stretch 3.5 Fig. 16 Strain field obtained from the image processing for uniaxial tension tests on material A at loading speed mm/sec

(a) Before deformation (b) At stretch 1.5 (c) At stretch (d) At stretch.5 (e) At stretch 3 (f) At stretch 3.5 Fig. 17 Strain field obtained from the image processing for strip-biaxial tension tests on material D at loading speed mm/sec

(a) Before deformation (b) At stretch 1. (c) At stretch 1.8 (d) At stretch. (e) At stretch.6 (f) At stretch 3 Fig. 18 Strain field obtained from the image processing for biaxial tension tests on material D at loading speed 0 mm/sec 6 Conclusion In this report, we firstly presented the outline of the biaxial loading device and discuss the correction of input data to the device. Secondly, we explained the outline of the specimen and the testing method employed in this working group. Finally, by interpolating the displacement values from those of the center point of white markers, we calculated a strain field.