On a geometrical approach in contact mechanics
|
|
- Kenneth Johns
- 5 years ago
- Views:
Transcription
1 Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb Karsruhe Te.: +49 (0) 721/ Fax: +49 (0) 721/ E-Mai: ifm@uni-karsruhe.de
2 On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof 2005 Abstract. The focus of the contribution is on a detaied discussion of the 2D formuation in contact which can be viewed either as a reduction of the genera 3D case to a specia cyindrica geometry, or as the contact of 2D bodies bounded by pane curves. In addition, typica frictiona characteristics, such as the yied surface and the update of the siding dispacements aow a geometrica interpretation in the chosen coordinate system on the contact surface. 1 Geometry and kinematics of contact In the iterature various contact descriptions adapted for an effective finite eement impementation are avaiabe, which can be characterized by the foowing: from 2D to 3D formuations, from non-frictiona to frictiona contact. One of the first contributions on finite eement soutions of 2D frictiona contact probems based on an easto-pastic anaogy has been made by Wriggers et. a. [1]. Genera overviews over contact conditions and contact agorithms which are nowadays used in practice, are covered by the books of Wriggers [2] and Laursen [3]. The covariant description, see Konyukhov and Schweizerhof [4] for the frictioness case and [5] for the frictiona contact, has been found as a universa penaty based approach for contact with various approximation of the surfaces. In the current contribution, we wi show the unity of 2D and 3D formuations, where the 2D case can be derived, from one hand, as a simpified case of the particuar 3D geometry of contact surfaces and, from the other hand, can be constructed separatey based on the differentia geometry of 2D pane curves. Considering a specia contact case contact between two cyindrica infinite bodies with pane strain deformations, see Fig. 1, eads to the definition of a 2D contact. In this case a generatrix GH of the first cyindrica body is a contact ine and corresponds to a contact ine G H which is aso a generatrix but of the second cyindrica body. Thus, 3D contact which can be seen as an interaction between two surfaces is reduced to an interaction between two boundary curves in the 2D case. One of both boundary curves resp. surfaces is chosen as the master contact curve resp. surface. A coordinate system is considered on the boundary, either for a surface in 3D or for a curve in 2D. On the pane we define a curviinear coordinate system associated with the curve by introducing two principa vectors as a basis: the tangent vector ρ ξ = 1
3 Figure 1: Two dimensiona contact as a specia case of three dimensiona contact ρ ξ and the unit norma vector ν. Then a save contact point S is found as r s (ξ, ζ) = ρ(ξ) + ζν(ξ). (1) The norma unit vector ν in the case of arbitrary Lagrangian parameterizations with ξ can be defined via a cross product in a Cartesian coordinate system as: ρ ξ = ρ ξ ; = ν = [k ρ ξ] ρξ ρ ξ, (2) where k is the third unit vector in this Cartesian coordinate system. According to this definition the traction vector R s is defined via the contravariant basis vectors: R s = T ρ 1 + Nρ 2 ρ ξ = T + Nν. (3) (ρ ξ ρ ξ ) The decomposition of the traction in eqn. (3) eads to the foowing contact integra δw c = (Nδζ + T δξ)d. (4) A coser ook reveas that the contact integra (4) contains the work of the contact tractions T and N defined on the master contact curve and is computed aong the save curve s. 2
4 2 Reguarization of the contact tractions For the norma traction N, the foowing reguarized equation in the cosed form is taken N = ɛ N ζ, if ζ 0, (5) where ɛ N is a penaty parameter for the norma interaction. As a reasonabe equation for the reguarization of the tangent traction T we choose a proportiona reation between the fu time derivative dt and the dt reative veocity vector expressed in covariant form on the tangent ine ζ = 0: D 1 T dt ρ1 = ɛ T ξρξ. (6) Expressing the covariant derivative in eqn. (6) the foowing equation is obtained dt dt = ɛ T (ρ ξ ρ ξ ) ξ + ρ ξξ ρ ξ (ρ ξ ρ ξ ) ξ h 11 a 11 ζ, (7) where ɛ T is a penaty parameter for the tangentia interaction, a 11 and h 11 are components of the metric tensor resp. of the curvature tensor for the cyindrica surface. If a path ength s = ξ is used for parameterization, then eqn. (7) is transformed as dt dt = ɛ T ṡ κ ζ, (8) where κ is a curvature of the master contact pane curve. 3 Linearization for a Newton type soution scheme The derivation of the contact matrices can be performed either according to the cyindrica geometry, or according to the geometry of the pane curves. The atter provides a straightforward geometrica expanation of each part of the contact matrix and can be used for the judgment of their necessity within the soution process. Thus, the inearization of the norma part in eqn. (4) given by the foowing integra: δwc N = Nδζd, (9) eads to the foowing. D(δW N c ) = ( dn dt ) dδζ δζ + N d = dt = ɛ N (δr s δρ) (ν ν)(v s v)d (10) 3
5 ( ɛ N ζ δτ (ν τ )(v s v) + (δr s δρ) (τ ν) τ t ) d (10 a) ɛ N ζκ(δr s δρ) (τ τ )(v s v)d. (10 b) The form written via the path ength aows a simpe geometrica interpretation of each part in eqn. (10) and even aows to estabish situations where some of the parts are zero. The first part eqn. (10) is caed main part and defines the constitutive reation for norma contact conditions. The second part eqn. (10 a) is caed rotationa part and defines the geometrica stiffness due to the rotation of the tangent vector of the master curve. It disappears when a master segment is moving in parae, because ony in this case the derivative of a unit vector τ becomes zero. The third part eqn. (10 b) is caed curvature part. This part disappears when the curvature κ of a master segment is zero, i.e. in the case of inear approximations of the master segment. The structure of the tangentia part for the sticking case is as foows: D v (δwc T ) = ( dt dδξ δξ + T dt dt ) d = ε T (δr s δρ) (τ τ )(v s v)d (11) [ T (δr s δρ) (τ τ ) τ ] t + δτ (τ τ ) (v s v) d (11 a) + κ(δr s δρ) (τ ν + ν τ ) (v s v)d. (11 b) Here T is a tria tangent traction computed from the discrete evoution equation (6) at each oad step, e.g. under the assumption that the save point obeys the eastic deformation aw. One can see that the symmetry is preserved for the fu sticking case for any curviinear geometry of the contact bodies. If siding is detected, i.e. if T > µ N, then the siding force is computed according to Couomb s friction aw within the return-mapping scheme, see [2] and [3]. We aso keep a covariant form: T s = µ N T tr T tr = µ N (ρ ξ ρ ξ ) 1/2 sgn(t tr ). (12) The inearized contact integra for the tangentia part in the case of siding has the same geometrica structure, but contains non-symmetric parts due to the non-associativity of the Couomb friction aw. 4
6 Detais for the impementation of the inear contact eement as we as for the computation according to the return-mapping scheme can be found in Konyukhov and Schweizerhof [7]. 4 Treatment of specia cases Some specia cases can appear when a direct appication of the return-mapping scheme can ead to improper resuts. The first probem is arising when the appied oad is not simpy modified proportionay. In this situation a tria oad can not be computed ony via the evoution equation, because the attraction point ξ 0 must be updated. Thus we have to extend the agorithm as is shown in the foowing. The second probem is arising when the projection point is crossing an eement boundary during the incrementa oading. In this case, the computation according to the incrementa evoution equation wi produce a jump, because the convective coordinate ξ beongs to different eements. 4.1 Update of the siding dispacements in the case of reversibe oading A geometrica interpretation of the tria step in the return-mapping scheme eads to the definition of the eastic region with an attraction point ξ 0 T (n) tr < µ N (n) a 11 = ɛ T ξ (n) ξ 0 < µ N (n) (13) If a point ξ (n+1) appears to be outside of the domain at oad step (n + 1), then its ony admissibe position is on the boundary of the domain. A siding force is appied then at the contact point. As ong as we have a motion of the contact point ony in one direction the sign function for the siding force sgn(t (n+1) tr ) = sgn( ξ (n+1) ) does not change and the computation wi be correct. However, when a reversibe oad is appied and it forces the contact point to move forward or backward, the attraction point ξ 0 must be updated in order to define the sign function for the siding force correcty. This update can be defined geometricay from Fig. 2: ξ (n+1) = ξ (n+1) µ N (n). (14) ɛ T Thus, computation of the tria force at the next oad step (n + 1) has to be made in accordance to the update procedure, see more in Konyukhov and Schweizerhof [7]. 5
7 Figure 2: Couomb friction. Updating of siding dispacements in convective coordinates. Motion of friction cone and center of attraction. 4.2 Crossing an eement boundary continuous integration scheme Consider two adjacent eements AB and BC, see Fig. 3. If we foow the computation of the tria force according to the formua expressed in convective coordinates as T (n+1) tr = ɛ T a 11 (ξ (n+1) ξ (n) ), (15) then this eads to a jump in the contact tractions. The maximum of the jump is computed as ( ) T jump = ɛ T a 11 im (ξ) ξ BC im (ξ) ξ AB = 2ɛ T a 11. (16) ξ 1+0 ξ +1 0 This jump appears ony due to the different approximation of the adjacent eements. In order to overcome this, we can compute the force in geometrica form. Foowing this procedure the incrementa coordinate ξ can be expressed as ξ = ( ρ(ξ (n+1) ) ξ BC ( ρ(ξ (n) ) + u(ξ (n) ) ) ) ξ AB, ρ (n+1) ξ a (n+1) 11 (17) where u is an incrementa dispacement vector. The evoution equation becomes then T (n+1) tr = T (n) ɛ T a (n+1) 11 ξ. (18) In the 2D case, the computation can be made via the ength parameter s eading to the continuous scheme as we, see Wriggers [2]. 5 Numerica exampes 5.1 Siding of a bock. Linear approximation of the contact surfaces. Reversibe oading process. We consider the siding of an eastic bock on the rigid base oaded with horizontay prescribed reversibe dispacements, see Fig. 4(a). The geometrica 6
8 Figure 3: Crossing an eement boundary within a oad increment. Typica case for the continuous integration scheme. and mechanica parameters are easticity moduus E = , Poisson ratio ν = 0.3, ength a = 20, height b = 5, Couomb friction coefficient µ = 0.3. The penaty parameters are chosen as ε N = ε T = The main point is to show the update procedure for siding dispacements. The oading is appied as prescribed dispacements at the top side of the deformabe bock. For contact the node-to-segment approach is taken. The hysteresis curve representing the computed horizonta dispacement at point D vs. the appied dispacement at point C is given in Fig. 4(b) G F6 F F4 F3 F2 F 5 1 F h u 2 u 1 B A v C B 1 C 1 D b Observed dispacement (pt. D) H O Appied dispacement (pt. C) O 1 O2 (a) Pane deformation of a bock. Appied dispacement oading at top of the bock. (b) Hysteresis curve. Observed horizonta dispacement at point D vs. appied horizonta dispacement at point C. 5.2 Drawing of an eastic strip into a channe with sharp corners. A particuar exampe in which appication of the continuous integration is absoutey necessary is a deep drawing of an eastic strip into a channe with sharp corners, see Fig. 4. The crucia point during the anaysis is the siding of a sharp corner C over the eement boundaries 1, 2, 3. A oad-dispacement curve computed for the oading point is chosen as the representative parameter to compare various contact approaches. We obtain, see aso [7], that the appication of the segment-to-segment approach, as described in [6], without the 7
9 continuous integration scheme aows to compute the force-dispacement curve ony for non-frictiona cases. For frictiona cases a more carefu transport of the history variabes is necessary as it is suggested here. 6 Concusions Figure 4: Drawing of an eastic strip into a channe with sharp corners. In this contribution a convective description is reconsidered for the 2D quasistatica frictiona contact probem. Specia attention is paid to the derivation of the necessary equations either as a reduction of the known 3D covariant formuation, or directy from the specia 2D cyindrica geometry of the contact surfaces. The agorithmic inearization in the covariant form aows to obtain the tangent matrices before the inearization process. Thus, an impementation can be easiy carried without providing any specia attention to the approximation of the contact surfaces. Further it is shown that specia agorithmic techniques have to be taken to provide robust answerers for oad reversion and for the crossing of eement boundaries. References [1] Wriggers P., Vu Van, Stein E. Finite eement formuation of arge deformation impact-contact probem with friction. Computers and Structures, 37, pp , [2] Wriggers P. Computationa Contact Mechanics. John Wiey & Sons, [3] Laursen T. A. Computationa contact and impact mechanics. Fundamentas of modeing interfacia phenomena in noninear finite eement anaysis. Springer, [4] Konyukhov A., Schweizerhof K. Contact formuation via a veocity description aowing efficiency improvements in frictioness contact anaysis. Computationa Mechanics, 33, pp ,
10 [5] Konyukhov A., Schweizerhof K. Covariant description for frictiona contact probems. Computationa mechanics, 35, 3, pp , [6] Harnau M, Konyukhov A, Schweizerhof K. Agorithmic aspects in arge deformation contact anaysis using Soid-She eements. Computers and Structures, 83, , [7] Konyukhov A., Schweizerhof K. A specia focus on 2D formuations for contact probems using a covariant description. submitted to Internationa Journa for Numerica Methods in Engineering. 9
Nonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationDYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE
3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationPost-buckling behaviour of a slender beam in a circular tube, under axial load
Computationa Metho and Experimenta Measurements XIII 547 Post-bucking behaviour of a sender beam in a circuar tube, under axia oad M. Gh. Munteanu & A. Barraco Transivania University of Brasov, Romania
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More informationUI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE
UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,
More informationСРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS
СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS
More informationPHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I
6 n terms of moment of inertia, equation (7.8) can be written as The vector form of the above equation is...(7.9 a)...(7.9 b) The anguar acceeration produced is aong the direction of appied externa torque.
More informationForces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment
Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This
More informationParallel-Axis Theorem
Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationVibrations of beams with a variable cross-section fixed on rotational rigid disks
1(13) 39 57 Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks Abstract The work is focused on the probem of vibrating beams with a variabe cross-section fixed on a rotationa
More informationUNCOMPLICATED TORSION AND BENDING THEORIES FOR MICROPOLAR ELASTIC BEAMS
11th Word Congress on Computationa Mechanics WCCM XI 5th European Conference on Computationa Mechanics ECCM V 6th European Conference on Computationa Fuid Dynamics ECFD VI E. Oñate J. Oiver and. Huerta
More information2.1. Cantilever The Hooke's law
.1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because
More informationElectromagnetism Spring 2018, NYU
Eectromagnetism Spring 08, NYU March 6, 08 Time-dependent fieds We now consider the two phenomena missing from the static fied case: Faraday s Law of induction and Maxwe s dispacement current. Faraday
More informationStrain Energy in Linear Elastic Solids
Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationSTABILITY ANALYSIS FOR 3D FRAMES USING MIXED COROTATIONAL FORMULATION
SDSS Rio 200 STABIITY AND DUCTIITY OF STEE STRUCTURES E. Batista, P. Veasco,. de ima (Eds.) Rio de Janeiro, Brazi, September 8-0, 200 STABIITY ANAYSIS FOR 3D FRAMES USING MIXED COROTATIONA FORMUATION Rabe
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationDislocations in the Spacetime Continuum: Framework for Quantum Physics
Issue 4 October PROGRESS IN PHYSICS Voume 11 15 Disocations in the Spacetime Continuum: Framework for Quantum Physics Pierre A. Miette PierreAMiette@aumni.uottawa.ca, Ottawa, Canada This paper provides
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationLecture 17 - The Secrets we have Swept Under the Rug
Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics
More informationTransverse Anisotropy in Softwoods
Transverse Anisotropy in Softwoods Modeing and Experiments CARL MODÉN Licenciate Thesis Stockhom, Sweden 2006 TRITA-AVE 2006:30 ISSN 1651-7660 ISBN 91-7178-385-7 KTH Engineering Sciences Department of
More informationKeywords: Functionally Graded Materials, Conical shell, Rayleigh-Ritz Method, Energy Functional, Vibration.
Journa of American Science, ;8(3) Comparison of wo Kinds of Functionay Graded Conica Shes with Various Gradient Index for Vibration Anaysis Amirhossein Nezhadi *, Rosan Abdu Rahman, Amran Ayob Facuty of
More informationModule 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur
odue 2 naysis of Staticay ndeterminate Structures by the atri Force ethod Version 2 E T, Kharagpur esson 12 The Three-oment Equations- Version 2 E T, Kharagpur nstructiona Objectives fter reading this
More informationTHE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE
THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationModeling Friction Phenomena in Flexible Multibody Dynamics
Modeing Friction Phenomena in Fexibe Mutibody Dynamics Oivier A. Bauchau and Changuan Ju Danie Guggenheim Schoo of Aerospace Engineering, Georgia Institute of Technoogy, 270 Ferst Dr., Atanta, GA 30332,
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationMECHANICAL ENGINEERING
1 SSC-JE SFF SELECION COMMISSION MECHNICL ENGINEERING SUDY MERIL Cassroom Posta Correspondence est-series16 Rights Reserved www.sscje.com C O N E N 1. SIMPLE SRESSES ND SRINS 3-3. PRINCIPL SRESS ND SRIN
More informationSIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT
8 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT P. Wang, N. Hamia *, P. Boisse Universite de Lyon, INSA-Lyon,
More informationDynamic equations for curved submerged floating tunnel
Appied Mathematics and Mechanics Engish Edition, 7, 8:99 38 c Editoria Committee of App. Math. Mech., ISSN 53-487 Dynamic equations for curved submerged foating tunne DONG Man-sheng, GE Fei, ZHANG Shuang-yin,
More informationThe Bending of Rectangular Deep Beams with Fixed at Both Ends under Uniform Load
Engineering,,, 8-9 doi:.6/eng..7 Pubised Onine December (ttp://.scirp.org/journa/eng) Te Bending of Rectanguar Deep Beams it Fied at Bot Ends under Uniform Load Abstract Ying-Jie Cen, Bao-Lian Fu, Gang
More informationBending Analysis of Continuous Castellated Beams
Bending Anaysis of Continuous Casteated Beams * Sahar Eaiwi 1), Boksun Kim ) and Long-yuan Li 3) 1), ), 3) Schoo of Engineering, Pymouth University, Drake Circus, Pymouth, UK PL4 8AA 1) sahar.eaiwi@pymouth.ac.uk
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationPhysicsAndMathsTutor.com
. Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationPhase Change Equation of State for FSI Applications
15 th Internationa LS-DYNA Users Conference FSI / ALE Phase Change Equation of State for FSI Appications Mhamed Soui, Ramzi Messahe Lie Uniersity France Cyri Regan, Camie Ruiuc Ingeiance Technoogies, Agence
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationHigh Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings
Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationReliability: Theory & Applications No.3, September 2006
REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationSECTION A. Question 1
SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationTheory and implementation behind: Universal surface creation - smallest unitcell
Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationPath planning with PH G2 splines in R2
Path panning with PH G2 spines in R2 Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru To cite this version: Laurent Gajny, Richard Béarée, Eric Nyiri, Oivier Gibaru. Path panning with PH G2 spines
More informationWork and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame
Work and energy method EI EI T x-axis Exercise 1 : Beam with a coupe Determine the rotation at the right support of the construction dispayed on the right, caused by the coupe T using Castigiano s nd theorem.
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationSlender Structures Load carrying principles
Sender Structures Load carrying principes Cabes and arches v018-1 ans Weeman 1 Content (preiminary schedue) Basic cases Extension, shear, torsion, cabe Bending (Euer-Bernoui) Combined systems - Parae systems
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationAbout the Torsional Constant for thin-walled rod with open. cross-section. Duan Jin1,a, Li Yun-gui1
Internationa Forum on Energy, Environment Science and Materias (IFEESM 17) bout the Torsiona Constant for thin-waed rod with open cross-section Duan Jin1,a, Li Yun-gui1 1 China State Construction Technica
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationTHE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia
More informationFlexural to shear and crushing failure transitions in RC beams by the bridged crack model
exura to shear and crushing faiure transitions in RC beams by the bridged crack mode A. Carpinteri & G. Ventura Poitecnico di Torino, Torino, Itay J. R. Carmona Universidad de Castia-La Mancha, Ciudad
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationNumerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process
ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More informationLobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z
Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationXI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com
XI PHYSICS M. Affan Khan LECTURER PHYSICS, AKHSS, K affan_414@ive.com https://promotephysics.wordpress.com [TORQUE, ANGULAR MOMENTUM & EQUILIBRIUM] CHAPTER NO. 5 Okay here we are going to discuss Rotationa
More informationPhysics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27
Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that
More informationLECTURE 10. The world of pendula
LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive
More informationMaryam Tabatabaei, Dy Le, Satya N. Atluri
Neary Exact and Highy Efficient Eastic-Pastic Homogenization and/or Direct Numerica Simuation of Low-Mass Metaic Systems with Architected Ceuar Microstructures Maryam Tabatabaei, Dy Le, Satya N. Aturi
More informationModal analysis of a multi-blade system undergoing rotational motion
Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationNumerical simulation of javelin best throwing angle based on biomechanical model
ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationFRIEZE GROUPS IN R 2
FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the
More informationAAPT UNITED STATES PHYSICS TEAM AIP 2012 DO NOT DISTRIBUTE THIS PAGE
2012 Semifina Exam 1 AAPT UNITED STATES PHYSICS TEAM AIP 2012 Semifina Exam DO NOT DISTRIBUTE THIS PAGE Important Instructions for the Exam Supervisor This examination consists of two parts. Part A has
More information9. EXERCISES ON THE FINITE-ELEMENT METHOD
9. EXERCISES O THE FIITE-ELEMET METHOD Exercise Thickness: t=; Pane strain proem (ν 0): Surface oad Voume oad; 4 p f ( x, ) ( x ) 0 E D 0 0 0 ( ) 4 p F( xy, ) Interna constrain: rigid rod etween D and
More information1 Equations of Motion 3: Equivalent System Method
8 Mechanica Vibrations Equations of Motion : Equivaent System Method In systems in which masses are joined by rigid ins, evers, or gears and in some distributed systems, various springs, dampers, and masses
More information830. Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation
8. Noninear dynamic characteristics of SMA simpy supported beam in axia stochastic excitation Zhi-Wen Zhu 1, Wen-Ya Xie, Jia Xu 1 Schoo of Mechanica Engineering, Tianjin University, 9 Weijin Road, Tianjin
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationAnalysis of Ball-Type Constant-Velocity Joints BasedonDynamics
736 Anaysis of Ba-Type Constant-Veocity Joints BasedonDynamics Kei KIMATA, Haruo NAGATANI and Masayuki IMOTO Universa joints, which transmit torque through the bas guided in such a manner that they aways
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationA PROCEDURE ON STUDYING A CRITICAL FORCE ON FEM AS PLASTIC DEFORMATIONS ARISE WITHIN STRAIGHT TRUSSES
11 th Nationa Congress on Theoretica and Appied Mechanics, -5 Sept. 009, Borovets, Bugaria A OCDU ON STUDYING A CITICAL FOC ON FM AS LASTIC DFOMATIONS AIS WITHIN STAIGHT TUSSS IVAN VAISILOV VSU L. Karaveov,
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationA 3-D hybrid finite-difference finite-element viscoelastic modelling of seismic wave motion
The definitive version is avaiabe at www.backwe-synergy.com http://www.interscience.wiey.com/journa/4007/abstract Geophys. J. Int. 008) 75, 5 84 doi: 0./j.65-46X.008.0866.x A -D hybrid finite-difference
More informationDisplacement Solution for a Static Vertical Rigid Movement of an Interior Circular Disc in a Transversely Isotropic Tri-Material Full-Space
Vo:8, No:, 4 Dispacement Soution for a Static Vertica Rigid Movement of an Interior Circuar Disc in a Transversey Isotropic Tri-Materia Fu-Space D. Mehdizadeh, M. Rahimian, M. Eskandari-Ghadi Internationa
More informationContribution to a formulation of integral type of the mechanics of Cosserat continua ( * ).
ontributo per una formuazione di tipo integrae dea meccanica dei continui di osserat Ann Mat pura ed App (976) 75-83 ontribution to a formuation of integra type of the mechanics of osserat continua ( *
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationAPPENDIX C FLEXING OF LENGTH BARS
Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More information