Computational Methods of Solid Mechanics. Project report
|
|
- David Booth
- 5 years ago
- Views:
Transcription
1 Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng
2 Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate the adhesive contact between a soft tip with ough suface and a igid substate, see Fig. (a). Since the tip size is much smalle than the substate, we take the substate to be semi-infinite. The pofile of the ough suface of tip in the efeence configuation is given by [ ( )] x 2 (x ) = x2 2R + A 2πx cos () λ whee the fist tem is a paabola shape with maximum cuvatue R and the second tem is a sinusoid undulation with wavelength λ and amplitude A. We take λ, A R to epesent the small scale oughness of the suface. The tip mateial is descibed by Neo-Hookean model as discussed in the class. We model the adhesive contact behavio by assuming that the tip and substates mateial points inteact with one anothe though the Lennad-Jones type inteaction potential [ ( σ ) 6 ( σ ) ] 2 V LJ = 4ɛ + (2) whee ɛ and σ ae constants and have units of enegy and length espectively, is the distance between two points. The inteaction foce between the tip and substate mateial point is obtained by diffeentiating the inteaction potential. Since the substate is semi-infinite lage and igid, the body foce acted on the tip mateial point and the wok of adhesion of the mateial can be deived analytically. Due to symmety, the only nonzeo component of body foce is [ ( ) (y) = 4πɛσ 2 σ ρ s ρ t (y) ( ) ] 4 σ (3) 5 y 2 2 y 2 as shown in Fig. (b), whee ρ t and ρ s ae the tip and substate densities. The wok of adhesion is w = 4 ( ) /3 2 πɛρ s ρ t σ 4 (4) 5 Note that the body foce is descibed in the defomed configuation. Rewite eq. 3 in the efeence configuation we have [ ( ) (x) = 4πɛσ 2 ρ s ρ t J σ (x) ( ) ] 4 σ (5) 5 x 2 + u 2 2 x 2 + u 2 whee u 2 is the displacement in ê 2 diection. 2 Govening equations The equilibium equation in the defomed configuation is σ ij y j + b i = in V σ ij n j = t i on S 2 (6) u i = u i on S whee σ ij is the Cauchy stess. Let η i be the test function with η i = u i on S. The weak fom is η i σ ij dv b i η i dv t i η i da = (7) y j V S 2 V 2
3 (a) e (b) 2 y 2 e.5 y Figue : (a) The schematic of adhesive contact between the soft indente and semi-infinite igid substate and (b) the Lennad-Jones type body foce. Map the integal ove the defomed configuation to the efeence configuation we have η i τ ij d b i η i d t i η i da = (8) y j whee τ ij = J σ ij is the Kichoff stess, b i = J b i, and t i is the taction mapped back to the efeence configuation. In the poject, the body foce is dependent on the position as descibed by eq. 5, and fo simplicity, the taction bounday condition is not consideed. Hence the thid tem on the left hand side of the weak fom is zeo. S 2 3 FE implementation Intoduce FE intepolation Substitute into the weak fom we obtain whee u i = N a (x)u a i, η i = N a (x)η a i (9) R a i [u b k] + B a i [u b k] = F a i () Ri a = τ ij [B kl ] V N a d y j Bi a = b i N a d Fi a = t i N a da S 2 Note that the foce vecto Bi a is now a function of displacement since the body foce is nonunifom distibuted given by [ ( ) b i = 4πɛσ 2 ρ s ρ σ t δ i2 ( ) ] 4 σ (2) No summation of i hee. The nonlinea equation will be solved using Newton-Raphson iteation. Stat with initial guess wi a and ty to coect the guess with coection dwa i. Lineaizing eq. in dwa i yields a system of linea equation (Kaibk R + Kaibk)dw F k b = Ri a Bi a + Fi a (3) () 3
4 to solve fo dwi a, whee the tem KR aibk is the usual one we deived in the class. The new tem KB aibk appeaing in the stiffness matix is Kaibk B = Ba i b i u b = k u b N a d (4) k Eveything is famila except the new tem we need to take cae of. The FE pocedue will be implemented with UEL in ABAQUS. Let us denote f(y i ) = 5 ( σ y i ) ( ) 4 σ (5) 2 y i To avoid f(y i ) goes to infinity as y i =, we can set the body foce to the following fom: { b δi2 c f(y i ), y i i = δ i2 c [h + k (y i )], y i < (6) whee c = 4πɛσ 2 ρ s ρ t, is a vey small positive paamete. Fo example, we can take =.5σ, and h = f yi= = k = df yi= = 432 dy i σ (7) Theefoe, fo y i, B a i [ ( ) σ = δ i2 c ( ) ] 4 σ N a d (8) and Fo y i <, and [ Kaibk B = δ 2k c B a i 2σ 4 (x i + u i ) 5 2σ (x i + u i ) ] N a N b d (9) = δ i2 c [h + k (y i )] N a d (2) Kaibk B = δ 2k c k N a N b d (2) 4 Benchmaks In this section, we veify the UEL suboutine by seveal simple test cases. In cases -3, the simple body foces (Fig. 2(b-d)) ae applied to a single element ( in size), which is specified with a igid body motion in the -ê 2 diection (Fig. 2(a)). The body foces have the fom Case : =., > Case 2: =., < <., <, > Case 3: =, 2 < < 3., < < 2., < (22) 4
5 The paamete is set to be.. The total eaction foce of the element, fo the thee cases, is shown in Fig. 3, which is the same as the exact solution. The total foce is negative since a positive body foce is applied. (a) 4 3 (b) v e 2 2 e (c) (d) 3 - Figue 2: Schematic of element test with thee simple body foces..8 case case 2 case 3 Total eaction foce Nomalized displacement Figue 3: The eaction foce F 2 of the element with igid body motion in test cases -3. In test case 4, we compae the Hetz contact simulation of the body foce method with that in ABAQUS/CAE. The contact of a elastic indente of adius 2.5 with igid substate is simulated. A indent depth of.2 is applied. The Young s modulus E and Poisson s atio µ ae set to be 2. and.3. Theoetically, the contacting suface cannot penetating each othe. Hypothetically, it can be viewed as a infinite foce is 5
6 applied on these suface. Theefoe, the body foce has the fom = {, > k, < (23) whee k is a lage numbe. (a) (b) (c) (d) Figue 4: The compaison of u 2 and σ 22 of ABAQUS/CAE (a,c) and body foce method fo k = 2 (b,d) in the Hetz contact simulation Hetz contact Body foce k= 2 Body foce k= 8 Body foce k= 2 Reaction foce 8 6 k Indentation depth Figue 5: The compaison of foce-depth cuves of Hetz contact of ABAQUS/CAE and body foce method. 6
7 Fig. 4 shows the compaison of u 2 and σ 22 plots calculated fom ABAQUS/CAE and the body foce method with k = 2. As we can see, the distibution of u 2 and σ 22 is the same, but the magnitude is diffeent. The esult of body foce method conveges to that of ABAQUS/CAE as the slope k deceases. This is eflected in the compaison of foce-depth cuves, as shown in Fig. 5. As the slope k deceases fom -2 to -2, the diffeence of foce-depth cuve of ABAQUS/CAE and body foce method becomes smalle and smalle. 5 Adhesive contact In this section, we study the adhesive contact though the method of Lennad-Jones type body foce, as descied in Sec.. As pedicted by JKR theoy, thee exist pull-in and pull-off instabilities duing the loading and unloading of the tip due to the adhesive inteactions. To ovecome the convegence poblem in the numeical simulation, a viscosity tem is added to the body foce, i.e., [ ( ) b σ i = δ i2 c ( ) ] 4 ( σ + η du ) i. (24) dt whee η is a small numbe. Accodingly, the foce vecto Bi a and stiffness matix Kaibk B become [ ( ) Bi a σ = δ i2 c ( ) ] 4 ( σ + η u ) i N a d (25) t and [ 2σ 4 (x i + u i ) 5 2σ (x i + u i ) [ ( ) σ δ 2k c ( ) ] 4 σ Kaibk B = δ 2k c ] ( + η u i t ) N a N b d η t N a N b d In the simulation, a 2D tip is moving in the -ê 2 diection and subjected to the body foce as descibed by Eq. 24. The paametes σ and η ae taken to be nm and., espectively. The paametes and k ae taken to be.5σ and. The mateial is descibed by Neo-Hookean model. The shea modulus is.5 MPa. The wok of adhesion of the mateial is mj/m 2. Equivalently, the constant c is taken to be The foce-displacement cuve duing the loading is shown in Fig. 6. As we can see, as the the foce deceases suddenly at a citical point when the pull-in instability occus. Then the foce inceases as the indent depth inceases. The contou plot of σ 22 afte the pull-in instability is shown in Fig. 7. Remak: The convegence is a big issue in this poject. The mesh should be vey efined to be able to esolve the body foce field. Although the pull-in instability is captued in the loading pocess, the computation does not convege duing the unloading pocess, hence the pull-off instability is not shown in the esult. A moe plausible way is to do the dynamic simulation by using VUEL suboutine. (26) 7
8 4 x Foce (un) Displacement (um) 2 Figue 6: The foce-displacement cuve of the adhesive contact. Figue 7: The contou plot of σ 22 afte the pull-in instability. 8
Numerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationFE FORMULATIONS FOR PLASTICITY
G These slides ae designed based on the book: Finite Elements in Plasticity Theoy and Pactice, D.R.J. Owen and E. Hinton, 970, Pineidge Pess Ltd., Swansea, UK. Couse Content: A INTRODUCTION AND OVERVIEW
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationyou of a spring. The potential energy for a spring is given by the parabola U( x)
Small oscillations The theoy of small oscillations is an extemely impotant topic in mechanics. Conside a system that has a potential enegy diagam as below: U B C A x Thee ae thee points of stable equilibium,
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationStress, Cauchy s equation and the Navier-Stokes equations
Chapte 3 Stess, Cauchy s equation and the Navie-Stokes equations 3. The concept of taction/stess Conside the volume of fluid shown in the left half of Fig. 3.. The volume of fluid is subjected to distibuted
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationSupplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in
Supplementay Figue 1. Cicula paallel lamellae gain size as a function of annealing time at 50 C. Eo bas epesent the σ uncetainty in the measued adii based on image pixilation and analysis uncetainty contibutions
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More information, the tangent line is an approximation of the curve (and easier to deal with than the curve).
114 Tangent Planes and Linea Appoimations Back in-dimensions, what was the equation of the tangent line of f ( ) at point (, ) f ( )? (, ) ( )( ) = f Linea Appoimation (Tangent Line Appoimation) of f at
More informationMomentum is conserved if no external force
Goals: Lectue 13 Chapte 9 v Employ consevation of momentum in 1 D & 2D v Examine foces ove time (aka Impulse) Chapte 10 v Undestand the elationship between motion and enegy Assignments: l HW5, due tomoow
More informationHW Solutions # MIT - Prof. Please study example 12.5 "from the earth to the moon". 2GmA v esc
HW Solutions # 11-8.01 MIT - Pof. Kowalski Univesal Gavity. 1) 12.23 Escaping Fom Asteoid Please study example 12.5 "fom the eath to the moon". a) The escape velocity deived in the example (fom enegy consevation)
More informationEM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1)
EM- Coulomb s law, electic field, potential field, supeposition q ' Electic field of a point chage ( ') E( ) kq, whee k / 4 () ' Foce of q on a test chage e at position is ee( ) Electic potential O kq
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More informationis the instantaneous position vector of any grid point or fluid
Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in
More informationThermo-Mechanical Model for Wheel Rail Contact using Coupled. Point Contact Elements
IM214 28-3 th July, ambidge, England hemo-mechanical Model fo heel Rail ontact using oupled Point ontact Elements *J. Neuhaus¹ and. Sexto 1 1 hai of Mechatonics and Dynamics, Univesity of Padebon, Pohlweg
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to
More informationPES 3950/PHYS 6950: Homework Assignment 6
PES 3950/PHYS 6950: Homewok Assignment 6 Handed out: Monday Apil 7 Due in: Wednesday May 6, at the stat of class at 3:05 pm shap Show all woking and easoning to eceive full points. Question 1 [5 points]
More informationLiquid gas interface under hydrostatic pressure
Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationStress Intensity Factor
S 47 Factue Mechanics http://imechanicaog/node/7448 Zhigang Suo Stess Intensity Facto We have modeled a body by using the linea elastic theoy We have modeled a cack in the body by a flat plane, and the
More informationFracture Mechanics: Linear Elastic Fracture Mechanics 2/2
Task 6 - Safety Review and Licensing On the Job Taining on Stess Analysis Factue Mechanics: Linea Elastic Factue Mechanics / Davide Mazzini Cio Santus Pisa (taly) June 15 July 14, 015 Table of content
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationME 210 Applied Mathematics for Mechanical Engineers
Tangent and Ac Length of a Cuve The tangent to a cuve C at a point A on it is defined as the limiting position of the staight line L though A and B, as B appoaches A along the cuve as illustated in the
More information1 Similarity Analysis
ME43A/538A/538B Axisymmetic Tubulent Jet 9 Novembe 28 Similaity Analysis. Intoduction Conside the sketch of an axisymmetic, tubulent jet in Figue. Assume that measuements of the downsteam aveage axial
More informationA dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids
APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing
More informationINFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract
INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS Mohammad Mohammadi, National Cente fo Physical Acoustics, Univesity of Mississippi, MS Caig J. Hicey, National Cente fo Physical Acoustics,
More information10.2 Parametric Calculus
10. Paametic Calculus Let s now tun ou attention to figuing out how to do all that good calculus stuff with a paametically defined function. As a woking eample, let s conside the cuve taced out by a point
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationLINEAR PLATE BENDING
LINEAR PLATE BENDING 1 Linea plate bending A plate is a body of which the mateial is located in a small egion aound a suface in the thee-dimensional space. A special suface is the mid-plane. Measued fom
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationLINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE
LINEAR AND NONLINEAR ANALYSES O A WIND-TUNNEL INTRODUCTION BALANCE R. Kakehabadi and R. D. Rhew NASA LaRC, Hampton, VA The NASA Langley Reseach Cente (LaRC) has been designing stain-gauge balances fo utilization
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More information2. Plane Elasticity Problems
S0 Solid Mechanics Fall 009. Plane lasticity Poblems Main Refeence: Theoy of lasticity by S.P. Timoshenko and J.N. Goodie McGaw-Hill New Yok. Chaptes 3..1 The plane-stess poblem A thin sheet of an isotopic
More informationChapter Introduction to Finite Element Methods
Chapte 1.4 Intoduction to Finite Element Methods Afte eading this chapte, you should e ale to: 1. Undestand the asics of finite element methods using a one-dimensional polem. In the last fifty yeas, the
More information17.1 Electric Potential Energy. Equipotential Lines. PE = energy associated with an arrangement of objects that exert forces on each other
Electic Potential Enegy, PE Units: Joules Electic Potential, Units: olts 17.1 Electic Potential Enegy Electic foce is a consevative foce and so we can assign an electic potential enegy (PE) to the system
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationMagnetic Dipoles Challenge Problem Solutions
Magnetic Dipoles Challenge Poblem Solutions Poblem 1: Cicle the coect answe. Conside a tiangula loop of wie with sides a and b. The loop caies a cuent I in the diection shown, and is placed in a unifom
More informationNumerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.
Numeical solution of diffusion mass tansfe model in adsoption systems Pof., D.Sc. Agenda Mass Tansfe Mechanisms Diffusion Mass Tansfe Models Solving Diffusion Mass Tansfe Models Paamete Estimation 2 Mass
More informationCHAPTER 25 ELECTRIC POTENTIAL
CHPTE 5 ELECTIC POTENTIL Potential Diffeence and Electic Potential Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic foce on the paticle given by F=E. When
More informationAH Mechanics Checklist (Unit 2) AH Mechanics Checklist (Unit 2) Circular Motion
AH Mechanics Checklist (Unit ) AH Mechanics Checklist (Unit ) Cicula Motion No. kill Done 1 Know that cicula motion efes to motion in a cicle of constant adius Know that cicula motion is conveniently descibed
More information7.2.1 Basic relations for Torsion of Circular Members
Section 7. 7. osion In this section, the geomety to be consideed is that of a long slende cicula ba and the load is one which twists the ba. Such poblems ae impotant in the analysis of twisting components,
More informationMolecular dynamics simulation of ultrafast laser ablation of fused silica
IOP Publishing Jounal of Physics: Confeence Seies 59 (27) 1 14 doi:1.188/1742-6596/59/1/22 Eighth Intenational Confeence on Lase Ablation Molecula dynamics simulation of ultafast lase ablation of fused
More informationGauss Law. Physics 231 Lecture 2-1
Gauss Law Physics 31 Lectue -1 lectic Field Lines The numbe of field lines, also known as lines of foce, ae elated to stength of the electic field Moe appopiately it is the numbe of field lines cossing
More information(read nabla or del) is defined by, k. (9.7.1*)
9.7 Gadient of a scala field. Diectional deivative Some of the vecto fields in applications can be obtained fom scala fields. This is vey advantageous because scala fields can be handled moe easily. The
More informationDesigning a Sine-Coil for Measurement of Plasma Displacements in IR-T1 Tokamak
Designing a Sine-Coil fo Measuement of Plasma Displacements in IR-T Tokamak Pejman Khoshid, M. Razavi, M. Ghoanneviss, M. Molaii, A. TalebiTahe, R. Avin, S. Mohammadi and A. NikMohammadi Dept. of Physics,
More informationTHERMODYNAMICS OF SURFACES AND INTERFACES
THERMODYNAMIC OF URFACE AND INTERFACE 1. Intoduction Eveything has to end somewhee. Fo solids, o liquids that "somewhee" is a suface, o an inteface between phases. Fo liquids, the inteface is between the
More informationThree dimensional flow analysis in Axial Flow Compressors
1 Thee dimensional flow analysis in Axial Flow Compessos 2 The ealie assumption on blade flow theoies that the flow inside the axial flow compesso annulus is two dimensional means that adial movement of
More informationMagneto-Elastic Analysis of an Annular FGM Plate Based on Classical Plate Theory Using GDQ Method
736 Magneto-Elastic Analysis of an Annula FGM Plate Based on Classical Plate Theoy Using GDQ Method Abstact Using GDQ method, the adial and cicumfeential stesses in an annula FGM plate with a unifom thickness
More informationPhysics 235 Chapter 5. Chapter 5 Gravitation
Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus
More information1) Consider an object of a parabolic shape with rotational symmetry z
Umeå Univesitet, Fysik 1 Vitaly Bychkov Pov i teknisk fysik, Fluid Mechanics (Stömningsläa), 01-06-01, kl 9.00-15.00 jälpmedel: Students may use any book including the tetbook Lectues on Fluid Dynamics.
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationVoltage ( = Electric Potential )
V-1 of 10 Voltage ( = lectic Potential ) An electic chage altes the space aound it. Thoughout the space aound evey chage is a vecto thing called the electic field. Also filling the space aound evey chage
More informationMONTE CARLO SIMULATION OF FLUID FLOW
MONTE CARLO SIMULATION OF FLUID FLOW M. Ragheb 3/7/3 INTRODUCTION We conside the situation of Fee Molecula Collisionless and Reflective Flow. Collisionless flows occu in the field of aefied gas dynamics.
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More informationPhysics 2212 GH Quiz #2 Solutions Spring 2016
Physics 2212 GH Quiz #2 Solutions Sping 216 I. 17 points) Thee point chages, each caying a chage Q = +6. nc, ae placed on an equilateal tiangle of side length = 3. mm. An additional point chage, caying
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationSTUDY ON 2-D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING
Study Rev. Adv. on -D Mate. shock Sci. wave 33 (13) pessue 111-118 model in mico scale lase shock peening 111 STUDY ON -D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING Y.J. Fan 1, J.Z. Zhou,
More informationNewton s Laws, Kepler s Laws, and Planetary Orbits
Newton s Laws, Keple s Laws, and Planetay Obits PROBLEM SET 4 DUE TUESDAY AT START OF LECTURE 28 Septembe 2017 ASTRONOMY 111 FALL 2017 1 Newton s & Keple s laws and planetay obits Unifom cicula motion
More informationMobility of atoms and diffusion. Einstein relation.
Mobility of atoms and diffusion. Einstein elation. In M simulation we can descibe the mobility of atoms though the mean squae displacement that can be calculated as N 1 MS ( t ( i ( t i ( 0 N The MS contains
More informationMSE 561, Atomic Modeling in Material Science Assignment 1
Depatment of Mateial Science and Engineeing, Univesity of Pennsylvania MSE 561, Atomic Modeling in Mateial Science Assignment 1 Yang Lu 1. Analytical Solution The close-packed two-dimensional stuctue is
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationA NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM
Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496
More informationTopic 4a Introduction to Root Finding & Bracketing Methods
/8/18 Couse Instucto D. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: cumpf@utep.edu Topic 4a Intoduction to Root Finding & Backeting Methods EE 4386/531 Computational Methods in EE Outline
More informationThis gives rise to the separable equation dr/r = 2 cot θ dθ which may be integrated to yield r(θ) = R sin 2 θ (3)
Physics 506 Winte 2008 Homewok Assignment #10 Solutions Textbook poblems: Ch. 12: 12.10, 12.13, 12.16, 12.19 12.10 A chaged paticle finds itself instantaneously in the equatoial plane of the eath s magnetic
More informationTwo Dimensional Inertial Flow of a Viscous Fluid in a Corner
Applied Mathematical Sciences, Vol., 207, no. 9, 407-424 HIKARI Ltd, www.m-hikai.com https://doi.og/0.2988/ams.207.62282 Two Dimensional Inetial Flow of a Viscous Fluid in a Cone A. Mahmood and A.M. Siddiqui
More informationDiffusion and Transport. 10. Friction and the Langevin Equation. Langevin Equation. f d. f ext. f () t f () t. Then Newton s second law is ma f f f t.
Diffusion and Tanspot 10. Fiction and the Langevin Equation Now let s elate the phenomena of ownian motion and diffusion to the concept of fiction, i.e., the esistance to movement that the paticle in the
More informationElectrostatics (Electric Charges and Field) #2 2010
Electic Field: The concept of electic field explains the action at a distance foce between two chaged paticles. Evey chage poduces a field aound it so that any othe chaged paticle expeiences a foce when
More information$ i. !((( dv vol. Physics 8.02 Quiz One Equations Fall q 1 q 2 r 2 C = 2 C! V 2 = Q 2 2C F = 4!" or. r ˆ = points from source q to observer
Physics 8.0 Quiz One Equations Fall 006 F = 1 4" o q 1 q = q q ˆ 3 4" o = E 4" o ˆ = points fom souce q to obseve 1 dq E = # ˆ 4" 0 V "## E "d A = Q inside closed suface o d A points fom inside to V =
More informationSolution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere
Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic
More informationDEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS
DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS TSOKOS LESSON 10-1 DESCRIBING FIELDS Essential Idea: Electic chages and masses each influence the space aound them and that influence can be epesented
More informationLevitation force analysis of ring and disk shaped permanent magnet-high temperature superconductor
Inn Jounal of Pue & Applied Physics Vol. 55, Apil 017, pp. 61-68 Levitation foce analysis of ing and disk shaped pemanent magnet-high tempeatue supeconducto Sinan Basaan & Selim Sivioglu* Depatment of
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.07: Electomagnetism II Septembe 5, 202 Pof. Alan Guth PROBLEM SET 2 DUE DATE: Monday, Septembe 24, 202. Eithe hand it in at the lectue,
More informationPhys-272 Lecture 17. Motional Electromotive Force (emf) Induced Electric Fields Displacement Currents Maxwell s Equations
Phys-7 Lectue 17 Motional Electomotive Foce (emf) Induced Electic Fields Displacement Cuents Maxwell s Equations Fom Faaday's Law to Displacement Cuent AC geneato Magnetic Levitation Tain Review of Souces
More informationFalls in the realm of a body force. Newton s law of gravitation is:
GRAVITATION Falls in the ealm of a body foce. Newton s law of avitation is: F GMm = Applies to '' masses M, (between thei centes) and m. is =. diectional distance between the two masses Let ˆ, thus F =
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationPhysics 181. Assignment 4
Physics 181 Assignment 4 Solutions 1. A sphee has within it a gavitational field given by g = g, whee g is constant and is the position vecto of the field point elative to the cente of the sphee. This
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationMATH 415, WEEK 3: Parameter-Dependence and Bifurcations
MATH 415, WEEK 3: Paamete-Dependence and Bifucations 1 A Note on Paamete Dependence We should pause to make a bief note about the ole played in the study of dynamical systems by the system s paametes.
More informationA scaling-up methodology for co-rotating twin-screw extruders
A scaling-up methodology fo co-otating twin-scew extudes A. Gaspa-Cunha, J. A. Covas Institute fo Polymes and Composites/I3N, Univesity of Minho, Guimaães 4800-058, Potugal Abstact. Scaling-up of co-otating
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationPoint defects in two-dimensional colloidal crystals: simulation vs. elasticity theory.
Point defects in two-dimensional colloidal cystals: simulation vs. elasticity theoy. Wolfgang Lechne Univesity of Vienna, Faculty of Physics, Boltzmanngasse 5, Vienna, Austia. E-mail: wolfgang.lechne@univie.ac.at
More informationCOUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION
ENOC 008 Saint Petesbug Russia June 30-July 4 008 COUPLED MODELS OF ROLLING SLIDING AND WHIRLING FRICTION Alexey Kieenkov Ins ti tu te fo P ob le ms in Me ch an ic s Ru ss ia n Ac ad em y of Sc ie nc es
More information4.1 The state of stress at a certain point in a body is given by:[ T ] =
CHARTER 4 4. The state of stess at a cetain point in a body is given by:[ T ] 4 5. 5 On each of the coodinate planes (with nomal in e, e, ediections), (a) what is the nomal stess and (b) what is the total
More informationMath 2263 Solutions for Spring 2003 Final Exam
Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate
More informationMotion in One Dimension
Motion in One Dimension Intoduction: In this lab, you will investigate the motion of a olling cat as it tavels in a staight line. Although this setup may seem ovesimplified, you will soon see that a detailed
More informationMathematical Analysis and Numerical Simulation of High Frequency Electromagnetic Field in Soft Contact Continuous Casting Mold
, pp. 974 981 Mathematical Analysis and Numeical Simulation of High Fequency Electomagnetic Field in Soft Contact Continuous Casting Mold Xianzhao NA, Xingzhong ZHANG and Yong GAN National Engineeing &
More informationPart V: Closed-form solutions to Loop Closure Equations
Pat V: Closed-fom solutions to Loop Closue Equations This section will eview the closed-fom solutions techniques fo loop closue equations. The following thee cases will be consideed. ) Two unknown angles
More informationCHAPTER 10 ELECTRIC POTENTIAL AND CAPACITANCE
CHAPTER 0 ELECTRIC POTENTIAL AND CAPACITANCE ELECTRIC POTENTIAL AND CAPACITANCE 7 0. ELECTRIC POTENTIAL ENERGY Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic
More informationTRANSILVANIA UNIVERSITY OF BRASOV MECHANICAL ENGINEERING FACULTY DEPARTMENT OF MECHANICAL ENGINEERING ONLY FOR STUDENTS
TNSILVNI UNIVSITY OF BSOV CHNICL NGINING FCULTY DPTNT OF CHNICL NGINING Couse 9 Cuved as 9.. Intoduction The eams with plane o spatial cuved longitudinal axes ae called cuved as. Thee ae consideed two
More informationStatic equilibrium requires a balance of forces and a balance of moments.
Static Equilibium Static equilibium equies a balance of foces and a balance of moments. ΣF 0 ΣF 0 ΣF 0 ΣM 0 ΣM 0 ΣM 0 Eample 1: painte stands on a ladde that leans against the wall of a house at an angle
More informationBlack Body Radiation and Radiometric Parameters:
Black Body Radiation and Radiometic Paametes: All mateials absob and emit adiation to some extent. A blackbody is an idealization of how mateials emit and absob adiation. It can be used as a efeence fo
More informationGraphs of Sine and Cosine Functions
Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the
More informationLab 10: Newton s Second Law in Rotation
Lab 10: Newton s Second Law in Rotation We can descibe the motion of objects that otate (i.e. spin on an axis, like a popelle o a doo) using the same definitions, adapted fo otational motion, that we have
More information