Theoretical and Experimental Investigation of Magnetostrictive Tagged Composite Beams

Size: px
Start display at page:

Download "Theoretical and Experimental Investigation of Magnetostrictive Tagged Composite Beams"

Transcription

1 Theoretical and Experimental Investigation of Magnetostrictive Tagged Composite Beams by Florin G. Jichi Bachelor of Science Technical University of Timisoara, Romania, 1995 Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Science in the Department of Mechanical Engineering College of Engineering & Information Technology University of South Carolina 2000 Department of Mechanical Engineering Director of Thesis Department of Mechanical Engineering Second Reader Dean of the Graduate School

2 ACKNOWLEDGEMENTS As any research endeavor cannot be done alone, but more likely is a team effort, there are many members involved. Therefore, I would like to express my gratitude and deeply thanks to Dr. Victor Giurgiutiu for all his guidance and help throughout my entire time spent here at USC. Without his support, none of what I accomplished would not have been possible. In addition, I would like to sincerely thank Dr Abdel-Moez E. Bayoumi for his time from his busy schedule. My gratitude also goes to, Adrian Cuc, Paulette Goodman, Jingjing Jack Bao, Greg Nall, Radu Pomirleanu, and Andrei Zagrai, my colleagues, who made my time spent here so pleasant and fun. I would like to express my forever gratitude and thanks to my mother, father, sister and brother-in-law for their trust, incommensurable support, love and encouragement.

3 ABSTRACT Among novel non-destructive evaluation techniques for structural health monitoring, the magnetostrictive-tagged fiber-reinforced composites stand out as especially suitable due to: (a) distributed sensory properties; (b) non-contact damage detection; and (c) straight forward manufacturing implementation. Experimental data and mathematical modeling of a magnetostrictive-tagged fiber reinforced composite specimen under bending (flexural) loading are presented. A brief review of the state of the art identifies previous work on axially loaded magnetostrictive composites, but finds no previous work on bending. Description of bending specimen design and fabrication is followed by the theoretical analysis and by the description of the experimental set-up and equipment used. Several analysis models were used. Test data, with and without magnetic annealing between loading cycles, is presented and results are discussed. Numerical values for the stress and strain versus magnetic flux density coefficients are given for both annealed and non-annealed cases. Piezomagnetic coefficients for the magnetostrictive composite are calculated. The correlation between the results developed in the presented paper for bending and previously published results for axial loading is found to be within 10% after correction factors depending on the quantity of the magnetostrictive material are applied. In conclusion, the usefulness of this method for structural health monitoring and further work are discussed.

4 TABLE OF CONTENTS ACKNOWLEDGEMENTS...II ABSTRACT... III TABLE OF CONTENTS... IV LIST OF FIGURES... VI LIST OF TABLES...XII 1 INTRODUCTION PREVIOUS WORK ON MAGNETOSTRICTIVE TAGGED COMPOSITES FOR STRUCTURAL HEALTH MONITORING (TENSILE EXPERIMENTS) CHARACTERISTICS OF ETREMA S TERFENOL-D MAGNETOSTRICTIVE MATERIAL PRESENT INVESTIGATION ANALYSIS OF A MAGNETOSTRICTIVE COMPOSITE BEAM STATIC ANALYSIS OF SIMPLY SUPPORTED COMPOSITE BEAM UNDER CENTRAL LOAD ANALYSIS OF A MAGNETOSTRICTIVE LAMINATED COMPOSITE Micromechanics analysis of the MS tagged composite specimen Lamination Analysis PLY BALANCED-ORTHOTROPIC MODEL PLY CROSS-PLY ANGLE MODEL PLY CROSS-PLY ANGLE MODEL PLY CROSS-PLY ANGLE MODEL CONVERGENCE ANALYSIS Convergence of the strain...26

5 2.7.2 Convergence of stress MAGNETOSTRICTIVE COMPOSITE BEAM EXPERIMENT SPECIMEN DIMENSIONAL DESIGN DESCRIPTION OF THE MS TAGGED COMPOSITE SPECIMEN SPECIMEN PREPARATION FOR THE EXPERIMENT EXPERIMENTAL DESIGN Description of the experiment Clamping fixture design Protective wood fixture design for the strengthen the gaussmeter probe List of the equipment used Equipment Calibration LabView Virtual Instrument Environment EXPERIMENTAL PROCEDURE AND RESULTS MAGNETIC ANNEALING TESTING PROCEDURE FOR MS TAGGED COMPOSITE BENDING EXPERIMENT RESULTS WITHOUT MAGNETIC ANNEALING Displacement and strain without magnetic annealing Magnetic flux density without magnetic annealing RESULTS WITH MAGNETIC ANNEALING BETWEEN LOADING-UNLOADING CYCLES Displacement and strain with magnetic annealing between loading-unloading cycles Magnetic flux density with magnetic annealing between loading-unloading cycles MAGNETIC FLUX DENSITY ANALYSIS PIEZOMAGNETIC STRESS AND STRAIN COEFFICIENTS Piezomagnetic stress coefficient, e 31 without magnetic annealing...73

6 4.6.2 Piezomagnetic stress coefficient, e 31, with magnetic annealing Piezomagnetic stress coefficient,e 31,using experimental magnetic flux density without magnetic annealing and the 56-ply strain results model Piezomagnetic stress coefficient, e 31 using experimental magnetic flux density with magnetic annealing between loading-unloading cycles and the 56-ply strain results model Piezomagnetic strain coefficient d 31 using experimental magnetic flux density without magnetic annealing and 56-ply stress model Piezomagnetic strain coefficient,d 31,for annealed specimen experiment and 56-ply stress model ANALYSIS AND DISCUSSION COMPARISON OF PIEZOMAGNETIC COEFFICIENTS FOR MAGNETOSTRICTIVE TAGGED COMPOSITES WITH PREVIOUSLY PUBLISHED DATA COMPARISON BETWEEN EXPERIMENTAL, DESIGN AND MODEL STRAIN FOR WITH AND WITHOUT MAGNETIC ANNEALING COMPARISON BETWEEN EXPERIMENTAL AND MODEL PIEZOMAGNETIC STRESS COEFFICIENT WITHOUT MAGNETIC ANNEALING COMPARISON BETWEEN EXPERIMENTAL AND MODEL PIEZOMAGNETIC STRESS COEFFICIENT WITH MAGNETIC ANNEALING CONCLUSIONS BIBLIOGRAPHY APPENDIX...91

7 LIST OF FIGURES Figure 1 Magnetic flux in axial and thickness direction during axial loading of magnetostrictive tagged composite specimens (White, 1999)... 3 Figure 2 Transverse (x-axis) magnetic flux density versus load: (a) neat resin;... 4 Figure 3 Stress/magnetic flux density cyclic testing results: (a) without annealing between cycles; (b) with annealing between cycles (Quattrone, Berman, and White, 1998)... 5 Figure 4 Strain versus magnetic strength (Butler, 1988)... 6 Figure 5 Relationships between a) magnetic flux density versus magnetic strength, b) stress versus magnetic flux density, c) stress versus magnetic strength... 6 Figure 6 Schematic of the static system... 9 Figure 7 Schematic of a generic lay-up Figure 8 Balanced-orthotropic lay-up model used in the analysis Figure 9 Cross-ply model used in the analysis Figure 10 Strain distribution of the 7, 14, 28 model with 0 /90 degree ply angle in longitudinal direction Figure 11 Convergence of the strain result of the models to the experimental data Figure 12 Strain distribution for the 7, 14, 28 mathematical models with 0 /90 degree ply angle in transversal direction

8 Figure 13 Stress thickness distribution of the 7-ply balanced-orthotropic model with 0 degree fiber angle Figure 14 Stress thickness distribution in longitudinal direction for the 14-ply 0 /90 degree fiber angle model Figure 15 Thickness distributions of longitudinal stress for 28-ply 0 /90 degree fiber angle model Figure 16 Thickness distributions of longitudinal stress for 56-ply 0 /90 degree fiber angle model Figure 17 Stress-thickness distribution for the 7, 14, 28-ply models with 0 /90 degree fiber angle Figure 18 Convergence of results on the specimen surface under maximum load condition (F max = 58.9 N, M max = Nm/m) as predicted by various models and measured by experiment for stress results Figure 19 Stress thickness distribution in longitudinal direction for 14-ply 90 /0 degree fiber angle model Figure 20 Stress - thickness distribution in longitudinal direction for the 28-ply 90 /0 degree fiber angle model Figure 21 Stress thickness distribution in longitudinal direction for 56-ply 90 /0 degree fiber angle model Figure 22 Stress-Thickness Distribution for the 7 ply 0 degree angle model and 14, 28- ply models with 90 /0 degree fiber angle... 38

9 Figure 23 General view of the experiment set-up Figure 24 Schematic of the experiment and data flow Figure 25 Position of the gaussmeter probe, strain gage, and displacement transducer47 Figure 26 Concept design of the protective wood fixture Figure 27 Overview of the assembly of magnetic sensor and carved wood Figure 28 Front panel of the program created for the MS composite beam experiment using the LabView virtual instrument Figure 29 The logic flow of the program created for the MS composite beam experiment using the LabView virtual instrument Figure 30 Mechanical data resulting from the experiments without magnetic annealing between loading cycles: displacement vs. load Figure 31 Mechanical data resulting from the experiments without magnetic annealing between loading cycles: strain vs. load Figure 32 Mechanical data resulting from the experiments without magnetic annealing: strain vs. displacement Figure 33 Magnetic flux density responses, without magnetic annealing, to bending strain in top surface of MS-tagged composite specimens: superposed 10 consecutive cycles Figure 34 Mean value and standard deviation of the magnetic field Figure 35 The resulted magnetic field as function of measured strain Figure 36 Mean values of magnetic flux density and standard deviation versus strain

