Anisotropic elastoplasticity for cloth, knit and hair frictional contact supplementary technical document
|
|
- Geoffrey Dalton
- 5 years ago
- Views:
Transcription
1 Anisotroic elastolasticity for cloth, knit and hair frictional contact sulementary technical document Chenfanfu Jiang, Theodore Gast, Joseh Teran 1 QR differentiation Here we discuss how to comute δq and δr from F and δf. This is used for comuting the linearized force in imlicit time integration. If F = QR where Q is orthonormal with Q T Q = I, and r 11 r 12 r 13 R = r 22 r 23 is uer triangular, and we define the function Ψ(F) = ˆΨ(R) then we have r 33 δf = QδR + δqr δψ(f) = δ ˆΨ(R) δq T Q + Q T δq = 0 δr is uer triangular Q T δf = δr + Q T δqr 0 ω 3 ω 2 Ω = Q T δq = ω 3 0 ω 1 ω 2 ω ω 3 ω 2 Q T δf = δr + ω 3 0 ω 1 R ω 2 ω 1 0 a b c 0 ω 3 ω 2 r 11 r 12 r 13 d e f = 0 + ω 3 0 ω 1 0 r 22 r 23 g h i 0 0 ω 2 ω r 33 This imlies d = ω 3 r 11 g = ω 2 r 11 h = ω 2 r 12 + ω 1 r 22 We can use these three equations to solve for ω 1, ω 2, ω 3. After that, we can construct δq = QΩ and then δr = Q T δf Q T δqr. 1
2 The corresonding 2D result is a = (QT δf) 21 δq = aq r 11 ( 0 ) δr = Q T δf Q T δqr 2 Comuting stress δψ = δ ˆΨ : δf = F R : δr : (QδR + δqr) = F R : δr : (QδR) + : (δqr) = : δr (1) F F R It can be shown that the above holds for any δq and δr that satisfy δq T Q + Q T δq = 0 and δr being uer triangular. Secifically, it can be shown that given arbitrary δq, δr with δq T Q + Q T δq = 0 and δr uer triangular, δf such that δr = R Q (F) : δf, and δq = (F) : δf. F F Using null[ Q F (F)] = san {N 1, N 2, N 3, N 4, N 5, N 6 }, for all uer triangular δr, δf Q null[ Q F (F)] such that δq = Q F (F) : δf Q = Q F (F) : δ(f Q + F R ) F R null[ Q F (F)]. Further, it can be shown that QδR + δqr δf Q null[ Q F (F)] and with δf R = QδR + δqr δf Q, δf = δf Q + δf R roduces If we choose δq = 0 in Equation (1), then δr = R Q (F) : δf, and δq = (F) : δf. F F : (QδR) = F R : δr (Q T ) : δr = F R : δr ) and R have the same uer triangular art. is lower triangular, it is easy to show (entry wise rovable) that Q T F RT and Recall δr is any uer triangular matrix, therefore (Q T F Further more, since R T R RT have the same uer triangular art. 2
3 If we choose δr = 0 in Equation (1), then F : (δqr) = 0 ( F RT ) : δq = 0 ( F RT ) : (δqq T Q) = 0 ( F RT Q T ) : (δqq T ) = 0 ( F FT ) : (δqq T ) = 0 Since δqq T is an arbitrary skew symmetric matrix (due to δq T Q + Q T δq = 0), for the above equation to hold, we know F FT has to be symmetric. This also roves that Kirchoff stress τ = F FT is symmetric without needing to use conservation of angular momentum. Now we have τ = τ T F FT = F( F )T F RT Q T = QR( F )T Q T F RT = R( F )T Q Q T F RT is symmetric. In summary, Q T F RT and R RT have the same uer triangular art. Q T F RT is symmetric. We further denote this tensor with A := Q T F RT = Q T F FT Q = Q T τq. Therefore A is just τ written in the Q basis. 3 A curve in 3D We can construct the uer triangular art of A with R RT, then fill the rest using symmetry of A. Assuming then ˆΨ = f(r 11 ) g(r r 2 13) + h(r 22, r 23, r 33 ), R = [ f g r T 0 ˆP [ where r = (r 12, r 13 ) T r11 r T ], R =, 0 ˆR ˆP = h ˆR. We can show f A = r 11 + g r T r g r T ˆRT g r T ˆRT ˆP ˆRT Here, we choose ˆP ˆR T to be the dilational art of the Kirchoff stress from the 2 2 Stvk Hencky Drucker- Prager. assuming we get some ˆτ after the return maing of dry sand from an inut F = ˆR, we relace the bottom right corner of A with I, where = tr(ˆτ)/2. ], 3
4 Under such a choice, f A = r 11 + g r T r g r T ˆRT g r T. ˆRT I Mohr friction criteria leads to the following maximization roblem: Maximize d T τn + c F n T τn over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length. Maximize d T QAQ T n + c F n T QAQ T n over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length. Maximize d T Añ + c F ñ T Añ over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length, where ñ = Q T n and d = Q T d. Since in our discretization, the fiber direction is q 1, therefore q T 1 n = 0, then we know ñ = (0, c, s) for some θ. Maximize d T Añ + c F ñ T Añ over all ossible n that is (0, c, s) over all θ, with all d erendicular to ñ and has unit length. Lagrangian multilier can solve this roblem and gives the maximum ˆRr g + c F. If we choose g(x) = x, then we have our yield surface ˆRr + c F < 0. The return maing is then simly a scaling on r so that the yield criteria is satisfied. 4 A surface in 3D Surface is very similar to curve in our framework. With our codimensional discretization, we ma the triangle back to x-y lane. If the inut 3D triangle at rest is Â, ˆB, Ĉ, we can define ˆD 1 = ˆB Â and ˆD 2 = Ĉ Â. This forms an imaginary ˆD s, whose QR decomosition gives us the rotation from x-y lane of this triangle ˆQ, as well as the to art of ˆR 3 2 being the 2 2 version D m. The third column of ˆQ is the rotated D 3 where D 3 = e 3. With these recomutations, for any triangle d 1, d 2 in world sace, we can construct the full F as D 1 m F = [d 1, d 2, d 3 ]. 1 In MPM, we have d n+1 3 = (I+ t v)d n 3, with d 0 3 = ˆq 3. Kee in mind that D m is uer triangular, therefore D 1 m is uer triangular. Now do the thin QR decomosition [d 1, d 2 ] = [q 1, q 2 ] R where R is 2 2 uer. If we construct q 3 = q 1 q 2, then R (hx, h [d 1, d 2, d 3 ] = [q 1, q 2, q 3 ] y ) T 0 T h z where Qh := d 3. Now we have constructed the (unique) QR decomosition of F as R (hx, h F = [q 1, q 2, q 3 ] y ) T D 1 0 T m := QR = Q r 11 r 12 r 13 r h z 1 22 r 23 r 33 We can see h = r 3 = Q T d 3. 4
5 The revious lemma in curve still holds: have the same uer triangular art. Q T F RT is symmetric. We further denote this tensor with A := Q T F RT = Q T F FT Q = Q T τq. Therefore A is just τ written in the Q basis. We can construct the uer triangular art of A with R RT, then fill the rest using symmetry of A. 4.1 Surface elastolasticity Recall D 3 = e 3, Q T F RT and R RT r 11 r 12 r 13 F = Q r 22 r 23 Physically, the to left is in lane (x-y) deformation of the triangle, since r 11 r 12 r 13 r 22 r 23 e 3 = r 3, r 33 we know r 3 reresents the deformation of D 3. From h = r 3 = Q T d 3, we can say r 3 is the length change, r r 2 23 is the shearing (or the deviation from being erendicular to the x-y triangle lane). Note that shearing also ernalizes length change. From these, we can define then r 33 ˆΨ = f(r 33 ) g(r r 2 23) + h(r 11, r 12, r 21, r 22 ), R = ˆP g r 0 T f, ˆR where r = (r 13, r 23 ) T r, R = 0 T, r ˆP = h 33 ˆR. We can show ˆP A = ˆRT + g rr T g r 33 r g r 33 r T f r 33 Mohr friction criteria leads to the following maximization roblem: Maximize d T τn + c F n T τn over all ossible n that is erendicular to the manifold lane, with all d erendicular to n and has unit length. Maximize d T QAQ T n + c F n T QAQ T n over all ossible n that is erendicular to the manifold lane, with all d erendicular to n and has unit length. Maximize d T Añ + c F ñ T Añ over all ossible n that is erendicular to the manifold lane, with all d erendicular to n and has unit length, where ñ = Q T n and d = Q T d. n erendicular to manifold lane means n = k(d 1 d 2 ) for some k (unit length constraint is extra). Recall [d 1, d 2 ] = [q 1, q 2 ] R, where R is uer triangular 2 2, this means k(d 1 d 2 ) = z(q 1 q 2 ) for some z. Therefore, n = ±q 3. Therefore ñ = Q T n = ±e 3. Maximize d T Añ + c F ñ T Añ over ñ = ±e 3, with all d = (c, s, 0) for some θ 5
6 The maximum is ±g r 33 r + c F f r 33 Assume f(x) = 1 3 k(1 x)3 for x 1, 0 otherwise. g(x) = γx. When r 33 > 1, f = 0, the maximum is γr 33 r. In this case the return maing is making r to be zero. When r 33 < 1, the maximum is ±γr 33 r c F k(r 33 1) 2 r 33. The yield surface is therefore max(± γ k r 33 r c F (r 33 1) 2 r 33 ) 0 If r 33 is negative (corresonding to inverted collision), the max should choose γ k. The return maing is setting r to be 0. Otherwise, try to scale r (when necessary) to satisfy 5 A curve in 2D γ k r c F (r 33 1) 2 0 2D derivation follows 3D curve derivation: Q T F RT and R RT have the same uer triangular art. Q T F RT is symmetric. We further denote this tensor with A := Q T F RT = Q T F FT Q = Q T τq. Therefore A is just τ written in the Q basis. The energy choice is then r11 r where R = 12, 0 r ˆP = h 22 ˆR. We can show ˆΨ = f(r 11 ) g(r2 12) + h(r 22 ), R = f g r 12 0 h, f A = r 11 + g r12 2 g r 12 r 22 g r 12 r 22 h r 22 Mohr friction criteria leads to the following maximization roblem: Maximize d T τn + c F n T τn over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length. Maximize d T QAQ T n + c F n T QAQ T n over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length. Maximize d T Añ + c F ñ T Añ over all ossible n that is erendicular to the fiber direction, with all d erendicular to n and has unit length, where ñ = Q T n and d = Q T d. Since in our discretization, the fiber direction is q 1, therefore q T 1 n = 0, then we know ñ = ±(0, 1) Maximize d T Añ + c F ñ T Añ over ñ = (0, ±1) and d = (±1, 0) 6
7 Maximize d T Añ + c F ñ T Añ over ñ = (0, ±1) and d = (±1, 0) The maximum is ±g r 12 r 22 + c F h r 22 Assume g(x) = γx, h(x) = 1 3 s(1 x)3 for x 1, 0 otherwise. When r 22 > 1, h = 0, the maximum is γr 12 r 22, the return maing is setting r 12 = 0. When r 22 < 1, the maximum is The yield surface is therefore ±γ r 12 r 22 c F s(1 r 22 ) 2 r 22 max(± γ s r 12 r 22 c F (1 r 22 ) 2 r 22 ) 0 If r 22 < 0, the max should choose γ s. The return aming is setting r 12 = 0. Otherwise, r 22 [0, 1], try to scale r 12 to satisfy γ s r 12 c F (1 r 22 ) Derivative of ˆF E In the aer we mention that ˆF E x i is a third order tensor, and does not deend on ˆx because ˆF E is linear in ˆx i. Here we give the derivation of comuting this derivative. We have by definition Plugging in ˆF E,α = i Differentiating we have ˆF E,α x = ˆd E,β(ˆx) = ˆx mesh(,β) ˆx mesh(,0) (2) (wi,mesh(,β) n wn i,mesh(,0) )x i,αd 1 (wi,mesh(,β) n wn i,mesh(,0) )δ αζd 1 ˆd E,β(ˆx) = ( ˆx),β (3) ˆx q = ˆx i wiq n i (4) ˆF E (ˆx) = ˆd E D 1 (5),β +,β + w n i x κ,κβ x i,αd 1,β w n i x κ,κβ δ αζd 1,β which doesn t deend on ˆx. Note we can think of the undeformed segment or triangle as being aligned with the x-axis or xy-lane resectively. This will mean that the initial F is no longer I, but rather the rotation which mas the axis aligned element to its initial osition in sace. This is not an issue as our elastic energy density ψ is world sace rotation invariant. However this allows us to assume that D is block diagonal of the form D D = for segments and D 11 D 12 0 D = 0 D for triangles. This saves memory for D, and simlifies comuting D 1. 7
8 7 Comutation of Force on the Grid Define Then i (ˆx) = V 0 (ˆx) = (ˆx) = (ˆx) = F : ˆF E (6) x i α f (ii) qζ (ˆx) = = =1 F ζ F ζ ˆF,α E (7) F α x ( (wi,mesh(,β) n wn i,mesh(,0) )D 1,β + (8) wi n ),κβ x D 1,β κ ( (wi,mesh(,β) n wn i,mesh(,0) )D 1,β + (10) F ζ F ζ (ˆx) = f (ii) qζ (ˆx)wn i I (ii) (ˆx) = f (ii) qζ (ˆx)wn i I (ii) wi n ),κβ x δ β κ (δ q,mesh(,β) δ q,mesh(,0) )D 1 (δ q,mesh(,β) δ q,mesh(,0) )D 1 =γ+1 F ζ F ζ κ=1 (9) (11),β (12),β (13) w n i x κ,κβ δ β (14) w n i x κ,κ (15) 8
9 8 Pseudocode Algorithm 1 Simulate 1: rocedure Time Ste 2: Transfer to grid 3: Grid ste 4: Transfer to articles 5: Udate article state 6: Plasticity 1: rocedure Transfer to grid 2: for all grid nodes i do 3: m n i wn i m 4: vi n 1 m n wn i m ( v n i + C n (x i x n ) ) 1: rocedure Grid ste 2: vi vn i + Force increment( FE,n ) 3: v n+1 i Grid collisions( vi ) 4: ṽ n+1 i Friction( v n+1 i, v n+1 i vi ) 1: rocedure Transfer to articles 2: for all articles of tye (i) and (ii) do 3: v n+1 i β=0 w n iṽn+1 i 4: for all articles of tye (iii) do 5: v n+1 1 γ vn+1 mesh(,β) 6: for all articles do 7: C n+1 i w n iṽn+1 i (x i x n ) T 1: rocedure Udate article state 2: for all articles of tye (i) do 3: x n+1 wi(x n n i + tv n+1 i ) i 4: v i v n+1 i ( wi) n T 5: ˆFE,n+1 (I + t v )F E,n 6: for all articles of tye (ii) do 7: x n+1 i w n i(x n i + tv n+1 i ) 8: for all articles of tye (iii) do 9: x n+1 1 γ xn+1 mesh(,β) β=0 10: v v n+1 i ( wi) n T i 11: for β = 1 to γ do 12: ˆdE,n+1,β x n+1 mesh(,β) xn+1 mesh(,0) 13: for β =γ + 1 to 3 do 14: ˆdE,n+1,β 15: ˆFE,n+1 (I + t v ) E,n+1 ˆd D 1 1: rocedure Plasticity 2: for all articles of tye (i) and (iii) do 3: F E,n+1 Return maing(ˆf E,n+1 ) 9
Supplementary Technical Document
Supplementary Technical Document Qi Guo, Xuchen Han, Chuyuan Fu, Theodore Gast, Rasmus Tamstorf, Joseph Teran 1 FEM Force computation We compute forces on the control points x p by fp KL = ΨS F KL,Etr
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationFind the equation of a plane perpendicular to the line x = 2t + 1, y = 3t + 4, z = t 1 and passing through the point (2, 1, 3).
CME 100 Midterm Solutions - Fall 004 1 CME 100 - Midterm Solutions - Fall 004 Problem 1 Find the equation of a lane erendicular to the line x = t + 1, y = 3t + 4, z = t 1 and assing through the oint (,
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationPROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2)
PROFIT MAXIMIZATION DEFINITION OF A NEOCLASSICAL FIRM A neoclassical firm is an organization that controls the transformation of inuts (resources it owns or urchases into oututs or roducts (valued roducts
More informationStrong Matching of Points with Geometric Shapes
Strong Matching of Points with Geometric Shaes Ahmad Biniaz Anil Maheshwari Michiel Smid School of Comuter Science, Carleton University, Ottawa, Canada December 9, 05 In memory of Ferran Hurtado. Abstract
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationLecture 6. 2 Recurrence/transience, harmonic functions and martingales
Lecture 6 Classification of states We have shown that all states of an irreducible countable state Markov chain must of the same tye. This gives rise to the following classification. Definition. [Classification
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationF(p) y + 3y + 2y = δ(t a) y(0) = 0 and y (0) = 0.
Page 5- Chater 5: Lalace Transforms The Lalace Transform is a useful tool that is used to solve many mathematical and alied roblems. In articular, the Lalace transform is a technique that can be used to
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationSimplifications to Conservation Equations
Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationFinding Shortest Hamiltonian Path is in P. Abstract
Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs
More informationExtension of Minimax to Infinite Matrices
Extension of Minimax to Infinite Matrices Chris Calabro June 21, 2004 Abstract Von Neumann s minimax theorem is tyically alied to a finite ayoff matrix A R m n. Here we show that (i) if m, n are both inite,
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationarxiv:math/ v4 [math.gn] 25 Nov 2006
arxiv:math/0607751v4 [math.gn] 25 Nov 2006 On the uniqueness of the coincidence index on orientable differentiable manifolds P. Christoher Staecker October 12, 2006 Abstract The fixed oint index of toological
More informationINTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL
Symosium on Ice (26) INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL Mohamed O. ElSeify and Thomas G. Brown University of Calgary, Calgary, Canada ABSTRACT Current ice rubble load
More informationP043 Anisotropic 2.5D - 3C Finite-difference Modeling
P04 Anisotroic.5D - C Finite-difference Modeling A. Kostyukevych* (esseral echnologies Inc.), N. Marmalevskyi (Ukrainian State Geological Prosecting Institute), Y. Roganov (Ukrainian State Geological Prosecting
More informationGeo-E2010 Advanced Soil Mechanics L Wojciech Sołowski. 07 March 2017
Geo-E2010 Advanced Soil Mechanics L Wojciech Sołowski 07 March 2017 Soil modeling: critical state soil mechanics and Modified Cam Clay model Outline 1. Refresh of the theory of lasticity 2. Critical state
More information19th Bay Area Mathematical Olympiad. Problems and Solutions. February 28, 2017
th Bay Area Mathematical Olymiad February, 07 Problems and Solutions BAMO- and BAMO- are each 5-question essay-roof exams, for middle- and high-school students, resectively. The roblems in each exam are
More informationSECTION 12: HOMOTOPY EXTENSION AND LIFTING PROPERTY
SECTION 12: HOMOTOPY EXTENSION AND LIFTING PROPERTY In the revious section, we exloited the interlay between (relative) CW comlexes and fibrations to construct the Postnikov and Whitehead towers aroximating
More informationStochastic integration II: the Itô integral
13 Stochastic integration II: the Itô integral We have seen in Lecture 6 how to integrate functions Φ : (, ) L (H, E) with resect to an H-cylindrical Brownian motion W H. In this lecture we address the
More informationSolutions to Problems in Goldstein, Classical Mechanics, Second Edition. Chapter 9
Solutions to Problems in Goldstein, Classical Mechanics, Second Edition Homer Reid October 29, 2002 Chater 9 Problem 9. One of the attemts at combining the two sets of Hamilton s equations into one tries
More informationGENERICITY OF INFINITE-ORDER ELEMENTS IN HYPERBOLIC GROUPS
GENERICITY OF INFINITE-ORDER ELEMENTS IN HYPERBOLIC GROUPS PALLAVI DANI 1. Introduction Let Γ be a finitely generated grou and let S be a finite set of generators for Γ. This determines a word metric on
More informationALGEBRAIC TOPOLOGY MASTERMATH (FALL 2014) Written exam, 21/01/2015, 3 hours Outline of solutions
ALGERAIC TOPOLOGY MASTERMATH FALL 014) Written exam, 1/01/015, 3 hours Outline of solutions Exercise 1. i) There are various definitions in the literature. ased on the discussion on. 5 of Lecture 3, as
More informationProfessor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x
Spring, 997 ME85 - Introduction to Continuum Mechanics Midterm Exam II roblem. (+ points) (a) Let ρ be the mass density, v be the velocity vector, be the Cauchy stress tensor, and b be the body force per
More informationIMPROVED BOUNDS IN THE SCALED ENFLO TYPE INEQUALITY FOR BANACH SPACES
IMPROVED BOUNDS IN THE SCALED ENFLO TYPE INEQUALITY FOR BANACH SPACES OHAD GILADI AND ASSAF NAOR Abstract. It is shown that if (, ) is a Banach sace with Rademacher tye 1 then for every n N there exists
More informationIntroduction to Group Theory Note 1
Introduction to Grou Theory Note July 7, 009 Contents INTRODUCTION. Examles OF Symmetry Grous in Physics................................. ELEMENT OF GROUP THEORY. De nition of Grou................................................
More informationWe collect some results that might be covered in a first course in algebraic number theory.
1 Aendices We collect some results that might be covered in a first course in algebraic number theory. A. uadratic Recirocity Via Gauss Sums A1. Introduction In this aendix, is an odd rime unless otherwise
More informationComputing Differential Modular Forms
Comuting Differential Modular Forms Chris Hurlburt Northern Illinois University hurlburt@math.niu.edu July 28, 2003 Let > 3 be a rime number. Let R be a comlete, discrete valuation ring whose maximal ideal
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationTopic 7: Using identity types
Toic 7: Using identity tyes June 10, 2014 Now we would like to learn how to use identity tyes and how to do some actual mathematics with them. By now we have essentially introduced all inference rules
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationThe Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule
The Grah Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule STEFAN D. BRUDA Deartment of Comuter Science Bisho s University Lennoxville, Quebec J1M 1Z7 CANADA bruda@cs.ubishos.ca
More informationChapter 7 Rational and Irrational Numbers
Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers
More informationHARMONIC EXTENSION ON NETWORKS
HARMONIC EXTENSION ON NETWORKS MING X. LI Abstract. We study the imlication of geometric roerties of the grah of a network in the extendibility of all γ-harmonic germs at an interior node. We rove that
More informationQuantization of the Photon Field QED
Quantization of the Photon Field QED 21.05.2012 0.1 Reminder: Classical Electrodynamics Before we start quantizing the hoton field, let us reflect on classical electrodynamics. The Hamiltonian is given
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationBEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH
BEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH DORIN BUCUR, ALESSANDRO GIACOMINI, AND PAOLA TREBESCHI Abstract For Ω R N oen bounded and with a Lischitz boundary, and
More informationSAT based Abstraction-Refinement using ILP and Machine Learning Techniques
SAT based Abstraction-Refinement using ILP and Machine Learning Techniques 1 SAT based Abstraction-Refinement using ILP and Machine Learning Techniques Edmund Clarke James Kukula Anubhav Guta Ofer Strichman
More informationREFINED STRAIN ENERGY OF THE SHELL
REFINED STRAIN ENERGY OF THE SHELL Ryszard A. Walentyński Deartment of Building Structures Theory, Silesian University of Technology, Gliwice, PL44-11, Poland ABSTRACT The aer rovides information on evaluation
More informationMTH234 Chapter 15 - Multiple Integrals Michigan State University
MTH24 Chater 15 - Multile Integrals Michigan State University 6 Surface Area Just as arc length is an alication of a single integral, surface area is an alication of double integrals. In 15.6 we comute
More informationRobustness of classifiers to uniform l p and Gaussian noise Supplementary material
Robustness of classifiers to uniform l and Gaussian noise Sulementary material Jean-Yves Franceschi Ecole Normale Suérieure de Lyon LIP UMR 5668 Omar Fawzi Ecole Normale Suérieure de Lyon LIP UMR 5668
More informationJournal of Mathematical Analysis and Applications
J. Math. Anal. Al. 44 (3) 3 38 Contents lists available at SciVerse ScienceDirect Journal of Mathematical Analysis and Alications journal homeage: www.elsevier.com/locate/jmaa Maximal surface area of a
More informationComplex Analysis Homework 1
Comlex Analysis Homework 1 Steve Clanton Sarah Crimi January 27, 2009 Problem Claim. If two integers can be exressed as the sum of two squares, then so can their roduct. Proof. Call the two squares that
More informationOPTIMAL DESIGNS FOR TWO-LEVEL FACTORIAL EXPERIMENTS WITH BINARY RESPONSE
Statistica Sinica: Sulement OPTIMAL DESIGNS FOR TWO-LEVEL FACTORIAL EXPERIMENTS WITH BINARY RESPONSE Jie Yang 1, Abhyuday Mandal, and Dibyen Majumdar 1 1 University of Illinois at Chicago and University
More informationMath 104B: Number Theory II (Winter 2012)
Math 104B: Number Theory II (Winter 01) Alina Bucur Contents 1 Review 11 Prime numbers 1 Euclidean algorithm 13 Multilicative functions 14 Linear diohantine equations 3 15 Congruences 3 Primes as sums
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationOn the capacity of the general trapdoor channel with feedback
On the caacity of the general tradoor channel with feedback Jui Wu and Achilleas Anastasooulos Electrical Engineering and Comuter Science Deartment University of Michigan Ann Arbor, MI, 48109-1 email:
More informationA Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars
Aalborg Universitet A Simle And Efficient FEM-Imlementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars Published in: Proceedings of the 9th Nordic Seminar on Comutational
More informationThe Euler Phi Function
The Euler Phi Function 7-3-2006 An arithmetic function takes ositive integers as inuts and roduces real or comlex numbers as oututs. If f is an arithmetic function, the divisor sum Dfn) is the sum of the
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationOn the Toppling of a Sand Pile
Discrete Mathematics and Theoretical Comuter Science Proceedings AA (DM-CCG), 2001, 275 286 On the Toling of a Sand Pile Jean-Christohe Novelli 1 and Dominique Rossin 2 1 CNRS, LIFL, Bâtiment M3, Université
More informationGeometric Constraint Solving in 2D
Universitat Politècnica de Catalunya Deartament de Llenguatges i Sistemes Informàtics PhD Thesis Geometric Constraint Solving in 2D November, 1998 Antoni Soto i Riera Advisor: Robert Joan i Arinyo Contents
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationOn Two-Finger Grasping of Deformable Planar Objects
2011 IEEE International Conference on Robotics and Automation Shanghai International Conference Center May 9-13, 2011, Shanghai, China On Two-Finger Grasing of Deformable Planar Objects Yan-Bin Jia Feng
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationNOTES. Hyperplane Sections of the n-dimensional Cube
NOTES Edited by Sergei Tabachnikov Hyerlane Sections of the n-dimensional Cube Rolfdieter Frank and Harald Riede Abstract. We deduce an elementary formula for the volume of arbitrary hyerlane sections
More informationThe. Consortium. Continuum Mechanics. Original notes by Professor Mike Gunn, South Bank University, London, UK Produced by the CRISP Consortium Ltd
The C R I S P Consortium Continuum Mechanics Original notes b Professor Mike Gunn, South Bank Universit, London, UK Produced b the CRISP Consortium Ltd THOR OF STRSSS In a three dimensional loaded bod,
More informationp-adic Measures and Bernoulli Numbers
-Adic Measures and Bernoulli Numbers Adam Bowers Introduction The constants B k in the Taylor series exansion t e t = t k B k k! k=0 are known as the Bernoulli numbers. The first few are,, 6, 0, 30, 0,
More informationBest approximation by linear combinations of characteristic functions of half-spaces
Best aroximation by linear combinations of characteristic functions of half-saces Paul C. Kainen Deartment of Mathematics Georgetown University Washington, D.C. 20057-1233, USA Věra Kůrková Institute of
More informationFAQ and Errata Sheet for The Global Positioning System and Inertial Navigation: Theory and Practice
FAQ and Errata Sheet for The Global Positioning System and Inertial Navigation: Theory and Practice Jay Farrell Deartment of Electrical Engineering University of California, Riverside jay.farrell@ucr.edu
More informationDifference of Convex Functions Programming for Reinforcement Learning (Supplementary File)
Difference of Convex Functions Programming for Reinforcement Learning Sulementary File Bilal Piot 1,, Matthieu Geist 1, Olivier Pietquin,3 1 MaLIS research grou SUPELEC - UMI 958 GeorgiaTech-CRS, France
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationarxiv: v2 [quant-ph] 2 Aug 2012
Qcomiler: quantum comilation with CSD method Y. G. Chen a, J. B. Wang a, a School of Physics, The University of Western Australia, Crawley WA 6009 arxiv:208.094v2 [quant-h] 2 Aug 202 Abstract In this aer,
More informationA Note on Guaranteed Sparse Recovery via l 1 -Minimization
A Note on Guaranteed Sarse Recovery via l -Minimization Simon Foucart, Université Pierre et Marie Curie Abstract It is roved that every s-sarse vector x C N can be recovered from the measurement vector
More informationPrincipal Components Analysis and Unsupervised Hebbian Learning
Princial Comonents Analysis and Unsuervised Hebbian Learning Robert Jacobs Deartment of Brain & Cognitive Sciences University of Rochester Rochester, NY 1467, USA August 8, 008 Reference: Much of the material
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationarxiv: v2 [math.na] 6 Apr 2016
Existence and otimality of strong stability reserving linear multiste methods: a duality-based aroach arxiv:504.03930v [math.na] 6 Ar 06 Adrián Németh January 9, 08 Abstract David I. Ketcheson We rove
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationIntrinsic Approximation on Cantor-like Sets, a Problem of Mahler
Intrinsic Aroximation on Cantor-like Sets, a Problem of Mahler Ryan Broderick, Lior Fishman, Asaf Reich and Barak Weiss July 200 Abstract In 984, Kurt Mahler osed the following fundamental question: How
More informationInference for Empirical Wasserstein Distances on Finite Spaces: Supplementary Material
Inference for Emirical Wasserstein Distances on Finite Saces: Sulementary Material Max Sommerfeld Axel Munk Keywords: otimal transort, Wasserstein distance, central limit theorem, directional Hadamard
More informationSeries Handout A. 1. Determine which of the following sums are geometric. If the sum is geometric, express the sum in closed form.
Series Handout A. Determine which of the following sums are geometric. If the sum is geometric, exress the sum in closed form. 70 a) k= ( k ) b) 50 k= ( k )2 c) 60 k= ( k )k d) 60 k= (.0)k/3 2. Find the
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationReal Analysis 1 Fall Homework 3. a n.
eal Analysis Fall 06 Homework 3. Let and consider the measure sace N, P, µ, where µ is counting measure. That is, if N, then µ equals the number of elements in if is finite; µ = otherwise. One usually
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationRobust Performance Design of PID Controllers with Inverse Multiplicative Uncertainty
American Control Conference on O'Farrell Street San Francisco CA USA June 9 - July Robust Performance Design of PID Controllers with Inverse Multilicative Uncertainty Tooran Emami John M Watkins Senior
More informationA construction of bent functions from plateaued functions
A construction of bent functions from lateaued functions Ayça Çeşmelioğlu, Wilfried Meidl Sabancı University, MDBF, Orhanlı, 34956 Tuzla, İstanbul, Turkey. Abstract In this resentation, a technique for
More informationThe Fekete Szegő theorem with splitting conditions: Part I
ACTA ARITHMETICA XCIII.2 (2000) The Fekete Szegő theorem with slitting conditions: Part I by Robert Rumely (Athens, GA) A classical theorem of Fekete and Szegő [4] says that if E is a comact set in the
More informationExact Solutions in Finite Compressible Elasticity via the Complementary Energy Function
Exact Solutions in Finite Comressible Elasticity via the Comlementary Energy Function Francis Rooney Deartment of Mathematics University of Wisconsin Madison, USA Sean Eberhard Mathematical Institute,
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationHOMEWORK # 4 MARIA SIMBIRSKY SANDY ROGERS MATTHEW WELSH
HOMEWORK # 4 MARIA SIMBIRSKY SANDY ROGERS MATTHEW WELSH 1. Section 2.1, Problems 5, 8, 28, and 48 Problem. 2.1.5 Write a single congruence that is equivalent to the air of congruences x 1 mod 4 and x 2
More informationPOINTS ON CONICS MODULO p
POINTS ON CONICS MODULO TEAM 2: JONGMIN BAEK, ANAND DEOPURKAR, AND KATHERINE REDFIELD Abstract. We comute the number of integer oints on conics modulo, where is an odd rime. We extend our results to conics
More informationEdinburgh Research Explorer
Edinburgh Research Exlorer The shar A) constant for weights in a reverse-holder class Citation for ublished version: Dindos M & Wall T 2009 'The shar A) constant for weights in a reverse-holder class'
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More informationThe Knuth-Yao Quadrangle-Inequality Speedup is a Consequence of Total-Monotonicity
The Knuth-Yao Quadrangle-Ineuality Seedu is a Conseuence of Total-Monotonicity Wolfgang W. Bein Mordecai J. Golin Lawrence L. Larmore Yan Zhang Abstract There exist several general techniues in the literature
More informationON FREIMAN S 2.4-THEOREM
ON FREIMAN S 2.4-THEOREM ØYSTEIN J. RØDSETH Abstract. Gregory Freiman s celebrated 2.4-Theorem says that if A is a set of residue classes modulo a rime satisfying 2A 2.4 A 3 and A < /35, then A is contained
More information1 Riesz Potential and Enbeddings Theorems
Riesz Potential and Enbeddings Theorems Given 0 < < and a function u L loc R, the Riesz otential of u is defined by u y I u x := R x y dy, x R We begin by finding an exonent such that I u L R c u L R for
More informationLandau Theory of the Fermi Liquid
Chater 5 Landau Theory of the Fermi Liquid 5. Adiabatic Continuity The results of the revious lectures, which are based on the hysics of noninteracting systems lus lowest orders in erturbation theory,
More informationOutline. CS21 Decidability and Tractability. Regular expressions and FA. Regular expressions and FA. Regular expressions and FA
Outline CS21 Decidability and Tractability Lecture 4 January 14, 2019 FA and Regular Exressions Non-regular languages: Puming Lemma Pushdown Automata Context-Free Grammars and Languages January 14, 2019
More informationFig. 4. Example of Predicted Workers. Fig. 3. A Framework for Tackling the MQA Problem.
217 IEEE 33rd International Conference on Data Engineering Prediction-Based Task Assignment in Satial Crowdsourcing Peng Cheng #, Xiang Lian, Lei Chen #, Cyrus Shahabi # Hong Kong University of Science
More informationEconometrica Supplementary Material
Econometrica Sulementary Material SUPPLEMENT TO WEAKLY BELIEF-FREE EQUILIBRIA IN REPEATED GAMES WITH PRIVATE MONITORING (Econometrica, Vol. 79, No. 3, May 2011, 877 892) BY KANDORI,MICHIHIRO IN THIS SUPPLEMENT,
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More information