Uncertainty on Fatigue Damage Accumulation for Composite Materials Toft, Henrik Stensgaard; Sørensen, John Dalsgaard

Similar documents
Pedersen, Ivar Chr. Bjerg; Hansen, Søren Mosegaard; Brincker, Rune; Aenlle, Manuel López

Risk-based Operation and Maintenance of Offshore Wind Turbines Sørensen, John Dalsgaard

Chapter 11: Simple Linear Regression and Correlation

Probabilistic Design of Wind Turbines

FUZZY FINITE ELEMENT METHOD

Composite Hypotheses testing

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

RELIABILITY ASSESSMENT

Simulation and Probability Distribution

Lecture Notes on Linear Regression

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Comparison of Regression Lines

Uncertainty as the Overlap of Alternate Conditional Distributions

DUE: WEDS FEB 21ST 2018

Negative Binomial Regression

8/25/17. Data Modeling. Data Modeling. Data Modeling. Patrice Koehl Department of Biological Sciences National University of Singapore

Linear Approximation with Regularization and Moving Least Squares

9.2 Seismic Loads Using ASCE Standard 7-93

Research Article Green s Theorem for Sign Data

Limited Dependent Variables

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

Fatigue Life Prediction Based on Variable Amplitude Tests

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

x i1 =1 for all i (the constant ).

Uncertainty and auto-correlation in. Measurement

MECHANICS OF MATERIALS

Optimal, risk-based operation and maintenance planning for offshore wind turbines

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

The Ordinary Least Squares (OLS) Estimator

Risk assessment of wind turbines close to highways

Second Order Analysis

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

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

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

Goodness of fit and Wilks theorem

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

The Study of Teaching-learning-based Optimization Algorithm

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

Laboratory 1c: Method of Least Squares

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

/ n ) are compared. The logic is: if the two

Laboratory 3: Method of Least Squares

The SAS program I used to obtain the analyses for my answers is given below.

SIMPLE LINEAR REGRESSION

Code_Aster. Identification of the model of Weibull

Credit Card Pricing and Impact of Adverse Selection

DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION

Indeterminate pin-jointed frames (trusses)

Chapter 9: Statistical Inference and the Relationship between Two Variables

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

On the Influential Points in the Functional Circular Relationship Models

Tilburg University. Note on Integer-Valued Bilinear Time Series Models Drost, Feico; van den Akker, Ramon; Werker, Bas. Publication date: 2007

Appendix B. The Finite Difference Scheme

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

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

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

SPANC -- SPlitpole ANalysis Code User Manual

Chapter 13: Multiple Regression

e i is a random error

Feb 14: Spatial analysis of data fields

Lecture 3 Stat102, Spring 2007

The Geometry of Logit and Probit

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH

Linear Regression Analysis: Terminology and Notation

An Application of Fuzzy Hypotheses Testing in Radar Detection

Introduction to the R Statistical Computing Environment R Programming

Resource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis

Interval Regression with Sample Selection

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Analytical Chemistry Calibration Curve Handout

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

STAT 405 BIOSTATISTICS (Fall 2016) Handout 15 Introduction to Logistic Regression

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

A Robust Method for Calculating the Correlation Coefficient

A new Approach for Solving Linear Ordinary Differential Equations

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Constitutive Modelling of Superplastic AA-5083

STAT 511 FINAL EXAM NAME Spring 2001

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

COMPONENT-LEVEL AND SYSTEM-LEVEL SENSITIVITY STUDY FOR EARTHQUAKE LOSS ESTIMATION

Statistical analysis using matlab. HY 439 Presented by: George Fortetsanakis

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

Uncertainty in measurements of power and energy on power networks

10) Activity analysis

Convexity preserving interpolation by splines of arbitrary degree

STAT 3008 Applied Regression Analysis

Generalized Linear Methods

Calibration and validation of model describing complete autotrophic nitrogen removal in granular sludge

Kernel Methods and SVMs Extension

Rockefeller College University at Albany

Transcription:

Aalborg Unverstet Uncertanty on Fatgue Damage Accumulaton for Composte Materals Toft, Henrk Stensgaard; Sørensen, John Dalsgaard Publshed n: Proceedngs of the Twenty Second Nordc Semnar on Computatonal Mechancs Publcaton date: 009 Document Verson Publsher's PDF, also known as Verson of record Lnk to publcaton from Aalborg Unversty Ctaton for publshed verson (APA): Toft, H. S., & Sørensen, J. D. (009). Uncertanty on Fatgue Damage Accumulaton for Composte Materals. In L. Damklde, L. Andersen, A. S. Krstensen, & E. Lund (Eds.), Proceedngs of the Twenty Second Nordc Semnar on Computatonal Mechancs (pp. 53-56). Department of Cvl Engneerng, Aalborg Unversty. DCE Techncal Memorandum, No. General rghts Copyrght and moral rghts for the publcatons made accessble n the publc portal are retaned by the authors and/or other copyrght owners and t s a condton of accessng publcatons that users recognse and abde by the legal requrements assocated wth these rghts.? Users may download and prnt one copy of any publcaton from the publc portal for the purpose of prvate study or research.? You may not further dstrbute the materal or use t for any proft-makng actvty or commercal gan? You may freely dstrbute the URL dentfyng the publcaton n the publc portal? Take down polcy If you beleve that ths document breaches copyrght please contact us at vbn@aub.aau.dk provdng detals, and we wll remove access to the work mmedately and nvestgate your clam. Downloaded from vbn.aau.dk on: jun 4, 08

