Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
|
|
- Bethany Daniel
- 5 years ago
- Views:
Transcription
1 Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at 0 ( 0 ) ; ( 0) o o o Defne the followng varables, Y= X; Y = X () and wrte EOM () as two frst-order Eqs. M Y + DY + KY = F() t & Y = Y (3) whch can be wrtten n matrx form as Y 0 Y 0 = Y + M K M D Y M F (4) Or, Y = AY+ b (5) Y = X 0 0 wth Y= ; ; Y X A= = M K M D b = M F Ths s known as the state-space formulaton. Eq. (5) s to be ntegrated numercally wth ntal condton vector Y = X V [ ] T o o o MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008
2 If the appled load s NOT a functon of tme, then an equlbrum state s defned after a very long tme as Y 0 YE = A bo (6) Computatonal software such as Mathcad, Mapple, Mathematca, Matlab, etc has bult-n functons or commands to perform the numercal ntegraton of equatons set n the form Y = AY+ b, even when system s nonlnear,.e. A=A(Y). A few words about numercal ntegraton methods Typcal numercal ntegraton methods nclude a) Euler (smple) method b) Fourth and Ffth-Order Runge-Kutta ntegrators, c) Rosenbrock Method see references on page, d) Adams Predctor Corrector Methods e) Average Acceleraton and Wlson-θ (Implct) Methods In most methods, the selecton of an adequate tme step s crucal for numercally stable and accurate results. (a)-(b) are favored by the young ntates nto numercal computng and because of ther ready avalablty n modern computatonal software. (c) (d) are more modern (mplct) methods wth automated ntermedate reszng of the tme step whle performng the ntegraton. Methods (e) have long been favored by structural mechanc analysts when ntegratng Multple DOF (lnear) systems All methods suffer from defcences when nonlneartes are apparent thus forcng extremely small tme steps and the ensung MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008
3 cost wth lots of numercal computng (tme). (Memory) Storage appears not to be an ssue anymore. State-space method for MDOF systems. Recall the EOMS for a lnear system are where MU+DU+KU () t =F() t (7) U,U, and U are the vectors of generalzed dsplacement, F s the vector of velocty and acceleraton, respectvely; and () t generalzed (external forces) actng on the system. M,D,K represent the matrces of nerta, vscous dampng and stffness coeffcents, respectvely. Defne the followng varables, Y=U ; Y=U and wrte EOM (7) as a set of n-frst-order Eqs. (t) (8) MY +DY +KY =F & Y =Y (9) whch can be wrtten n matrx form as Y 0 I Y 0 = Y -M K -M D Y -M F (0) Or, Y = AY+ b () U 0 I 0 wth Y= ; A= ; b= U -M K -M D M F () The matrces are square wth n-rows = n columns, whle the vectors are n- rows. MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 3
4 and ntal condtons Y( t= 0) = Uo U o. A s a n x n matrx. I s the nxn dentty matrx, and 0 s a nxn matrx full of zeroes. Condtons for a good numercal ntegrator In general a numercal ntegraton scheme should a) reproduce EOM as tme step 0 b) provde, as wth physcal model, bounded solutons for any sze of tme step,.e. method should be stable c) reproduce the physcal response wth fdelty and accuracy. The numercal ntegraton reles n representng tme dervatves of a functon wth an algebrac approxmaton, for example T dx Δ x x = = lm dt Δ t 0 xt+ xt x x x ~ Eq. above s exact only f 0 + (3) Numercal ntegraton methods are usually dvded nto two categores, mplct and explct. x = f x, t (4) Consder the ODE In an explct numercal scheme, the ODE s represented n terms of known values at a pror tme step,.e. x = x +Δ + t f x, t, (5) MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 4
5 whle n an mplct numercal scheme x, + x t f x+ t = +Δ (6) Explct numercal schemes are condtonally stable. That s, they provde bounded numercal solutons for (very) small tme steps. For example, Tn τ crt = (7) π π where T and K n= ωn = M are the natural perod and natural ωn frequency of the system, respectvely. The restrcton on the tme step s too severe when analyzng stff systems,.e. those wth large natural frequences. Implct numercal schemes are uncondtonally stable,.e. do not mpose a restrcton on the sze of the tme step Δ t. (However, accuracy may be compromsed f too large tme steps are used). ANALOGY between numercal schemes and a flter A few words of wsdom released n class MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 5
6 The average acceleraton method for numercal ntegraton of SDOF equaton: M X + DX + K X = F() t () Consder a change of thnkng frame by defnng Arthmetc ~ Contnuum (functon) X X t, F F t,, X X t F F t (7) Wrte Eq. () at two tmes, t = t M X + DX + K X = F (8a) at t = t M X + DX + K X = F (8b) at Subtract (b) from (a) to obtan: where M Δ X + DΔ X + KΔ X =ΔF (9) + Δ X = X+ X Δ F= F+ F, Δ X = X X, etc. Note that known quanttes at t=t are{ X, X, X} (0) Ths numercal method s extremely popular among the structural dynamcs communty. Its extenson to MDOF systems wll be shown later. The other favorte method, Wlson-θ scheme, wll also be gven n later lectures (MDOF systems). MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 6
7 Now, assume the acceleraton s constant wthn the tme nterval Δ t = t t + ( τ) for 0 () X = a < τ set as an average value a= X + + X. The velocty and dsplacement follow from ntegraton of Eq. (6) wthn the tme nterval,.e. ( τ ) = + X X a τ (a) ( τ ) = + τ + τ X X X a (b) Acceleraton velocty X + X X ( τ ) = a X + X ( τ ) = + X X a τ τ t t + t τ t t + t dsplacement X + X τ = X + X τ + a τ X τ t t + t MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 7
8 At the end of the tme nterval, the velocty and dsplacement equal X = X Δ t = + X + a Δ t (3a) X = X Δ t = X + + X Δ t + a Δ t (3b) And the dfferences n velocty and dsplacement re ( + ) ( + ) ( X X X + ) t ( X X ) t Δ X = X X = a Δ t = X + X = + Δ = = +Δ Δ Δ X = X X = X Δ t + a + ( + ) = X Δ t + Δ t X + X X + X 4 = X Δ t + Δ X + X 4 (4a) (5b) 4 (6a) from (5b), Δ X = X+ ( ΔX X) and nto (4a) ( ) Δ X = X +ΔX 4 = X X + ( ΔX X ) MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 8
9 Δ X = ΔX X (6b) Note that n Eqs (6), { ΔX, ΔX } obtaned at the pror tme step,.e. { X, X } Δ X. Thus, replace { ΔX, ΔX }, depend on the known values and the unknown nto the dfference equaton (9), M Δ X + DΔ X + KΔ X =ΔF 4 M X + ( ΔX X Δ t ) + D ΔX X + KΔ X =ΔF Rearrangng terms leads to where K Δ X =Δ F (7) * * * 4 K = K+ D+ M (8a) * 4 Δ F =Δ F+ M X + D+ M X (8b) are known as pseudo dynamc stffness and dynamc force, respectvely The recpe for the numercal ntegraton to fnd the system tme response s MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 9
10 At tme t, known varables are { X, X } (current state) () fnd from EOM: X = M ( F DX K X) () form pseudo stffness and forcng functons, ( K *, F) Eqs. (8), * (3) Calculate Δ X = K Δ F, and (4) X + = X+Δ X, X + = X +ΔX at t + (5) Increase tme to t + and return to step () Δ from Δ X = ΔX X ; The average acceleraton method s an mplct method,.e. numercally stable and consstent. The dsadvantage s that t requres memory 3 to store X, X, ΔX, ΔX, ΔF. Average acceleraton methods for numercal ntegraton of a nonlnear system Consder the system wth EOM where g( X, X) ( 3 ) o o 3 M X + g X, X = F( t) (30) s a nonlnear functon, for example g X, X = g + k X + k X + F sgn( X ) As wth the lnear system, evaluate Eq. (30) at two tmes (closely spaced): at t = t M X + g X, X = F (3a) μ 3 A non-ssue n the st century MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 0
11 at t = t M X + g X, X = F (3b) Subtract (b) from (a) to obtan: where (3) M Δ X + g g =ΔF + Δ X = X X Δ F = F F, + + (, ), (, ) g = g X X g = g X X (33) A Taylor seres expanson of the nonlnear functon at t gves g g g = g + Δ X + Δ X + O ΔX ΔX (, ) + X X X, X = 0 X (34) defne local lnearzed stffness and dampng coeffcents as Hence, K g = ; D = X g X, X X X, X g g K Δ + X + D Δ X (35) and the dfference Eq. (3) becomes lnear M Δ X + D Δ X + K Δ X =ΔF (36) Eq. (35) s formally dentcal to the one devsed for a lnear system. Thus, the numercal treatment s smlar, except that at each tme MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008
12 step, lnearzed stffness and dampng coeffcents need be calculated. The recpe s thus dentcal; however wth the apparent nonlnearty, the method does not guarantee stablty for (too) large tme steps. The recpe for the numercal ntegraton to fnd the system tme response s (current state) At tme t, known varables are { X, X } () fnd from EOM: X = M F g( X, X ) (a) fnd local (lnearzed) stffness and dampng coeffcents, (K, D ) from eq. (35) K *, Δ F from (b) form pseudo stffness and forcng functons, ( ) * 4 K = K+ D + M * (3) Calculate ; * 4 Δ F =Δ F + M X + D + M X Δ X = K Δ F, and Δ X = ΔX X ; X = X +Δ X, X = X +ΔX at t + (4) Set + + (5) Increase tme to t + and return to step () References The followng are a must! Press, W.H., Flannery, B.P., Teukolsky, S.A., and Vetterlng, W.T., 986 ( st edton), Numercal Recpes, The Art of Scentfc Computng, Cambrdge Unversty Press, UK. Bathe, J-K, 98 ( st ed), 007 latest, Fnte Element Procedures, Prentce Hall. Pche, R., and P. Nevalanen, 999, Varable Step Rosenbrock Algorthm for Transent Response of Damped Structures, Proc. IMechE, Vol. 3, part C, Paper C MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008
13 Rosenbrock, H.H., Some General Implct Processes for the Numercal Soluton of Dfferental Equatons, Appendx A 4 Numercal Integraton Usng Modfed Euler s Method It s often dffcult to solve (exactly) a nonlnear dfferental equaton. Numercal ntegraton s then employed to obtan results (predctons of moton). The Modfed Euler Method s one type of numercal ntegraton scheme. Solve for q(t) governed by M q + Cq + Kq= Q( q, q, t) (A.) wth ntal condtons q0and q 0at t = 0. In Eqn. (A.), Qq (, qt, ) may contan terms that are nonlnear n Let V dq dt = (A.) and wrte Eq. (A.) as a system of TWO frst order dfferental equatons,.e. C K Q( q, V, t) V = V q +, M M M q = V (A.3) 4 Ths Appendx s ncluded because most young engneerng students have learned only about Euler s method. Hence, the appendx complements ther educaton.. However, I encourage you to abandon the usage of ths poor method. MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 3
14 Defne K C ωn = andζ = (A.4) M K M as the natural frequency and vscous dampng rato, respectvely. Wth the noted defntons Eqs. (A.3) become dv Q( q, V, t) V = = F ( q, V, t) = ζωnv ωn q + dt M dq q = == F ( q, V, t) = V dt (A.5) Note that F and F are the slopes of the (V vs t) and (q vs t) curves, respectvely. q(t) V(t) q+ V q V+ t t t+ t t+ t Let the numercal approxmatons (arthmetc values) be,and V V t q q t (A.6) t = ; = and Δ t s a sutably small tme step for numercal ntegraton. where 0,,... F q+ In Euler s numercal scheme, approxmate the tme dervatves (or slopes) as: q V F V+ MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 4
15 V V +Δ t f, f = F q, V, t, + q q +Δ t f, f = F q, V, t + for =,.. (A.7) wth ntal condtons q0and q 0at t 0 = 0 Eq. (7) offer a recursve relaton to calculate the numercal (arthmetc) values of the varables V and q at successve tmes t). The regular Euler method s frst order wth a truncaton error of order ( ). A modfed Euler method (second order accurate) wth error order ( 3 ) follows: (a) Compute prelmnary estmates of ( V, q ) (,, ) V V +F q, V, t, + q q +F q V t as (A8.a) (b) Use these prelmnary estmates to obtan mproved slopes as + + ( ) ( ) f = F q, V, t, f = F q, V, t, (A8.b) (c) Defne average slopes as per f = f + f f = f + f + + (A8.c) MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 5
16 (c) and obtan new estmatons usng the averaged slopes,.e. V V + f + q q + f +, (A8.d) (d) Repeat steps (a) through (d) at each tme step, t= Δ t ; = 0,,..., and startng wth the ntal condtons q and q at t = The sze of the tme step s very mportant to obtan accurate and numercally stable results. If s too large, then numercal predctons wll be n error and very lkely show unstable (oscllatng) results. If s too small, then the numercal method wll be slow and costly snce the number of operatons ncreases accordngly. In practce, a tme step of the order =T n /60, where T n s the natural perod of moton gven as (π/ω n ). Euler s method s most tmes NOT a good choce to perform the numercal ntegraton of lnear or nonlnear ODES. Alas, t s wdely used by rooke engneers and engneerng students because t s easy to mplement. Often enough, however, predctons can be wrong and msleadng. I call Euler s method a BRUTE FORCE approach, Often regarded as an art, numercal computng s n actualty a well establshed branch of mathematcs. MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 6
17 MEEN 67 HD 6 Numercal Integraton for Tme Response: SDOF system L. San Andrés 008 7
CHAPTER 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 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 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 informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationE91: Dynamics. E91: Dynamics. Numerical Integration & State Space Representation
E91: Dnamcs Numercal Integraton & State Space Representaton The Algorthm In steps of Δ f ( new ) f ( old ) df ( d old ) Δ Numercal Integraton of ODEs d d f() h Δ Intal value problem: Gven the ntal state
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.
Handou # 6 (MEEN 67) Numercal Inegraon o Fnd Tme Response of SDOF mechancal sysem Sae Space Mehod The EOM for a lnear sysem s M X DX K X F() () X X X X V wh nal condons, a 0 0 ; 0 Defne he followng varables,
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 informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
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 informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
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 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 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 informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
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 information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
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 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 information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationMEEM 3700 Mechanical Vibrations
MEEM 700 Mechancal Vbratons Mohan D. Rao Chuck Van Karsen Mechancal Engneerng-Engneerng Mechancs Mchgan echnologcal Unversty Copyrght 00 Lecture & MEEM 700 Multple Degree of Freedom Systems (ext: S.S.
More informationSolution of Equilibrium Equation in Dynamic Analysis. Mode Superposition. Dominik Hauswirth Method of Finite Elements II Page 1
Soluton of Equlbrum Equaton n Dynamc Analyss Mode Superposton Domnk Hauswrth..7 Method of Fnte Elements II Page Contents. Mode Superposton: Idea and Equatons. Example 9.7 3. Modes 4. Include Dampng 5.
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationAs it can be observed from Fig. a) and b), applying Newton s Law for the tangential force component results into: mg sin mat
PHY4HF Exercse : Numercal ntegraton methods The Pendulum Startng wth small angles of oscllaton, you wll get expermental data on a smple pendulum and wll wrte a Python program to solve the equaton of moton.
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 informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
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 informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
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 informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy
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 informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
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 informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationNUMERICAL SCHEME INFLUENCE ON PSEUDO DYNAMIC TESTS RESULTS
October 1-17, 8, Bejng, Chna NUMERICAL SCHEME INFLUENCE ON PSEUDO DYNAMIC TESTS RESULTS ABSTRACT : G. Amato 1, L. Cavaler 1 PhD Student, Dpartmento d Ingegnera Strutturale e Geotecnca, Palermo. Italy Assocate
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationCALCULUS CLASSROOM CAPSULES
CALCULUS CLASSROOM CAPSULES SESSION S86 Dr. Sham Alfred Rartan Valley Communty College salfred@rartanval.edu 38th AMATYC Annual Conference Jacksonvlle, Florda November 8-, 202 2 Calculus Classroom Capsules
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
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 informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
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 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 informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More information6.3.7 Example with Runga Kutta 4 th order method
6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
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 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 informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationTracking with Kalman Filter
Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,
More informationHidden Markov Models
Hdden Markov Models Namrata Vaswan, Iowa State Unversty Aprl 24, 204 Hdden Markov Model Defntons and Examples Defntons:. A hdden Markov model (HMM) refers to a set of hdden states X 0, X,..., X t,...,
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
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 informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationarxiv: v1 [math.ho] 18 May 2008
Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc
More informationAssessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion
Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
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 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 informationStart with the equation of motion for a linear multi-degree of freedom system with base ground excitation:
SE 80 Earthquake Enneern November 3, 00 STEP-BY-STEP PROCEDURE FOR SETTING UP A SPREADSHEET FOR USING NEWMARK S METHOD AND MODAL ANALYSIS TO SOLVE FOR THE RESPONSE OF A MULTI-DEGREE OF FREEDOM (MDOF) SYSTEM
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 informationModified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls [Verson: 2000/03/22 v.0] Modfed Mass Matrces and Postvty Preservaton for Hyperbolc and
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationWeek 5: Neural Networks
Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationG4023 Mid-Term Exam #1 Solutions
Exam1Solutons.nb 1 G03 Md-Term Exam #1 Solutons 1-Oct-0, 1:10 p.m to :5 p.m n 1 Pupn Ths exam s open-book, open-notes. You may also use prnt-outs of the homework solutons and a calculator. 1 (30 ponts,
More informationAdjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.
th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Adjont Methods of Senstvty Analyss for Lyapunov Equaton Bopng Wang, Kun Yan Department of Mechancal and Aerospace
More informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
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 informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
More informationSOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Kim Gaik Universiti Tun Hussein Onn Malaysia
SOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Km Gak Unverst Tun Hussen Onn Malaysa Kek Se Long Unverst Tun Hussen Onn Malaysa Rosmla Abdul-Kahar
More informationDiscretization. Consistency. Exact Solution Convergence x, t --> 0. Figure 5.1: Relation between consistency, stability, and convergence.
Chapter 5 Theory The numercal smulaton of PDE s requres careful consderatons of propertes of the approxmate soluton. A necessary condton for the scheme used to model a physcal problem s the consstency
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationAE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Comptatonal Fld Dynamcs (CFD) Comptatonal Fld Dynamcs (AE/ME 339) Implct form of dfference eqaton In the prevos explct method, the solton at tme level n,,n, depended only on the known vales of,
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationPHYS 705: Classical Mechanics. Hamilton-Jacobi Equation
1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More information