Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

Similar documents
Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Finite Element Modelling of truss/cable structures

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Constitutive Modelling of Superplastic AA-5083

Please initial the statement below to show that you have read it

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

Effect of loading frequency on the settlement of granular layer

Simulation analysis for the ultimate behavior of full-scale lead-rubber seismic isolation bearings

In this section is given an overview of the common elasticity models.

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Effects of Polymer Concentration and Molecular Weight on the Dynamics of Visco-Elasto- Capillary Breakup

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Nordic Insulation Symposium - Nord-IS 13 - Trondheim, Norway, June 9-12, 2013

CHAPTER 9 CONCLUSIONS

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

Second Order Analysis

Analytical and Numerical Analysis of Free Bulge Tube Hydroforming

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

SIMPLIFIED PREDICTION METHOD FOR SEISMIC RESPONSE OF ROCKING STRUCTURAL SYSTEMS WITH YIELDING BASE PLATES

THE IDENTIFICATION OF MATERIAL PARAMETERS IN NONLINEAR DEFORMATION MODELS OF METALLIC-PLASTIC CYLINDRICAL SHELLS UNDER PULSED LOADING

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Determinationn of gasket load drop at welding neck flange joints operating at medium high temperature

2. PROBLEM STATEMENT AND SOLUTION STRATEGIES. L q. Suppose that we have a structure with known geometry (b, h, and L) and material properties (EA).

SIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

Lecture Note 3. Eshelby s Inclusion II

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

LAYING METHOD OF GEOGRID MATERIAL AND REINFORCEMENT EFFECT OF SOFT GROUND

The Finite Element Method

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

DETERMINATION OF COMPOSITE EFFECTIVE CHARACTERISTICS AND COMPLEX MODULES FOR SHELL STRUCTURES FROM THE SOLUTION OF INVERSE PROBLEM

Dynamic Contact Problem for Slide Hinge

Modal Identification of the Elastic Properties in Composite Sandwich Structures

Numerical Nonlinear Analysis with the Boundary Element Method

DEFORMATION ANALYSIS OF AN INFLATED CYLINDRICAL MEMBRANE OF COMPOSITE WITH RUBBER MATRIX REINFORCED BY TEXTILE CORDS

Modeling of Dynamic Systems

PRATICAL STATIC CALCULATION METHOD FOR ESTIMATING ELASTO-PLASTIC DYNAMIC RESPONSES OF SPACE FRAMES

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

SEISMIC ANALYSIS AND DESIGN OF BUILDING STRUCTURES WITH SUPPLEMENTAL LEAD DAMPERS

OPTIMAL DESIGN OF VISCOELASTIC COMPOSITES WITH PERIODIC MICROSTRUCTURES

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model

829. An adaptive method for inertia force identification in cantilever under moving mass

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

An influence line shows how the force in a particular member changes as a concentrated load is moved along the structure.

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

Note 10. Modeling and Simulation of Dynamic Systems

Influential Factors Affecting Inherent Deformation during Plate Forming by Line Heating (Report 1)

Lecture 8 Modal Analysis

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

Parameter Estimation for Dynamic System using Unscented Kalman filter

The Analysis of Convection Experiment

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT

Structural Dynamics and Earthquake Engineering

MECHANICS OF MATERIALS

Effects of internal=external pressure on the global buckling of pipelines

THEORETICAL AND EXPERIMENTAL ANALISYS OF ELASTOMERS

Modal Identification of Non-Linear Structures and the Use of Modal Model in Structural Dynamic Analysis

STUDY ON SEISMIC BEHAVIOR OF RC COMPOSITE CORE WALLS WITH CONCEALED STEEL TRUSS SUBJECTED TO COMBINED ACTION

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Finite Element Solution Algorithm for Nonlinear Elasticity Problems by Domain Decomposition Method

Numerical Simulation of ph-sensitive Hydrogel Response in Different Conditions

Determination of the response distributions of cantilever beam under sinusoidal base excitation

Indeterminate pin-jointed frames (trusses)

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

Topology Optimization of Elastic Material Microstructures with Classic Models of Micromechanics

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

Numerical Transient Heat Conduction Experiment

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

PREDICTIVE CONTROL BY DISTRIBUTED PARAMETER SYSTEMS BLOCKSET FOR MATLAB & SIMULINK

Fast Simulation of Pyroshock Responses of a Conical Structure Using Rotation-Superposition Method

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Peridynamic Modeling of plain concrete structures under monotonic loading Jiezhi Lu1, a, Yaoting Zhang1, b, Zhijun Chen1

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

