EVALUATION OF THERMAL CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PLASTICS

Similar documents
Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate

Journal of Solid Mechanics and Materials Engineering

Long-term Life Prediction of CFRP Structures Based on MMF/ATM Method

A method to determine relative stroke detection efficiencies from multiplicity distributions

ID 1163 CURE MONITORING OF THERMOSETTING RESIN COMPOSITES BY LACOMTECH DIELECTROMETRY

SHAPE IDENTIFICATION USING DISTRIBUTED STRAIN DATA FROM EMBEDDED OPTICAL FIBER SENSORS

STUDY OF THERMAL DIFFUSIVITY IN HEAT-INSULATING MATERIALS

Predictions of Storage Modulus of Glass Bead-Filled Low-Density-Polyethylene Composites

Laminated Composite Plates

Density. Density 22/10/2015. Rules of Mixture for Elastic Properties. Materi ke-6

Determination of Cold Spot Location for Conduction-Heated Canned Foods Using an Inverse Approach

Estimating flow properties of porous media with a model for dynamic diffusion

ESTIMATION OF THE VISCOELASTIC PARAMETERS OF LAMINATED COMPOSITES. PART I. ANALYTICAL CONSIDERATIONS

NUMERICAL SIMULATION AND PARAMETRIC STUDY OF LAMINAR MIXED CONVECTION NANOFLUID FLOW IN FLAT TUBES USING TWO PHASE MIXTURE MODEL

SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD

A NEW APPROACH TO DYNAMIC BUCKLING LOAD ESTIMATION FOR PLATE STRUCTURES

Critical Analysis of FRP-NSM double-shear pull-out tests. Abstract

High Temperature Chemical Processes, Stadiou str., Platani, Patras, Greece

Supplement of Detailed characterizations of the new Mines Douai comparative reactivity method instrument via laboratory experiments and modeling

Effect of fiber volume fraction on fracture mechanics in continuously reinforced fiber composite materials

Supplementary Materials for

DIRECT NUMERICAL SIMULATION OF DAMAGE PROGRESSION IN LAMINATED COMPOSITE PLATES USING MULTI-SCALE MODELLING

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles

Application of Eddy Current Damping Effect to Design a Novel Magnetic Damper

Determination of fiber/matrix interface debond growth parameters from cyclic loading of single fiber composites

Adsorption of 2,4-Xylidine on Granulated Activated Carbon in a Packed Bed Absorber. f f

Chapter 3 Micromechanical Analysis of a Lamina Ultimate Strengths of a Unidirectional Lamina

Strain Rate and Temperature Effects on the Nonlinear Behavior of Woven Composites

Micromechanical models for stiffness prediction of alumina trihydrate (ATH)

REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING MEDIA

Modelling of the Through-air Bonding Process

EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE

Dual porosity DRM formulation for flow and transport through fractured porous media

Part I: How Dense Is It? Fundamental Question: What is matter, and how do we identify it?

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Simplified Explicit Model to Measure Transient Heat Transfer in Foamcrete Panel System

EVALUATION OF THERMOELASTIC PROPERTIES OF CARBON NANOTUBE-BASED COMPOSITES USING FINITE ELEMENT METHOD

Accuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates

Experimental studies on equivalent thermal properties of particle-reinforced flexible mould materials

Analytical investigation of unsteady CuO nanofluid flow, heat and mass transfer between two parallel disks

Determination of the Young's modulus of an aluminium specimen

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

Multi-Scale Statistical Approach of the Elastic and Thermal Behavior of a Thermoplastic Polyamid-Glass Fiber Composite

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

A Fiber Optic Phase Modulator with an Optical Frequency Shift of up to 20 GHz

Revealing spatially heterogeneous relaxation in a model. nanocomposite

Functions: Review of Algebra and Trigonometry

CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION

1 Analysis of heat transfer in a single-phase transformer

Modeling and Analysis of Thermal Bimorph Using COMSOL

Transfer from Arbitrary Pure State to Target Mixed State for Quantum Systems

Farid Samara 1, Dominic Groulx 1 and Pascal H. Biwole 2 1

SIMULATION OF SOIL TEMPERATURE VARIATION FOR GEOTHERMAL APPLICATIONS

Modelling of damage in composite materials using interface elements

EFFECT OF SHEAR STRESS ON FATIGUE LIFE OF COMPOSITE LAMINATES

Materials and Design

Structural optimization of an automobile transmission case to minimize radiation noise using the model reduction technique

Comparison Studies on Dynamic Packaging Properties of Corrugated Paperboard Pads