10 for the specimen without annealing Figure 37 Displacement results versus applied force for annealed specimen Figure 38 Strain results for the MS composite specimen with magnetic annealing between cycles Figure 39 Magnetic flux density results for specimen with magnetic annealing between cycles Figure 40 Mean value and standard deviation of magnetic flux density results vs. load for the specimen with magnetic annealing between cycles Figure 41 Experimental magnetic flux density results versus strain with magnetic annealing of the specimen between cycles Figure 42 Mean value and standard deviation of magnetic flux density results vs. strain for the specimen with magnetic annealing between cycles Figure 43 Magnetic flux density in relationship with model longitudinal stresses on the surface of the composite specimen model Figure 44 Standard deviation of the stress from the mean value versus magnetic flux density Figure 45 Magnetic flux density versus strain on the top surface of the model material considered and 0 /90 degree fiber angle without magnetic annealing Figure 46 Plot of the average magnetic flux density vs. strain showing the trend lines in determination of the piezomagnetic stress coefficient e 31 for non -annealed specimen case 74

11 Figure 47 Plot of the average magnetic flux density vs. strain showing the trend lines in determination of the piezomagnetic stress coefficient e 31 of magnetic annealing of the specimen Figure 48 Plot of the average magnetic flux density with no magnetic annealing specimen vs. 56-ply model strain showing the trend lines in determination of the piezomagnetic stress coefficient e Figure 49 Plot of the average magnetic flux density with magnetic annealing vs. 56-ply model strain showing the trend lines in determination of the piezomagnetic stress coefficient e Figure 50 Trend line for calculating the piezomagnetic strain coefficient Figure 51 Trend line for computing the piezomagnetic strain coefficient d 31 for the case of magnetic annealing of the specimen between cycles... 79

12 LIST OF TABLES Table 1 Dimensions of the MS tagged composite specimen... 9 Table 2 Mechanical properties considered for the magnetostrictive-tagged woven composite material specimen Table 3 Lay-up of the composite models for 14, 28 and 56-ply models Table 4 Comparison of the longitudinal strain on the specimen surface under maximum load condition (F max = 58.9 N, M max = Nm/m) as predicted by various models and measured by experiment Table 5 Stress results on the surface for 7 balanced-orthotropic, 14, 28-ply cross-ply models 34 Table 6 Stress results on the surface for 7, 14, and 28 ply 0 and 90 /0 cross-ply models 38 Table 7 List of the equipment used in the magnetostrictive composite beam experiment Table 8 The displacement data without magnetic annealing Table 9 Strain data without magnetic annealing Table 10 Magnetic flux density results without magnetic annealing Table 11 Displacement results with annealing between loading-unloading cycles Table 12 Strain results with annealing between loading-unloading cycles... 65

13 Table 13 Experimental magnetic flux density with magnetic annealing between loading-unloading cycles Table 14 Summary of magnetostrictive coefficient for MS tagged composites, as determined in the present work and by previous investigators Table 15 Comparison of the design, model and measured strain INTRODUCTION In recent years, numerous applications in civil engineering construction used composite materials. The increase of the usage of the composite materials in civil construction imposed problems of evaluating the in-service composite civil engineering structures. The evaluation of the composite is a wide area in which engineers and researchers have proposed several technologies, for in-service Non-destructive evaluation (NDE). Conventional NDE methods, initially developed for metallic structures, have been shown to be less effective in monitoring composite structures due to the micromechanical complexity of the composite material. New NDE technologies are required. The NDE technology analyzed in this paper is proposed for the inspection of advanced composite by using the magnetostrictive particle tagging technique. Magnetostrictive-tagged composites permit: (a) distributed sensory properties; (b) noncontact damage detection; and (c) straight forward manufacturing implementation. 1.1 PREVIOUS WORK ON MAGNETOSTRICTIVE TAGGED COMPOSITES FOR

14 STRUCTURAL HEALTH MONITORING (TENSILE EXPERIMENTS) Terfenol-D is a magnetic anisotropy-compensated alloy Tb x Dy 1-x Fe 2 that shows a strong magnetostrictive behavior. The name Terfenol-D represents the composition of the material and the original name of the Navy Laboratory at which the work begun. Ter represents Terbium, Fe from iron, nol Naval Ordnance Laboratory and D Dysprosium. A magnetostrictive (MS) material, such as Terfenol-D, produces magnetic field when subjected to mechanical strain. This phenomenon is known as the converse magnetostrictive effect. White and collaborators (White, Albers and Quattrone, 1996; White and Brouwers, 1998) did extensive work on magnetostrictivetagged composites under axial loading. White (1999) reviewed the magnetostrictive tagging methodology of composites for structural health monitoring and gave an update on recent results. An experiment in which MS-tagged composite specimens were subjected to uniaxial tension in a testing machine was presented. Neat resin specimens tagged 2.24% by volume with magnetostrictive Terfenol-D powder were used. The resulting magnetic field was measured in both the axial and thickness directions (Figure 1). Trovillion et al. (1999) studied the magnetic characteristics of neat resin and glassfiber-reinforced magnetostrictive composites subjected to axial load. The fiber reinforced polymer composite (FRP) specimens consisted of 4 layers of continuous strand glass mat fibers embedded in a polyester resin. Trovillion et al. (1999) studied the magnetic characteristics of neat resin and glass-fiber reinforced magnetostrictive composites subjected to axial load. The fiber reinforced polymer composite (FRP) specimens consisted of 4 layers of continuous strand glass mat fibers embedded in a polyester resin.

15 Figure 1 Magnetic flux in axial and thickness direction during axial loading of magnetostrictive tagged composite specimens (White, 1999) The top lamina of the composite was impregnated with Terfenol-D powder at a volume fraction of 2.24% for that lamina. The specimens were subjected to uniaxial loading under load control at a rate of 0.02 kn/s. Hall-effect device were used to measure the magnetic field response to subsequent loading and unloading. As seen in Figure 2, both measuring devices gave similar results. Magnetic annealing (i.e. rearranging the magnetic dipoles chains of the magnetostrictive molecules through application of strong magnetic field) was performed by applying a magnetic field through the thickness of the specimen using a pair of 800 Gauss permanent magnets. Nersessian N. and Carman G.P. presented at ASME conference in Orlando results of tests for five different volume fraction composites namely 10 %, 20%, 30 %, 40 % and 50 %. The composite were tested under constant magnetic field with varying the mechanical load and constant mechanical load with varying magnetic field conditions. Results for the constant magnetic filed test indicated that modulus generally increases with volume fraction and increasing H/H max. For low fields, initial dip is noticed in modulus attributed to domains becoming more mobile at lower magnetic field levels. Results presented for the constant load test show a strong dependence of strain output on applied pre-stress.

16 Magnetic Flux Density (gauss) 2.5 Gauss Probe Hall Effect Chip Magnetic Flux Density (gauss) Gauss Probe Hall Effect Chip (a) Load (kn) (b) Load (kn) Figure 2 Transverse (x-axis) magnetic flux density versus load: (a) neat resin; (b) composite sample (Trovillion et al, 1999) Quattrone, Berman, and White (1998) studied the magnetic response repeatability of MS tagged composites under cyclic loading. Terfenol-D active-tagged composites were subjected to uniaxial tension and the magnetic response in the axial direction under repeated loading and unloading was measured. Two types of experiments were performed: without magnetic annealing between loading cycles (Figure 3a) and with annealing (Figure 3b).

17 Slope 45/10 mg/mpa Slope 140/10 mg/mpa (a) (b) Figure 3 Stress/magnetic flux density cyclic testing results: (a) without annealing between cycles; (b) with annealing between cycles (Quattrone, Berman, and White, 1998) Krishnamurthy, Anjanappa, and Wang (1999) consider health-monitoring detection of delaminations in composite materials using an excitation coil and a sensing coil. The open-circuit voltage induced in the sensing coil is proportional to the stress generated in the magnetostrictive layer by the presence of the delamination. 1.2 CHARACTERISTICS OF ETREMA S TERFENOL-D MAGNETOSTRICTIVE MATERIAL An application manual for designing the magnetostrictive transducers using Etrema Terfenol-D was prepared by J. Butler (1988). A physical description of the magnetostrictive material, with emphasis on Etrema Terfenol-D, was given. The nominal Etrema Terfenol-D strain versus magnetic strength, H, as illustrated in this manual, is shown in Figure 4.

18 Figure 4 Strain versus magnetic strength (Butler, 1988) The theoretical stress versus magnetic flux density, stress versus magnetic strength and magnetic flux density versus magnetic strength curves are presented in Figure 5. Butler (1988) also provides an analysis of the properties of Etrema Terfenol-D with emphasis on strain, stress and magnetic field relationships. Some design considerations and fundamental concepts for designing magnetostrictive transducers are also presented. Figure 5 Relationships between a) magnetic flux density versus magnetic strength, b) stress versus magnetic flux density, c) stress versus magnetic strength 1.3 PRESENT INVESTIGATION In the present investigation, studies of the theoretical and experimental of stressstrain and magnetic flux density response of a woven tagged composite subjected to bending loading were performed. Several two-dimensional simulation models for strain and stress prediction in the composite were investigated. The analysis was performed under two separate assumptions regarding woven composites modeling: (a) balancedorthotropic equivalent and (b) cross-ply equivalent. The experimental investigation considers the bending moment created by the loading and unloading of a simply supported MS-tagged composite beam using several weights. At each loading and unloading step, the magnetic flux density induced by the converse magnetostrictive effect

