Nonlinear correction to the bending stiffness of a damaged composite beam
|
|
- Sherman Horn
- 6 years ago
- Views:
Transcription
1 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Nonlinear correction to te bending stiffness of a damaged composite beam W. Van Paepegem a*, R. Decaene a and J. Degrieck a a Professor, Dept. of Mecanical Construction and Production, Sint-Pietersnieuwstraat 41, 9 Gent, Belgium Abstract Wen studying damage in composite materials, te classic beam teory is sometimes used in a modified manner to calculate te bending response of a damaged composite laminate. Te longitudinal stiffness E is ten replaced by a field variable E(x,y) = E [1-D(x,y)]. Te damage distribution D(x,y) affects te calculation of stresses and strains, and requires a modified calculation of te neutral fibre y (x) and te bending stiffness EI(x). It seems obvious to postulate tat te bending stiffness EI(x) is also multiplied wit [1-D(x,y)], as was proposed by, amongst oters, Sidoroff and Subagio (Sidoroff, F. and Subagio, B. (1987). Fatigue damage modelling of composite materials from bending tests. In : Proceedings of ICCM-VI & ECCM-II, pp ). It is sown ere tat tis formulation is not correct and tat a nonlinear correction term to te bending stiffness degradation is necessary. Tis term is especially important wen te damage D(x,y) is reacing iger values. Te results of a semi-analytical simulation for cantilever beams is sown to support te findings. Keywords: beam teory, damage, composite, bending, residual stiffness 1 Introduction Today, a lot of researc is dedicated to te understanding of damage in fibre-reinforced composites. Simple and more advanced simulation tools are used to model te mecanical response of damaged composites. One of te straigtforward calculation tools is te classic beam teory. In 1987, Sidoroff and Subagio [1] proposed a local damage model to predict te damage growt in composites subject to bending fatigue tests. Tis local damage model was combined wit an adapted formulation of te classic beam teory to calculate te stresses and strains in te damaged composite beam under fatigue loading. In tis paper it is demonstrated tat te linear relation between te degradation of te bending stiffness and te damage is not correct and a nonlinear correction term sould be added. Te relevant equations are derived and te formulations are implemented in a semi-analytical simulation of a bending fatigue test. Review of classic beam teory for damaged beams In tis paragrap, te existing approac for damage modelling wit te classic beam teory is discussed. Te work of Sidoroff and Subagio is taken as a guideline as tey worked out te equations and applied te teory to te bending fatigue of a composite beam. * Autor to wom correspondence sould be addressed ( Wim.VanPaepegem@UGent.be ).
2 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Te local damage model, proposed by Sidoroff and Subagio [1], was a one-dimensional damage model based on te Continuum Damage Mecanics teory [-8]. Te well-known relation between stress and strain was: σ = E (1 D) ε (1) wit E te virgin elastic modulus and D te scalar damage ( D 1). Tis relation evolves from te concept of strain equivalence, introduced by Lemaitre [9], wic states tat a damaged volume of material under te nominal stress σ sows te same strain response as a comparable undamaged volume under te effective stress σ ~ (= σ/(1-d)). Te damage growt under fatigue loading was establised as follows [1]: dd dn ( Δε) c A in tension b = (1 D) () in compression were : - D : local damage variable - N : number of cycles - Δ ε : strain amplitude of te applied cyclic loading - A, b and c : tree material constants Te relevant equations for te bending response of te damaged composite beam were derived from classic beam teory as follows: wen te axial force is supposed to remain zero and wen only a bending moment exists, te position of te neutral fibre y (x) at eac moment of time is calculated as: + [ 1 D(x, y) ] + 1 y dy y (x) = (3) D(x, y) dy were : - x : coordinate along te beam lengt - y : tickness-coordinate, wit y = in te middle of te specimen s tickness - : total tickness of te specimen Te degraded bending stiffness EI(x) becomes (wit b te specimen s widt): + [ 1 D(x, y) ] EI (x) = b E y dy (4) In te next paragrap it will be sown tat tis latter equation is not correct and sould include a nonlinear correction term for te bending stiffness EI(x).
3 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. 3 Derivation of te nonlinear correction term for te bending stiffness EI(x) For te derivation of te nonlinear correction term for te bending stiffness EI(x), it is again assumed tat: σ ε (5) = E (1 D) were D is a field variable D(x,y). Te x-coordinate is again along te beam lengt, wile te y-coordinate represents te beam tickness. Te following displacement field is proposed for te bending beam: u u y x (x, y) = u(x) (x, y) = u (x) + u' y (6) were u(x) is te vertical displacement of te bending beam, u (x) is te longitudinal displacement of te midplane and y is te tickness coordinate lying in te interval [-/, +/]. Te widt b is assumed to be unity for te moment. Te strain ε can ten be calculated as: And te corresponding stress σ equals: u x du ε = = + u'' y = ε + u'' y (7) x dx σ = E ε = E (1 D) ( ε + u'' y) (8) Te normal force N is assumed to be zero in pure bending: N = σ dy = E ε (1 D) dy + E u' ' (1 D) y dy = (9) From Equation (9), te strain ε can be calculated as: ε (1 D) y dy = u '' (1) (1 D) dy Ten te position of te neutral fibre y is determined from te equations (7) and (1):
4 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. ε (1 D) y dy = u' ' y = u' '(y y) (1 D) dy (11) Te normal stress σ can ten be written as: σ = (1 D) u'' (y y ) (1) E Te bending moment M(x) is calculated from te expression: M (x) = σ y dy (13) Taking into account Equation (1), te bending moment M(x) can be written as: M (x) = E u' ' (1 D) y dy y (1 D) y dy (14) Te bending stiffness EI(x) is ten defined from te equation: Hence: M(x) = EI(x) u'' (15) (1 D) y dy = EI (x) E (1 D) y dy (16) (1 D) dy If Equation (16) is compared wit Equation (4) proposed by Sidoroff and Subagio [1], it is clear tat a nonlinear correction term appears in te expression for te bending stiffness EI(x). From Equation (11) and (15), it follows directly tat:
5 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. M(x) (y y) ε = (17) EI(x) 4 Simulations To sow te importance of te nonlinear correction term in Equation (16), a semi-analytical simulation is performed for cantilever bending fatigue tests on a composite specimen. In Figure 1, te conventions for te simulations are sown. Te composite specimen is modelled as a cantilever beam, and te moving clamp (rigt) is represented by a rigid rod. Te prescribed displacement in te displacement-controlled bending fatigue test is u max = u(c). L a y A x B u(b) u(c) C α(b) = α(c) M [N.mm] F x [mm] - F.(L + a) Figure 1 Bending of te cantilever composite beam following te conventions of te classic beam teory. Te force necessary to impose te prescribed displacement u max = u(c), can be calculated from te following equations: u' ' M(x) L + a x = = F EI(x) EI(x) (18) x L + a x = F dx EI(x) (19) u' (x) u' (B) L L + a x = F dx () EI(x)
6 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. L L + a x' u (B) = u' (x)dx = F dx dx' (1) EI(x' ) L x u (C) = u'(b) a + u(b) () L L x L + a x L + a x' u = = max u(c) F a dx + dx dx' (3) EI(x) EI(x' ) Te equation (3) provides te relation between te force F and te prescribed displacement u max. Te local damage model is very similar to te one used by Sidoroff and Subagio and as te following form: dd(x, y) dn Δσ(x,y) A σts = in tension (4) ( 1 D(x, y) ) b c in compression Here, te stress amplitude Δσ(x,y) in te damage law (4) corresponds wit te maximum stress in eac material point during te simulated fatigue loading cycle. Te material parameters E and σ TS of te (unidirectional) composite laminate were respectively 4.57 GPa and 39.7 MPa [1]. Te constants A, b and c were respectively [1/cycle],.45 [-] and 6.5 [-] [1]. Unlike te local damage model by Sidoroff and Subagio [1], tis model includes an elastic limit of te stress and can be more easily related wit S-N curves by use of te stress amplitude instead of te strain amplitude. Of course, if one and te same damage model is applied for bot te linear and nonlinear bending stiffness degradation, it does not matter really wic damage model is cosen to demonstrate te effect of te nonlinear correction term. Tis constitutive model can be easily implemented in a matematical software package suc as Matcad. Te numerical integration formulae must be cosen suc tat te second degree polynomials are exactly integrated. Tis is te case for te Simpson s rule, wic is a Newton-Cotes quadrature formula. Because te increase of te damage variable D during one cycle is so small, te integrations must be exact indeed, oterwise te relative error on te calculation of te bending stiffness EI may be larger tan te increase of te damage variable itself. Terefore te conventional first-order trapezium metod is not suited for tis purpose. First te distribution of te bending moment along te lengt of te specimen is determined. Secondly te stresses and strains in eac integration point can be calculated. Te damage law is applied and a new cycle is evaluated. From equation (3) te necessary force to impose te displacement wit amplitude u max can be calculated for eac cycle. Tis algoritm must be solved for eac calculated loading cycle, and as a consequence, te Matcad TM workseet must be called several times. Tis is only possible witin te Matconnex TM environment wic allows to repetitively call te same Matcad TM workseet wit different starting values.
7 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. 5 Simulations In te following semi-analytical simulation, te displacement-controlled bending fatigue test, sown in Figure 1, is simulated for about 3, loading cycles. Te displacement amplitude u max was cosen to be 3.4 mm. Te simulations wit linear and nonlinear bending stiffness degradation EI(x) were done wit te same material parameters E and σ TS and te same constants A, b and c for te damage law (4). Te only difference is due to te nonlinear correction term EI(x). Figure sows te simulated force-cycle istory wit linear and nonlinear bending stiffness degradation EI(x). Te initial force F = N is te force amplitude necessary to bend te specimen into its maximum deformed state (u = u max ) during te first loading cycle. Ten, due to te damage growt, te force amplitude is decreasing during fatigue life. It is clear tat te simulated degradation of te composite specimen is muc larger in te case of te nonlinear correction term, altoug te applied fatigue damage law (4) and its constants are te same in bot cases Force versus number of cycles for composite specimen u max = 3.4 mm 9 8 Force [N] Linear EI Nonlinear EI No. of cycles [-] Figure Simulated force-cycle istory wit linear and nonlinear bending stiffness degradation EI(x). Of course, te difference in force-cycle istory sould be reflected in a difference in damage distribution as well. Tis is confirmed by Figure 3 and Figure 4 wic sow te simulated damage distribution in te clamped cross-section (x = ) wit linear and nonlinear bending stiffness degradation respectively.
8 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Damage distribution in te clamped cross-section wit linear EI 1..6 eigt y [mm] damage [-] cycle 38 cycle,76 cycle 388,3 Figure 3 Simulated damage distribution in te clamped cross-section at several stages during fatigue life wit linear EI. Damage distribution in te clamped cross-section wit nonlinear EI 1..6 eigt y [mm] damage [-] cycle 95 cycle,7 cycle 36,35 Figure 4 Simulated damage distribution in te clamped cross-section at several stages during fatigue life wit nonlinear EI. Figure 5 and Figure 6 sow te corresponding simulated strain distribution in te clamped cross-section wit linear and nonlinear bending stiffness degradation. For te simulation witout te correction term, te strain increases at te tensile side of te specimen, but decreases at te compressive side. Tis is not realistic as due to te increasing damage at te tensile side, stress redistribution occurs and te strain sould increase at te compressive side. Tis is confirmed by te simulation in Figure 6 wit te nonlinear correction term.
9 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Strain distribution in te clamped cross-section wit linear EI cycle 38 cycle,76 cycle 388,3 eigt y [mm] strain [-] Figure 5 Simulated strain distribution in te clamped cross-section at several stages during fatigue life wit linear EI. Strain distribution in te clamped cross-section wit nonlinear EI eigt y [mm] strain [-] cycle 95 cycle,7 cycle 36,35 Figure 6 Simulated strain distribution in te clamped cross-section at several stages during fatigue life wit nonlinear EI. Finally, Figure 7 and Figure 8 sow te corresponding simulated stress distribution in te clamped cross-section at several stages during fatigue life, wit linear and nonlinear bending stiffness degradation. Again te linear simulation in Figure 7 predicts a decrease of te compressive stress wit increasing damage at te tensile side and tat is not realistic.
10 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Stress distribution in te clamped cross-section wit linear EI cycle 38 cycle,76 cycle 388,3 eigt y [mm] stress [MPa] Figure 7 Simulated stress distribution in te clamped cross-section at several stages during fatigue life wit linear EI. Stress distribution in te clamped cross-section wit nonlinear EI cycle 95 cycle,7 cycle 36,35 eigt y [mm] stress [MPa] Figure 8 Simulated stress distribution in te clamped cross-section at several stages during fatigue life wit nonlinear EI. It could be argued tat te final damage values are very ig and would not be tolerated in fatigue design for composites. Tat is true, but tese simulations serve te only purpose of clearly demonstrating te effect of te nonlinear correction term. Of course te larger te damage values, te larger te effect of te correction term. 6 Conclusions
11 Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Te classic beam teory can be used in a modified manner to calculate te damage growt in fibre-reinforced composites during fatigue. It is reasonable to assume tat if te longitudinal stiffness decreases in a linear relation wit te damage, te bending stiffness will do as well. However in tis paper it is demonstrated tat a nonlinear correction term for te bending stiffness degradation is necessary to predict reliable results. Semi-analytical simulations of a displacement-controlled bending fatigue test are sown to support te validity of te correction. It is evident tat more sopisticated calculation tools like finite elements are widely used to simulate damage growt in composite materials, but if one wants to use a simplified tool like te modified classic beam teory, it is important to be aware of te necessity of te nonlinear correction term for te degradation of te bending stiffness. Acknowledgements Te autor W. Van Paepegem gratefully acknowledges is finance troug a grant of te Fund for Scientific Researc Flanders (F.W.O.). References [1] Sidoroff, F. and Subagio, B. (1987). Fatigue damage modelling of composite materials from bending tests. In : Mattews, F.L., Buskell, N.C.R., Hodgkinson, J.M. and Morton, J. (eds.). Sixt International Conference on Composite Materials (ICCM-VI) & Second European Conference on Composite Materials (ECCM-II) : Volume 4. Proceedings, -4 July 1987, London, UK, Elsevier, pp [] Kacanov, L.M. (1958). On creep rupture time. Izv. Acad. Nauk SSSR, Otd. Tecn. Nauk, No.8, [3] Kacanov, L.M. (1986). Introduction to continuum damage mecanics. Dordrect, Martinus Nijoff Publisers, 135 pp. [4] Krajcinovic, D. and Lemaitre, J. (eds.) (1987). Continuum damage mecanics: teory and applications. Wien, Springer - Verlag, 94 pp. [5] Caboce, J.L. (1988). Continuum damage mecanics : part I - General concepts. Journal of Applied Mecanics, 55, [6] Caboce, J.L. (1988). Continuum damage mecanics : part II - Damage growt, crack initiation and crack growt. Journal of Applied Mecanics, 55, [7] Krajcinovic, D. (1985). Continuous damage mecanics revisited : basic concepts and definitions. Journal of Applied Mecanics, 5, [8] Sidoroff, F. (1984). Damage mecanics and its application to composite materials. In : Cardon, A.H. and Vercery, G. (eds.). Mecanical caracterisation of load bearing fibre composite laminates. Proceedings of te European Mecanics Colloquium 18, 9-31 August 1984, Brussels, Belgium, Elsevier, pp [9] Lemaitre, J. (1971). Evaluation of dissipation and damage in metals, submitted to dynamic loading. Proceedings I.C.M. I, Kyoto, Japan. [1] Van Paepegem, W. (). Development and finite element implementation of a damage model for fatigue of fibre-reinforced polymers. P.D. tesis. Gent, Belgium, Gent University Arcitectural and Engineering Press (ISBN ), 43 p.
Calculation of Damage-dependent Directional Failure Indices from the Tsai-Wu Static Failure Criterion
Van Paepegem, W. and Degrieck, J. (3. alculation of Damage-dependent Directional Failure Indices from the sai-wu Static Failure riterion. omposites Science and echnology, 63(, 35-3. alculation of Damage-dependent
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationExperimental setup for and numerical modelling of bending fatigue experiments on plain woven glass/epoxy composites
Experimental setup for and numerical modelling of bending fatigue experiments on plain woven glass/epoxy composites W. Van Paepegem a* and J. Degrieck b a Research Assistant of the Fund for Scientific
More informationModelling 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 information3. Using your answers to the two previous questions, evaluate the Mratio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationOptimization of the thin-walled rod with an open profile
(1) DOI: 1.151/ matecconf/181 IPICSE-1 Optimization of te tin-walled rod wit an open profile Vladimir Andreev 1,* Elena Barmenkova 1, 1 Moccow State University of Civil Engineering, Yaroslavskoye s., Moscow
More informationCALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY
CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations
More informationFinite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading
nternational Journal of ntegrated Engineering, Vol. 9 No. 2 (207) p. -8 Finite Element Analysis of J-ntegral for Surface Cracks in Round Bars under Combined Mode Loading A.E smail, A.K Ariffin 2, S. Abdulla
More informationStability of Smart Beams with Varying Properties Based on the First Order Shear Deformation Theory Located on a Continuous Elastic Foundation
Australian Journal of Basic and Applied Sciences, 5(7): 743-747, ISSN 99-878 Stability of Smart Beams wit Varying Properties Based on te First Order Sear Deformation Teory ocated on a Continuous Elastic
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationPath to static failure of machine components
Pat to static failure of macine components Load Stress Discussed last week (w) Ductile material Yield Strain Brittle material Fracture Fracture Dr. P. Buyung Kosasi,Spring 008 Name some of ductile and
More informationA Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment
American Journal of Mecanical and Industrial Engineering 28; 3(4): 47-54 ttp://www.sciencepublisinggroup.com/j/ajmie doi:.648/j.ajmie.2834.2 ISSN: 2575-679 (Print); ISSN: 2575-66 (Online) A Multiaxial
More informationAnalysis of Stress and Deflection about Steel-Concrete Composite Girders Considering Slippage and Shrink & Creep Under Bending
Send Orders for Reprints to reprints@bentamscience.ae Te Open Civil Engineering Journal 9 7-7 7 Open Access Analysis of Stress and Deflection about Steel-Concrete Composite Girders Considering Slippage
More informationEuler-Bernoulli Beam Theory in the Presence of Fiber Bending Stiffness
IOSR Journal of Matematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 13, Issue 3 Ver. V (May - June 017), PP 10-17 www.iosrjournals.org Euler-Bernoulli Beam Teory in te Presence of Fiber Bending
More informationTheoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump
TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing
More informationNCCI: Simple methods for second order effects in portal frames
NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames Tis NCC presents information
More informationDEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING
Asia-Pacific Conference on FRP in Structures (APFIS 2007) S.T. Smit (ed) 2007 International Institute for FRP in Construction DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT:
More informationFlapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method
Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationFabric Evolution and Its Effect on Strain Localization in Sand
Fabric Evolution and Its Effect on Strain Localization in Sand Ziwei Gao and Jidong Zao Abstract Fabric anisotropy affects importantly te overall beaviour of sand including its strengt and deformation
More informationSample Problems for Exam II
Sample Problems for Exam 1. Te saft below as lengt L, Torsional stiffness GJ and torque T is applied at point C, wic is at a distance of 0.6L from te left (point ). Use Castigliano teorem to Calculate
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationStatic Response Analysis of a FGM Timoshenko s Beam Subjected to Uniformly Distributed Loading Condition
Static Response Analysis of a FGM Timoseno s Beam Subjected to Uniformly Distributed Loading Condition 8 Aas Roy Department of Mecanical Engineering National Institute of Tecnology Durgapur Durgapur, Maatma
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationEmpirical models for estimating liquefaction-induced lateral spread displacement
Empirical models for estimating liquefaction-induced lateral spread displacement J.J. Zang and J.X. Zao Institute of Geological & Nuclear Sciences Ltd, Lower Hutt, New Zealand. 2004 NZSEE Conference ABSTRACT:
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationComment on Experimental observations of saltwater up-coning
1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:
More informationPoisson s ratio as a sensitive indicator of (fatigue) damage in fibre-reinforced plastics
10.1111/j.1460-2695.2007.01095.x Poisson s ratio as a sensitive indicator of (fatigue) damage in fibre-reinforced plastics W. VAN PAEPEGEM, I. DE BAERE, E. LAMKANFI and J. DEGRIECK Ghent University, Department
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationLehrstuhl für Technische Mechanik
Ä Preprint 02-2007 Rur University Bocum Lerstul für Tecnisce Mecanik A novel variational algoritmic formulation for wrinkling at finite strains based on energy minimization: application to mes adaption
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationModel development for the beveling of quartz crystal blanks
9t International Congress on Modelling and Simulation, Pert, Australia, 6 December 0 ttp://mssanz.org.au/modsim0 Model development for te beveling of quartz crystal blanks C. Dong a a Department of Mecanical
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationVidmantas Jokūbaitis a, Linas Juknevičius b, *, Remigijus Šalna c
Availale online at www.sciencedirect.com Procedia Engineering 57 ( 203 ) 466 472 t International Conference on Modern Building Materials, Structures and Tecniques, MBMST 203 Conditions for Failure of Normal
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationDimensioning of Isolated Footing Submitted to the under Biaxial Bending Considering the Low Concrete Consumption
International Journal Materials Engineering 2017, 7(1): 1-11 DOI: 10.5923/j.ijme.20170701.01 Dimensioning Isolated Footing Submitted to te under Biaxial Bending Considering te Low Concrete Consumption
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationCritical control in transcritical shallow-water flow over two obstacles
Lougboroug University Institutional Repository Critical control in transcritical sallow-water flow over two obstacles Tis item was submitted to Lougboroug University's Institutional Repository by te/an
More informationEFFECTS OF LINE AND PASSIVATION GEOMETRY ON CURVATURE EVOLUTION DURING PROCESSING AND THERMAL CYCLING IN COPPER INTERCONNECT LINES
Acta mater. 48 (000) 3169±3175 www.elsevier.com/locate/actamat EFFECTS OF LINE AND PASSIVATION GEOMETRY ON CURVATURE EVOLUTION DURING PROCESSING AND THERMAL CYCLING IN COPPER INTERCONNECT LINES T.-S. PARK
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationProblem Solving. Problem Solving Process
Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and
More informationFriction Coefficient s Numerical Determination for Hot Flat Steel Rolling at Low Strain Rate
American ournal of Applied Matematics 05; 3(5: -8 Publised online September 6, 05 (ttp://www.sciencepublisinggroup.com/j/ajam doi: 0.648/j.ajam.050305.3 ISS: 330-0043 (Print; ISS: 330-006X (Online riction
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationWeek #15 - Word Problems & Differential Equations Section 8.2
Week #1 - Word Problems & Differential Equations Section 8. From Calculus, Single Variable by Huges-Hallett, Gleason, McCallum et. al. Copyrigt 00 by Jon Wiley & Sons, Inc. Tis material is used by permission
More informationDevelopment of new and validation of existing convection correlations for rooms with displacement ventilation systems
Energy and Buildings 38 (2006) 163 173 www.elsevier.com/locate/enbuild Development of new and validation of existing convection correlations for rooms wit displacement ventilation systems Atila Novoselac
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E. Publised in: IEA Annex 41 working meeting, Kyoto, Japan Publised:
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationPECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER
Materials Pysics and Mecanics (5) 5- Received: Marc 7, 5 PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER Т.I. Belyankova *, V.V. Kalincuk Soutern Scientific Center of Russian
More informationADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS
ADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS Yasuo NIHEI (1) and Takeiro SAKAI (2) (1) Department of Civil Engineering, Tokyo University of Science, 2641 Yamazaki,
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationStress analysis of laminated glass with different interlayer materials
Alexandria Engineering Journal (01) 51, 61 67 Alexandria University Alexandria Engineering Journal www.elsevier.com/locate/aej www.sciencedirect.com Stress analysis of laminated glass wit different interlayer
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationUNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics
UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Annex 41 Kyoto, April 3 rd to 5 t, 2006 Adam Neale 1, Dominique Derome 1, Bert Blocken 2 and Jan Carmeliet 2,3
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationA general articulation angle stability model for non-slewing articulated mobile cranes on slopes *
tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of
More informationThe application of thin-wall shell model in plate problem of syntactic foams
Available online www.ocpr.com Journal of Cemical and Parmaceutical Researc, 14, 6(6:1748-1755 Researc Article ISSN : 975-7384 CODEN(USA : JCPRC5 Te application of tin-wall sell model in plate problem of
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationAPPENDIXES. Let the following constants be established for those using the active Mathcad
3 APPENDIXES Let te following constants be establised for tose using te active Matcad form of tis book: m.. e 9.09389700 0 3 kg Electron rest mass. q.. o.6077330 0 9 coul Electron quantum carge. µ... o.5663706
More informationPolitecnico di Torino. Porto Institutional Repository
Politecnico di Torino Porto Institutional Repository [Article] Coesive Crack Model Description of Ductile to Brittle Size-scale Transition: Dimensional Analysis vs. Renormalization Group Teory Original
More informationAvailable online at ScienceDirect. Procedia Engineering 125 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 15 (015 ) 1065 1069 Te 5t International Conference of Euro Asia Civil Engineering Forum (EACEF-5) Identification of aerodynamic
More informationCLOSED CONVEX SHELLS Meunargia T. Mixed forms of stress-strain relations are given in the form. λ + 2µ θ + 1
Seminar of I. Vekua Institute of Applied Matematics REPORTS, Vol. 43, 207 CLOSED CONVEX SHELLS Meunargia T. Abstract. If Ω is a closed convex sell, ten S : x 3 = 0 is an ovaloid. It is proved tat in tis
More informationJune : 2016 (CBCS) Body. Load
Engineering Mecanics st Semester : Common to all rances Note : Max. marks : 6 (i) ttempt an five questions (ii) ll questions carr equal marks. (iii) nswer sould be precise and to te point onl (iv) ssume
More informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
More informationNotes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.
Capter 6. Fluid Mecanics Notes: Most of te material in tis capter is taken from Young and Freedman, Cap. 12. 6.1 Fluid Statics Fluids, i.e., substances tat can flow, are te subjects of tis capter. But
More informationWhy gravity is not an entropic force
Wy gravity is not an entropic force San Gao Unit for History and Pilosopy of Science & Centre for Time, SOPHI, University of Sydney Email: sgao7319@uni.sydney.edu.au Te remarkable connections between gravity
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationHT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008
Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 August 10-14, 2008, Jacksonville, Florida USA Proceedings of HT2008 2008 ASME Summer Heat Transfer Conference August 10-14, 2008, Jacksonville,
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationarxiv: v1 [math.pr] 28 Dec 2018
Approximating Sepp s constants for te Slepian process Jack Noonan a, Anatoly Zigljavsky a, a Scool of Matematics, Cardiff University, Cardiff, CF4 4AG, UK arxiv:8.0v [mat.pr] 8 Dec 08 Abstract Slepian
More information