PROBLEM 7.2 1/3. (b) The local convection coefficient, Eq. 7.23, and heat flux at x = L are 1/2 1/3

MODELING OF SUBSURFACE BIOBARRIER FORMATION

HEAT CONDUCTION IN ISOTROPIC HETEROGENEOUS MEDIA

CHAPTER 5 THEORETICAL MODELING OF ELASTIC PROPERTIES

Dynamic buckling of laminated composite stringer stiffened CFRP panels under axial compression

CFD SIMULATION OF PRESSURE LOSS IN HVDC TRANSFORMER WINDING

School of Aerospace Engineering Equilibrium Diagrams and Saturated Liquid/Vapor Systems

Chapter 3 Micromechanical Analysis of a Lamina Volume Fractions, Weight Fractions, Density, and Void Content

A Pulley System Apparatus for a Laboratory Experience in Dynamics

TIP TRAJECTORY TRACKING OF FLEXIBLE-JOINT MANIPULATORS

CFD SIMULATION OF A MEMBRANE DISTILLATION MODULE CHANNEL

DEPENDENCE OF PISTON RING PARAMETERS DETERMINATION ACCURACY ON THE METHOD OF CALCULATION

The Development of a Generalised Jet Mixing Model for the Derivation of Exclusion Zones in Hazardous Area Classification

Excitability of guided waves in composites with PWAS transducers

L2 IntroSurvey. EIEN20 Design of Electrical Machines, IEA, Today s goal. L2: Variety of electrical machines.

An Optimal Analysis of Flow and Heat Transfer over a Slender Permeable Elastic Sheet with Variable Fluid Properties

Solidification of Porous Material under Natural Convection by Three Phases Modeling

DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION

Archimedes Principle

PULSED POWER ELECTROMECHANICS - PERMANENT MAGNETS VERSUS COPPER COILS

Fatigue Life Prediction in Composite Materials

Micro-Crack Initiation and Propagation in Fiber Reinforced Composites. Andrejs Pupurs

H S = Hashin-Shtrikman lower bound of bulk modulus, N/mm 2 K U G L G U. H = Hashin-Shtrikman lower bound of shear modulus, N/mm 2

Physics 41 Homework #2 Chapter 16. fa. Here v is the speed of the wave. 16. The speed of a wave on a massless string would be infinite!

2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL

Arithmetic Unit for Complex Number Processing

Block failure in connections - including effets of eccentric loads

Numerical Modeling of Self-Compacting Mortar Flow Using Discrete Element Method

S. Srinivasan, Technip Offshore, Inc., Houston, TX

WIRELESS DRIVE OF PIEZOELECTRIC COMPONENTS SATYANARAYAN BHUYAN 2011

HALF-LIFE TIME FOR VOC EMISSION AND SORPTION OF POROUS BUILDING MATERIALS

MHD Flow of Thixotropic Fluid with Variable Thermal Conductivity and Thermal Radiation

Data-Driven Imaging in Anisotropic Media

Title. Author(s)Izumida, Yuki; Okuda, Koji. CitationPhysical review E, 80(2): Issue Date Doc URL. Rights. Type.

Characterizing and Modeling Mechanical Properties of Nanocomposites- Review and Evaluation

Heat Transfer Enhancement in channel with obstacles

Analytical Study of AC Magnetic Susceptibility of (Bi, Pb) Sr-Ca-Cu-O Superconducting Systems Using Bean Critical State Model

Random Vibration Fatigue Analysis with LS-DYNA

Optimization of Flywheel Weight using Genetic Algorithm

National 5 Summary Notes

Optimum Design of Assembled Cavity Dies for Precision Forging Process

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

Transcription:

16 TH INTERNATIONA CONFERENCE ON COMPOSITE MATERIAS EVAUATION OF THERMA CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PASTICS Shinji Ogihara*, Makoto Yaaguchi**, Takahito Chiba**, Junichi Shiizu****, Yoji Okabe**** and Nobuo Takeda****, *Departent o Mechanical Engineering, Tokyo University o Science, Chiba, Japan **Graduate Student, Tokyo University o Science, Chiba, Japan ***Nippon Steel Materials Copany iited **** Departent o Aeronautics and Astronautics, The University o Tokyo Keywords: Theral conductivity, Theral resistance, CFRP, Unidirectional coposite Abstract Theral conductivities in longitudinal and transverse directions o a CFRP unidirectional coposite are evaluated experientally. aser lash ethod is eployed to easure theral conductivity. To discuss the experiental results, rules o ixture predictions and two-diensional steady-state and transient heat transer siulations are conducted. In addition, a heat sink is ade with the CFRP. The heat radiation perorance under natural convection (the theral resistance) is copared with heat sinks ade o aluinu with dierent surace treatent. Eect o heat radiation is discussed both experientally and analytically. Three-diensional steady-state heat transer siulations are perored to copare with the experient result. 1 Introduction Carbon iber reinorced plastics (CFRP) are used in various ields, especially in aerospace engineering, due to their high speciic odulus and strength. Because pitch-based carbon ibers have high theral conductivity, their coposites are expected to be used as heat control aterials [1]. Considering a carbon iber unidirectional coposite, it ay have anisotropy in theral conductivity, that is, it is high in iber direction whereas low in the direction perpendicular to the ibers. Moreover, it is expected that we can develop a aterial with higher theral conductivity. In the present study, theral conductivities in pitch-based carbon iber (YS9A) reinorced epoxy unidirectional coposites with dierent iber volue ractions are experientally evaluated by using the laser lash ethod. To discuss the experiental results, rules o ixture prediction and two-diensional steady-state and transient heat transer siulations are conducted. In addition, a heat sink is ade with the CFRP. The heat radiation perorance under natural convection (the theral resistance) is copared with heat sinks ade o aluinu with dierent surace treatent, that is, with and without aluite treatent. To discuss the experiental results, three-diensional steady-state heat transer siulations are conducted. Experient.1 Materials A pitch-based carbon iber, YS9A (Nippon Graphite Fiber td.), is used. An epoxy resin is used as the atrix aterial. Unidirectional coposites with dierent iber volue ractions (55%, 4%, % and % (neat resin)) are abricated or the theral conductivity easureent by the laser lash ethod. For theral resistance easureent, CFRP heat sink o 55% volue raction is used. Both speciens whose iber directions are parallel and perpendicular to the heater plane are used. Properties o iber and atrix in the present study are listed in Table 1. Table 1 Properties o iber and atrix Theral Speciic heat Conductivity capacity (W/K) (J/gK) Density (g/c3) Fiber 5.837 1.75 Epoxy.6 1.1 1.3 1

SHINJI OGIHARA, Makoto Yaaguchi, Takahito Chiba, Junichi Shiizu Yoji Okabe and Nobuo Takeda. Measureent o theral conductivity by the laser lash ethod Theral conductivities in the unidirectional coposites are easured by using the laser lash ethod. In the ethod, a laser bea is lashed on the top surace o a cylindrical specien (Fig.1). The teperature change at the botto surace o the specien is onitored. Based on the easured teperature history, the haltie t 1/ (tie or the teperature rise up to the hal o the axiu teperature) is deterined, which in turn used to calculate specien theral diusivity, α, by using the ollowing equation. α =.1388 (1) 1/ t where is specien length. The diaeter and the length o the speciens are 1 and 3, respectively. Based on the theral diusivity, theral conductivity can be calculated by using ollowing equation. λ = αcρ () where C is speciic heat capacity and ρ is density. Specien aser pulse Inrared sensing device α =.1388 t α:theral diusivity : Specien length t / :Haltie T ax T ax Fig.1 Principle o the aser Flash Method.3 Heat sinks A heat sink is ade with a CFRP prepreg ade o a carbon iber, YS9A (Nippon Graphite Fiber td.), and Epoxy resin. Prepreg sheets are stacked such that the iber direction is parallel to the in height. The heat sink is ade with integrated olding. The size o the heat sink is shown in Fig.. Coercially available heat sinks o the sae size ade o aluinu alloy (Ryosan Inc.) are also used. Both aluinu heat sinks with and without aluite treatent are used. Properties o each heat sink are shown in Table. 1/ T t 1/ t Fig. Size o heat sink Table Properties o heat sinks Al Al(A663-T5)- (A663-T5) aluite treatent CFRP (YS9A) Density (g/c 3 ).7.7 1.75 Theral conductivity (W/K) 9 9 7 Mass (g) 44.6 44.8 3.8.4 Measureent o theral resistance Heat radiation perorance was evaluated by easuring theral resistance under the condition o natural convection. The theral resistance R is deined by ( T1 T ) R( / W ) = (3) Q where T 1 is the highest teperature at the botto surace o the heat sink, T is the abient teperature and Q is the heat lux which into the heat sink. The experiental setup used in the present study is shown in Fig.3. A rubber heater was used as a heat source. Heat lux sensors were used to easure the heat lux which lows in to the heat sink, Q. therocouples o type T were used to the teperatures T 1 and T. The easureents were conducted until the syste reaches the steady state. Therocouple logger Notebook coputer DC power source Heat sink Radiation sheet Styrooa Copper Heater Thero-couple (T-type) Heat lux sensor Fig.3 Experiental setup or theral resistance easureent A V