19 was measured. Comparison between the experimental and analysis results are then presented. Based on the comparison of experimental and analysis data, piezomagnetic stress and strain coefficients under bending are calculated. Comparison of these coefficients with the results obtained by other investigators for axial loading is performed. It is shown that, although exact comparison is not fully possible since bending and axial stress distributions are essentially different, the results are similar. The piezomagnetic stress and strain coefficients determined in this paper for MS tagged composite under bending can be used in industry as reference data for design of novel NDE devices based in the magnetostrictive effect in MS tagged composites.

20 2 ANALYSIS OF A MAGNETOSTRICTIVE COMPOSITE BEAM The analytical work performed on the MS tagged woven composite specimen consisted of: (a) loading analysis; (b) micromechanical analysis; (c) lamination analysis. One balanced-orthotropic and several cross-ply lamination models were used. Converge analysis of the strain and stress predictions versus number of plies in the lamination model were also performed. A magnetostrictive composite stress-strain model has been developed for simulating the behavior of the woven magnetostrictive composite. Two dimensional mathematical model study was performed for the 7 ply balanced-orthotropic model, and 14, 28 and 56-ply multidimensional composite model with 0 /90 and 90 /0 degree. Influence of the Terfenol-D was considered for the first and last ply of the 7-ply model, in first 2 plies and last 2 plies for the 14-ply model, for the first 4 plies and last 4 plies of the 28-ply model and the 8 plies for the 56-ply model. 2.1 STATIC ANALYSIS OF SIMPLY SUPPORTED COMPOSITE BEAM UNDER CENTRAL LOAD To determine the bending moment generated by the centrally-placed load, a static analysis of simply supported composite beam was performed. In this analysis, the dimensions and mechanical properties of the woven composite presented in Table 1 were used. A schematic of the system and how the load is applied is shown in Figure 6. The

21 loading and unloading of the specimen was assumed to be performed using 2-kg weights. Table 1 Dimensions of the MS tagged composite specimen Characteristics Dimension [mm] Span, L 600 Width, w 100 Thickness, t 6.5 Thickness of the lamina Force given by the bricks Strain Gages 2-Kg bricks Composite Material Plate w L Figure 6 Schematic of the static system The force exerted by the one brick is calculated considering the weight of 2 kg and gravitational acceleration of 9.81 m/s 2. Thus the force of the brick is The moment per unit width exerted by the weights is: F = = N (1) FL M 4w (2) Where L is the actual length of the beam, w- the width of the beam, F is the forced applied by one weight. During the loading cycle, up to three weights were sequentially put on the specimen and then removed. The maximum force value, corresponding to

22 three weights, was F max = 58.9 N. The corresponding moment per unit width was M max = Nm/m. 2.2 ANALYSIS OF A MAGNETOSTRICTIVE LAMINATED COMPOSITE Micromechanics analysis of the MS tagged composite specimen Since a woven layer has fibers in two directions, warp and fill, the analysis of a woven composite layer cannot be directly performed through the Classical Lamination Theory (CLT), which assumes fibers aligned with just one direction. An equivalence principle needs to be applied. Tsai (1992) suggested that the predictions of elastic constants and strength of woven composite could be made using classical micro and macro mechanics with appropriate empirical correction factors (Tsai, 1992, page 7-14). One approach is to replace a woven composite layer with two equivalent conventional layers representing the weave and the warp of the original fabric. The micro-mechanics stiffness and strength formulas are then applied to the equivalent plies. One shortcoming of this approach is that the order of 0 and 90 plies may strongly affect the results. Another approach is the replace the balanced woven composite layer by an orthotropic layer with averaged properties. For a woven composite having fibers aligned with the loading axes, this hypothesis yields balanced-orthotropic behavior. In our analysis, both approaches were taken. In either case, a conventional micromechanics analysis (Jones, 1999) was first applied to determine the basic MS composite properties from the properties of its constituents, given in Table 2.

23 Table 2 Mechanical properties considered for the magnetostrictive-tagged woven composite material specimen Mechanical properties Dimension 1 Young s modulus of the fiber, E f 72.4 GPa 1 Young s modulus of the matrix, E m 3.25 GPa 2 Young s modulus for the Terfenol-D, E MS 30 GPa 1 Poisson ratio of the fiber, f Poisson ratio of the matrix, m 0.3 Fiber weight fraction, w f Fiber density, f 2.54 g/ml 1 Matrix density, m 1.18 g/ml 3 Terfenol-D density, MS kg/m 3 1 Thermal expansion coefficient for the fiber, α f m/m per ºC 1 Thermal expansion coefficient for the matrix, α m m/m per ºC Terfenol-D weight fraction 15 Terfenol-D particle size μm Note: 1 Malik, 1992; 2 Butler, 1988; 3 De Laicheisserie, 1993, Trovillion et al., The fiber volume fraction v f is calculated as (Malik, 1992, page 81, equation 2.5): v f wf f wf 1 w f m f (3) For the values given in Table 2, the fiber volume fraction values yield v f = Since the sum of the volume fractions is equal to 1, the volume fraction for the matrix is v m = 1-v f for the plies which does not include the Terfenol-D, and v m = 1- v f - v MS for the plies that contain MS material. The basic moduli of elasticity in the material axes were calculated following (Jones) 1999 as:

24 E E v E v E v 1 f f m m MS MS E [(E ) v (E ) v (E ) v ] f f m m MS MS (4) The ply Poisson s ratio was calculated as: v v v 12 f f m m MS MS E E1 (5) Shear modulus for the fiber, matrix and Terfenol-D for each ply were defined as: G f E f m MS Gm GMS 2 1 2, 1 2, 1 f E m E MS (6) Following Jones (1999, page 134), we calculate the shear modulus for the plies as: G 12 Gf GmGMS v G G v G G v G G f m MS m _ MS f MS MS f m (7) Thermal expansion coefficients for the ply were: E v E v E v 1 E v E v E v f1 f f m1 m m MS MS MS f f m m MS MS (8) 1 v 1 v 1 v 2 f f1 f m m1 m MS MS MS 1 12 (9) and 1 [ ] 2 (10)

25 2.2.2 Lamination Analysis The elastic behavior of multidirectional plies can be described in terms of the stiffness matrix, the compliance matrix, (Jones, 1975). The goal is to determine the strain-stress distribution in the material. The plane stress stiffness matrix [Q] of a orthotropic composite ply is computed as: E1 21E E E G12 (11) [Q] 0 Where E 1 and E 2 are given by Equation (4). The composite is assumed to consist of a layup of multiple plies. The lay-up schematic for generic composite is presented in Figure 7 where x-y are the loading axes. In our case, the angle of the plies is either 0 or 90-degree. θ 1 θ k y x θ n z Figure 7 Schematic of a generic lay-up Considering the transformation matrix:

26 2 2 cos sin 2cossin 2 2 cossin cossin cos sin (12) 2 2 [T] sin cos 2 cos sin The stiffness matrix in loading axes [Q] is: T 1 1 [Q] [T] [Q] [T] (13) Knowing the stiffness matrix, one calculates the extensional stiffness matrix [A], coupling stiffness matrix [B], bending stiffness matrix [D] for the composite laminate: [A] [Q] (z z ) 1 k k1 k k 2 2 [B] [ [Q] k1(z k z k1)] 2 k [D] [ [Q] k1(z k z k1)] 3 k (14) The variable z k represents the distance from the midplane to the bottom of the kth lamina; z k-1 represents the distance from the midplane to the top of the kth lamina. These distances were calculate using the formula: N layer zk h (k ) 2 (15) In matrix form the general equation is: 0 N A B M B D (16) Where {N} is the load vector, {M} is the moment vector, {ε 0 } is the mid-surface strain vector, and {κ} is the curvature vector. For the design purposes, knowing the loading, the

27 state of deformation general solution is: 0 1 A B N B D M (17) Next, the effect of thermal conditions is analyzed by considering only the thermal effect, i.e. without loading F=0. The loads derived from the above-mentioned conditions are: N layer [N] [Q] (z z ) T T k1 k1 k k1 k1 k1 N 2 2 k 1 T k1 layer z z [M] [Q] k1( ) T, k1 2 (18) Here, the term k1represents the coefficient of thermal expansion in the loading axes, which is calculated using the transformation matrix and the coefficients of thermal expansion in the material axes: { } [T] { } T k 1 k 1 k1 (19) The loads and moments due to thermal effect: N layer {N} [Q] { } (z z ) T T k1 k1 k k1 k1 k1 (20) Nlayer 2 2 z k z k1 T k1 k1 k1 k1 2 {M} [Q] { } ( ) T (21) The load and moment vector are composed as: N NN MMM T T (22) A general loading vector it is defined as:

28 N0 N1 N2 {NM} M0 M1 M 2 (23) Solving the state of deformation Equation (17) yields the loading axes strain on the top and bottom of each ply is calculated. { } { } z { } top k 1 0 k 1 { } { } z { } botk1 0 k (24) Using the strains in loading axes, the strains in material axes are calculated for the top and the bottom ply: 1 {} [T] { k 1} T (25) Consequently, the stresses in the longitudinal (L) and transversal (T) directions at the top and the bottom of each ply are: { } [Q] { } topk1 k1 top k 1 (26) { } [Q] { } bot k1 bot k1 k1 (27)

29 2.3 7-PLY BALANCED-ORTHOTROPIC MODEL Ply 1 Ply k Ply n Z X Y Figure 8 Balanced-orthotropic lay-up model used in the analysis The balanced-orthotropic model, presented in Figure 8, uses the assumption that a balanced woven composite can be represented by a composite having averaged properties. Using the longitudinal and transverse moduli for a uni-directional fiber composite layer we calculated the averaged properties: E1 E 1 2 (EL E T) 2 (28) The MS Terfenol-D volume fraction and Poisson s ratio were v MS = 2.24% and ν MS = 0.3 for the plies with Terfenol-D, and zero for the other plies. Two indexes were defined in order to make a distinction between the laminas with and without Terfenol-D, N layer =7, N layers =6.. Performing the CLT analysis yields the balanced-orthotropic (BO) matrices:

