Finite element method for structural dynamic and stability analyses
|
|
- Barnard Davidson
- 5 years ago
- Views:
Transcription
1 Finie elemen mehod for srucural dynamic and sabiliy analyses Module- Nonlinear FE Models Lecure-39 Toal and updaed Lagrangian formulaions Prof C Manohar Deparmen of Civil Engineering IIc, Bangalore 56 India
2 Linear analysis Displacemens and srains are small ress-srain relaions are linear Equillibrium equaions are derived based on undeformed geomery Principle of superposiion holds Nonlinear analysis Displacemen and srains need no be small; special effors needed o characerize roaions The geomery of he obec, sress-srain relaions, and BCs could change during he process of deformaion. The definiions of sress and srain and he formulaion of he governing equaions need o ake ino accoun hese changes. Necessiaes inroducion of newer measures of sress and srain Principle of superposiion does no hold
3 rain measures Infiniesimal srains: for body under rigid body roaions, he srains would no be zero. New measures needed: Rigid body moions imply zero srains For small srains, he infiniesimal srain definiions are o be resored. 3
4 Green - Lagrange srain measure dx FdX ds d dx dx dx dx dx F FdX dx dx dx F F I dx dx EdX E F F I Green-Lagrange srain measure ds MF E N N & AB cos E M N d Almansi - Hamel (Eulerian) srain dx F dx ds d dx dx dx dx dx dx dx F F dx dx I F F dx dx edx e I F F 4
5 rain measure Engineering Green-Lagrange L L L Lagarihmic L L ln L Exensional raio L L Engineering Green-Lagrange logarhmic Exension raio L L MF L L/L) 5
6 ress measures Cauchy sress ensor: defined wih respec o deformed geomery. This would no be known in advance. Two alernaives: ress as a measure which conugaes wih a measure of srain o produce inernal energy As a quaniy which produces a racion vecor in conuncion wih a normal vecor defined wih respec o a surface elemen 6
7 Cauchy - Euler sress ress a P : n T n n lim a da da JF da N f a Firs Piola - Kirchoff sress (P) Based on deformed configuraion Based on undeformed configuraion sress vecor acing on elemen da wih ouward normal N which produces force = df. T n N df da T da N PN PdA da J F da P J F P is no symmeric; i has 9 independen componens econd Piola - Kirchoff sress () Based on undeformed configuraion Inroduce a pseudo-force vecor dpˆ F df ˆ n dp F df F da F nda JF nf da nda JF F econd Piola-Kirchoff sress ensor is symmeric & JF F PF 7
8 Noaions for configuraions and deformaions x, x, x x, x, x C Q ds P u u u Q x, x, x ds P C Q ds P C Configuraion Coordinaes of a poin olumes Areas Densiy Toal displacemen i Q iq quaniy Q is measured in C wih respec o is value in C i???? Check Noaionson u C C C x x x A A A u u u J N Reddy, 4, An inroducion o nonlinear FEA, Oxford Univ Press, NY, 4 K J Bahe, 996, Finie elemen procedures, Prenice Hall India, New Delhi. 8
9 Deformaion: A paricle in C a x x x x 3 moves o Moion from C o C : u x x x x x x3 ui xi xi i C a ;,,3 and o Moion from C o C : u x x C a x x x x3 ui xi xi ; i,,3 Conservaion of mass : xi xi d Jd d d d x x x x x x3 x x x d Jd wih J F x x x3 J & J x x x x x x3 9
10 rain ensors for C and C configuraions Green Lagrange srain ensor ui uk uk Ei x xi xi x u u u i uk u k Ei x xi xi x Green Lagrange incremenal srain ensor id xid x ds ds ds ds ds ds Ei Ei d xid x ei i d xid x Ei Ei ei i Ei Ei ei i where i
11 d x d x E E d x d x i i i i i i i i i k k k k ei x xi xi x xi x k i xi u uk (nonlinear in incremen u i ) x e d x d x u u u u u u (linear in incremen u ) i
12 Updaed Green Lagrange srain ensor u u i uk u k Ei is helpful in oal Lagrangian formulaions. x xi xi x To faciliae updaed Lagrangian formulaion we inroduce i i d xd x ds ds updaed Green-Lagrange srain ensor. Using he following resuls i k k i i i x xi x k k i i i & ki xi xi we ge i i x x x d x F d x d x ds d x d x u x x wih e x ui xi u u uk u k e xi xi x i u u u u & i k k i i xi xi xi x i,
13 Euler srain ensor Consider ha he body has reached configuraion C from C in several incremens. The deformaion from C o C could be large. We now wish o move o configuraion C. The incremen from C o C is aken o be small and we could refer o srains wih respec o C. id xid x ds ds k k id xid x ds ds i d x id x xi x x k i i xi We have x x k x u x x x u k k k k k k k k x x x u u u u i k k i x xi xi x x x Euler srain ensor u 3
14 u u u u i k k i x xi xi x The linear par of he Euler srain ensor is given by u u i ei ei x x i This quaniy is called he infiniesimal srain ensor. These srain componens conugae wih Cauchy sress componens o produce he expressions for inernal energy sored. i 4
15 ress ensors Cauchy sress: inernal force reckoned in he deformed configuraion and area reckoned in he deformed configuraion. Configuraion C : Configuraion C : i i i i econd Piola - Kirchoff sress ensor: force in C ransformed o C and area reckoned in C x ˆ n d A d f F d f d f x nˆ uni normal o d A in C. Updaed Kirchoff sress ensor useful in updaed Lagrangian formulaions Consider he poin P x, x, x in C. 3 Cauchy sress ensor a P in C is denoed by i 5
16 i i i Kirchoff sress incremen ensor i i i C C C i i i i i Updaed Kirchoff ress 6
17 i i i Updaed Kirchoff ress i i i ince Kirchoff sress incremen ensor x x i i i J mn xm xn x i i J mn x m x n i x J, we also ge x x i i xm xn x x i i x m x n mn mn Relaion beween Cauchy sress in C and updaed Kirchoff sress xi i x p x q x x i i pq x p x q PK- sress in differen configuraions i x x x i x p x q x pq Relaions beween incremenal sresses xi i x p x q x x i i x p x q x i i x p x pq pq pq pq x q
18 Consiuive relaions Aenion is limied o linear relaions beween conugae sress-srain pairs. Maerial behaviour is elasic: consiuive behaviour is funcion of curren sae of deformaion. Relaion beween second Piola - Kirchoff sress and Green - Lagrange srain C ikl E kl i where C maerial elasiciy ensor ress - srain relaions in incremenal form Kirchoff sress incremen C i ikl kl Updaed Kirchoff sress incremen C i ikl kl C in differen configuraions x i & Green-Lagrange srain incremen i & Green-Lagrange srain incremen i xi xk xl xi xk xl Cikl C ikl & Cikl Cikl xp xq xr xs xp xq xr xs x i 8
19 Principle of virual displacemens um of virual exernal work done on a body and he virual work sored in he body should be zero. Consider configuraion C W : e d f u d u d W i ei d fi ui d i ui d 9
20 Remarks We canno use his equaion direcly since he descripion of configuraion C would no be known. This descripion keeps changing as he deformaion evolves. The assumpion made in he linear analysis ha he configuraion of he body does no change so ha he equaions can be formulaed based on undeformed geomery is no valid in nonlinear analysis. This calls for inroducion of measures of sress and srain which ake ino accoun he changes in configuraion during deformaion. This enables he evaluaion of inegrals in he expression for he inernal work done over known configuraions. We will use ress: he second Piola Kirchoff sress ensor rain: he Green-Lagrange srain ensor
21 Toal Lagrangian approach Curren configuraion Base & reference configuraions coincide (Taken o remain fixed)
22 Updaed Lagrangian approach Base configuraion Curren configuraion Reference configuraion updaed a each incremen while solving he equilibrium equaion
23 Toal Lagrangian approach All quaniies reckoned wih respec o undeformed configuraion (C ). i i i i i i i i i i i i i e d E d f u d f u d u d u d E i d f u d u d i i i i 3
24 E d f u d u d i i i i i i Recall: Consider Recall: i i i E i Ei Ei ei i Ei ei Kirchoff sress incremen ensor E E e i i i i E i i E does no depend on unknown displacemens i 4
25 Consider Recall: i k k k k ei x xi xi x xi x i k k k k ei x xi xi x xi x Consider Recall: e i i u uk x k i xi u u u u u u u u u u u u k k k k i xi x xi x u u u u 5
26 i i i i i i i i i i i E d R Consider he erm R f u d u d E d d i i i ei i d e d i i This represens he virual inernal energy sored in he body in configuraion C. By applying virual work principle o body in configuraion C, one ges R e d f u d u d i i i i i i 6
27 i i d i i d R R i i d d i i Change in virual srain energy due o virual incremenal displacemens u beween configuraions C and C. i irual work done by forces due o iniial sress. i This arises due o change in geomery beween he wo configuraions. Consiuive relaion C i ikl kl ikl kl i i i C d d R R
28 i i i i Assuming he displacemens i ikl kl e i Cikl ekl ei d i i d R R o be small reasonable assumpion for small load incremens d d R R C e C E i ikl kl i u i 8
29 Ineria forces i i i i i i u u d u u d u u Equaion of moion i i irs rs i i i R R u u d C e e d d 9
30 Toal Lagrangian formulaion for a D elemen All quaniies measured wr C x, x, x C Q ds P u u u Q ds P Q ds P x, x, x C C x x u u u u x y u v u v ikl kl i i i C e e d d R R
31 Weak form ikl kl i i i e C e e d d R R u u u v v x x x x x e xx v u u v v eyy y y y y y e xy u v u u v v u u v v y x x y x y y x y x u v x x x x x u v u D y y y y y v u u v v y x y x x y y x x y Du u 3
32 e D D u ikl kl i u u D Du C D Du d C C C orhoropically elasic u C e e d e C e d C C C 66 3
33 i i d d u u v v x x x x xx u u v v yy y y y y xy u u v v x y x y u v x x x x u v u y y y y v u u v v y x x y y x x y D u u 33
34 D u D u u u v v x x x x xx u u v v y y y y xy u u u u v v v v x y x y x y x y yy u u C C E xx E yy E xy xx yy C C xy C66 u u v x x x E xx v u v Eyy y y y E xy u v u v u v y x x y y x E 34
35 The resuling forma of i i d d is no sill in a form suiable for developmen of FE soluion. Noe ha he erms are nonlinear in he incremenal displacemen u. We re-arrange he erms as follows: [ xx d u u v v x x x x xx u u v v yy d y y y y xy u u u u v v v v x y x y x y x y u u v v x x x x u u v v yy y y y y u u u u v v v v x y x y x y x y xy ] d 35
36 u xx xy u y xy yy y d v xx xy v x x xy yy u D D u d u u x x v v y y x xx xy xy yy y u & D xx xy v x xy yy y We are now ready o launch he FE formulaion. d 36
37 Incremenal displacemen: Toal displacemen: u u n n u v u v u v u v u v u n u x u n v v x n u x u n v v x n n n v n 37
38 u u u u BL C BL d BL D Du u u NL NL D Du C D Du d D D C D D d wih. e C e d d D D d D D d B B d wih BNL D 38
39 i i R e d e d B L i i i i f d d R f u d u d u d i 39
40 ikl kl i i i C e e d d R R KL KNL F F KL BL C BL d K NL BNL BNL d F f d i x x fy y F BL d ; f ; u d f Remarks L NL K K K K & C are symmeric This is an incremenal formulaion. The siffness marix here is he angen siffness marix. F K NL For a linear analysis,,,& 4
41 u v BL B L B L B L n x x x B L n x x x u u u n u x x x x x x B L u u u n x x x x x x v v v n v x x x x x B x L v v v n x x x x x x 4
42 B NL n x x x n y y y n x x x n y y y KL KNL F F L N L K K K u F F L L N K K K v F F 4
43 Updaed Lagrangian approach All quaniies reckoned wih respec o he laes known configuraion (C ). i i i i i i i i i i i i i e d d f u d f u d u d u d updaed Green-Lagrange srain ensor body force referred o in C. i surface racion referred o in C. i f i d R i i i i i R f u d u d 43
44 We have e i k k i i i xi xi xi x i i xi u u u u e k k k k i xi x xi x i Consider d R Recall u i u x i i i i Updaed Kirchoff ress u u u u i i i i i d R d R d e d R i i i i i 44
45 i i i i d R d e d R i i i i i i i i i Considering he equillibrium of body in C d d R R R e d R f u d u d i i i We have he consiuive relaion ikl kl i i i i i i ikl kl C d d R e d As we did in he oal Lagrangian approach we ake C e & e i ikl kl i i i C 45
46 This leads o he required weak form ikl kl i i i i i The Cauchy sress kl C e e d d R e d Almansi srain componens d x d x i i are evaluaed using C i i ikl kl ds ds ui uk uk i x xi xi x u 46
47 FE model for D elemen based on updaed Lagrangian formulaion K K C e e d d R e d ikl kl i i i i i L NL F F KL B L C B L d K NL BNL BNL d F f d x F B L d ; f ; d f x fy y 47
48 n x x x B L n x x x B NL n x x x n y y y n x x x n y y y 48
49 C C xx yy xy xx yy C C xy C66 u u v x x x xx v u v yy y y y xy u v u v u v y x x y y x 49
50 Follow up Maerial nonlineariy abiliy analysis: inclusion of nonlineariy a differen levels Hybrid esing Bayesian filering Uncerainy modelling and fem Thermal loads: fire Anisoropy 5
Finite Element Analysis of Structures
KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using
More informationFinite Element Formulation for Large Deformation Problem
Finie Elemen Formulaion for Large Deformaion Problem Sushan Verma 1*, Sandeep Kumar Singh 1, Anuj Gupa 1, Anwer Ahmed 1 1 Mechanical Engineering Deparmen, G.L Bajaj Insiue of Technology and Managemen,
More informationLecture 10: Wave equation, solution by spherical means
Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationFinite Element Formulation. for Large Displacement Analysis
Finie Elemen Formulaion for Large Displacemen Analysis Mario Gomez Boero Deparmen of Mechanical Engineering EAFIT Universiy March 16, 009 Conens 1 Inroducion o he Variaional Formulaion 9 1.1 Minimum Toal
More informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationCH.5. BALANCE PRINCIPLES. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.5. BALANCE PRINCIPLES Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Balance Principles Convecive Flux or Flux by Mass Transpor Local and Maerial Derivaive of a olume Inegral Conservaion
More informationv A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3
The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped
More informationIntroduction to continuum mechanics-ii
Ioan R. IONESCU LSPM, Universiy Paris 3 (Sorbonne-Paris-Cié) July 3, 27 Coninuum kinemaics B uneforme/reference configuraion, B eforme/acual configuraion Moion: x(, X ) = X + u(, X ) wih x(, B ) = B. Maerial/Lagrangian
More information4.1.1 Mindlin plates: Bending theory and variational formulation
Chaper 4 soropic fla shell elemens n his chaper, fia shell elemens are formulaed hrough he assembly of membrane and plae elemens. The exac soluion of a shell approximaed by fia faces compared o he exac
More informationFinite element method for structural dynamic and stability analyses
Finie elemen mehod for srucural dynamic and sabiliy analyses Modules- & 3 Finie elemen analysis of dynamics of planar russes and frames. Analysis of equaions of moion. Lecure-7: FE modeling of planar srucures:
More informationSummary of shear rate kinematics (part 1)
InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial
More informationFrom Particles to Rigid Bodies
Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and
More informationANALYSIS OF REINFORCED CONCRETE BUILDINGS IN FIRE
ANALYSIS OF REINFORCED CONCRETE BUILDINGS IN FIRE Dr Zhaohui Huang Universiy of Sheffield 6 May 2005 1 VULCAN layered slab elemens: connecion o beam elemens Plae Elemen Slab nodes y x Reference Plane h
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationG. =, etc.
Maerial Models υ υ3 0 0 0 υ υ 3 0 0 0 υ3 υ3 0 0 0 = 0 0 0 0 0 0 0 0 0 0 3 0 0 0 0 0 3 l (9..4) he subscris denoe he maerial axes, i.e., υ = υ and = (9..5) i j xi xj ii xi Since l is symmeric υ υ =, ec.
More informationStochastic Structural Dynamics. Lecture-12
Sochasic Srucural Dynamics Lecure-1 Random vibraions of sdof sysems-4 Dr C S Manohar Deparmen of Civil Engineering Professor of Srucural Engineering Indian Insiue of Science Bangalore 56 1 India manohar@civil.iisc.erne.in
More informationOptimal Path Planning for Flexible Redundant Robot Manipulators
25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering
More informationHomework Set 2 Physics 319 Classical Mechanics
Homewor Se Physics 19 Classical Mechanics Problem.7 a) The roce velociy equaion (no graviy) is m v v ln m Afer wo minues he velociy is m/sec ln = 79 m/sec. b) The rae a which mass is ejeced is ( 1 6-1
More informationTwo- and Three Dimensional Solid Elements; Plane Stress, Plane Strain, and Axisymmetric Conditions
Topic 7 Two- and Three Dimensional Solid Elemens; Plane Sress, Plane Srain, and Axisymmeric Condiions Conens: soparameric inerpolaions of coordinaes and displacemens Consisency beween coordinae and displacemen
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationTR/06/83 September 1983 LOCAL AND GLOBAL BIFURCATION PHENOMENA IN PLANE STRAIN FINITE ELASTICITY R. W. OGDEN
TR/06/8 Sepember 98 LOCAL AND GLOBAL BIFURCATION PHENOMENA IN PLANE STRAIN FINITE ELASTICITY by R W OGDEN w95987x Local and global bifurcaion phenomena in plane srain finie elasiciy By RW OGDEN Deparmen
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationCH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear
More informationProgram: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member
Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam
More informationCh1: Introduction and Review
//6 Ch: Inroducion and Review. Soli and flui; Coninuum hypohesis; Transpor phenomena (i) Solid vs. Fluid No exernal force : An elemen of solid has a preferred shape; fluid does no. Under he acion of a
More informationMECHANICS OF MATERIALS Poisson s Ratio
Poisson s Raio For a slender bar subjeced o axial loading: ε x x y 0 The elongaion in he x-direcion i is accompanied by a conracion in he oher direcions. Assuming ha he maerial is isoropic (no direcional
More informationShells with membrane behavior
Chaper 3 Shells wih membrane behavior In he presen Chaper he sress saic response of membrane shells will be addressed. In Secion 3.1 an inroducory example emphasizing he difference beween bending and membrane
More informationNEW ALGORITHMS FOR THE SIMULATION OF THE SEQUENTIAL TUNNEL EXCAVATION WITH THE BOUNDARY ELEMENT METHOD (EURO:TUN 2007)
ECCOMAS Themaic Conference on Compuaional Mehods in Tunnelling (EURO:TUN 2007) J. Eberhardseiner e.al. (eds.) Vienna, Ausria, Augus 27-29, 2007 NEW ALGORITHMS FOR THE SIMULATION OF THE SEQUENTIAL TUNNEL
More informationNew effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification
IOP Conference Series: Maerials Science and Engineering PAPE OPEN ACCESS New effecive moduli of isoropic viscoelasic composies. Par I. Theoreical jusificaion To cie his aricle: A A Sveashkov and A A akurov
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationA Fixed Grid method for hyperelastic models in large strain analysis using the Level Set Method
A Fixed Grid mehod for hyperelasic models in large srain analysis using he Level Se Mehod William Ramírez Beniez Manuel J Garcia. Grupo de Invesigación Mecánica Aplicada, universidad EAFIT, Medellín, Colombia
More informationPrediction of Concrete Fracture Mechanics Behavior and Size Effect using Cohesive Zone Modeling
Predicion of Concree Fracure Mechanics Behavior and Size Effec using Cohesive Zone Modeling Kyoungsoo Park, Glaucio H. Paulino, Jeffery R. Roesler Deparmen of Civil and Environmenal Engineering Universiy
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More information2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance
Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion
More informationDynamic Analysis of Loads Moving Over Structures
h Inernaional ongress of roaian ociey of echanics epember, 18-, 3 Bizovac, roaia ynamic nalysis of Loads oving Over rucures Ivica Kožar, Ivana Šimac Keywords: moving load, direc acceleraion mehod 1. Inroducion
More informationBasilio Bona ROBOTICA 03CFIOR 1
Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationFlow-Induced Vibration Analysis of Supported Pipes with a Crack
Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall
More informationψ(t) = V x (0)V x (t)
.93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in
More informationFinite element method for structural dynamic and stability analyses
Finie elemen mehod for srucural dynamic and sabiliy analyses Module-1 Closure Lecure-40 Closure and a few opics for furher sudy Prof C S Manohar Deparmen of Civil Engineering IISc, Bangalore 560 01 India
More informationAnalytic nonlinear elasto-viscosity of two types of BN and PI rubbers at large deformations
Bulgarian Chemical Communicaions, Volume 48, Special Issue E (pp. 59-64) 016 Analyic nonlinear elaso-viscosiy of wo ypes of BN and PI rubbers a large deformaions K. B. Hadjov, A. S. Aleksandrov, M. P.
More informationa. Show that these lines intersect by finding the point of intersection. b. Find an equation for the plane containing these lines.
Mah A Final Eam Problems for onsideraion. Show all work for credi. Be sure o show wha you know. Given poins A(,,, B(,,, (,, 4 and (,,, find he volume of he parallelepiped wih adjacen edges AB, A, and A.
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationElements of Computer Graphics
CS580: Compuer Graphics Min H. Kim KAIST School of Compuing Elemens of Compuer Graphics Geomery Maerial model Ligh Rendering Virual phoography 2 Foundaions of Compuer Graphics A PINHOLE CAMERA IN 3D 3
More informationProfessor Terje Haukaas University of British Columbia, Vancouver terje.civil.ubc.ca. St. Venant Torsion
S. Venan Torsion Torque in srucural members is carried by S. Venan orsion, and possible warping orsion. In S. Venan orsion he orque is carried by shear sresses; in warping orsion he orque is carried by
More informationLie Derivatives operator vector field flow push back Lie derivative of
Lie Derivaives The Lie derivaive is a mehod of compuing he direcional derivaive of a vecor field wih respec o anoher vecor field We already know how o make sense of a direcional derivaive of real valued
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationCh.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationDiffusion & Viscosity: Navier-Stokes Equation
4/5/018 Diffusion & Viscosiy: Navier-Sokes Equaion 1 4/5/018 Diffusion Equaion Imagine a quaniy C(x,) represening a local propery in a fluid, eg. - hermal energy densiy - concenraion of a polluan - densiy
More informationOn the applicability and stress update formulations for corotational stress rate hypoelasticity constitutive models
Available online a www.sciencedirec.com Finie Elemens in Analysis and Design 39 (003) 783 816 www.elsevier.com/locae/nel On he applicabiliy and sress updae formulaions for coroaional sress rae hypoelasiciy
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationStochastic Structural Dynamics. Lecture-6
Sochasic Srucural Dynamics Lecure-6 Random processes- Dr C S Manohar Deparmen of Civil Engineering Professor of Srucural Engineering Indian Insiue of Science Bangalore 560 0 India manohar@civil.iisc.erne.in
More informationThe motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationEight-Node Solid Element for Thick Shell Simulations
Eigh-Node Solid Elemen for Thick Shell Simulaions Yong Guo Livermore Sofware Technology Corporaion Livermore, California, USA Absrac An eigh-node hexahedral solid elemen is incorporaed ino LS-DYNA o simulae
More informationNon-uniform circular motion *
OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform
More informationComparison between the Discrete and Continuous Time Models
Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationNumerical Analysis of Cable Structures
Paper 240 Numerical Analysis of Cable Srucures Civil-Comp Press, 2012 Proceedings of he Elevenh Inernaional Conference on Compuaional Srucures Technology, B.H.V. Topping, (Edior), Civil-Comp Press, Sirlingshire,
More information2 Some Property of Exponential Map of Matrix
Soluion Se for Exercise Session No8 Course: Mahemaical Aspecs of Symmeries in Physics, ICFP Maser Program for M 22nd, January 205, a Room 235A Lecure by Amir-Kian Kashani-Poor email: kashani@lpensfr Exercise
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationth World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3256
11111 1 h World Conference on Earhquake Engineering Vancouver, B.C., Canada Augus 1-6, 2004 Paper No. 256 TOUCHING ANALYSIS OF TWO BUILDINGS USING FINITE ELEMENT METHOD Mircea IEREMIA 1, Silviu GINJU 1,
More informationFinite element method for structural dynamic and stability analyses
Finie elemen mehod for srucural dynamic and sabiliy analyses Module-3 Analysis of equaions of moion Lecure-8: FRF-s and Damping models Prof C S Manohar Deparmen of Civil Engineering IISc, Bangalore 56
More informationAP Calculus BC Chapter 10 Part 1 AP Exam Problems
AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a
More informationMechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,
Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas
More informationTests of Nonlinear Resonse Theory. We compare the results of direct NEMD simulation against Kawasaki and TTCF for 2- particle colour conductivity.
ess of Nonlinear Resonse heory We compare he resuls of direc NEMD simulaion agains Kawasaki and CF for 2- paricle colour conduciviy. 0.0100 0.0080 0.0060 Direc CF, BK and RK J x 0.0040 0.0020 DIR KAW CF
More informationThe mathematical shell model underlying general shell elements
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING In. J. Numer. Meh. Engng. 2000; 48:289 313 The mahemaical shell model underlying general shell elemens Dominique Chapelle 1; ; and Klaus-Jurgen
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationModified one dimensional elastic wave equations that retain time synchronization under spatial coordinate transformations
Modified one dimensional elasic wave equaions ha reain ime synchronizaion under spaial coordinae ransformaions Ruiwen Yao, Zhihai Xiang * AML, Deparmen of Engineering Mechanics, Tsinghua Universiy, Beijing
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationMath 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm
Mah Spring 6 Differenial Equaions Final Exam Due Dae: Tuesday, June 6, 5pm Your name (please prin): Insrucions: This is an open book, open noes exam. You are free o use a calculaor or compuer o check your
More informationCash Flow Valuation Mode Lin Discrete Time
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationEE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:
EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he
More information2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS
Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described
More informationTurbulent Flows. Computational Modelling of Turbulent Flows. Overview. Turbulent Eddies and Scales
School of Mechanical Aerospace and Civil Engineering Turbulen Flows As noed above, using he mehods described in earlier lecures, he Navier-Sokes equaions can be discreized and solved numerically on complex
More informationMTH Feburary 2012 Final term PAPER SOLVED TODAY s Paper
MTH401 7 Feburary 01 Final erm PAPER SOLVED TODAY s Paper Toal Quesion: 5 Mcqz: 40 Subjecive quesion: 1 4 q of 5 marks 4 q of 3 marks 4 q of marks Guidelines: Prepare his file as I included all pas papers
More informationPolymer Engineering (MM3POE)
Polymer Engineering (MM3POE) VISCOELASTICITY hp://www.noingham.ac.uk/~eazacl/mm3poe Viscoelasiciy 1 Conens Wha is viscoelasiciy? Fundamenals Creep & creep recovery Sress relaxaion Modelling viscoelasic
More informationPhysics for Scientists & Engineers 2
Direc Curren Physics for Scieniss & Engineers 2 Spring Semeser 2005 Lecure 16 This week we will sudy charges in moion Elecric charge moving from one region o anoher is called elecric curren Curren is all
More information( ) = Q 0. ( ) R = R dq. ( t) = I t
ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationL p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity
ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.
More informationComputer-Aided Analysis of Electronic Circuits Course Notes 3
Gheorghe Asachi Technical Universiy of Iasi Faculy of Elecronics, Telecommunicaions and Informaion Technologies Compuer-Aided Analysis of Elecronic Circuis Course Noes 3 Bachelor: Telecommunicaion Technologies
More informationFormulation of the Stress Distribution Due to a Concentrated Force Acting on the Boundary of Viscoelastic Half-Space
Formulaion of he Sress Disribuion Due o a Concenraed Force Acing on he Boundar of Viscoelasic Half-Space Yun eng and Debao Zhou Deparmen of Mechanical and Indusrial Engineering Universi of Minnesoa, Duluh
More informationParameter identification of computational fracture models
Parameer idenificaion of compuaional fracure models C. Iacono & L. J. Sluys Faculy of Civil Engineering and Geosciences, Delf Universiy of Technology. Delf, The Neherlands. J. G. M. van Mier ETH Hönggerberg,
More informationCombined Bending with Induced or Applied Torsion of FRP I-Section Beams
Combined Bending wih Induced or Applied Torsion of FRP I-Secion Beams MOJTABA B. SIRJANI School of Science and Technology Norfolk Sae Universiy Norfolk, Virginia 34504 USA sirjani@nsu.edu STEA B. BONDI
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationEECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits
EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive
More informationKinematics and kinematic functions
ROBOTICS 01PEEQW Basilio Bona DAUIN Poliecnico di Torino Kinemaic funcions Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions(called kinemaic funcions or KFs) ha mahemaically
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
More information