Proceedngs of the Twenty Second Nordc Semnar on Computatonal Mechancs Aalborg Unversty 009 ISSN 90-778 DCE Techncal Memorandum No. Uncertanty on Fatgue Damage Accumulaton for Composte Materals Henrk Stensgaard Toft* and John Dalsgaard Sørensen Department of Cvl Engneerng Aalborg Unversty, Aalborg, Denmark emal: hst@cvl.aau.dk Summary In the present paper stochastc models for fatgue damage accumulaton for composte materals are presented based on publc avalable constant and varable ampltude fatgue tests. The methods used for estmatng the SN-curve and accumulated fatgue damage are presented. Introducton Damage accumulaton models for composte materals exposed to fatgue loadng have been wdely consdered n the lterature, see e.g. [] for a revew. However, even though new emprcal and physcal models for the accumulaton of damage are proposed, these models do not seem to perform much better than the lnear damage accumulaton proposed by Mner []. For ths reason the accumulated damage s normally determned by Mners rule as recommended n [3] for wnd turbne blades. The uncertantes n damage accumulaton based on Mners rule can n general be dvded nto three parts: Physcal uncertanty on the SN-curves Statstcal uncertanty on the SN-curves due to a lmted number of tests Model uncertanty on Mners rule The physcal uncertanty on the SN-curves s due to the natural nherent uncertanty n the materal and can not be reduced. The statstcal uncertanty can be reduced by performng addtonal fatgue tests and the model uncertanty can n prncple be reduced by adoptng a better model. In the present paper s the physcal and statstcal uncertanty on SN-curves for composte materal determned based on a number of constant ampltude fatgue tests performed wth dfferent mean stresses. Based on varable ampltude fatgue tests for the same materal usng a standard load spectrum and the estmated SN-curves based on constant ampltude fatgue tests s the model uncertanty on Mners rule determned. In most standards and regulatons ncludng [3] are fatgue desgn performed by usng a determnstc desgn approach. However, calculaton of the accumulated fatgue damage ncludes sgnfcant uncertantes for whch reason a probablstc desgn approach should be adopted n order to take the ndvdual uncertantes nto account n a ratonal manner. The stochastc models determned n the present paper forms the bass for a probablstc modelng of the consdered class of composte materals n fatgue loadng. Constant ampltude fatgue tests The constant ampltude and varable ampltude fatgue tests used n the present paper are gven n the OptDAT database [4] for geometry R04 MD (MultDrectonal lamnate). Ths geometry has been selected due to the many fatgue tests performed wth ths geometry. For composte materals the mean stress can have a sgnfcant nfluence on the fatgue propertes whch can be taken nto 53

account by calculaton of an SN-curve for dfferent R-ratos and arrangng these n a constant lfe dagram. The R-rato s defned by: mn R () max where mn and max are the mnmum and maxmum stress n a stress cycle respectvely. Dfferent types of SN-curves have been used for composte materal, but no specfc SN-curve have been recommended n [3]. In the present study s a log-log SN-curve used: log N log K mlog () where N s the number of cycles to falure, s the stress range and s parameter whch model the lack of ft and s assumed normal dstrbuted wth mean value zero and standard devaton. The constants K and m are materal parameters. By assumng that the resduals are normal dstrbuted on a log-log scale the lkelhood functon n case of n constant ampltude tests and n 0 run-outs s gven by: n 0 nn log N log K mlog log N log K mlog Llog K, exp (3) n where N and s the number of cycles to falure and stress range for test specmen number respectvely. The parameter m s determned by least square method and the parameters log K and are estmated usng the Maxmum-Lkelhood Method where the optmzaton problem max L log K, s solved usng a standard nonlnear optmzer, e.g. the NLPQL algorthm [5]. log K, In ths paper s m assumed fxed determned by the least square method, but ths parameter could also be ncluded n the optmzaton. Snce the parameters log K and are estmated by the Maxmum-Lkelhood technque they become asymptotcally Normal dstrbuted stochastc varables wth expected values equal to the Maxmum-Lkelhood estmators and covarance equal to, see e.g. [6]: log K log K, log K Clog K, H log K, (4) log K, log K where Hlog K, s the Hessan matrx wth second order dervatves of the log-lkelhood functon. log K and denote the standard devaton on log K and respectvely. log K, s the correlaton coeffcent between log K and. The Hessan matrx s estmated by numercal dfferentaton. In table are the estmated parameters shown and the SN-curves are ftted usng all vald constant ampltude fatgue tests for the partcular R-rato and runouts are taken nto account. The parameters log K and can be assumed uncorrelated. It s noted that represents the physcal uncertanty and that and represents the statstcal uncertanty. In table are the statc log K 54