30 (N/m) (29) 6 7 [A] BO (N) (30) [B] BO [D] BO (N m) (31) The thermal effect was included for complete analysis procedure but since the cure temperature was the same with environmental temperature, Tk 1 0, then all the stresses due to thermal effects are zero. vector is: Considering only the bending moment, the load vector is zero and the moment M {M vector} 0 0 (Nm/m) (32) The values for this loading vector are: (N/m) NM (N) 0 0 (33)

31 The strain and curvature vector is determined then using the general state of deformation equation mentioned above. Thus the strain and curvatures are: and (1/m) (34) The values for the stresses in longitudinal direction for the top and bottom layer are: L_top (MPa) L _ bot (MPa) (35) The strains in the longitudinal direction for the top and bottom layer are: top (με) bot (με) (36)

32 PLY CROSS-PLY ANGLE MODEL Ply 1 Ply k Ply n Z X Y Figure 9 Cross-ply model used in the analysis In the cross-ply analysis (Figure 9), each woven composite layer is approximated by a couple of unidirectional cross-ply layers. The cross-ply layer can be either 0 /90 or 90 /0. Since the choice between 0 /90 and 90 /0 directly affects the accuracy of the results on the specimen surfaces, we initially considered both cases in our analysis. Additionally, the number of the layers was gradually increased from 14 to 28 and 56, and a convergence study was performed. During this convergence study, as the number of layers were doubled, the thickness of the layer was correspondingly halved, such that the overall thickness of the specimen was maintained. The difference in 14, 28 and 56-ply mathematical models procedure from the balanced-orthotropic model is that the Young s modulus of elasticity in longitudinal and transversal direction were not averaged The 28- ply and 56-ply models were obtained by subsequent subdivision of the 0/90 layers and application of symmetry principles. The ply angle distributions for the 14, 28 and 56-ply

33 models are presented in Table 3: Table 3 Lay-up of the composite models for 14, 28 and 56-ply models Number of Layers Lay-up 14 [(0/90) 3 /0] S 28 [(0/90) 7 /90] S 56 [(0/90) 14 ] S Terfenol-D is assumed present in the first two layers and last two layers and the for the Young s modulus values in the longitudinal and transversal directions are: E 1 = GPa, E 2 = 3.94 GPa. The Young s modulus of elasticity values for the layers that does not contain magnetostrictive material are: E 1 = GPa, E 2 = 3.84 GPa. The A, B and D matrices for the 14-ply model were: [A] (N/m) (37) (N) (38) [B] [D] (N m) (39) Stress and strain results in longitudinal direction for the 14-ply 0 /90 cross-ply laminate, at the top and bottom of the plies, are:

34 L_top (MPa), L_bot (MPa), top (με), bot (με) (40) For the 90 /0 cross-ply laminate the stress and strain results are: L_top (MPa), L_bot (MPa), top (με), bot (με) (41)

35 PLY CROSS-PLY ANGLE MODEL The 28-ply model is computed in a similar manner with the 14-ply model. The changes occur in the indices, the ply angle, and the thickness of the lamina. The A, B and D matrices were: [A] (N/m) (42) [B] (N) (43) [D] (N m) (44) The stress and strain results are presented next:

36 L_top L _ bot top bot (45)

37 PLY CROSS-PLY ANGLE MODEL The 56-ply model was obtained by subsequent subdivision of the 0 and 90 layers and application of symmetry principles. The changes appear in setting the indexes and the thickness of the lamina. The A, B and D matrices for the 56-ply model were: [A] 56 = (N/m) (46) [B] 56= (N) (47) [D] 56 = (N m) (48) The maximum strain on the surface of the specimen model was found to be: = 1296 (49) The corresponding maximum stress was: = 19.6 MPa (50) 2.7 CONVERGENCE ANALYSIS The convergence of the strain and stress distribution for various models was studied. Attention was focused on the longitudinal stress and strain under maximum load

38 conditions (F max = 58.9 N, M max = Nm/m). The strain values predicted on the top and bottom surfaces were also compared with the experimental value measured during the tests described in Chapter Convergence of the strain Figure 10 shows the distribution of the longitudinal strain for the 7-ply balancedorthotropic model and the 14-ply and 28-ply 0 /90 cross-ply models. (The 56-ply results were not plotted on Figure 10 to avoid cluttering the drawing). The sign convention used in Figure 10 is negative to the left and positive to the right. 14 plies 0/90 deg 28 plies 0/90 deg 7-ply balancedorthotropic model Figure 10 Strain distribution of the 7, 14, 28 model with 0 /90 degree ply angle in longitudinal direction The study shows that the largest strain value was predicted by the 7-ply balancedorthotropic model. This value is also very close to the experimental value (1347 με vs με). The convergence of the strain distribution for the cross-ply models is apparent as the models are becoming more refined i.e. strain distribution is getting closer to the experimental data as the number of plies employed in the model is increased. The strain

39 values predicted by each model, at the surface of the composite, and comparison with experimental data, are presented in Table 4. Table 4 Comparison of the longitudinal strain on the specimen surface under maximum load condition (F max = 58.9 N, M max = Nm/m) as predicted by various models and measured by experiment 7-ply balancedorthotropic model 14-ply 0 /90 model Strain [με] 28-ply 0 /90 model 56-ply 0 /90 model Experimental The reason for the difference between the balanced-orthotropic model and the other models is that the balanced-orthotropic model considers each of the plies at the same angle. Thus, the strain that could be carried is larger than the more realistic models in which the strains are lower because the angle of ply is alternating from 0 to 90 degree. The results show that the plies with transverse fiber to the longitudinal direction of measuring (90 degree fiber angle) shrink and oppose the longitudinal stress direction. The convergence of the strain data of the models is presented in Figure 11.

40 Strain [micro strain] Strain response of the models for surface Experimental Baseline Strain Strain response for quasi-isotropic model Number of the plies Figure 11 Convergence of the strain result of the models to the experimental data The figure shows that the maximum strain given by the 7-ply balanced-orthotropic model, is approximately the same as the experimentally obtained strain. The graph also shows that the strain prediction by the cross-ply models improves with increasing number of layers.

41 14-ply 0/90 deg model Balanced orthotropic model 28-ply 0/90 deg model Figure 12 Strain distribution for the 7, 14, 28 mathematical models with 0 /90 degree ply angle in transversal direction. The strain distribution for each model in transversal direction is presented in Figure 12.The strains on the left side of the symmetry line represent positive strains and on the right side represent negative strains. This case shows the strain distribution for the 14 and 28 ply models have close results. It is also observed that for the 7-ply model the strains are larger.

42 2.7.2 Convergence of stress For the 7-ply balanced-orthotropic model, the longitudinal stress distribution is presented in the Figure 11. It can be noticed the stress distribution is rather linear, in accordance with the balanced-orthotropic assumption. The slight discontinuity between the outer first and last layers and the other layers is due to the outer layers being stiffer due to presence of MS tagging in the modified layers, layer 0 and 6. Stresses -Thickness Distribution in Longitudinal Direction for 7 plies Figure 13 Stress thickness distribution of the 7-ply balanced-orthotropic model with 0 degree fiber angle. 15. The stress distribution for the 14-ply 0 /90 cross-ply model is presented in Figure

43 Stress - Thickness Distribution for 14 plies 0/90 degree fiber angle model Figure 14 Stress thickness distribution in longitudinal direction for the 14-ply 0 /90 degree fiber angle model In this model, the stress distribution is different then the stress distribution for the 7-ply balanced-orthotropic model.it can be noticed the stress in the layers with 90-degree fiber angle is lower than the stress in the layers with 0-degree fiber angle. For the layers that contain the magnetostrictive material an increased stress is noticed. The cause of this is the increased stiffness due to Terfenol-D. Similar model behavior is observed for 28 ply and 56-ply 0 /90 cross-ply models as can be shown in Figure 15 and 16.

44 Stress - Thickness Distribution for 28 plies 0/90 degree fiber angle model Figure 15 Thickness distributions of longitudinal stress for 28-ply 0 /90 degree fiber angle model. Stress - Thickness Distribution for 7, 14, 28, 56 plies models Figure 16 Thickness distributions of longitudinal stress for 56-ply 0 /90 degree fiber angle model.

45 The longitudinal stress distribution for the 14, 28-ply 0 /90 cross-ply model is presented in Figure 17. In the cross-ply models, the stress distribution is substantially different then the stress distribution from the 7-ply balanced-orthotropic model. The stresses in the layers with 90-degree fiber angle are substantially smaller. This produce alternating changes from high stress to low stress as the layer stiffness change from the high to low in accordance with the 0 and 90 orientations of the cross-ply laminate. 7 plies balanced-orthotropic model 28 plies 0/90 deg model 14 plies 0/90 deg model Figure 17 Stress-thickness distribution for the 7, 14, 28-ply models with 0 /90 degree fiber angle The 56 ply 0 /90 degree fiber angle model is presented separately since the graph would have been hard to visualize and the stress distribution for each model would have been too hard to follow if they were superimposed on Figure 17. The results reveal that the layers that contain Terfenol-D are stiffer than the other layers. The 7 ply balanced-orthotropic model stress distribution shows that the model does not accurately simulate the real woven composite material because larger stress values were observed than in the cross-ply. This statement is made from a stress distribution