EVAUATION OF THERMA CONDUCTIVITY IN PITCH-BASED CARBON FIBER REINFORCED PASTICS 3 Analysis 3.1 Rules o ixture predictions Theral conductivity o a unidirectional coposite in the iber direction, λ, can be predicted by the rules o ixture as λ λ V + λ V = (4) conducted ro the next expression by using the value o this heat lux. q λ = (6) T / d Where ΔT is the teperature dierence, that is 3-1=, and d is the odel thickness. λ is theral conductivity, V is volue raction and subscripts and denote iber and atrix, respectively. Theral conductivity o a transverse coposite in the iber direction, λ T, can be predicted by the rules o ixture as 3.5μ 1 λ T V = λ V + λ (5) Fiber direction Analytical odels o the rules o ixture are shown in Fig.4. 7μ λ λ Transverse direction iber atrix V 1-V Fig.5 Models o two-diensional steady-state heat transer siulation (FEM) λ λ 1-V V :iber :atrix Fig.4 Models o the rules o ixture prediction 3. Calculation o theral conductivity by twodiensional steady-state heat transer siulation Coposite theral conductivity is calculated by using the two-diensional steady-state heat transer siulation. Finite eleent ethod is applied. Analytical odels are shown in Fig.5. The top edge o the odel is ixed at 3 and the botto edge is ixed at 1, and the right and let edges are therally isolated. Heat lux q was calculated by suing the node heat lux at the botto surace []. The calculation o theral conductivity λ was 3.3 Calculation o theral conductivity by twodiensional transient heat transer siulation Coposite eective theral conductivity is also evaluated by using a two-diensional transient heat transer siulation. In the analysis, the laser lash ethod is siulated. Analytical odels are shown in Fig.6. Finite eleent ethod is applied. To siulate laser lash ethod, the top layer o 3μ thick is set 1 degrees hotter as the initial condition, and the right and let edge are therally isolated. We onitor the teperature rise at the botto edge o the odel which enables us to deterine the haltie t 1/. Based on the haltie calculated, we can deterine the coposite eective theral diusivity. 3

SHINJI OGIHARA, Makoto Yaaguchi, Takahito Chiba, Junichi Shiizu Yoji Okabe and Nobuo Takeda 3.5μ 3.5μ Fiber direction Fig.6 Models o Two-diensional transient heat transer siulation (FEM) To calculate the theral conductivity, coposite eective speciic heat capacity C c and density ρ c are deterined by using the ollowing rules o ixture as C ρ V ρ V C = C + C (7) ρv + ρv ρv + ρv ρ ρ V + ρ V c = (8) 7μ Transverse direction where λ is theral conductivity, V is volue raction and subscripts and denote iber and atrix, respectively. 3.4 Calculation o theral resistance by threediensional steady-state thero-luid analysis The heat radiation perorances o each heat sink are calculated by steady-state thero-luid analysis. To calculate the theral resistance, coposite eective speciic heat capacity C c and density ρ c are deterined by Eq. 7 and Eq. 8. Analytical odels are shown in Fig.7. To siulate the experiental setup, all the coponents, the heater, the heat radiation sheets, the copper and the heat sink, and the styrene oa are odeled. Only the hal part o the experiental setup is odeled due to syetry. iber atrix Fig.7 Model o three-diensional steady-state thero-luid analysis As the initial condition, the odel is set to be 4. The eect o radiation is considered. The eissivities o CFRP, Aluinu alloy with and without aluite treatent are.94,.8, and.13, respectively. The gravitational acceleration was set to be 9.8 /s. The theral resistance is calculated ro Eq. 3. 4 Results and discussion 4.1 aser lash ethod and analytical result Table 3 and Fig.8 shows the experiental results o the unidirectional coposite theral conductivity in the iber direction and results o the predictions o the analysis as a unction o iber volue raction. In the case o iber direction, the theral conductivity increases linearly with increasing iber volue raction. Experiental results o theral conductivity o CFRP agree with the rules o ixture predictions and the predictions o both two-diensional steady-state and transient heat transer analysis. Table 3 Results o theral conductivity in iber direction easureents and analytical predictions Volue raction (%) 55 4 Experient (W/K) 7.44 19.1 91.63 Rules o ixture prediction (W/K) 75.1.16 1.1 Steady-state siulation (W/K) 75.1.16 1.1 Transient siulation (W/K) 71.55 197.44 97.47 4