tenson and compresson strength gven. In fgure (left) s the constant lfe dagram contanng the SN-curves and statc strengths shown. Table : SN-curves for dfferent R-ratos for geometry R04 MD. R-rato Tests n Runouts n 0 m log K logk 0.5 5 0 0.54 7.768 0.358 0.09 0.065 0. 45 9.508 7.9 0.59 0.039 0.07-0.4 8 0 7.58 3.398 0.435 0.08 0.058 -.0 84 3 6.79.359 0.878 0.094 0.068 -.5 0.983 35.3 0.633 0.97 0.43 0.0 34 0. 58.664 0.644 0.0 0.078.0 6 3 9.686 73.780 0.354 0.43 0.03 Table : Statc tenson strength (STT) and statc compresson strength (STC) for geometry R04 MD. Test-type Tests n Mean [MPa] Std. [MPa] STT 66 556.5 64. STC 55-458.6 33. amp [MPa] 300 00 N=0 3 N=0 4 N=0 5 N=0 6 N=0 7 N=0 8 N=0 9 R=.0 R=0.0 R=-.5 R=-.0 R=-0.4 R=0. R=0.5 00 0-600 -400-00 0 00 400 600 mean [MPa] Fgure. Left: Constant lfe dagram for geometry R04 MD. Rght: Lnear nterpolaton. Varable ampltude fatgue tests Varable ampltude fatgue tests are also performed wth geometry R04 MD. The load spectrum used s the wsper and wsperx spectra developed for representng the flap bendng moment of a wnd turbne blade. In order to calculate the accumulated damage D Mners rule for lnear damage accumulaton s used: n D (5) N In order to calculate the accumulated damage at falure the number of cycles to falure N has to be determned for an arbtrary R-rato. Ths s n the present paper done by lnear nterpolaton n the constant lfe dagram usng the procedure gven n the followng, see also fgure (rght). The stress cycle P s located n the constant lfe dagram Draw a lne a from the orgn through and beyond the pont P Identfy the constant lfe lnes closest to P, denoted n and n Calculate the length a on lne a between the two constant lfe lnes n and n 55

Calculate the length a on lne a between pont P and the constant lfe lne n Fnd the R-rato closest to P and calculate the length b between n and n ba Calculate b a Determne the stress ampltude CLD correspondng to pont Q Determne the expected number of cycles to falure N usng the SN-curve for the R-rato In table 3 are the conducted varable fatgue tests lsted together wth the mean and standard devaton on the accumulated damage at falure usng the materal parameters gven n table and. It s noted that COV D represents the model uncertanty. Table 3: Mean and standard devaton for estmated damage at falure for varable ampltude tests. Spectrum Tests N Mean Std. COV D Wsper 0 0.9000 0.5355 0.595 Wsperx 3 0.763 0.98 0.77 Reverse Wsper 0.00 - - Reverse Wsperx 0 0.379 0.576 0.496 Wsper, Wsperx 3 0.5474 0.4886 0.893 All 35 0.46 0.496 0.908 From table 3 t s seen that except for the wsper spectrum the estmated accumulated damage at falure s sgnfcantly below one. The uncertanty for fatgue damage accumulaton are often modelled by a lognormal dstrbuton n order two avod negatve values of Mners rule whch are physcal mpossble. The mean and standard devatons gven n table 3 can be used n a lognormal dstrbuton. Concludng remarks In the present paper are stochastc models for the uncertanty related to fatgue damage accumulaton for composte materals presented. The stochastc models are based on publc avalable constant and varable ampltude fatgue tests and the procedure used for estmatng the stochastc models are presented. Acknowledgement The work presented n ths paper s part of the project Probablstc desgn of wnd turbnes supported by the Dansh Research Agency, grant no. 04-05-0075 and the project Improvement of methods for fatgue assessment supported by the Dansh Energy Authorty, EFP007 grant no. 33033-0077. The fnancal support s greatly apprecated. References [] Post NL, Case SW, Lesko JJ. Modelng the varable ampltude fatgue of composte materals: A revew and evaluaton of the state of the art for spectrum loadng. Internatonal Journal of Fatgue 008 Dec;30():064-86. [] Mner MA. Cumulatve Damage n Fatgue. Journal of Appled Mechancs-Transactons of the Asme 945;(3):A59-A64. [3] IEC 6400-. Wnd turbnes - Part: Desgn requrements. 3 rd edton. 005. [4] OPTIMAT BLADES. 006. www.kc-wmc.nl. [5] Schttkowsk K. NLPQL: A FORTRAN Subroutne Solvng Non-Lnear Programmng Problems. Annals of Operatons Research 986. [6] Lndley DV. Introducton to Probablty and Statstcs from a Baysan Vewpont. Cambrdge Unversty Press, Cambrdge; 976. 56