46 point of view. Depending on the fiber orientation and the stress amplitude, the stress carried is consequently proportional. For the 0 degree fiber orientation the stress is larger. For the fiber oriented to a 90-degree angle to the longitudinal direction, the stress amplitude is small. However, this shows that even the transversal fiber carry load. For analysis purposes, Table 4 presents the largest stresses on the surface for each model. Table 5 Stress results on the surface for 7 balanced-orthotropic, 14, 28-ply crossply models 7 ply balancedorthotropic model [MPa] 14 ply 0 /90 model [MPa] 28 ply 0 /90 model [MPa] 56 ply 0 /90 model [MPa] Table 5 shows that the longitudinal stresses obtained by the 14, 28, and 56-ply models are larger than the stresses obtained with the balanced-orthotropic model. These results are due to increased stiffness in the longitudinal direction and less stiffness in transverse direction typical of cross-ply models. In the case of the balanced-orthotropic model the stiffness in both direction are the same, but of lower value than the longitudinal cross-ply stiffness.. Figure 18 presents the results of Table 5 in graphical form. It can be appreciated that the balanced-orthotropic model grossly underestimates the stresses, whereas the cross-ply models shown a definite convergence. It can be estimated that, by further increasing the number of layers in the model, convergence towards an asymptotic value of around 20 MPa would be obtained.

47 Longitudinal stress [MPa] balanced-orthotropic model cross-ply models Number of the layers Figure 18 Convergence of results on the specimen surface under maximum load condition (F max = 58.9 N, M max = Nm/m) as predicted by various models and measured by experiment for stress results Stress analysis for 90 /0 cross-ply models was also performed. As before, the maximum load case (F max = 58.9 N, M max = Nm/m) was considered. The next graphs refer to the mathematical models with 90 /0 degree fiber angle. In Figure 19, 20 and 21 the 14, 28 and 56-ply models with 90 /0 degree fiber angle are presented.

48 Figure 19 Stress thickness distribution in longitudinal direction for 14-ply 90 /0 degree fiber angle model Figure 20 Stress - thickness distribution in longitudinal direction for the 28-ply 90 /0

49 degree fiber angle model Stress - Thickness Distribution in Longitudinal Direction for 56 plies 90/0 Degree Fiber Angle Model Figure 21 Stress thickness distribution in longitudinal direction for 56-ply 90 /0 degree fiber angle model In Figure 22 the stress distributions for the 14 and 28 cross-ply models with 90 /0 degree fiber angle are presented. 56-ply model is not included to not complicate the graph and make it hard to visualize the stress thickness distribution. In this graph, the 7-ply balanced-orthotropic model is also represented. The thickness distribution for the balanced-orthotropic model is similar to 0 /90 models stress distribution. Here is visible, in addition, the influence of the Terfenol-D in increasing the stiffness of the material with results in larger stress

50 7-ply balanced-orthotropic model 14-ply 90/0 deg model 28-0ply90/0 deg model Figure 22 Stress-Thickness Distribution for the 7 ply 0 degree angle model and 14, 28-ply models with 90 /0 degree fiber angle For analysis purposes Table 6 presents the largest stresses on the surface for each model. In this case, the comparison shows that the stress carried in the transverse direction, when the top ply has 90 degree ply angle, is smaller than the stress given by the balanced-orthotropic model. Table 6 Stress results on the surface for 7, 14, and 28 ply 0 and 90 /0 cross-ply models 7 ply balancedorthotropic model [MPa] 14 ply 90 /0 model [MPa] 28 ply 90 /0 model [MPa] 56 ply 90 /0 model [MPa]

51 3 MAGNETOSTRICTIVE COMPOSITE BEAM EXPERIMENT The experimental results alluded to the previous section were obtained during a carefully conducted experiments, as described next. 3.1 SPECIMEN DIMENSIONAL DESIGN For pre-design purpose, the mechanical properties of the specimen were evaluated with a very simplified theory that took into account the following data: 1 Volume fraction of the fiber: v f = Volume fraction of the matrix v m = 1- v f = Young s modulus of the E-glass fiber the: Ef Pa 4 Young s modulus of the resin Em Pa. We assumed an average Young s modulus for random fiber composite calculated with (Malik, 1992, equation 3.42, page 130): 3 5 E E E (51) E 11 and E 22 are calculated using the rule of mixture for longitudinal modulus, E 11, and transverse modulus E 22 (Malik, 1992, equations 3.26,3.27, page 123).

52 E11 Ef vf Em vm (52) E 22 Ef Em E v E v f m m f (53) The result is E Pa. The shear modulus for the composite (Malik, 1992, equation 3.43, page 130) material is assumed to be: 1 1 G E E (54) The result is G Pa. Simple statics was used to determine the bending moment of a simply supported composite beam centrally loaded by a concentrated force. To dimension the specimen for further computations, some mechanical and dimensional quantities were assumed. The equivalent Young s modulus, E, calculated above was considered. The width was measured on the composite specimen, w = 100 mm. The maximum deflection was considered not to exceed the value of q max = 36 mm based work previously done in the area of materials testing. calculated: To determine the thickness of the composite, the length, L to thickness, t, ratio was k 1 3 2wEq max L F k max, t (55) The maximum force was calculated based on the force of one weight and the

53 number of weights to be used in the experiment. The maximum force to be used in the experiment was considered to be the force exerted by ten weights, F max = N. The thickness is then 3F max k t E w a (56) t = mm. The length of the specimen was calculated as: L t k (57) with L = mm, thus t = 6.5 mm. Furthermore the length of the specimen is adjusted to L = mm. Rounding up, the final length for following calculations is L = 600 mm. The design dimensions of the specimen we choose to be: Length 1000 (mm) (58) Width 100 (mm) (59) Thickness 6.75 (mm) (60) The next step was to calculate the moment of inertia for determining the deflection. The moment of inertia was calculated using the well-known formula for a rectangle: w t I 12 3 (61) The result is: I mm. Then the deflection is

54 q F L 3 F E I 24 (62) The result is q F [mm] for F 98.1 (N) (63) At this point, the maximum stress and strain were calculated. Because of symmetry, the stress and strain have same absolute value when equally distanced from mid-plane. The values for the strain are: (MPa) ( ) (64) The maximum stress was max t M 2 I and the maximum strain was max E max. Assuming the flexural strength to be SF Pa and a safety factor of SF=1.5, admissible tensile stress is calculated as

55 SF a SF (65) The result is a Pa. The admissible strain value is a a E (66) ε a = 5617 με. (67) Comparison of the result with the strain predicted above shows that adequate safety is built into the specimen. These simple calculations were used to ensure that the specimen was adequate for experimentation. 3.2 DESCRIPTION OF THE MS TAGGED COMPOSITE SPECIMEN An MS tagged composite specimen was fabricated at Reichhold Chemicals (Raleigh, N.C.) by binding 7 layers of fiberglass woven roving 36oz./sq.yd. and Atlac Urethane-modified Vinyl Ester Resins. The specimen was 1000 mm long, and had a 100 mm by 6.5 mm cross section. MS Terfenol-D tagging powder was used in the two outside layers in the middle 500 mm of the span. Out of 1000 grams of resin, 250 grams had the MS powder. The resin was cured with 1 % MEKP (methyl ethyl ketone peroxide) at room temperature for approximately 90 minutes. The target weight fraction of the glass fibers in the composite was w f = 30 %. The weight fraction of MS Terfenol-D tagging in the resin was 15 %.

56 3.3 SPECIMEN PREPARATION FOR THE EXPERIMENT The specimen was instrumented with strain gages in the mid-span section on the upper and lower surfaces. The specimen surface was prepared for bonding the strain gages using the procedures given by Measurements Group (Raleigh, N.C.). The bonding area was made planar using abrasive sandpaper. Next, the surface was cleaned using Measurements Group surface conditioner M-Prep Conditioner A (a water based acidic surface cleaner). Measurements Group strain gages types CEA UT-120 were applied with M-Bond AE-10 Kit gage adhesive. To protect the strain gages and reduce their exposure to the environment, the strain gages were covered with one sided Scotch tape.

57 3.4 EXPERIMENTAL DESIGN Description of the experiment Lakeshore Gaussmeter Strain Indicator Magnetostrictive Composite Solartron Displacement Transducer National Instrument Amplifier Figure 23 General view of the experiment set-up The overall configuration of the experiment can be seen in Figure 23. The experimental setup permitted the simultaneous measurement of beam deflection, mechanical strains, and magnetic response of the magnetostrictive (MS) tagged beam. Data flow through the data acquisition and processing modules is presented schematically in Figure 24. The MS-tagged composite beam was supported on concrete blocks (500 mm equivalent span) and loaded gradually with an incremental number of clay weights (2 kgf = 19.6 N each). Strain gages were placed on both the upper and lower surfaces of the specimen at the mid-span and were connected in a half bridge configuration to the strain indicator. The magnetic flux density produced from the magnetostrictive particles was measured by the gaussmeter.

58 Strain Gages Strain Indicator LVDT Displacement Transducer SCXI Unit DAQ PC Magnetic Sensor Gaussmeter Figure 24 Schematic of the experiment and data flow For correlation purposes, the mid-span displacement was also measured. An LVDT displacement transducer and a non-magnetic (aluminum and brass) clamping fixture were used. Details of the mid-span instrumentation are given in Figure 25. Initial trials showed that the strain gauge and LVDT electromagnetic fields did not influence the Gaussmeter reading of the magnetostrictively induced magnetic field. The MS tagged composite beam, supported on concrete blocks, was loaded gradually with 2-kg weights. In order to avoid excessive strain in the composite material, the number of loading weights was limited to three. The test procedure was as follow: load the specimen, acquire data and repeat until all the three weights were on the specimen. Next, reverse the procedure and perform unloading. The strain gages were connected to a strain indicator, which provided the actual strain value at top and bottom surfaces. At the same time, the displacements were measured with an LVDT displacement transducer. Hence, the relationship between loads, displacements and strains could be established. Simultaneously, the magnetic field developed by the MS

