IMPULSE-BASED SIMULATION OF INEXTENSIBLE CLOTH
|
|
- May Carter
- 6 years ago
- Views:
Transcription
1 IMPULSE-BASED SIMULATION OF INEXTENSIBLE CLOTH Jn Bender nd Dniel Byer Institut für Betries- und Dilogssysteme Universität Krlsrue Am Fsnengrten Krlsrue Germny ABSTRACT In tis pper n impulse-sed metod for clot simultion is presented. Te simultion of clot is required in different ppliction res like computer nimtion, virtul relity or computer gmes. Simultion metods often ssume tt clot is n elstic mteril. Wit tis ssumption te simultion cn e performed very efficiently using spring forces. Te prolem is tt mny textiles cnnot e stretced significntly. A relistic simultion of tese textiles wit spring forces leds to stiff differentil equtions wic cuse deteriortion of performnce. Te impulse-sed metod descried in tis pper solves tis prolem nd llows te relistic simultion of inelstic textiles. KEYWORDS Clot simultion, pysiclly-sed modelling, impulse-sed simultion, inelstic textiles 1. INTRODUCTION Te simultion of clot nd frics is n import re of reserc nd s mny pplictions like computer nimtion, gmes nd virtul environments. Terefore tis re s long istory of reserc. House nd Breen (2000) nd Mgnent-Tlmnn nd Volino (2005) give good survey of te reserc done in clot nimtions. Current reserc prolems in cloting simultion re summrized in Coi nd Ko (2005). Clot is often simulted s n elstic mteril due to performnce resons ut mny textiles re not noticely stretcle. Mny pproces use spring forces for te simultion of stretcle clot. In fct tody's cloting simultions re minly sed on spring-mss systems of prticles (see, e.g. Coi nd Ko (2002)). Te relistic simultion of n inextensile piece of clot wit suc n pproc would require lrge spring constnts. Tis leds to stiff differentil equtions wic re rd to integrte nd decrese te numericl stility (Hut et l. (2003)). Te integrtion cn only e performed y reducing te time step size significntly or y using specil metods. In ot cses te consequence is low performnce. Anoter pproc is to use constrints to enforce te conditions insted of integrting te spring forces directly (see, e.g. Furmnn et l. (2003)). Especilly if smll tolernces re demnded tis pproc provides good lterntive. A system for constrined-sed simultion of inextensile clot sed on te lgrngin mecnics is given in Goldentl et l. (2007). In tis pper new constrint-sed pproc is presented tt uses impulses for te simultion of clot. It is sown tt te use of impulses insted of spring forces llows n ccurte solution of te constrints. 2. CLOTH SIMULATION A piece of clot is represented y mes of prticles tt re linked y distnce constrints. Te following sections descrie te simultion of prticles nd constrints nd te ndling of teir dependencies in mes.
2 2.1 Prticle simultion A prticle is ody tt s mss ut no volume. Since prticle s no dimension, it s just trnsltionl degrees of freedom nd no rottionl ones. In te simultion te stte of prticle is descried y its mss m its position c nd its liner velocity v. In generl te mss of prticle is constnt during te simultion. A time step of prticle is performed y integrting its velocity nd its position over time. It is ssumed tt te sum of ll externl forces cting on prticle, like e.g. grvity, is constnt during te time step. In tis cse te velocity nd position fter step of size cn e computed directly y te following equtions: Fext Fext v( t0 + ) = v( + dt = v( t ) m 0 + m 0 Fext 1 Fext 2 c( t0 + ) = c( + v( + t dt = c( + v( + m 2 m 0 were Fext is te sum of ll externl forces nd t0 is te strt time of te simultion step. In te simultion it is differentited etween dynmic nd sttic prticles. In contrst to dynmic prticles sttic ones ve fix position nd no velocity. 2.2 Distnce constrints A distnce constrint for two prticles nd consists of two prts: position constrint nd constrint for te velocities of te prticles. Te position constrint is defined s follows: d( d(0) = 0 wit d( = c( c (. Tis mens tt te difference of teir ctul distnce nd teir distnce t te eginning of te simultion must e zero to stisfy tis constrint. Hence teir distnce must sty constnt over time. In te simultion tis constrint is stisfied y impulses. At time t tese impulses re computed y using preview of te simultion step. For distnce constrint te distnce of te corresponding prticles d(t+) fter te next simultion step of size is determined y integrting te prticle positions (see section 2.1). Te difference etween te vlue d(t+) nd te distnce d(0) t te eginning of te simultion is exctly te error tt would occur, if te simultion step would e performed witout regrding te position constrint. Tis error e = d( t + ) d(0) cn e eliminted y computing impulses nd pplying tem t te eginning of te simultion step t time t. Since n unconstrined prticle s liner motion, te two prticles must cnge teir reltive velocity y e / in direction of te constrint in order to cnge teir distnce y e in time step of size. Te direction of te constrint is given y te vector from prticle to prticle. In order to stisfy te conservtion of momentum two impulses p nd p of te sme mgnitude nd opposite directions re pplied to perform tis velocity cnge. Te impulse p cn e determined y solving te eqution: e Δ v ( p, Δv ( p, = ( c ( c ( ) were Δv ( p, is te velocity cnge of prticle, wen te impulse p is pplied. Te differentition etween dynmic nd sttic prticles is done using te following vlue for computing te velocity cnge: k 1 / m, if prticle is dynmic = 0, oterwise. Te resulting eqution for te impulse is: Δ v ( p, Δv ( p, = k p k( p) = ( k + k) p = ( c ( c ( ) e
3 wic cn e solved, if t lest one of te prticles is dynmic nd te prticles ve different positions. It is ssumed tt te constrint is stisfied t te eginning of simultion step, so te prticles cnnot ve te sme positions. After te determintion of te impulse p it must e pplied to te prticles in positive nd negtive direction respectively. Due to te resulting velocity cnge, te position constrint will e stisfied fter te simultion step. Te second prt of te distnce constrint is constrint for te velocities of te prticles: ( v () t v () t )( c ( c ( ) = 0. Tis mens tt te reltive velocity of te prticles in direction of te constrint must e zero. In contrst to te position constrint, ere te required impulse cn e determined directly witout preview ecuse n impulse cuses n instntneous velocity cnge. Te impulse tt stisfies te constrint is computed y: ( v () t v ( )( c () t c () t ) p =. k + k Te computtion of n impulse in order to correct te velocities is not solutely necessry for te simultion of distnce constrint. A iger degree of ccurcy cn e cieved y regrding te velocity constrint in te simultion ut t te cost of performnce. 2.3 Mes simultion In te simultion clot is represented y mes. In ec vertex of tis mes one prticle is creted. A simple pproc to simulte stretcle clot is to introduce dmped spring on ec edge of te mes. Te elsticity of te clot is controlled y te spring constnt nd te dmping coefficient. However te simultion of inextensile clot using spring forces leds to stiff differentil equtions wic cuse significnt decrese of performnce. In order to simulte clot using impulses distnce constrint s descried ove is defined for ec pir of prticles tt re connected y n edge in te mes. Te constrint is stisfied, if te distnce etween te positions of te two corresponding prticles equls te lengt of te edge t te eginning of te simultion. If ec constrint is exctly stisfied, te simulted mes is inextensile. In suc mes prticle is linked y multiple distnce constrints to oter prticles. Te impulses of constrints wic ve common prticle influence ec oter. Te resulting dependencies cn e resolved, if te constrints re ndled in n itertive process. In n itertion step n impulse is computed for ec distnce constrint s descried in section 2.2. Since te impulses influence ec oter te constrints re not stisfied witin one step ut te errors re reduced. Te itertive process stops, wen ll constrints re stisfied witin predefined tolernce. Scmitt et l. (2005) proofed tt tis itertive metod converges to te pysiclly correct solution. Te itertive metod s mny dvntges. It is very simple to implement. Collision nd contct ndling wit friction cn e esily integrted in te itertive process (Bender et l., 2006). Te mximl llowed extensiility of te simulted clot is directly defined y te used tolernce vlue. Te itertive process cn e interrupted t ny time to get preliminry result. Even if te process is interrupted nd te tolernce vlue is not reced, te simultion stys stle. All constrints wic ve no common prticle nd terefore no direct dependency cn e solved in prllel. Te disdvntge of te itertive pproc is tt te simultion of complex mes requires mny itertions, if smll tolernce vlue is used. Te numer of itertions cn e reduced significntly y sudividing te mes into strips of constrints. Te distnce constrints of strip must e cyclic. In tis cse ec strip cn e simulted y te liner-time metod of Jn Bender (2007). Te dependencies etween te single strips re still resolved y n itertive process ut since only few strips re required for mes, now tis process converges very fst. 3. RESULTS Te presented metod ws implemented in C++. For te ccurte ndling of collisions nd contcts wit friction te constrint-sed collision response metod of Bender et l. (2006) ws integrted in te itertive process descried ove.
4 Figure 1 sows te results of two different simultions wit te impulse-sed metod descried in tis pper. Te tolernce vlue used in te itertive process defines te mximl extensiility of te simulted clot. In ot simultions te cosen vlue prevented te clot from stretcing more tn 0.01 percent. Te introduced metod lso works wit muc smller tolernce vlues ut t te cost of performnce. By integrting te liner-time metod of Jn Bender (2007) s mentioned ove mes of 900 prticles nd 1740 distnce constrints cn e simulted out tree times fster tn rel-time on PC wit 2.4 GHz Intel Core 2 Qud processor. Figure 1. Results of te impulse-sed clot simultion 4. CONCLUSION Te presented metod for clot simultion llows te ccurte ndling of inextensile textiles. Tese textiles re simulted y using mes of prticles linked y distnce constrints wic re stisfied y te computtion of impulses. An impulse is determined y using preview of te corresponding constrint stte. Te constrint dependencies in mes re resolved eiter itertively or in comintion wit liner-time metod to increse te performnce. Te extensiility of te clot cn e controlled y tolernce vlue. One of te min fetures of te new metod is its stility. Since constrint is directly resolved y using preview, even totlly destroyed models cn e repired. Hence in ec step n erly result cn e otined t ny time y interrupting te itertive process witout te loss of stility. Tis cn e used to increse te performnce. Oter fetures re tt te metod is simple to implement nd tt collision nd contct ndling wit friction cn e esily integrted in te presented simultion process. REFERENCES Bender, J. nd Scmitt, A., 2006, Constrint-sed collision nd contct ndling using impulses. Proceedings of te 19t interntionl conference on computer nimtion & socil gents, Genev, Switzerlnd, pp Bender, J., 2007, Impulse-sed dynmic simultion in liner time. In Journl of Computer Animtion nd Virtul Worlds, Vol. 18, No. 4-5, pp Coi, K.-J. nd Ko, H.-S, Stle ut responsive clot. ACM Trns. on Grpics, Vol. 21, No. 3, pp Coi, K.-J. nd Ko, H.-S, Reserc prolems in cloting simultion. Computer-Aided Design, Vol. 37, No. 6, pp Furmnn, A. et l., Interctive nimtion of clot including self collision detection. WSCG 03, pp Goldentl, Rony et l., Efficent Simultion of Inextensile Clot. ACM Trnsctions on Grpics (Proceedings of SIGGRAPH 2007), Vol. 26, No. 3. Hut, M. et l., Anlysis of numericl metods for te simultion of deformle models. Te Visul Computer, Vol. 19, No. 7-8, pp
5 House, D.H. nd Breen, D.E., Clot modeling nd nimtion. A.K. Peters, Ltd., Ntick, MA, USA. Mgnent-Tlmnn, N. nd Volino, P., From erly drping to ute couture models: 20 yers of reserc. Te Visul Computer, Vol. 21, No. 8-10, pp Scmitt, A. et l., 2005, On te Convergence nd Correctness of Impulse-Bsed Dynmic Simultion. Tecnicl Report University of Krlsrue, Germny.
Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),
More informationKey words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.
O SPECTRA ACCURACY OF QUADRATURE FORMUAE BASED O PIECEWISE POYOMIA ITERPOATIO A KURGAOV AD S TSYKOV Abstrct It is well-known tt te trpezoidl rule, wile being only second-order ccurte in generl, improves
More informationForces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc
Physics 170 Summry of Results from Lecture Kinemticl Vribles The position vector ~r(t) cn be resolved into its Crtesin components: ~r(t) =x(t)^i + y(t)^j + z(t)^k. Rtes of Chnge Velocity ~v(t) = d~r(t)=
More informationThe practical version
Roerto s Notes on Integrl Clculus Chpter 4: Definite integrls nd the FTC Section 7 The Fundmentl Theorem of Clculus: The prcticl version Wht you need to know lredy: The theoreticl version of the FTC. Wht
More information4 The dynamical FRW universe
4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationSection 4: Integration ECO4112F 2011
Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic
More information12 Basic Integration in R
14.102, Mt for Economists Fll 2004 Lecture Notes, 10/14/2004 Tese notes re primrily bsed on tose written by Andrei Bremzen for 14.102 in 2002/3, nd by Mrek Pyci for te MIT Mt Cmp in 2003/4. I ve mde only
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationI1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3
2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is
More informationEcon 401A Three extra questions John Riley. Homework 3 Due Tuesday, Nov 28
Econ 40 ree etr uestions Jon Riley Homework Due uesdy, Nov 8 Finncil engineering in coconut economy ere re two risky ssets Plnttion s gross stte contingent return of z (60,80) e mrket vlue of tis lnttion
More informationMath Week 5 concepts and homework, due Friday February 10
Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationLogarithms and Exponential Functions. Gerda de Vries & John S. Macnab. match as necessary, or to work these results into other lessons.
Logritms nd Eponentil Functions Gerd de Vries & Jon S. Mcn It is epected tt students re lred fmilir wit tis mteril. We include it ere for completeness. Te tree lessons given ere re ver sort. Te tecer is
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationDerivatives, Finding Roots, Interpolating and Integration
Derivtives, Finding Roots, Interpolting nd Integrtion Doing Numericl Derivtives Altoug tking derivtives is very esy opertion for ll of you, tere re situtions were it must e done numericlly. For exmple,
More informationTopic 6b Finite Difference Approximations
/8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationSection 6.1 Definite Integral
Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationSuppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.
Mth 43 Section 6. Section 6.: Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationPhysicsAndMathsTutor.com
1. A uniform circulr disc hs mss m, centre O nd rdius. It is free to rotte bout fixed smooth horizontl xis L which lies in the sme plne s the disc nd which is tngentil to the disc t the point A. The disc
More informationA P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES
A POWERS OF TEN AND SCIENTIFIC NOTATION In science, very lrge nd very smll deciml numbers re conveniently expressed in terms of powers of ten, some of wic re listed below: 0 3 0 0 0 000 0 3 0 0 0 0.00
More informationParse trees, ambiguity, and Chomsky normal form
Prse trees, miguity, nd Chomsky norml form In this lecture we will discuss few importnt notions connected with contextfree grmmrs, including prse trees, miguity, nd specil form for context-free grmmrs
More informationList all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.
Mth Anlysis CP WS 4.X- Section 4.-4.4 Review Complete ech question without the use of grphing clcultor.. Compre the mening of the words: roots, zeros nd fctors.. Determine whether - is root of 0. Show
More informationFarey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University
U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationPKK Energy chains. PKK 320 / 200 x 2340 / 200; 5Pz, 1PT55
Energy cins Te order for plstic energy cin sould contin te following dt: Type / dius x Lengt / Widt "Arrngement"; Seprtor rrngement 1 3 Type selection is in ccordnce wit dimeter nd quntity of te lines
More informationChapter 2 Differentiation
Cpter Differentition. Introduction In its initil stges differentition is lrgely mtter of finding limiting vlues, wen te vribles ( δ ) pproces zero, nd to begin tis cpter few emples will be tken. Emple..:
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationChapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS
S Cpter Numericl Integrtion lso clled qudrture Te gol of numericl integrtion is to pproximte numericlly. f(x)dx Tis is useful for difficult integrls like sin(x) ; sin(x ); x + x 4 Or worse still for multiple-dimensionl
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationPLK VICWOOD K.T. CHONG SIXTH FORM COLLEGE Form Six AL Physics Optical instruments
AL Pysics/pticl instruments/p.1 PLK VICW K.T. CHNG SIXTH FRM CLLEGE Form Six AL Pysics pticl Instruments pticl instruments Mgniying glss Microscope Rercting telescope Grting spectrometer Qulittive understnding
More informationSection 4.7 Inverse Trigonometric Functions
Section 7 Inverse Trigonometric Functions 89 9 Domin: 0, q Rnge: -q, q Zeros t n, n nonnegtive integer 9 Domin: -q, 0 0, q Rnge: -q, q Zeros t, n non-zero integer Note: te gr lso suggests n te end-bevior
More informationMethod of Localisation and Controlled Ejection of Swarms of Likely Charged Particles
Method of Loclistion nd Controlled Ejection of Swrms of Likely Chrged Prticles I. N. Tukev July 3, 17 Astrct This work considers Coulom forces cting on chrged point prticle locted etween the two coxil,
More informationA Discussion on Formulas of Seismic Hydrodynamic Pressure
Interntionl Forum on Energy Environment Science nd Mterils (IFEESM 2017) A Discussion on Formuls of Seismic Hydrodynmic Pressure Liu Himing1 To Xixin2 1 Cin Mercnts Congqing Communiction Reserc & Design
More information1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation
1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationSpecial Relativity solved examples using an Electrical Analog Circuit
1-1-15 Specil Reltivity solved exmples using n Electricl Anlog Circuit Mourici Shchter mourici@gmil.com mourici@wll.co.il ISRAE, HOON 54-54855 Introduction In this pper, I develop simple nlog electricl
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationJURONG JUNIOR COLLEGE
JURONG JUNIOR COLLEGE 2010 JC1 H1 8866 hysics utoril : Dynmics Lerning Outcomes Sub-topic utoril Questions Newton's lws of motion 1 1 st Lw, b, e f 2 nd Lw, including drwing FBDs nd solving problems by
More information12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION
12. DYNAMIC ANALYSIS Force Equilibrium is Fundmentl in the Dynmic Anlysis of Structures 12.1 INTRODUCTION { XE "Newton's Second Lw" }All rel physicl structures behve dynmiclly when subjected to lods or
More information5: The Definite Integral
5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce
More information13.4 Work done by Constant Forces
13.4 Work done by Constnt Forces We will begin our discussion of the concept of work by nlyzing the motion of n object in one dimension cted on by constnt forces. Let s consider the following exmple: push
More informationANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS
ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS A. N. Anoshkin *, V. Yu. Zuiko, A.V.Glezmn Perm Ntionl Reserch Polytechnic University, 29, Komsomolski Ave., Perm, 614990, Russi
More informationAN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir
RGMIA Reserch Report Collection, Vol., No., 999 http://sci.vu.edu.u/ rgmi AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS I. Fedotov nd S. S. Drgomir Astrct. An
More informationSOLUTIONS TO CONCEPTS CHAPTER
1. m = kg S = 10m Let, ccelertion =, Initil velocity u = 0. S= ut + 1/ t 10 = ½ ( ) 10 = = 5 m/s orce: = = 5 = 10N (ns) SOLUIONS O CONCEPS CHPE 5 40000. u = 40 km/hr = = 11.11 m/s. 3600 m = 000 kg ; v
More informationStudy on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method
Avilble online t www.sciencedirect.com Physics Procedi 4 () 9 97 Interntionl Conference on Applied Physics nd Industril Engineering Study on the Clcultion of Mgnetic Force Bsed on the Equivlent Mgnetic
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description
More informationSection 2.1 Special Right Triangles
Se..1 Speil Rigt Tringles 49 Te --90 Tringle Setion.1 Speil Rigt Tringles Te --90 tringle (or just 0-60-90) is so nme euse of its ngle mesures. Te lengts of te sies, toug, ve very speifi pttern to tem
More informationNumerical Integration Problems
Integrtion COS 33 Nuericl Integrtion Proles Bsic D nuericl integrtion Given ility to evlute or ny, ind Gol: est ccurcy wit ewest sples Clssic prole even nlytic unctions not necessrily integrle in closed
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationCOSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature
COSC 336 Numericl Anlysis I Numericl Integrtion nd Dierentition III - Guss Qudrture nd Adptive Qudrture Edgr Griel Fll 5 COSC 336 Numericl Anlysis I Edgr Griel Summry o the lst lecture I For pproximting
More informationOn the application of explicit spatial filtering to the variables or fluxes of linear equations
Journl of Computtionl Physics 225 (27) 2 27 www.elsevier.com/locte/jcp Short Note On the ppliction of explicit sptil filtering to the vriles or fluxes of liner equtions Christophe Bogey *, Christophe Billy
More informationPatrice Cassagnard Université Montesquieu Bordeaux IV LAREefi. Abstract
A useful grpicl metod under Cournot competition Ptrice Cssgnrd Université Montesquieu Bordeux IV LAEefi Astrct Tis note proposes grpicl pproc useful in gme teor. Tis metod consists in representing incentives
More informationMath 20C Multivariable Calculus Lecture 5 1. Lines and planes. Equations of lines (Vector, parametric, and symmetric eqs.). Equations of lines
Mt 2C Multivrible Clculus Lecture 5 1 Lines nd plnes Slide 1 Equtions of lines (Vector, prmetric, nd symmetric eqs.). Equtions of plnes. Distnce from point to plne. Equtions of lines Slide 2 Definition
More informationChapter 4: Techniques of Circuit Analysis. Chapter 4: Techniques of Circuit Analysis
Chpter 4: Techniques of Circuit Anlysis Terminology Node-Voltge Method Introduction Dependent Sources Specil Cses Mesh-Current Method Introduction Dependent Sources Specil Cses Comprison of Methods Source
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationClassical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011
Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,
More informationSolution to Fredholm Fuzzy Integral Equations with Degenerate Kernel
Int. J. Contemp. Mth. Sciences, Vol. 6, 2011, no. 11, 535-543 Solution to Fredholm Fuzzy Integrl Equtions with Degenerte Kernel M. M. Shmivnd, A. Shhsvrn nd S. M. Tri Fculty of Science, Islmic Azd University
More information1 Nondeterministic Finite Automata
1 Nondeterministic Finite Automt Suppose in life, whenever you hd choice, you could try oth possiilities nd live your life. At the end, you would go ck nd choose the one tht worked out the est. Then you
More informationLecture 2: January 27
CS 684: Algorithmic Gme Theory Spring 217 Lecturer: Év Trdos Lecture 2: Jnury 27 Scrie: Alert Julius Liu 2.1 Logistics Scrie notes must e sumitted within 24 hours of the corresponding lecture for full
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGI OIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. escription
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationNew Fuzzy Numerical Methods for Solving Cauchy Problems
Article New Fuzzy Numericl Metods for Solving Cucy Problems Hussein ALKssbe, *, Irin Perfiliev, Mummd Zini Amd nd Zinor Ridzun Yy Institute of Engineering Mtemtics, Universiti Mlysi Perlis, Kmpus Tetp
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationAverage Rate of Change (AROC) The average rate of change of y over an interval is equal to change in
Averge Rte o Cnge AROC Te verge rte o cnge o y over n intervl is equl to b b y y cngein y cnge in. Emple: Find te verge rte o cnge o te unction wit rule 5 s cnges rom to 5. 4 4 6 5 4 0 0 5 5 5 5 & 4 5
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 17
CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking
More informationQuantum transport (Read Kittel, 8th ed., pp )
Quntum trnsport (Red Kittel, 8t ed., pp. 533-554) Wen we ve structure in wic mny collisions tke plce s crriers trnsport cross it, te quntum mecnicl pse of te electron wvefunctions is essentilly rndomized,
More informationThomas Whitham Sixth Form
Thoms Whithm Sith Form Pure Mthemtics Unit C Alger Trigonometry Geometry Clculus Vectors Trigonometry Compound ngle formule sin sin cos cos Pge A B sin Acos B cos Asin B A B sin Acos B cos Asin B A B cos
More informationRGMIA Research Report Collection, Vol. 1, No. 1, SOME OSTROWSKI TYPE INEQUALITIES FOR N-TIME DIFFERENTIA
ttp//sci.vut.edu.u/rgmi/reports.tml SOME OSTROWSKI TYPE INEQUALITIES FOR N-TIME DIFFERENTIABLE MAPPINGS AND APPLICATIONS P. CERONE, S.S. DRAGOMIR AND J. ROUMELIOTIS Astrct. Some generliztions of te Ostrowski
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More information1 Online Learning and Regret Minimization
2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationA5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s
4. Cosmic Dynmics: The Friedmnn Eqution Reding: Chpter 4 Newtonin Derivtion of the Friedmnn Eqution Consider n isolted sphere of rdius R s nd mss M s, in uniform, isotropic expnsion (Hubble flow). The
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationHYPERBOLA. AIEEE Syllabus. Total No. of questions in Ellipse are: Solved examples Level # Level # Level # 3..
HYPERBOLA AIEEE Sllus. Stndrd eqution nd definitions. Conjugte Hperol. Prmetric eqution of te Hperol. Position of point P(, ) wit respect to Hperol 5. Line nd Hperol 6. Eqution of te Tngent Totl No. of
More informationTopics Covered AP Calculus AB
Topics Covered AP Clculus AB ) Elementry Functions ) Properties of Functions i) A function f is defined s set of ll ordered pirs (, y), such tht for ech element, there corresponds ectly one element y.
More informationLinear Momentum and Collisions
chpter 9 Liner Momentum nd Collisions 9.1 Liner Momentum 9. Anlysis Model: Isolted ystem (Momentum) 9.3 Anlysis Model: Nonisolted ystem (Momentum) 9.4 Collisions in One Dimension 9.5 Collisions in Two
More informationHomework Solution - Set 5 Due: Friday 10/03/08
CE 96 Introduction to the Theory of Computtion ll 2008 Homework olution - et 5 Due: ridy 10/0/08 1. Textook, Pge 86, Exercise 1.21. () 1 2 Add new strt stte nd finl stte. Mke originl finl stte non-finl.
More informationSOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS
ELM Numericl Anlysis Dr Muhrrem Mercimek SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Some of the contents re dopted from Lurene V. Fusett, Applied Numericl Anlysis using MATLAB.
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationLinear Systems with Constant Coefficients
Liner Systems with Constnt Coefficients 4-3-05 Here is system of n differentil equtions in n unknowns: x x + + n x n, x x + + n x n, x n n x + + nn x n This is constnt coefficient liner homogeneous system
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationPatch Antennas. Chapter Resonant Cavity Analysis
Chpter 4 Ptch Antenns A ptch ntenn is low-profile ntenn consisting of metl lyer over dielectric sustrte nd ground plne. Typiclly, ptch ntenn is fed y microstrip trnsmission line, ut other feed lines such
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More informationMathematics. Area under Curve.
Mthemtics Are under Curve www.testprepkrt.com Tle of Content 1. Introduction.. Procedure of Curve Sketching. 3. Sketching of Some common Curves. 4. Are of Bounded Regions. 5. Sign convention for finding
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationu( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.
nlyzing Dmped Oscilltions Prolem (Medor, exmple 2-18, pp 44-48): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $
More informationDETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE
Determintion RevAdvMterSci of mechnicl 0(009) -7 properties of nnostructures with complex crystl lttice using DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationPolynomials and Division Theory
Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the
More informationStudy Guide Final Exam. Part A: Kinetic Theory, First Law of Thermodynamics, Heat Engines
Msschusetts Institute of Technology Deprtment of Physics 8.0T Fll 004 Study Guide Finl Exm The finl exm will consist of two sections. Section : multiple choice concept questions. There my be few concept
More informationCalculus AB Section I Part A A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION
lculus Section I Prt LULTOR MY NOT US ON THIS PRT OF TH XMINTION In this test: Unless otherwise specified, the domin of function f is ssumed to e the set of ll rel numers for which f () is rel numer..
More informationSTEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.
STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t
More information2. General formula for Runge-Kutta methods
NOLTA, IEICE Pper Equivlent circuits for implicit unge-kutt metods in circuit simultors for nonliner circuits Ysuiko Toym 1), Ttusy Kuwzki 1, nd Jun Sirtki 1 Mkiko Okumur 1 1 Kngw Institute of Tecnology,
More informationThe Periodically Forced Harmonic Oscillator
The Periodiclly Forced Hrmonic Oscilltor S. F. Ellermeyer Kennesw Stte University July 15, 003 Abstrct We study the differentil eqution dt + pdy + qy = A cos (t θ) dt which models periodiclly forced hrmonic
More information