F8. Prototype models for plastic response; integration of plastic flow; structural analys of beam problem.
|
|
- Britney Lee
- 5 years ago
- Views:
Transcription
1 F8. Prototype models for plastic response; integration of plastic flow; structural analys of beam problem. Prototype model: Perfect plasticity Consider model beavior: Figure 1 stablis reological model: spring combined wit slider! Figure 2 Free energy formulated in internal plastic strain ε p fl κ =ε p ψ@ε, ε p ; θd = 1 2 Hε εp αθl 2 σ= ψ ε = Hε εp αθl, K = ψ ε p =σ Inelastic deformation in slider spontaneously as stress exceeds yield stress fl rate independent beavior Plastic micro stress K same as continuum stress σ, cf. Maxwell model for linear visco-elasticity! stablis reduced dissipation
2 2 Notes_F8.nb = K ε p 0OK proper formulation of ε p =? stablis evolution of plastic deformation in deformation in slider: Introduce: yield function: φ@σd =» σ» σ y, σ y = yield stress Yield criterion: φ@σd =» σ» σ y = 0 Figure 3 Postulate evolution of deformation in slider as il ε p = 0ifφ<0 Helastic loadingl iil ε p =λ sign@σd wit λ 0 if φ=0 Hplastic loadingl Reconsider dissipation Figure 4 il = K ε p def = 0OK! iil = K ε p =σsign@σd λ =» σ» λ 0 OK since λ 0 Plastic multiplier λ 0 φ =»σ» σ = sign@σd=gradient of yield function fl ε p =λ φ wen φ@σd = 0 Plastic evolution ε p =λ φ controlled by loading (or Kun-Tucker) conditions φ 0, λ 0,φλ = 0 ü Consequences of loading conditions - Continuum tangent relation Assume α θ=0 fl σ=hε ε p L σ = Hε ε p L Distinguis 1L φ@σd < 0 σ = ε = σ e 2L φ@σd = 0 σ = Hε ε p L = σ e Conclusion: lastic loading HL λ sign@σd Plastic loading HPL
3 Notes_F8.nb 3 σ = σ e λ sign@σd φ 0, λ 0,φλ = 0 = Consider tangent stiffness beavior! Assume φ@σd = 0, and increase loading so tat φ@σ +σ dtd =φ@σd dt dt 0 Beavior controlled by linearized loading conditions, i.e. σ = σ e λ sign@σd = φ 0, λ 0,φ λ = 0 In tis case, we ave φ = sign@σd σ = sign@σd Hσ e σ = σ e λ sign@σd =...= 0 HPL Plastic loading (P) for φ@σd = 0 occurs if λ > 0 for ε > 0 if sign@σd = 1 λ > 0 for ε < 0 if sign@σd = 1 Stress state remains on te yield surface! λ sign@σdl = 0 λ = sign@σd σ e = sign@σd ε Prototype model: Isotropic ardening plasticity Consider model beavior Figure 5 stablis reological model: spring combined wit slider+spring!
4 4 Notes_F8.nb Figure 6 Free energy formulated in internal plastic strain ε p and internal ardening variable κ fl ψ@ε, ε p, κ; θd = 1 2 Hε εp αθl H κ2 σ= ψ ε = Hε εp αθl Reduced dissipation wit =σ p ε p + K κ =σε p + K κ 0OK proper formulation of ε p =? σ p =σ= ψ ε p = Hε εp αθl, K = ψ κ = H κ Inelastic deformation in slider spontaneously as stress exceeds yield stress fl rate independent beavior Plastic micro-stress K σ! xtra contribution K κ in dissipation due to micro stress K in slider! Caracterization fl H > 0 Hardening material beavior H = 0 Perfectly plastic material beavior H < 0 softening material beavior Hassociated wit fracturel Formulate evolution of plastic deformation in slider: Introduce: yield function: φ@σ, KD :=» σ» K σ y wit σ y = initial yield stress Postulate evolution ε p for lastic loading () and plastic loading (P): HL ε p = 0ifφ@σ, KD < 0 HPL ε p =λ φ =λ sign@σd κ =λ φ K = λ = wit λ 0 for φ=0 stablis dissipation il =σ p ε p + K κ =σε p H κ κ def = 0OK!
5 Notes_F8.nb 5 iil = Hσ sign@σd KL λ = H»σ» KL λ =σ y λ 0 OK since λ 0 Plastic multiplier λ 0 φ =»σ» σ = sign@σd=gradient of yield function fl ε p =λ φ wen φ@σd = 0 Plastic evolution ε p =λ φ φ 0, λ 0,φλ = 0 controlled by loading (or Kun-Tucker) conditions Prototype model: Non-linear creep (Norton's creep law) Solution procedures - Structural analysis à Numerical integration flow rules - Backward uler Metod ü Isotropic ardening plasticity model Consider time interval I td subdivided into a finite number of timesteps mtim, i.e I = 80, t 1,..., t n, t n+1 < n=1,mtim Apply B finite differences to rate formulated state equation ε σ = +ε p B ε = σ + εp ε = σ σ n + ε p σ=σ n + ε ε p =σ e ε p wit σ =σ σ n and te elastic "trial" stress σ e :=σ n + ε. lastic loading () fl ε p := 0 fl σ =σ e if φ@σ e,k n D 0. Plastic loading (P) fl Apply B finite differences to elastic-plastic "flow rule": ε p φ =λ =λ sign@σd B ε p = λsign@σd κ φ =λ K = λ B κ = λ K = K n H κ = K n + H λ xplicit expression for integrated stress: σ =σ e λ sign@σd fl sign@σd = sign@σ e D» σ» =» σ e» λ φ@»σ», KD =» σ e» λ HK n + H λl σ y =» σ e» λ H λ σ y K n = 0 λ =» σe» σ y K n + H fl Integrated stress = φ el + H wit φ el = φ@»σ e», K n D σ=σ e ε p =σ e λ sign@σ e D =σ e σ=σ e» σ e» σ y K n + H» σ e» σ y K n + H σ e» σ e» =σe H σy+kn»σ e» + H sign@σ e D (1) (2)
6 6 Notes_F8.nb Tangent relationsip dσ dε =? wit dσ dε = d»σ e» dε + H sign@σe D = d»σ e» dε + H sign@σe D = + H = H + H (3) d» σ e» dε = sign@σ e D (4) ü Linear creep law ü Non-linear Norton creep law à Analysis of beam cross section - solution strategy for isotropic ardening plasticity ü General Consider beam cross section Figure 9 ü Problem Apply bending moment about y-axis and normal force along x-axis of beam. Compute (for given bending moment and axial force) uniaxial stress distribution along cross section eigt ü Analysis Kinematics: Assume uler-beroulli teory fl Plane beam sections remain plane after deformation:
7 Notes_F8.nb 7 Figure 10 Described via 1) curvature κ =κ@xd = w and 2) axial elongation ξ =ξ@xd = u Local uniaxial strain of "beam fiber" at given value of x along considered beam: ε@zd =κz +ξ=h1, zl J ξ κ N Constitutive law: Specify te relation σ = f Hε@z, tdl =?; Assume elasticity σ=ε@zd Stress resultants (section forces): Specify te definition of te bending moment and te axial force! fl N = σ A A M = σ z A A Out-of-balance force stablis out-of-balance force vector defined as g@ pd = i N Ÿ A σ A y κ j k M Ÿ A σ z A z := 0 wit p = J ξ N { Note! From kinematics we ave uniaxial fiber strain ε@zd = H1, zl J ξ ξ N ε=h1, zl J κ κ N Consider as generic case: isotropic ardening plasticity σ=σ e» σ e» σ y K n sign@σ + H e D σ e relates to current time step t = n t + t ; In context of beam cross section we obtain: σ e =σ n + ε wit ε@zd := H ξ + κ z L ü Main solution procedure - Non-linear solution strategy 1) For eac loadstep at time t = n t + t wit 0 t = 0
8 8 Notes_F8.nb 2) stablis p i D 3) Solve for g@ pd := 0; consider Newton procedure D i ξ i+1 = g i wit D = g p 4) Update improvement of solution = tangent stiffness matrix p i+1 = p i + ξ i+1 5) Ceck convergence, Compute error e =» g i» e tol go to step 2 e > tol got step 6 6) Set incremental soltion p = p i,σ = n σ; goto step 1. 7) nd of procedure ü Spatial integration - out-of-balance force Problem: determine generalized stresses fl consider "layerwise" numerical integration: Subdivide integration interval I D ( = is te cross section eigt) into finite segments 8 z k = ê Nint< k=1,nint Carry out integration as: -Normal force, integration interval subdivided into layers Nint N = σ b z = k=1 σ k b k z k wit σ k = Jσ e -Bending moment, integration interval subdivided into layers» σ e» σ y K n + H sign@σ e DN k Nint M = σ zb z = k=1 σ k b k z k z k ü Tangent stiffness matrix Consider linearization: we ave dg = D d p wit D = g p dg = J dn dm N = J 1 N dσ b x = 9dσ = z D = a J 1 z z z 2 N b x were (sow as omework) a = ε = tangent stiffness moduli =...= H ε d ε = a H1, zl J d ξ d κ N= + H
9 Notes_F8.nb 9 ü xample: lastic-plastic response of symmetric beam cross section Consider te symmetric beam cross section in te figure below. Figure 11 Te material is assumed to be isotropiic elastic-plastic wit linear ardening, wereby te integrated stress (according to te discussion above) becomes σ=σ e ε p =σ e λ sign@σ e D = (7) 9 λ = φ el + H =» σ =σe e» σ y K n sign@σ + H e D In te present case we assume total deformation teory wic means tat te (total) strain is applied in one single step, i.e. ε ε, λ λ = σ= ε» σe» σ y + H + H σ = ε = 1 + H H» ε» σ yl sign@εd + H Te absolute value of te uniaxial stress is tus given by:» ε»» σ» = 9 +H HH» ε» +σ yl H» ε» σ yl H» ε» σ yl = HH» ε» +σ yl + H HL HPL In tis case te strain of te cross section is given by ε@zd =κz sign@εd sign@zd for κ>0 Te elastic range of te cross section is given by σ y =ε y@z y D =κz y =κ el 2 z y 2 I κ el κ Hence, te integrated stress of te cross section is given by wit σ=» σ» sign@εd M (8) (9) (10) (11) σ=9 κ z +H HH κ z +σ z y Te bending moment is defined by»z» L HL HPL wit σ=» σ» sign@εd (12)
10 10 Notes_F8.nb 2 z y M = 2 σ zb z = 2 ba σ z z Tis is evaluated as: z y 2 ba κ z 2 z H z y 2 σ z z = z y 2 JH κ z 2 +σ y z 2» z» N z (13) M@κD = b 3 12 κ if κ κ el = 2 σ y M@κD = H + H wic can be plotted as: b 3 12 J 1 + m H I κ el κ M2 + 3 σ y H 1 κ J1 I κ el κ M2 NN κ if κ>κ el M_@κD Responses wit H 20 and H κ Figure 12
a) Give an example of a case when an (s,s) policy is not the same as an (R,Q) policy. (2p)
roblem a) Give an example of a case wen an (s,s) policy is not te same as an (R,) policy. (p) b) Consider exponential smooting wit te smooting constant α and moving average over N periods. Ten, tese two
More informationIn the following I do the whole derivative in one step, but you are welcome to split it up into multiple steps. 3x + 3h 5x 2 10xh 5h 2 3x + 5x 2
Mat 160 - Assignment 3 Solutions - Summer 2012 - BSU - Jaimos F Skriletz 1 1. Limit Definition of te Derivative f( + ) f() Use te limit definition of te derivative, lim, to find te derivatives of te following
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationHIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR 1D PARABOLIC EQUATIONS. Ahmet İzmirlioğlu. BS, University of Pittsburgh, 2004
HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR D PARABOLIC EQUATIONS by Ahmet İzmirlioğlu BS, University of Pittsburgh, 24 Submitted to the Graduate Faculty of Art and Sciences in partial fulfillment of
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationSurface Curvatures. (Com S 477/577 Notes) Yan-Bin Jia. Oct 23, Curve on a Surface: Normal and Geodesic Curvatures
Surface Curvatures (Com S 477/577 Notes) Yan-Bin Jia Oct 23, 2017 1 Curve on a Surface: Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves
More informationFinite difference method for the Black and Scholes PDE (TP-1)
Numerical metods for PDE in Finance - ENSTA - S1-1/MMMEF Finite difference metod for te Black and Scoles PDE (TP-1) November 2015 1 Te Euler Forward sceme We look for a numerical approximation of te European
More informationONERA Fatigue Model. Z-set group. March 14, Mines ParisTech, CNRS UMR 7633 Centre des Matériaux BP 87, Evry cedex, France
ONERA Fatigue Model Z-set group Mines ParisTech, CNRS UMR 7633 Centre des Matériaux BP 87, 91003 Evry cedex, France March 14, 2013 Plan 1 ONERA Fatigue Model 2 Basic Tools Multiaxial stress amplitude (SEH)
More informationACC 471 Practice Problem Set # 4 Fall Suggested Solutions
ACC 471 Practice Problem Set # 4 Fall 2002 Suggested Solutions 1. Text Problems: 17-3 a. From put-call parity, C P S 0 X 1 r T f 4 50 50 1 10 1 4 $5 18. b. Sell a straddle, i.e. sell a call and a put to
More informationCurves on a Surface. (Com S 477/577 Notes) Yan-Bin Jia. Oct 24, 2017
Curves on a Surface (Com S 477/577 Notes) Yan-Bin Jia Oct 24, 2017 1 Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves on the surface. Let
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationHeterogeneous Firm, Financial Market Integration and International Risk Sharing
Heterogeneous Firm, Financial Market Integration and International Risk Sharing Ming-Jen Chang, Shikuan Chen and Yen-Chen Wu National DongHwa University Thursday 22 nd November 2018 Department of Economics,
More informationWhy Surplus Consumption in the Habit Model May be Less Pe. May be Less Persistent than You Think
Why Surplus Consumption in the Habit Model May be Less Persistent than You Think October 19th, 2009 Introduction: Habit Preferences Habit preferences: can generate a higher equity premium for a given curvature
More informationPrincipal Curvatures
Principal Curvatures Com S 477/577 Notes Yan-Bin Jia Oct 26, 2017 1 Definition To furtheranalyze thenormal curvatureκ n, we make useof the firstandsecond fundamental forms: Edu 2 +2Fdudv +Gdv 2 and Ldu
More informationThe study guide does not look exactly like the exam but it will help you to focus your study efforts.
Mat 0 Eam Study Guide Solutions Te study guide does not look eactly like te eam but it will elp you to focus your study efforts. Here is part of te list of items under How to Succeed in Mat 0 tat is on
More informationPractice Exam 1. Use the limit laws from class compute the following limit. Show all your work and cite all rules used explicitly. xf(x) + 5x.
Practice Exam 1 Tese problems are meant to approximate wat Exam 1 will be like. You can expect tat problems on te exam will be of similar difficulty. Te actual exam will ave problems from sections 11.1
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationState of Stress in Three Dimensions
State of Stress in Three Dimensions Theories of failure Introduction: Due to large numbers of examples of compound stresses met with in engineering practice, the cause of failure or permanent set under
More informationUNIVERSITY OF CALIFORNIA College of Engineering Departments of Mechanical Engineering and Material Science & Engineering
Fall 006 UNIVERSITY OF CALIFORNIA College of Engineering Departments of Mechanical Engineering and Material Science & Engineering MSEc113/MEc14 Mechanical Behavior of Materials Midterm #1 September 19
More informationA Model of Financial Intermediation
A Model of Financial Intermediation Jesús Fernández-Villaverde University of Pennsylvania December 25, 2012 Jesús Fernández-Villaverde (PENN) A Model of Financial Intermediation December 25, 2012 1 / 43
More informationAK and reduced-form AK models. Consumption taxation. Distributive politics
Chapter 11 AK and reduced-form AK models. Consumption taxation. Distributive politics The simplest model featuring fully-endogenous exponential per capita growth is what is known as the AK model. Jones
More informationRiemannian Geometry, Key to Homework #1
Riemannian Geometry Key to Homework # Let σu v sin u cos v sin u sin v cos u < u < π < v < π be a parametrization of the unit sphere S {x y z R 3 x + y + z } Fix an angle < θ < π and consider the parallel
More informationAK and reduced-form AK models. Consumption taxation.
Chapter 11 AK and reduced-form AK models. Consumption taxation. In his Chapter 11 Acemoglu discusses simple fully-endogenous growth models in the form of Ramsey-style AK and reduced-form AK models, respectively.
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationPartitioned Analysis of Coupled Systems
Partitioned Analysis of Coupled Systems Hermann G. Matthies, Rainer Niekamp, Jan Steindorf Technische Universität Braunschweig Brunswick, Germany wire@tu-bs.de http://www.wire.tu-bs.de Coupled Problems
More informationChapter 8. Introduction to Endogenous Policy Theory. In this chapter we begin our development of endogenous policy theory: the explicit
Capter 8 Introduction to Endogenous Policy Teory In tis capter we begin our development of endogenous policy teory: te explicit incorporation of a model of politics in a model of te economy, permitting
More informationSENSITIVITY ESTIMATES FROM CHARACTERISTIC FUNCTIONS
Proceedings of te 27 Winter Simulation Conference S. G. Henderson, B. Biller, M.-H. Hsie, J. Sortle, J. D. Tew, and R. R. Barton, eds. SENSITIVITY ESTIMATES FROM CHARACTERISTIC FUNCTIONS Paul Glasserman
More informationPricing American Options Using a Space-time Adaptive Finite Difference Method
Pricing American Options Using a Space-time Adaptive Finite Difference Method Jonas Persson Abstract American options are priced numerically using a space- and timeadaptive finite difference method. The
More informationGrowth and Distributional Effects of Inflation with Progressive Taxation
MPRA Munich Personal RePEc Archive Growth and Distributional Effects of Inflation with Progressive Taxation Fujisaki Seiya and Mino Kazuo Institute of Economic Research, Kyoto University 20. October 2010
More informationOptimum thresholding using mean and conditional mean square error
Optimum tresolding using mean and conditional mean square error José E. Figueroa-López and Cecilia Mancini August 7, 7 Abstract We consider a univariate semimartingale model for te logaritm of an asset
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationUsing Fiber Reinforced Polymer to Restore Deteriorated Structural Members
International Journal of Material and Mechanical Engineering, 01, 1: 1-7 - 1 - Published Online April 01 http://www.ijm-me.org Using Fiber Reinforced Polymer to Restore Deteriorated Structural Members
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
More informationSentiments and Aggregate Fluctuations
Sentiments and Aggregate Fluctuations Jess Benhabib Pengfei Wang Yi Wen June 15, 2012 Jess Benhabib Pengfei Wang Yi Wen () Sentiments and Aggregate Fluctuations June 15, 2012 1 / 59 Introduction We construct
More informationDistorted Trade Barriers: A Dissection of Trade Costs in a Distorted Gravity Model
Distorted Trade Barriers: A Dissection of Trade Costs in a Distorted Gravity Model Tibor Besedeš Georgia Institute of Tecnology Mattew T. Cole California Polytecnic State University October 26, 2015 Abstract
More informationMATH 220: DUHAMEL S PRINCIPLE
MATH 22: DUHAMEL S PRINCIPLE Although we have lved only the homogeneous heat equation on R n, the same method employed there al lves the inhomogeneous PDE. As an application of these methods, let s lve
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationOn a Tractable Small Open Economy Model with Endogenous Monetary-policy Trade-offs
On a Tractable Small Open Economy Model wit Endogenous Monetary-policy Trade-offs Jaime Alonso-Carrera and Timoty Kam Tis version: November 30, 2009 ABSTRACT We propose a new igly tractable small open
More informationExercises in Growth Theory and Empirics
Exercises in Growth Theory and Empirics Carl-Johan Dalgaard University of Copenhagen and EPRU May 22, 2003 Exercise 6: Productive government investments and exogenous growth Consider the following growth
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationMarket shares and multinationals investment: a microeconomic foundation for FDI gravity equations
Market sares and multinationals investment: a microeconomic foundation for FDI gravity equations Gaetano Alfredo Minerva November 22, 2006 Abstract In tis paper I explore te implications of te teoretical
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #1 March 3, 2014 Name: SOLUTION ID#: PS.:
More informationThe coupling of electrical eddy current heat production and air cooling
Söllerhaus 2004 Institute of Applied Analysis and Numerical Simulation The coupling of electrical eddy current heat production and air cooling Institut of Applied Analysis and Numerical Simulation D-70569
More informationCredit and hiring. Vincenzo Quadrini University of Southern California, visiting EIEF Qi Sun University of Southern California.
Credit and hiring Vincenzo Quadrini University of Southern California, visiting EIEF Qi Sun University of Southern California November 14, 2013 CREDIT AND EMPLOYMENT LINKS When credit is tight, employers
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann
More informationCalculus I Homework: Four Ways to Represent a Function Page 1. where h 0 and f(x) = x x 2.
Calculus I Homework: Four Ways to Represent a Function Page 1 Questions Example Find f(2 + ), f(x + ), and f(x + ) f(x) were 0 and f(x) = x x 2. Example Find te domain and sketc te grap of te function
More informationPDE Approach to Credit Derivatives
PDE Approach to Credit Derivatives Marek Rutkowski School of Mathematics and Statistics University of New South Wales Joint work with T. Bielecki, M. Jeanblanc and K. Yousiph Seminar 26 September, 2007
More informationCombined mode I stress intensity factors of slanted cracks
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Combined mode I stress intensity factors of slanted cracks To cite this article: A.E Ismail et al 2017 IOP Conf. Ser.: Mater.
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationAdditional material D Descriptive statistics on interest rate spreads Figure 4 shows the time series of the liquidity premium LP in equation (1. Figure 5 provides time series plots of all spreads along
More informationEcon 582 Nonlinear Regression
Econ 582 Nonlinear Regression Eric Zivot June 3, 2013 Nonlinear Regression In linear regression models = x 0 β (1 )( 1) + [ x ]=0 [ x = x] =x 0 β = [ x = x] [ x = x] x = β it is assumed that the regression
More informationA General Welfare Decomposition for CGE Models
urdue University urdue e-ubs GTA Tecnical apers Agricultural Economics 1-1-2000 A General Welfare Decomposition for CGE Models Kevin J Hanslow roductivity Commission, Australia Follow tis and additional
More informationSentiments and Aggregate Fluctuations
Sentiments and Aggregate Fluctuations Jess Benhabib Pengfei Wang Yi Wen March 15, 2013 Jess Benhabib Pengfei Wang Yi Wen () Sentiments and Aggregate Fluctuations March 15, 2013 1 / 60 Introduction The
More informationIntroduction to Computable General Equilibrium Model (CGE)
Introduction to Computable General Equilibrium Model (CGE Dazn Gillig & ruce. McCarl Department of gricultural Economics Texas &M University Course Outline Overview of CGE n Introduction to te Structure
More informationECON 200 EXERCISES (1,1) (d) Use your answer to show that (b) is not the equilibrium price vector if. that must be satisfied?
ECON 00 EXERCISES 4 EXCHNGE ECONOMY 4 Equilibrium in an ecange economy Tere are two consumers and wit te same utility function U ( ) ln H {, } Te aggregate endowment is tat prices sum to Tat is ( p, p)
More informationSAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax:
ProSINTAP - A Probabilistic Program for Safety Evaluation Peter Dillström SAQ / SINTAP / 09 SAQ KONTROLL AB Box 49306, 100 29 STOCKHOLM, Sweden Tel: +46 8 617 40 00; Fax: +46 8 651 70 43 June 1999 Page
More informationECON 815. A Basic New Keynesian Model II
ECON 815 A Basic New Keynesian Model II Winter 2015 Queen s University ECON 815 1 Unemployment vs. Inflation 12 10 Unemployment 8 6 4 2 0 1 1.5 2 2.5 3 3.5 4 4.5 5 Core Inflation 14 12 10 Unemployment
More informationPRICE INDEX AGGREGATION: PLUTOCRATIC WEIGHTS, DEMOCRATIC WEIGHTS, AND VALUE JUDGMENTS
Revised June 10, 2003 PRICE INDEX AGGREGATION: PLUTOCRATIC WEIGHTS, DEMOCRATIC WEIGHTS, AND VALUE JUDGMENTS Franklin M. Fiser Jane Berkowitz Carlton and Dennis William Carlton Professor of Economics Massacusetts
More informationDENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE
DENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE Marcos Salvai FaMAF, Ciudad Universitaria, 5000 Córdoba, Argentina. e-mail: salvai@mate.uncor.edu Abstract Let S
More informationUnit M2.2 (All About) Stress
Unit M. (All About) Stress Readings: CDL 4., 4.3, 4.4 16.001/00 -- Unified Engineering Department of Aeronautics and Astronautics Massachusetts Institute of Technology LEARNING OBJECTIVES FOR UNIT M. Through
More informationMECHANICS OF MATERIALS
CHAPTER 7 Transformations MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech Universit of Stress and Strain 006 The McGraw-Hill Companies,
More informationAustralian Journal of Basic and Applied Sciences. Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model
AENSI Journals Australian Journal of Basic and Applied Sciences Journal home page: wwwajbaswebcom Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model Khawla Mustafa Sadiq University
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 2, Mikhail Zaslavsky 2 University of Michigan, Ann
More informationState-Dependent Pricing and the Paradox of Flexibility
State-Dependent Pricing and the Paradox of Flexibility Luca Dedola and Anton Nakov ECB and CEPR May 24 Dedola and Nakov (ECB and CEPR) SDP and the Paradox of Flexibility 5/4 / 28 Policy rates in major
More informationEstimating Maximum Smoothness and Maximum. Flatness Forward Rate Curve
Estimating Maximum Smoothness and Maximum Flatness Forward Rate Curve Lim Kian Guan & Qin Xiao 1 January 21, 22 1 Both authors are from the National University of Singapore, Centre for Financial Engineering.
More informationWhat Can Rational Investors Do About Excessive Volatility and Sentiment Fluctuations?
What Can Rational Investors Do About Excessive Volatility and Sentiment Fluctuations? Bernard Dumas INSEAD, Wharton, CEPR, NBER Alexander Kurshev London Business School Raman Uppal London Business School,
More informationan example Giacomo Boffi
an example http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 2, 27 vt P + x u(x, t L A uniform beam, (unit mass
More informationModeling Co-movements and Tail Dependency in the International Stock Market via Copulae
Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae Katja Ignatieva, Eckhard Platen Bachelier Finance Society World Congress 22-26 June 2010, Toronto K. Ignatieva, E.
More informationGeometry Factors pt. II, Error Analysis, and Optimizations
Geometry Factors pt. II, Error Analysis, and Optimizations Christopher Coppola November 11, 2014 Simulation Method Error Analysis Optimizations Simulation Goals -Calculate geometry factors -Optimize pressure
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationLeader or Follower? A Payoff Analysis in Quadratic Utility Harsanyi Economy
Leader or Follower? A Payoff Analysis in Quadratic Utility Harsanyi Economy Sai Ma New York University Oct. 0, 015 Model Agents and Belief There are two players, called agent i {1, }. Each agent i chooses
More informationarxiv: v1 [math.pr] 15 Dec 2011
Parameter Estimation of Fiber Lay down in Nonwoven Production An Occupation Time Approach Wolfgang Bock, Thomas Götz, Uditha Prabhath Liyanage arxiv:2.355v [math.pr] 5 Dec 2 Dept. of Mathematics, University
More informationHousehold Debt, Financial Intermediation, and Monetary Policy
Household Debt, Financial Intermediation, and Monetary Policy Shutao Cao 1 Yahong Zhang 2 1 Bank of Canada 2 Western University October 21, 2014 Motivation The US experience suggests that the collapse
More informationSupplementary Appendix (Not for Publication) Supplementary Appendix: Additional Proofs
Supplementary Appendix Not for Publication B Supplementary Appendix: Additional Proofs B.1 Two-Generations Model: Convexity of Welfare Gain Proposition 3. Applying Definition 1 to equation 4 gives the
More informationIntroduction to Probability Theory and Stochastic Processes for Finance Lecture Notes
Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,
More informationSolving Asset-Pricing Models with Recursive Preferences
Solving Asset-Pricing Models with Recursive Preferences Walter Pohl University of Zurich Karl Schmedders University of Zurich and Swiss Finance Institute Ole Wilms University of Zurich July 5, Abstract
More informationUnderstanding Predictability (JPE, 2004)
Understanding Predictability (JPE, 2004) Lior Menzly, Tano Santos, and Pietro Veronesi Presented by Peter Gross NYU October 19, 2009 Presented by Peter Gross (NYU) Understanding Predictability October
More informationFinancial Econometrics
Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value
More informationMAFS Computational Methods for Pricing Structured Products
MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )
More information3.1 THE 2 2 EXCHANGE ECONOMY
Essential Microeconomics -1-3.1 THE 2 2 EXCHANGE ECONOMY Private goods economy 2 Pareto efficient allocations 3 Edgewort box analysis 6 Market clearing prices and Walras Law 14 Walrasian Equilibrium 16
More informationarxiv: v1 [q-fin.cp] 1 Nov 2016
Essentially high-order compact schemes with application to stochastic volatility models on non-uniform grids arxiv:1611.00316v1 [q-fin.cp] 1 Nov 016 Bertram Düring Christof Heuer November, 016 Abstract
More informationA Term-Structure Model for Dividends and Interest Rates
A Term-Structure Model for Dividends and Interest Rates Sander Willems Joint work with Damir Filipović School and Workshop on Dynamical Models in Finance May 24th 2017 Sander Willems (SFI@EPFL) A TSM for
More informationInformation Processing and Limited Liability
Information Processing and Limited Liability Bartosz Maćkowiak European Central Bank and CEPR Mirko Wiederholt Northwestern University January 2012 Abstract Decision-makers often face limited liability
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More information5.1 Gauss Remarkable Theorem
5.1 Gauss Remarkable Theorem Recall that, for a surface M, its Gauss curvature is defined by K = κ 1 κ 2 where κ 1, κ 2 are the principal curvatures of M. The principal curvatures are the eigenvalues of
More informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More informationPrincipal Stresses: Interpreting Principal Stresses
Principal Stresses: This document talks about the principal stresses in a gross sense than its theoretical explanations. And, here I go... Most of the times we do wonder what use are these principal stresses
More informationSupplemantary material to: Leverage causes fat tails and clustered volatility
Supplemantary material to: Leverage causes fat tails and clustered volatility Stefan Turner a,b J. Doyne Farmer b,c Jon Geanakoplos d,b a Complex Systems Researc Group, Medical University of Vienna, Wäringer
More informationChapter 8: CAPM. 1. Single Index Model. 2. Adding a Riskless Asset. 3. The Capital Market Line 4. CAPM. 5. The One-Fund Theorem
Chapter 8: CAPM 1. Single Index Model 2. Adding a Riskless Asset 3. The Capital Market Line 4. CAPM 5. The One-Fund Theorem 6. The Characteristic Line 7. The Pricing Model Single Index Model 1 1. Covariance
More informationIntroductory Microeconomics
Prof. Wolfram Elsner Faculty of Business Studies and Economics iino Institute of Institutional and Innovation Economics Introductory Microeconomics More Formal Concepts of Game Theory and Evolutionary
More informationMATH60082 Example Sheet 6 Explicit Finite Difference
MATH68 Example Sheet 6 Explicit Finite Difference Dr P Johnson Initial Setup For the explicit method we shall need: All parameters for the option, such as X and S etc. The number of divisions in stock,
More informationOn Existence of Equilibria. Bayesian Allocation-Mechanisms
On Existence of Equilibria in Bayesian Allocation Mechanisms Northwestern University April 23, 2014 Bayesian Allocation Mechanisms In allocation mechanisms, agents choose messages. The messages determine
More informationFinal Exam Key, JDEP 384H, Spring 2006
Final Exam Key, JDEP 384H, Spring 2006 Due Date for Exam: Thursday, May 4, 12:00 noon. Instructions: Show your work and give reasons for your answers. Write out your solutions neatly and completely. There
More informationBluff Your Way Through Black-Scholes
Bluff our Way Through Black-Scholes Saurav Sen December 000 Contents What is Black-Scholes?.............................. 1 The Classical Black-Scholes Model....................... 1 Some Useful Background
More informationA Model of Central Bank Liquidity Provision
A of Central Bank Liquidity Provision James T.E. Chapman 1 Jonathan Chiu 1 Miguel Molico 1 1 Bank of Canada Bank of Canada 19 February 2009 A of Central Bank Liquidity Provision Policy Questions When a
More informationA Macro-Finance Model of the Term Structure: the Case for a Quadratic Yield Model
Title page Outline A Macro-Finance Model of the Term Structure: the Case for a 21, June Czech National Bank Structure of the presentation Title page Outline Structure of the presentation: Model Formulation
More informationGrowth transmission. Econ 307. Assume. How much borrowing should be done? Implications for growth A B A B
Growt transmission Econ 307 Lecture 5 GDP levels differ dramatically across countries Wy does tis not open up uge gains from trade? According to te most simple model, very low GDP countries sould ave very
More information