59 active tagging particles was detected using the gaussmeter. The gaussmeter provided the value of the magnetic flux density read by the probe on the surface of the MS tagged composite material. The data was collected using National Instruments LabView software and associate hardware consisting of a SCXI signal-conditioning module, and the Gateway computer. The information was processed using National Instruments LabView software. The data collection was developed for simultaneous on-line operation with the strain gages, displacement transducer and the gaussmeter. At each loading step, complete data collection was performed. The configuration of the strain gages, displacement transducer and position of the gaussmeter probe are presented in Figure 25. Strain Gage The rod of the Displacement Transducer Gaussmeter Probe Clamping Fixture Figure 25 Position of the gaussmeter probe, strain gage, and displacement transducer Clamping fixture design A clamping device was designed in order to achieve a better position of the

60 displacement transducer on the magnetostrictive composite material.the fixture is composed of 2 aluminum parts (see Figure 25), which are kept attached to the composite material with a brass bolt and nut. A brass bolt was chosen for use in order to diminish the magnetic influence of all the components to the gaussmeter probe Protective wood fixture design for the strengthen the gaussmeter probe A special fixture had to be constructed to ensure proper and consistent alignment of the gaussmeter probe with respect to the composite surface. Its design is presented in Figure 26: 2 mm Sensor Tip Sensor Rod Wood Figure 26 Concept design of the protective wood fixture Considering the fragility of the gaussmeter probe and the advice from Lakeshore Cryotonics, a wood fixture for the protection of the probe tip was designed and fabricated. The distance from the tip of the sensor to the upper surface of the specimen was 2 mm. The open area was covered with clear tape. The sensor rod was attached to the wood fixture using Pro Seal Blue RTV Silicone made by Pacer Technology (Rancho Cucamonga, CA.). A view of this apparatus is given in the Figure 27.

61 Sensor Tip Sensor Rod Machine Carved Wood Figure 27 Overview of the assembly of magnetic sensor and carved wood List of the equipment used In this experiment, many different pieces of equipment were used to achieve the simultaneous measurements of beam deflection, mechanical strains, and magnetic response of the MS tagging particles. The equipment list is presented in Table 7.

TABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA

TABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA TABLE OF CONTENTS 1. INTRODUCTION TO COMPOSITE MATERIALS 1.1 Introduction... 1.2 Classification... 1.2.1

More information

Bending of Simply Supported Isotropic and Composite Laminate Plates

Bending of Simply Supported Isotropic and Composite Laminate Plates Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,

More information

DEVELOPMENT OF THERMOELASTIC STRESS ANALYSIS AS A NON-DESTRUCTIVE EVALUATION TOOL

DEVELOPMENT OF THERMOELASTIC STRESS ANALYSIS AS A NON-DESTRUCTIVE EVALUATION TOOL DEVELOPMENT OF THERMOELASTIC STRESS ANALYSIS AS A NON-DESTRUCTIVE EVALUATION TOOL S. Quinn*, R.K. Fruehmann and J.M. Dulieu-Barton School of Engineering Sciences University of Southampton Southampton SO17

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

Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test

Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test M. Praveen Kumar 1 and V. Balakrishna Murthy 2* 1 Mechanical Engineering Department, P.V.P. Siddhartha Institute of Technology,

More information

Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost

Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost 1 Mihir A. Mehta, 2 Satyen D. Ramani 1 PG Student, Department

More information

Hygrothermal stresses in laminates

Hygrothermal stresses in laminates Hygrothermal stresses in laminates Changing environment conditions (temperature and moisture) have an important effect on the properties which are matrix dominated. Change in temperaturet and moisture

More information

QUESTION BANK Composite Materials

QUESTION BANK Composite Materials QUESTION BANK Composite Materials 1. Define composite material. 2. What is the need for composite material? 3. Mention important characterits of composite material 4. Give examples for fiber material 5.

More information

Composite models 30 and 131: Ply types 0 and 8 calibration

Composite models 30 and 131: Ply types 0 and 8 calibration Model calibration Composite Bi-Phase models 30 and 3 for elastic, damage and failure PAM-CRASH material model 30 is for solid and 3 for multi-layered shell elements. Within these models different ply types

More information

Open-hole compressive strength prediction of CFRP composite laminates

Open-hole compressive strength prediction of CFRP composite laminates Open-hole compressive strength prediction of CFRP composite laminates O. İnal 1, A. Ataş 2,* 1 Department of Mechanical Engineering, Balikesir University, Balikesir, 10145, Turkey, inal@balikesir.edu.tr

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

Non-conventional Glass fiber NCF composites with thermoset and thermoplastic matrices. F Talence, France Le Cheylard, France

Non-conventional Glass fiber NCF composites with thermoset and thermoplastic matrices. F Talence, France Le Cheylard, France 20 th International Conference on Composite Materials Copenhagen, 19-24th July 2015 Non-conventional Glass fiber NCF composites with thermoset and thermoplastic matrices. Thierry Lorriot 1, Jalal El Yagoubi

More information

Composite Structures. Indian Institute of Technology Kanpur

Composite Structures. Indian Institute of Technology Kanpur Mechanics of Laminated Composite Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 23 Analysis of an Orthotropic Ply Lecture Overview Introduction Engineering constants for an 2

More information

SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Mechanical Testing and Diagnosis ISSN 2247 9635, 2012 (II), Volume 3, 79-85 SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Valeriu DULGHERU, Viorel BOSTAN, Marin GUŢU Technical University

More information

Finite element modelling of infinitely wide Angle-ply FRP. laminates

Finite element modelling of infinitely wide Angle-ply FRP. laminates www.ijaser.com 2012 by the authors Licensee IJASER- Under Creative Commons License 3.0 editorial@ijaser.com Research article ISSN 2277 9442 Finite element modelling of infinitely wide Angle-ply FRP laminates

More information

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC. BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally

More information

Lecture 19. Measurement of Solid-Mechanical Quantities (Chapter 8) Measuring Strain Measuring Displacement Measuring Linear Velocity

Lecture 19. Measurement of Solid-Mechanical Quantities (Chapter 8) Measuring Strain Measuring Displacement Measuring Linear Velocity MECH 373 Instrumentation and Measurements Lecture 19 Measurement of Solid-Mechanical Quantities (Chapter 8) Measuring Strain Measuring Displacement Measuring Linear Velocity Measuring Accepleration and

More information

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic

More information

ON THE NUMERICAL ANALYSIS OF COMPOSITE MATERIAL

ON THE NUMERICAL ANALYSIS OF COMPOSITE MATERIAL The 4th International Conference Advanced Composite Materials Engineering COMAT 8- October, Brasov, Romania O THE UMERICAL AALYSIS OF COMPOSITE MATERIAL D.D. icoara Transilvania University, City Brasov,

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

Passive Damping Characteristics of Carbon Epoxy Composite Plates

Passive Damping Characteristics of Carbon Epoxy Composite Plates Journal of Materials Science and Engineering A 6 (-) 35-4 doi:.765/6-63/6.-.5 D DAVID PUBLISHING Passive Damping Characteristics of Carbon Epoxy Composite Plates Dileep Kumar K * and V V Subba Rao Faculty

More information

Getting Started with Composites Modeling and Analysis

Getting Started with Composites Modeling and Analysis Getting Started with Composites Modeling and Analysis IN THIS WEBINAR: PRESENTED BY: Orthotropic materials and how to define them Composite Laminate properties and modeling Composite failure theories and

More information

DRAPING SIMULATION. Recent achievements and future trends. Dr. Sylvain Bel LGCIE University Lyon 1

DRAPING SIMULATION. Recent achievements and future trends. Dr. Sylvain Bel LGCIE University Lyon 1 DRAPING SIMULATION Recent achievements and future trends 1 Dr. Sylvain Bel LGCIE University Lyon 1 2 DRAPING SIMULATION Why? How? What? DRAPING SIMULATION WHY? Clamps Punch Fabric Die 1 2 Resin 3 4 Fig.

More information

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

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric

More information

Composite Structural Mechanics using MATLAB

Composite Structural Mechanics using MATLAB Session 2520 Composite Structural Mechanics using MATLAB Oscar Barton, Jr., Jacob B. Wallace United States Naval Academy Annapolis, Md 21402 Abstract In this paper MATLAB is adopted as the programming

More information

MECHANICAL FAILURE OF A COMPOSITE HELICOPTER STRUCTURE UNDER STATIC LOADING

MECHANICAL FAILURE OF A COMPOSITE HELICOPTER STRUCTURE UNDER STATIC LOADING MECHANICAL FAILURE OF A COMPOSITE HELICOPTER STRUCTURE UNDER STATIC LOADING Steven Roy, Larry Lessard Dept. of Mechanical Engineering, McGill University, Montreal, Québec, Canada ABSTRACT The design and

More information

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011 Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.

More information

Laboratory 4 Bending Test of Materials

Laboratory 4 Bending Test of Materials Department of Materials and Metallurgical Engineering Bangladesh University of Engineering Technology, Dhaka MME 222 Materials Testing Sessional.50 Credits Laboratory 4 Bending Test of Materials. Objective

More information

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of Civil Engineering, IIT Bombay, Powai, Mumbai 400076, India. Tel.022-25767328 email: cejnm@civil.iitb.ac.in Module-13 LECTURE-

More information

THE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF ORIENTED STRANDBOARD: A FINITE ELEMENT SIMULATION