EVAUATION OF THERMA CONDUCTIVITY IN PITCH-BASED CARBON FIBER REINFORCED PASTICS Theral conductivity (W/K) 4 3 1 experient (YS9A) rules o ixture prediction (YS9A) steady-state siulation (YA9A) transient siulation (YS9A) 1 3 4 5 6 Fiver volue raction (%) Theral conductivity (W/K) 1.5 1.5 experient (YS9A) rules o iwture prediction (YS9A) steady-state siulation (YS9A) transient siulation (YS9A) 1 3 4 5 6 Fiver volue raction (%) Fig.8 Theral conductivity o CFRP in iber direction as a unction o iber volue raction or CFRP using YS9A.Coparison with the rules o ixture predictions and heat transer siulations (FEM) Table 4 and Fig.9 shows the experiental results o the unidirectional coposite theral conductivity in the transverse direction and results o the predictions o the analysis as a unction o iber volue raction. In the case o transverse direction, the increase in the theral conductivity due to iber volue raction is very sall. The rules o ixture gave lower prediction than the experiental results but the FEM predictions show a air agreeent with the experiental results. Table 4 Results o theral conductivity in transverse direction easureents and analytical predictions Volue raction (%) 55 4 Experient (W/K) 1.17.73.36 Rules o ixture prediction (W/K).58.43.3 Steady-state siulation (W/K).81.59.38 Transient siulation (W/K).8.59.38 Fig.9 Theral conductivity o CFRP in transverse direction as a unction o iber volue raction or CFRP using YS9A.Coparison with the rules o ixture predictions and heat transer siulation (FEM) The anisotropy o the theral conductivity in pitch-based carbon iber reinorced epoxy unidirectional coposites is conired. It is understood that the laser lash ethod is eective to the easureent o the theral conductivity in the direction o the iber. It is also understood the theral conductivity can be predicted by the rules o ixtures, and two-diensional transient heat transer siulation. 4. Experiental and analytical results o theral resistance Fig.1 shows the experiental results o teperatures o T 1. And Fig.11 shows the experiental results o the theral resistance under natural convection using each heat sink and results o the predictions o steady-state thero-luid analysis. As or the CFRP heat sink, the heat radiation perorance is higher than the aluinu heat sink without aluite treatent, and it can be understood that it is the sae heat radiation perorance as the aluinu heat sink with aluite treatent under the natural convection. Moreover, experiental and analytical results show a good agreeent. 5

SHINJI OGIHARA, Makoto Yaaguchi, Takahito Chiba, Junichi Shiizu Yoji Okabe and Nobuo Takeda Teperature ( ) Theral resistance ( /W) 7 6 5 4 3 Al Al-aluite tretent CFRP 1 5 1 15 5 3 Tie (in.) 3 5 15 1 5 Fig.1 Teperature o T 1 Experient Analysis The heat radiation perorance under natural convection is copared with heat sinks ade o CFRP and aluinu alloy with dierent surace treatent, that is, with and without aluite treatent. To discuss the experiental results, threediensional steady-state heat transer siulations are conducted. As a result, the heat radiation perorance is higher than the aluinu alloy heat sink without Aluite treatent, and it is understood that heat sink ade o CFRP is the sae heat radiation perorance as the heat sink ade o aluinu alloy with aluite treatent under the natural convection. Moreover, experiental and analytical results show a good agreeent. Reerences [1] J.Zier, Extended Abstr Progra Bienn Con Carbon.,pp.39-391, 1991. [] Yokota H., Ohta H. and Shibata H. Coputer Siulation or Estiating the Eective Theral Conductivity o Coposite Materials. Japan Journal o Therophysical Properties, Vol.13, pp.4-45, 1999. Al Al-Aluite YS9A Fig.11 Experiental and analytical results o theral resistance 5 Conclusion Unidirectional coposites consisting o a pitchbased carbon iber and epoxy atrix is abricated and theral conductivity in the iber direction and transverse direction with various iber volue ractions is evaluated experientally by using the laser lash ethod. To discuss the experiental results, the rules o ixture prediction and twodiensional steady-state and transient heat transer siulations are conducted. As a result, the experient result o unidirectional coposites is alost corresponding to results o each analysis. Fro this, it is understood the theral conductivity can be predicted by the rules o ixture, and twodiensional transient heat transer siulation. Moreover, the anisotropy o the theral conductivity o unidirectional coposites can be conired. 6