ES 240 Solid Mechanics Z. Suo. Vibration

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

LAB 4: Modulus of elasticity

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

USING THE INCREMENTAL HOLE-DRILLING TECHNIQUE FOR MEASURING RESIDUAL STRESSES IN FIBRE-REINFORCED POLYMER COMPOSITES

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

DYNAMIC BEHAVIOR OF PILE GROUP CONSIDERING SOIL-PILE-CAP INTERACTION

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

An Interactive Optimisation Tool for Allocation Problems

title: SEM Charging Effect Model for Chromium/Quartz Photolithography Masks authors: Adam Seeger, Alessandro Duci, Horst Haussecker

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Numerical Modeling of Woven Carbon Composite Failure

ANALYSIS OF CONTACT PROBLEM USING IMPROVED FAST MULTIPOLE BEM WITH VARIABLE ELEMENTS LENGTH THEORY

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Field computation with finite element method applied for diagnosis eccentricity fault in induction machine

Transcription:

ECCMR, Prague, Czech Republc; September 3 th, 2015 Lfetme predcton of EP and NBR rubber seal by thermos-vscoelastc model Kotaro KOBAYASHI, Takahro ISOZAKI, Akhro MATSUDA Unversty of Tsukuba, Japan Yoshnobu MIZUTANI, Yasuhko HORI Central Research Insttute of Electrcal Power Industry 1

Contents Background Objectves Proposed materal model Baxal loadng test Stress relaxaton test Numercal smulaton Concluson 2

Background Gas Insulated Swtchgear (GIS) s used n electrcal substatons Insulaton gas was contanment nsde of GIS Rubber O-rng s appled to GIS for sealng materal GIS has engneerng and envronmental problems caused by degradaton of seal materal Leakage accdent of nsulaton gas Insulaton gas has strong greenhouse effect Emsson of nsulated gas s regulated nternatonally Replacement of rubber O-rng s requred dsmantlng of GIS 3

Rubber O-rng used n GIS The O-rng s made of Ethylene Propylene Rubber (EPR) or Ntrle Butadene Rubber (NBR) Deteroraton of rubber O-rng s caused manly by heat of GIS Rubber O-rng s placed on ld groove n flange nsde GIS gas Restorng force Rubber O-rng Sealng mechansm of rubber O-rng 4

Seal performance of rubber O-rng Seal performance of rubber O-rng s evaluated by the compresson set compresson set s permanent stran remanng after compresson Threshold value for seal performance s 80% Cross secton of O-rng CS D D 0 1 0 D D 2 100(%) CS : compresson set Before use After use Compresson set occurs n after cross secton D 0 D 1 D 2 : thckness before compresson : thckness after compresson : heght of groove 5

Objectves Propose a thermos-vscoelastc model to predct sealng performance of EPR and NBR O-rng Introduce contact condton between O-rng and ld groove to boundary condton of FEM smulaton Evaluate applcablty of the proposed model We conduct 1, Baxal loadng tests and stress relaxaton tests 2, Numercal smulaton based on the model 6

Proposed materal model Proposed materal model consst of parallel connectons of the Maxwell model Proposed materal model show stress relaxaton of rubber Sprng elements of the Maxwell model s hyperelastcty model Dashpot elements of the Maxwell model s vscosty model Dashpot elements of the Maxwell model show temperature dependence Sprng [hyperelastcty] Dashpot [vscosty, temperature dependence] 7