THE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF ORIENTED STRANDBOARD: A FINITE ELEMENT SIMULATION Proceedings of McMat5: 5 Joint ASME/ASCE/SES Conference on Mechanics and Materials June 1-3, 5, Baton Rouge, Louisiana, USA 255 THE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF

More information

Bending Load & Calibration Module

Bending Load & Calibration Module Bending Load & Calibration Module Objectives After completing this module, students shall be able to: 1) Conduct laboratory work to validate beam bending stress equations. 2) Develop an understanding of

More information

Mechanical Behavior of Circular Composite Springs with Extended Flat Contact Surfaces

Mechanical Behavior of Circular Composite Springs with Extended Flat Contact Surfaces Mechanical Behavior of Circular Composite Springs with Extended Flat Contact Surfaces Ping-Cheung Tse epartment of Mechanical Engineering, The Hong Kong Polytechnic University, Hunghom, Kowloon, Hong Kong

More information

Module 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites

Module 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites Module 4: Behaviour of a Laminae-II Learning Unit 1: M1 M4.1 Mechanics of Composites M4.1.1 Introduction to Mechanics of Composites The relation between ply uniaxial strengths and constituent properties

More information

IJSER 1. INTRODUCTION. M.Elhadary

IJSER 1. INTRODUCTION. M.Elhadary 1591 A new failure criterion for GRP composite materials subjected to in-phase and out-of-phase biaxial fatigue loading under different stress ratios M.Elhadary Abstract this studying the fatigue behavior

More information

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms. NASA Langley Research Center. Hampton, VA *University of Denver

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms. NASA Langley Research Center. Hampton, VA *University of Denver Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms W. H. Prosser, M. A. Hamstad + *, J. Gary +, and A. O Gallagher + NASA Langley Research Center Hampton, VA 23681-1 *University of

More information

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES L.V. Smith 1 *, M. Salavatian 1 1 School of Mechanical and Materials

More information

Effects of Resin and Fabric Structure

Effects of Resin and Fabric Structure Fatigue of Wind Blade Laminates: Effects of Resin and Fabric Structure Details David Miller, Daniel D. Samborsky and John F. Mandell Montana State t University it MCARE 2012 Outline Overview of MSU Fatigue

More information

Flexural properties of polymers

Flexural properties of polymers A2 _EN BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS FACULTY OF MECHANICAL ENGINEERING DEPARTMENT OF POLYMER ENGINEERING Flexural properties of polymers BENDING TEST OF CHECK THE VALIDITY OF NOTE ON

More information

SANDWICH COMPOSITE BEAMS for STRUCTURAL APPLICATIONS

SANDWICH COMPOSITE BEAMS for STRUCTURAL APPLICATIONS SANDWICH COMPOSITE BEAMS for STRUCTURAL APPLICATIONS de Aguiar, José M., josemaguiar@gmail.com Faculdade de Tecnologia de São Paulo, FATEC-SP Centro Estadual de Educação Tecnológica Paula Souza. CEETEPS

More information

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Sandhyarani Biswas Department of Mechanical Engineering, N.I.T Rourkela, INDIA Abstract: Now-a-days,

More information

Effect of Thermal Stresses on the Failure Criteria of Fiber Composites

Effect of Thermal Stresses on the Failure Criteria of Fiber Composites Effect of Thermal Stresses on the Failure Criteria of Fiber Composites Martin Leong * Institute of Mechanical Engineering Aalborg University, Aalborg, Denmark Bhavani V. Sankar Department of Mechanical

More information

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,

More information

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,

More information

TESTING AND ANALYSIS OF COMPOSITE SANDWICH BEAMS

TESTING AND ANALYSIS OF COMPOSITE SANDWICH BEAMS TESTING AND ANALYSIS OF COMPOSITE SANDWICH BEAMS I. M. Daniel, J. L. Abot, and K. A. Wang Walter P. Murphy Professor, Departments of Civil and Mechanical Engineering, Robert R. McCormick School of Engineering

More information

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

Evaluation of in-plane orthotropic elastic constants of paper and paperboard Evaluation of in-plane orthotropic elastic constants of paper and paperboard T. Yokoyama and K. Nakai Department of Mechanical Engineering, Okayama University of Science - Ridai-cho, Okayama 7-5, Japan

More information

Chapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.

Chapter 2 - Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2. Chapter - Macromechanical Analysis of a Lamina Exercise Set. The number of independent elastic constants in three dimensions are: Anisotropic Monoclinic 3 Orthotropic 9 Transversely Orthotropic 5 Isotropic.

More information

VALIDATION of CoDA SOFTWARE for COMPOSITES SYNTHESIS AND PRELIMINARY DESIGN (or GETTING COMPOSITES USED - PART 2 )

VALIDATION of CoDA SOFTWARE for COMPOSITES SYNTHESIS AND PRELIMINARY DESIGN (or GETTING COMPOSITES USED - PART 2 ) VALIDATION of CoDA SOFTWARE for COMPOSITES SYNTHESIS AND PRELIMINARY DESIGN (or GETTING COMPOSITES USED - PART 2 ) Graham D Sims and William R Broughton Composites Design Data and Methods, Centre for Materials

More information

MESH MODELING OF ANGLE-PLY LAMINATED COMPOSITE PLATES FOR DNS AND IPSAP

MESH MODELING OF ANGLE-PLY LAMINATED COMPOSITE PLATES FOR DNS AND IPSAP 16 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS MESH MODELING OF ANGLE-PLY LAMINATED COMPOSITE PLATES FOR DNS AND IPSAP Wanil Byun*, Seung Jo Kim*, Joris Wismans** *Seoul National University, Republic

More information

VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS

VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP

More information

Most of the material in this package is based on a recently published book. This is:

Most of the material in this package is based on a recently published book. This is: Mechanics of Composite Materials Version 2.1 Bill Clyne, University of Cambridge Boban Tanovic, MATTER Assumed Pre-knowledge It is assumed that the student is familiar with simple concepts of mechanical

More information

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations W. Van Paepegem *, I. De Baere and J. Degrieck Ghent

More information

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP) KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP) 1 University of Science & Technology Beijing, China, niukm@ustb.edu.cn 2 Tsinghua University, Department of Engineering Mechanics, Beijing, China,

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

American Society for Testing and Materials (ASTM) Standards. Mechanical Testing of Composites and their Constituents

American Society for Testing and Materials (ASTM) Standards. Mechanical Testing of Composites and their Constituents Mechanical Testing of Composites and their Constituents American Society for Testing and Materials (ASTM) Standards Tests done to determine intrinsic material properties such as modulus and strength for

More information

Overview of Hygrothermally Stable Laminates with Improved Extension-twist Coupling

Overview of Hygrothermally Stable Laminates with Improved Extension-twist Coupling Overview of Hygrothermally Stable Laminates with Improved Extension-twist Coupling R. Haynes epartment of erospace Engineering, Georgia Institute of Technology 27 Ferst rive W, tlanta, G 3332-5, US robert.haynes@gatech.edu

More information

HIGH VELOCITY IMPACT ON TEXTILE REINFORCED COMPOSITES

HIGH VELOCITY IMPACT ON TEXTILE REINFORCED COMPOSITES 16 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS HIGH VELOCITY IMPACT ON TEXTILE REINFORCED COMPOSITES Warnet L., Akkerman R., Ravensberg M. University of Twente, Faculty of Engineering Technology,

More information

LAMINATED COMPOSITE PLATES

LAMINATED COMPOSITE PLATES LAMINATED COMPOSITE PLATES David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 February 10, 2000 Introduction This document is intended

More information

Module-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates

Module-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination

More information

INVESTIGATION OF THE PROCESSING PARAMETERS OF A 3D WOVEN REINFORCEMENT

INVESTIGATION OF THE PROCESSING PARAMETERS OF A 3D WOVEN REINFORCEMENT INVESTIGATION OF THE PROCESSING PARAMETERS OF A 3D WOVEN REINFORCEMENT Andreas Endruweit, Dhiren K. Modi and Andrew C. Long School of Mechanical, Materials and Manufacturing Engineering, University of

More information

Assessment Methods of Mechanical Properties of Composite Materials

Assessment Methods of Mechanical Properties of Composite Materials Mechanics and Mechanical Engineering Vol. 21, No. 4 (2017) 1005 1018 c Lodz University of Technology Assessment Methods of Mechanical Properties of Composite Materials Monika Kamocka Radoslaw J. Mania

More information

THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH

THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH M. R. Wisnom University of Bristol Advanced Composites Centre for Innovation and Science University Walk, Bristol BS8 1TR, UK M.Wisnom@bristol.ac.uk

More information

PLY LEVEL UNCERTAINTY EFFECTS ON FAILURE OF COMPOSITE

PLY LEVEL UNCERTAINTY EFFECTS ON FAILURE OF COMPOSITE 7th European Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France More Info at Open Access Database www.ndt.net/?id=17206 PLY LEVEL UNCERTAINTY EFFECTS ON FAILURE OF COMPOSITE

More information

Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA

Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA What programs are in PROMAL? Master Menu The master menu screen with five separate applications from

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 2, Issue 4, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 2, Issue 4, July 2013 Delamination Studies in Fibre-Reinforced Polymer Composites K.Kantha Rao, Dr P. Shailesh, K. Vijay Kumar 1 Associate Professor, Narasimha Reddy Engineering College Hyderabad. 2 Professor, St. Peter s Engineering

More information

Computational Analysis for Composites

Computational Analysis for Composites Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics

More information

MECHANICS OF MATERIALS

MECHANICS OF MATERIALS Third E CHAPTER 2 Stress MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech University and Strain Axial Loading Contents Stress & Strain:

More information

Materials and Structures. Indian Institute of Technology Kanpur

