A Modified Neuber Method Avoiding Artefacts Under Random Loads
|
|
- Jonah Bates
- 5 years ago
- Views:
Transcription
1 A Modfed Neuber Method Avodng Artefacts Under Random Loads T. Herbland a,b, G. Calletaud a, J. L. Chaboche c, S. Qulc a, F. Gallerneau c a Mnes Pars Pars Tech, CNRS UMR 7633, P 87, vry cedex, France b CTIM, Pôle DSL, Senls, France c ONRA, DMS, P72, Châtllon cedex, France Abstract Surprsngly, Neuber's method s n error when appled to non regular cyclc loads. Some examples are frst shown, n whch the predcted result by the classcal formulaton produces a non physcal local ratchettng n the notch, for a global load made of a large symmetrc cycle and of a small cycle nsde the large one. A modfcaton of the orgnal method s then propopsed to avod the occurence of ths artefact. It s successfully tested aganst fnte element computatons on several notched specmens. 1 INTRODUCTION Many engneerng components subjected to cyclc loadng contan notches lke grooves, holes, keyways, welds, etc. When such a component s loaded, a stress concentraton appears at the notch root. Fatgue can then be reduced snce, even f the component remans globally elastc, a plastc zone develops at the notch root, and an early crack ntaton can be observed. In order to predct fatgue lfe of such components, engneers have to compute the local stresses and strans. Such type of calculaton may be qute long, especally when non-lnear consttutve equatons are used. Smplfed methods are then appled, as an alternatve to FM computaton. Although many solutons have been nvestgated, lke Molsk-Glnka [1], Glnka [2,3], llyn and Kujawsk [4] Neuber s rule [5] remans the most currently used for ndustral applcatons. In the case of unaxal constant ampltude loadngs, t appears to be a good approxmaton; however, a ratchettng phenomenon may appear when the method s appled to varable-ampltude loadngs. In ths paper, a modfcaton to Neuber s rule s then proposed, and tested for several loadngs spectra. 2 PRSNTATION OF NUR S MTHOD Ths method makes t possble to get an approxmaton of the elasto-plastc stress and stran at a notch root, from an elastc computaton, as represented n Fgure 1. σ nom σ nom (σ e = K t σ nom, ε e ) (σ N,ε N ) Fgure 1 : Illustraton of Neuber s rule Result of an elastoplastc computaton ; local stress and stran accordng to equaton (1) If ths notched component, loaded by a nomnal tensle stress σ nom, would reman purely elastc (Fgure 1 ), ths load would nduce a local stress σ e and stran ε e at the notch root. Defnng K t as the stress concentraton factor, σ e can be expressed as σ e = K t σ nom. σ e and ε e can be reached through a smple elastc FM
2 computaton. To calculate the elasto-plastc stress σ N and stran ε N, Neuber [5] postulates that, for a local plastcty, the product of stress and stran at the notch root does not depend on the plastc flow : N N σ ε = e e σ ε (1) Fgure 2 s a graphcal representaton of Neuber s rule n a stress-stran space. quaton (1) mples that the area of the rectangle defned by σ e and ε e s equal to the one defned by σ N and ε N. One can also remark that the elastc pont (σ e, ε e ) and the elasto-plastc pont belong to a hyperbola defned by σ ε = σ e ε e Fgure 2 : Classcal Neuber's method Classcal extenson to cyclc loadng Ths method can be extended to cyclc loadngs. The classcal approach conssts n transferrng the classcal constructon on a cyclc dagram, as llustrated n Fgure 2, usng stress and stran ranges nstead of stress and stran n equaton (1). The subsequent pont s found on the cyclc curves, followng equaton (2) : N σ ε N e e = σ ε (2) Ths soluton s not satsfactory, snce t cannot take nto account mean stresses that would be present n global loadngs. Ths s why more precse strateges have been proposed [6,7] (Fgure 3), based on an update of the stress-stran curves for each reversal. For a gven loadng branch (), the current pont always refers to the last couple (σ N -1, ε N -1) reached n the prevous loadng step. The resultng computaton scheme s gven by equaton (3). For nstance, Fgure 3 llustrates the case of the frst reversal. σ ε 1 ε 2 σ Fgure 3 : xtenson of Neuber s rule to cyclc loadngs wth an updatng strategy
3 N N N N e e e e ( σ σ 1 )( ε ε 1 ) = ( σ σ 1 )( ε ε 1) (3) Once determned the frst maxmum stress pont 1 by the classcal method, one can smply compute the stress and stran for the second extremum 2 by nvertng the (σ,ε) axes and choosng pont 1 as the new orgn. In some case, the local stress tensor s not unaxal. A seres of possble extensons can also be found n the lterature, for nstance, usng Von Mses nvarants nstead of components (equaton (4)) N N J ( σ ) J ( ε ) = J ( σ ) J ( ε ) (4) Ths s stll an open problem, wth no unque soluton for plane stran or axsymmetrc problems. It wll not be treated n ths paper, whch s restrcted to plane-stress cases. N N J ( σ N σ N 1 ) J ( ε ε 1) = J ( σ σ 1) J ( ε ε 1 ) (5) 3 RATCHTTING FFCT FOR VARIAL-AMPLITUD LOADING A 2D plane-stress notched test specmen s chosen, so that the normal stress s null at the notch root, and an unaxal stress state appears wth only σ When applyng Neuber s method to varable ampltude loadngs, we observed a ratchettng phenomenon, even for a lnear knematc hardenng. Ratchettng was not observed when computng the same notch subjected to the same loadng n smulatons wth the FM code ZeuLoN. Moreover, such a localzed ratchet s unrealstc n plastc confnement condtons. Fgure 4 llustrates the loadng geometry. Only one quarter of the test specmen s modelled, to respect symmetry condtons, dsplacements are fxed n drecton 1 on the specmen axs, and fxed n drecton 2 on the bottom. A tensle stress σ nom s appled at the top of the specmen. Fgure 5 presents the response at the notch root under the loadng of Fgure 4. The local elastoplastc stress σ N s drawn wth respect to the local total elastoplastc stran ε N. On Fgure 5 these values are calculated followng the classcal Neuber algorthm, on Fgure 5 they are calculated followng the new Neuber algorthm, that wll be explaned n the next part. σ nom σ Ε (MPa) T t Fgure 4 : Mesh and boundary condtons appled to the specmen ; Loadng hstory, σ e = K t σ nom
4 Table 1 : Materal propertes (lnear knematc hardenng) R 0 (MPa) C (MPa) d (GPa) ν ,3 Fgure 5 : Ratchettng response under the loadng defned n Fgure 4 Intal algorthm Modfed algorthm 4 NUR S MTHOD MODIFICATION Fgure 6 shows an elastc loadng that wll be used to explan the prncple of the new strategy. The problem s that the cyclc rule must take nto account the pont where a small loop reaches the boundary left by the last large loop. The constructon s llustrated n Fgure 7. Steps A and are made accordng to the prevous defned algorthm. The proposed modfcaton s made durng step 3: accordng to classcal plastcty, when the reloadng branch reaches pont C, the subsequent curve must follow the ntal branch, wth ts orgn at pont O, whch produces the evoluton after pont C shown n step 4. Keepng as the orgn of the reloadng branch would produce an extra-hardenng, as shown n step 4. Fgure 6 : response of an elastc computaton σ 1 ε 2 A σ 3 C A O ε 1 σ 2 ε 3 Step 1 Step 2 Step 3
5 σ 4 D σ 4 D Step 4 ε 4 O ε 4 Step 4 Fgure 7 : Neuber's classcal (Steps 1, 2, 3, 4) and modfed (Steps 1, 2, 3, 4 ) algorthm applcaton to loadng of Fgure 6 Ths algorthm can be compared to the ranflow technque: the classcal Neuber method s appled untl pont C, where σ c = σ a, Once the elasto-plastc loop s closed, t s erased from the loadng hstory, and the elasto-plastc soluton for next ponts lke D s calculated by choosng O as orgn. The algorthm s fully detaled n Table 2. It s executed at each stress ncrement. n s the number of reversals to be erased at the begnnng of the algorthm, before computng the soluton (σ Ν, ε Ν ). The value of n used for ncrement s determned at the end of the algorthm for ncrement k-1. Pe represents the ampltude of the reversal whch ends at pont, accordng to Neuber s formula. The resultng expresson s then: Pe J ( σ σ 1 ) J ( ε ε 1). = ( Table 2 : The modfed algorthm (*) For each ncrement k endng at a pont of the loadng hstory: rase the n reversals prevous to the current ncrement from the loadng hstory (n s determned at ncrement k-1) The new orgn for Neuber s calculaton s now the latest peak before to the erased reversals In case of unloadng: The new orgn for Neuber s calculaton s now pont Re-ntalzaton: n = 0 (**) If loadng hstory contans less than 3 reversals: go to (***) If Pe > Pe 1 n n = n + 2 : the two last reversals wll be suppressed at the next ncrement go to (**) (***) Calculaton of σ N, ε Ν accordng to Neuber s rule and go to (*) As shown n Fgure 5, the present method allows the cycle to reman symmetrc for the loadng defned n Fgure 4. Fgure 8 shows a complex hstory computed wth the present algorthm. Reversal suppressons occur at ponts 1, 2 and 3. The result obtaned wth a knematc hardenng rule s shown n Fgure 9. The materal has the same propertes than n Table 1. A non lnear knematc hardenng rule p p X = C ε dx ε, materal propertes n Table 3) gves the result presented n Fgure 10.
6 σ e (MPa) A O C D 1 F G H 3 I 2 K J Tme Fgure 8 : Test loadng hstory showng asymmetrc reversals Fgure 9 : lastoplastc response to loadng of Fgure 8 wth a lnear knematc hardenng followng : Classcal Neuber s algorthm New Neuber algorthm Table 3 : Materal propertes (non-lnear knematc hardenng) R 0 (MPa) C (MPa) d (GPa) ν ,3 Fgure 10 : lastoplastc response to loadng of Fgure 8 wth a non-lnear knematc hardenng followng: Classcal Neuber s algorthm New Neuber algorthm
7 5 CONCLUSIONS Local stress and stran hstory was determned followng classcal Neuber method at the notch root of a notched component under varable-ampltude loadng. The elastoplastc soluton was compared to a soluton reached by FM computaton. A ratchettng effect that dd not appear by usng FM, was observed by usng Neuber s rule. Consderng that ths ratchet was due to asymmetrc cycles, a modfed Neuber method was proposed. Lke the classcal one, ths new algorthm can only be appled to unaxal stress states. It s also based on Neuber equvalence rule, but a change of orgn at certan moments of the hstory makes t possble to avod the shft due to asymmetrc cycles. Ths new algorthm was tested successfully under several complex load spectra, for lnear and non-lnear knematc hardenng rules. 6 RFRNCS [1] K. Molsk and Glnka G, Materal and Scence ngrg, 50, 1981, [2] G. Glnka, ngrg. Fract. Mech., 22, 1985, [3] G. Glnka, ngrg. Fract. Mech., 22, 1985, [4] F. llyn, D. Kujawsk, ngrg. Fract. Mech., 1989, 32, [5] H. Neuber, J. Appl. Mech. (Trans. ASM), 1961, 28, [6] J. Lematre, J.L. Chaboche, Mécanque des matéraux soldes, nd edton. [7] M. Chaudonneret, J.P. Culé, La Rech. Aérosp., nglsh Transl., 1985, , 33-40
χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationNumerical Nonlinear Analysis with the Boundary Element Method
Blucher Mechancal Engneerng Proceedngs May 2014, vol. 1, num. 1 www.proceedngs.blucher.com.br/evento/10wccm Numercal Nonlnear Analyss wth the Boundary Element Method E. Pneda 1, I. Vllaseñor 1 and J. Zapata
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationAn Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors
An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationGraph Reconstruction by Permutations
Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationChapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods
Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary
More informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More information12 MATH 101A: ALGEBRA I, PART C: MULTILINEAR ALGEBRA. 4. Tensor product
12 MATH 101A: ALGEBRA I, PART C: MULTILINEAR ALGEBRA Here s an outlne of what I dd: (1) categorcal defnton (2) constructon (3) lst of basc propertes (4) dstrbutve property (5) rght exactness (6) localzaton
More informationCalculation of time complexity (3%)
Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationConstitutive Modelling of Superplastic AA-5083
TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy
More informationCode_Aster. Identification of the model of Weibull
Verson Ttre : Identfcaton du modèle de Webull Date : 2/09/2009 Page : /8 Responsable : PARROT Aurore Clé : R70209 Révson : Identfcaton of the model of Webull Summary One tackles here the problem of the
More informationSIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS
SIMPLIFIED APPROACH TO THE NON-LINEAR BEHAVIOR OF RC MEMBERS Shahd NASIR 1, Supratc GUPTA 2 And Hdetaka UMEHARA 3 SUMMARY In ths paper, a smplfed one-dmensonal analytcal tool based on fnte dfference technque
More informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationCode_Aster. Identification of the Summarized
Verson Ttre : Identfcaton du modèle de Webull Date : 2/09/2009 Page : /8 Responsable : Aurore PARROT Clé : R70209 Révson : 609 Identfcaton of the Summarzed Webull model One tackles here the problem of
More informationKinematic Hardening Parameters Identification with Respect to Objective Function
Vol:8, No:4, 04 Knematc Hardenng Parameters Identfcaton wth Respect to Objectve Functon Marna Franulovc, Robert Basan, Bozdar Krzan Internatonal Scence Index, Mathematcal and Computatonal Scences Vol:8,
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationBallot Paths Avoiding Depth Zero Patterns
Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,
More informationGEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE
GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of cvl engneerng, IIT Bombay, Powa, Mumba 400076, Inda. Tel.022-25767328 emal: cejnm@cvl.tb.ac.n Module - 9 LECTURE - 48
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationAn Interactive Optimisation Tool for Allocation Problems
An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationPARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationUSE OF FATIGUE ANALYSIS SOFTWARE IN THE OPTIMIZATION OF STRUCTURAL
USE OF FAIGUE ANALYSIS SOFWARE IN HE OPIMIZAION OF SRUCURAL APPROVAL PROCESS FOR VEHICLE COMPONENS Scozzafave C.C., Ercsson L. G. S., Santana F.A., Wldmann U., Spnell D. M. Mercedes-Benz do Brasl E-mals:
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationBe true to your work, your word, and your friend.
Chemstry 13 NT Be true to your work, your word, and your frend. Henry Davd Thoreau 1 Chem 13 NT Chemcal Equlbrum Module Usng the Equlbrum Constant Interpretng the Equlbrum Constant Predctng the Drecton
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationNotes on Frequency Estimation in Data Streams
Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to
More informationFUZZY FINITE ELEMENT METHOD
FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationLNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS
CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS CONTENTS 1. Method for determnng transferred energy durng cargo transfer. Calculatng the transferred energy.1 Calculatng the gross transferred energy.1.1
More informationMAE140 - Linear Circuits - Fall 10 Midterm, October 28
M140 - Lnear rcuts - Fall 10 Mdterm, October 28 nstructons () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator
More informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationLifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model
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
More informationTOPICS MULTIPLIERLESS FILTER DESIGN ELEMENTARY SCHOOL ALGORITHM MULTIPLICATION
1 2 MULTIPLIERLESS FILTER DESIGN Realzaton of flters wthout full-fledged multplers Some sldes based on support materal by W. Wolf for hs book Modern VLSI Desgn, 3 rd edton. Partly based on followng papers:
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationA comprehensive study: Boundary conditions for representative volume elements (RVE) of composites
Insttute of Structural Mechancs A comprehensve study: Boundary condtons for representatve volume elements (RVE) of compostes Srhar Kurukur A techncal report on homogenzaton technques A comprehensve study:
More informationLecture 17 : Stochastic Processes II
: Stochastc Processes II 1 Contnuous-tme stochastc process So far we have studed dscrete-tme stochastc processes. We studed the concept of Makov chans and martngales, tme seres analyss, and regresson analyss
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
More informationMaximizing the number of nonnegative subsets
Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum
More informationLearning Theory: Lecture Notes
Learnng Theory: Lecture Notes Lecturer: Kamalka Chaudhur Scrbe: Qush Wang October 27, 2012 1 The Agnostc PAC Model Recall that one of the constrants of the PAC model s that the data dstrbuton has to be
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationMATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)
1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons
More informationOFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES
ICAMS 204 5 th Internatonal Conference on Advanced Materals and Systems OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES VLAD LUPĂŞTEANU, NICOLAE ŢĂRANU, RALUCA HOHAN, PAUL CIOBANU Gh. Asach Techncal Unversty
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More information