8 Hyperelastcty model 3 3 3 I 5 4 3 2 1 2 I C exp C ) I ( C ) I ( C C W Relatonshp between stress and stran s nonlnear n hyperelastcty model Stress s calculated from partal dfferental of elastc potental functon wth stran C S W 2 Materal parameters are dentfed by baxal loadng test S : Second Pola-Krchhoff stress W : Elastc potental functon C : Rght Cauchy-Green deformaton tensor I : 1 st and 2 nd nvarant of volume-preserved C

Vscosty model Stress relaxaton behavor s determned by vscosty model Stress relaxaton s represented from superposton of deferent dashpot Stress relaxaton s represented by S S S t S0 N d t 0 Nd t exp 1 : Stress at tme=t : Stress at tme=0 : Number of dashpot t : Stffness rato : Tme : Relaxaton tme : temperature Nd t exp 1 Vscous characterstc are dentfed by stress relaxaton test S S t 0 9

Identfcaton of parameters Baxal loadng tests were conducted for EPR and NBR Stress relaxaton tests were conducted for EPR and NBR We wll talk about the baxal loadng test and stress relaxaton test of NBR n ths presentaton 10

Baxal loadng test Rubber sheet specmen (80mm*80mm*1mm) One sde was loaded to 200% stretch Other sde was kept ntal length Two stress σ 1 and σ 2 were approxmated by hyperelastcty model Fx 5 4 σ1(test) σ2(test) σ1(smulaton) σ2(smulaton) Tensle Tensle Stress [MPa] 3 2 1 Fx 0 1 1.2 1.4 1.6 1.8 2 Stretch 11

Stress relaxaton test Rubber O-rng specmen (cross sectonal dameter ; 5.7mm, nternal dameter ; 94.6mm) O-rng was compressed at compresson rate of 27.7% Test temperature s 100C o, 110C o and 120C o Restorng force was measured by load cell durng the test O-rng Thermostat chamber Orgnal jg Load cell 12

Nomnalze stress S 0 /S t 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 Stress relaxaton curve Stress relaxaton curve at 120 C o was selected as master curve Master curve was approxmated by 100 110 120 0.2 0 2 10 6 4 10 6 6 10 6 8 10 6 1 10 7 S S t 0 S t S 0 N d t Nd t exp 1 : Stress at tme=t : Stress at tme=0 : Number of dashpot : Stffness rato : Tme : Relaxaton tme : temperature Degradaton perod [s] 13

Temperature dependence Master curve was multpled by acceleraton factor to ft other curves Acceleraton factor were arranged usng the Arrhenus plot Applcablty of the Arrhenus plot to relaxaton tme was confrmed Logarthms of acceleraton foctor 0.2 0 Acceleraton foctor -0.2-0.4-0.6-0.8-1 -1.2 0.00254 0.00256 0.00258 0.0026 0.00262 0.00264 0.00266 0.00268 Recprocal of absolute tempreture [1/K] Relaxaton tme at any temperature can be expressed as follows /( A* exp( B / )) A,B 0 0 : Relaxaton tme at temperature : Relaxaton tme at temperature of master curve : A constant obtaned from lnear approxmaton of Arrhenus plot 14

Fnte element method We developed an orgnal code wth proposed model 3-dmentonal sold element Dsplacement / pressure mxed method Newton-Raphson method for teraton Coupled analyss of thermal conductvty and mechancal property The half of cross secton of O-rng 15

Boundary condton Contact wth groove when O-rng s compressed f f y y 0 Load s added by penalty method y y 0 : Load : Penalty number : Pre-deformaton coordnate of node outsde deformable range : Post-deformaton coordnate of node outsde deformable range Groove Load Groove O-rng can not be deformed outsde grooves 16

Procedure of smulaton Smulaton conssts 4 step Step 1; Add load n vertcal drecton Step 2; Add load n horzontal drecton Step 3; Stress relaxaton analyss was conducted Step 4; Cancel the constrants of node Contact surface Load correspondng to each node Constraned node Contact surface Relaxaton of nternal stress Groove surface 17

Compresson set [%] 100 90 Numercal smulaton of NBR Compresson set of smulatons were compared wth and long compresson deteroraton set perod of deterorated rubber O-rng Good agreement wth expermental data n hgh temperature Change of stress relaxaton behavor wth temperature decreasng s suggested 80 70 60 50 40 30 20 10 0 0 1000 2000 3000 4000 5000 Tme [hr] 70 (smulaton) Lfetme predcton n short 70 (Experment) degradaton perod s 90 (smulaton) requre further nvestgaton 90 (Experment) Acquston of tress 110 (smulaton) relaxaton behavor n more 110 (Experment) low temperature would mprove predcton accuracy 18

Concluson Thermos-vscoelastc model was proposed to predct seal performance of EPR and NBR O-rng The proposed model could reproduce nonlnear vscoelastc of rubber dependng on temperature Boundary condton consderng contact between O- rng and ld groove s ntroduced to FEM smulaton Analyss on actual use condton has become possble Applcablty of proposed model was confrmed by comparson wth deterorated rubber O-rng Effectveness of lfetme predcton of NBR O-rng was shown n hgh temperature and long term degradaton perod. Stress relaxaton behavor n more low temperature would be requred to mprove accuracy 19

THANK YOU FOR YOUR ATTENTION 20