Materials and Structures. Indian Institute of Technology Kanpur Introduction to Composite Materials and Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 34 Thermal Stresses in Plates Lecture Overview Introduction Mechanical and Thermal Strains

More information

Experiment Five (5) Principal of Stress and Strain

Experiment Five (5) Principal of Stress and Strain Experiment Five (5) Principal of Stress and Strain Introduction Objective: To determine principal stresses and strains in a beam made of aluminum and loaded as a cantilever, and compare them with theoretical

More information

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

ENSC387: Introduction to Electromechanical Sensors and Actuators LAB 3: USING STRAIN GAUGES TO FIND POISSON S RATIO AND YOUNG S MODULUS ENSC387: Introduction to Electromechanical Sensors and Actuators LAB 3: USING STRAIN GAUGES TO FIND POISSON S RATIO AND YOUNG S MODULUS 1 Introduction... 3 2 Objective... 3 3 Supplies... 3 4 Theory...

More information

The science of elasticity

The science of elasticity The science of elasticity In 1676 Hooke realized that 1.Every kind of solid changes shape when a mechanical force acts on it. 2.It is this change of shape which enables the solid to supply the reaction

More information

Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load

Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Sriram Chintapalli 1, S.Srilakshmi 1 1 Dept. of Mech. Engg., P. V. P. Siddhartha Institute of Technology.

More information

In Situ Ultrasonic NDT of Fracture and Fatigue in Composites

In Situ Ultrasonic NDT of Fracture and Fatigue in Composites ECNDT 26 - Mo.2.6.5 In Situ Ultrasonic NDT of Fracture and Fatigue in Composites I. SOLODOV, K. PFLEIDERER, and G. BUSSE Institute for Polymer Testing and Polymer Science (IKP), Non-destructive Testing

More information

GATE SOLUTIONS E N G I N E E R I N G

GATE SOLUTIONS E N G I N E E R I N G GATE SOLUTIONS C I V I L E N G I N E E R I N G From (1987-018) Office : F-16, (Lower Basement), Katwaria Sarai, New Delhi-110016 Phone : 011-65064 Mobile : 81309090, 9711853908 E-mail: info@iesmasterpublications.com,

More information

A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites.

A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites. A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites. P.F.Giddings, C.R.Bowen, H.A.Kim University of Bath, UK Dept. Mech. ng, University

More information

GB/T / ISO 527-1:1993

GB/T / ISO 527-1:1993 Translated English of Chinese Standard: GB/T1040.1-2006 www.chinesestandard.net Sales@ChineseStandard.net GB NATIONAL STANDARD OF THE PEOPLE S REPUBLIC OF CHINA ICS 83.080.01 G 31 GB/T 1040.1-2006 / ISO

More information

Module 7: Micromechanics Lecture 25: Strength of Materials Approach. Introduction. The Lecture Contains. Effective Transverse Modulus

Module 7: Micromechanics Lecture 25: Strength of Materials Approach. Introduction. The Lecture Contains. Effective Transverse Modulus Introduction In the previous lecture we have introduced the concept of Representative Volume Element or Unit Cell. This is the basic building block in a micromechanical study. Further, we explained the

More information

ME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2)

ME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2) ME 582 Advanced Materials Science Chapter 2 Macromechanical Analysis of a Lamina (Part 2) Laboratory for Composite Materials Research Department of Mechanical Engineering University of South Alabama, Mobile,

More information

1. Demonstrate that the minimum cation-to-anion radius ratio for a coordination number of 8 is

1. Demonstrate that the minimum cation-to-anion radius ratio for a coordination number of 8 is 1. Demonstrate that the minimum cation-to-anion radius ratio for a coordination number of 8 is 0.732. This problem asks us to show that the minimum cation-to-anion radius ratio for a coordination number

More information

Strength of GRP-laminates with multiple fragment damages

Strength of GRP-laminates with multiple fragment damages Strength of GRP-laminates with multiple fragment damages S. Kazemahvazi, J. Kiele, D. Zenkert Kungliga Tekniska Högskolan, KTH 100 44 Stockholm, Sweden sohrabk@kth.se SUMMARY The strength of glass fibre

More information

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques 9 th International LS-DYNA Users Conference Optimization Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques Stefano Magistrali*, Marco

More information

WORKBOOK MECHANICS OF MATERIALS AND ELEMENTS OF ENGINEERING STRUCTURES

WORKBOOK MECHANICS OF MATERIALS AND ELEMENTS OF ENGINEERING STRUCTURES WORKBOOK MECHANICS OF MATERIALS AND ELEMENTS OF ENGINEERING STRUCTURES LUBLIN 014 Authors: Sylwester Samborski, Andrzej Teter and Marcin Bocheński Desktop publishing: Sylwester Samborski, Andrzej Teter

More information

Stress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina

Stress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina Chapter 2 Macromechanical Analysis of a Lamina Stress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina Islamic

More information

6. Bending CHAPTER OBJECTIVES

6. Bending CHAPTER OBJECTIVES CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where

More information

NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS

NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS 5 th European LS-DYNA Users Conference Composites NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS Kevin Brown 1, Richard Brooks, Nicholas Warrior School of Mechanical, Materials and

More information

AERSYS KNOWLEDGE UNIT

AERSYS KNOWLEDGE UNIT -7016 1. INTRODUCTION The scope of this document is to provide a clarification and a deeper understanding of the two different ways to move the mid plane of the element out of the nodal plane. Although

More information

Module 5: Laminate Theory Lecture 17: Laminate Constitutive Relations. The Lecture Contains: Laminate Constitutive Relations

Module 5: Laminate Theory Lecture 17: Laminate Constitutive Relations. The Lecture Contains: Laminate Constitutive Relations Lecture 17: Laminate Constitutive Relations The Lecture Contains: Laminate Constitutive Relations Classification of Laminates Cross-Ply Laminates Specially Orthotropic Laminates Examples Homework References

More information

Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics

Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Dr. A. Johnson DLR Dr. A. K. Pickett ESI GmbH EURO-PAM 99 Impact and Crash Modelling of Composite Structures: A Challenge

More information

A SELF-INDICATING MODE I INTERLAMINAR TOUGHNESS TEST

A SELF-INDICATING MODE I INTERLAMINAR TOUGHNESS TEST A SELF-INDICATING MODE I INTERLAMINAR TOUGHNESS TEST P. Robinson The Composites Centre, Department of Aeronautics, Imperial College London South Kensington, London, SW7 2AZ, UK p.robinson@imperial.ac.uk

More information

Massachusetts Institute of Technology Department of Aeronautics and Astronautics Cambridge, MA Problem Set 14

Massachusetts Institute of Technology Department of Aeronautics and Astronautics Cambridge, MA Problem Set 14 Massachusetts Institute of Technology Department of Aeronautics and Astronautics Cambridge, MA 02139 16.01/16.02 Unified Engineering I, II Fall 2003 Problem Set 14 Name: Due Date: 12/9/03 F18 F19 F20 M19

More information

MULTI-SCALE MODELLING OF FIBRE BUNDLES

MULTI-SCALE MODELLING OF FIBRE BUNDLES THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS MULTI-SCALE MODELLING OF FIBRE BUNDLES N. D. Chakladar 1, P. Mandal 1 *, P. Potluri 2 1 School of Mechanical, Aerospace and Civil Engineering,

More information

Multi Disciplinary Delamination Studies In Frp Composites Using 3d Finite Element Analysis Mohan Rentala

Multi Disciplinary Delamination Studies In Frp Composites Using 3d Finite Element Analysis Mohan Rentala Multi Disciplinary Delamination Studies In Frp Composites Using 3d Finite Element Analysis Mohan Rentala Abstract: FRP laminated composites have been extensively used in Aerospace and allied industries

More information

Composites Design and Analysis. Stress Strain Relationship

Composites Design and Analysis. Stress Strain Relationship Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines

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

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

Advanced Structural Analysis EGF Section Properties and Bending

Advanced Structural Analysis EGF Section Properties and Bending Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear

More information

LAMINATED COMPOSITE PLATES

LAMINATED COMPOSITE PLATES LAMINATED COMPOSITE PLATES David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 February 10, 2000 Introduction This document is intended

More information

A STRUCTURE DESIGN OF CFRP REAR PRESSURE BULKHEAD WITHOUT STIFFENERS

A STRUCTURE DESIGN OF CFRP REAR PRESSURE BULKHEAD WITHOUT STIFFENERS Xi an, 2-25 th August 217 A STRUCTURE DESIGN OF CFRP REAR PRESSURE BULKHEAD WITHOUT STIFFENERS LI Zhongyang 1, LI Dong 2 Mailbox72-35, Yanliang District, Xian, China, Email: zhongyangli@live.com Keywords:

More information

Benchmarking study of steel-composite structures in CAE crash applications. Master s thesis in Applied Mechanics MADELEINE ANDERSSON EMMA LARSSON

Benchmarking study of steel-composite structures in CAE crash applications. Master s thesis in Applied Mechanics MADELEINE ANDERSSON EMMA LARSSON Benchmarking study of steel-composite structures in CAE crash applications Master s thesis in Applied Mechanics MADELEINE ANDERSSON EMMA LARSSON Department of Applied Mechanics CHALMERS UNIVERSITY OF TECHNOLOGY

More information

STRENGTH OF MATERIALS-I. Unit-1. Simple stresses and strains

STRENGTH OF MATERIALS-I. Unit-1. Simple stresses and strains STRENGTH OF MATERIALS-I Unit-1 Simple stresses and strains 1. What is the Principle of surveying 2. Define Magnetic, True & Arbitrary Meridians. 3. Mention different types of chains 4. Differentiate between

More information