The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

Similar documents
2C09 Design for seismic and climate changes

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Notes The Incremental Motion Model:

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Systems of Particles: Angular Momentum and Work Energy Principle

Balancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a)

Chimica Inorganica 3

ECE 308 Discrete-Time Signals and Systems

The Space Redundant Robotic Manipulator Chaotic Motion Dynamics Control Algorithm

Two or more points can be used to describe a rigid body. This will eliminate the need to define rotational coordinates for the body!

Analytical mechanics

Principle Of Superposition

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium

Course 4: Preparation for Calculus Unit 1: Families of Functions

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

G r a d e 1 1 P r e - C a l c u l u s M a t h e m a t i c s ( 3 0 S )

CHAPTER 8 SYSTEMS OF PARTICLES

Warped, Chirp Z-Transform: Radar Signal Processing

TESTING OF THE FORCES IN CABLE OF SUSPENSION STRUCTURE AND BRIDGES

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

Analysis of composites with multiple rigid-line reinforcements by the BEM

3 Balance equations ME338A CONTINUUM MECHANICS

Vibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

Signal Processing in Mechatronics

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

EXPERIMENT OF SIMPLE VIBRATION

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

Building an NMR Quantum Computer

A new formulation of internal forces for non-linear hypoelastic constitutive models verifying conservation laws

Today. Homework 4 due (usual box) Center of Mass Momentum

Diploma Programme. Mathematics HL guide. First examinations 2014


Physics Supplement to my class. Kinetic Theory

NEWTON'S SECOND LAW AND CURVILINEAR MOTION OF A MATERIAL PARTICLE. Yuriy Zhivotov Yuzhnoye State Design Office, Dniepropetrovsk, Ukraine

Appendix: The Laplace Transform

EDEXCEL NATIONAL CERTIFICATE UNIT 4 MATHEMATICS FOR TECHNICIANS OUTCOME 4 - CALCULUS

SOME RELATIONS ON HERMITE MATRIX POLYNOMIALS. Levent Kargin and Veli Kurt

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

ELEG 4603/5173L Digital Signal Processing Ch. 1 Discrete-Time Signals and Systems

Classical Mechanics Qualifying Exam Solutions Problem 1.

Chapter 5 Vibrational Motion

Chapter 7 z-transform

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

EXAM-3 MATH 261: Elementary Differential Equations MATH 261 FALL 2006 EXAMINATION COVER PAGE Professor Moseley

True Nature of Potential Energy of a Hydrogen Atom

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Machine Learning for Data Science (CS 4786)

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

Internal Information Representation and Processing

In number theory we will generally be working with integers, though occasionally fractions and irrationals will come into play.

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

SOLID MECHANICS TUTORIAL BALANCING OF RECIPROCATING MACHINERY

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains

Analytical solutions for multi-wave transfer matrices in layered structures

R is a scalar defined as follows:

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

Using Spreadsheets as a Computational Tool in Teaching Mechanical. Engineering

Additional Notes on Power Series

A NOTE ON THE TOTAL LEAST SQUARES FIT TO COPLANAR POINTS

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

DISTRIBUTION LAW Okunev I.V.

Estimation of Multibody Kinematics Using Position Measurements

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

Chapter 7: The z-transform. Chih-Wei Liu

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Teaching Mathematics Concepts via Computer Algebra Systems

, then cv V. Differential Equations Elements of Lineaer Algebra Name: Consider the differential equation. and y2 cos( kx)

ROBUST ATTITUDE DETERMINATION AND CONTROL SYSTEM OF MICROSATELLITE. Principle of the control system construction

6.003 Homework #3 Solutions

NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS

Streamfunction-Vorticity Formulation

ARISTOTELIAN PHYSICS

Inelastic spherical caps with defects

Math 113 Exam 3 Practice

a for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α

Limits for the Use of Auto-Balancing

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square

CALCULUS BASIC SUMMER REVIEW

Free Surface Hydrodynamics

Using An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically

DESIGN, PRODUCTION, AND APPLICATION OF A STAND FOR TESTING FRICTION OF THE BEARINGS

PILOT STUDY ON THE HORIZONTAL SHEAR BEHAVIOUR OF FRP RUBBER ISOLATORS

PROPOSING INPUT-DEPENDENT MODE CONTRIBUTION FACTORS FOR SIMPLIFIED SEISMIC RESPONSE ANALYSIS OF BUILDING SYSTEMS

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

2.004 Dynamics and Control II Spring 2008

Finite Element Modeling of Seismic Response of Field Fabricated Liquefied Natural Gas (LNG) Spherical Storage Vessels

Math 113 Exam 3 Practice

Mathematics Extension 1

Frequency Domain Filtering

Mechanics Physics 151

Introduction to Machine Learning DIS10

Transcription:

The Mathematical Model ad the Simulatio Modellig Algoritm of the Multitiered Mechaical System Demi Aatoliy, Kovalev Iva Dept. of Optical Digital Systems ad Techologies, The St. Petersburg Natioal Research Uiversity of Iformatio Techologies, Sait Petersburg, RUSSIAN FEDERATION (LOMO PLC) Moder mechaical systems are usually have complex spatial ad multitiered structure. The importat stage of the desig of such systems is to estimate its behavior i the dyamics which ca be obtaied by computer simulatio. Figure 1 below shows a multi-mechaical cotrol system of the supervisory system, icludig a mobile base (dec of a ship) ad the observatio device i gimbal suspesio. The feature of this system is that the mass ceters of compoets of the system are separated i space. Fig. 1 Multitiered mechaical cotrol system of the supervisory system. Asia Busiess Cosortium ABC-JAR Page 44

Figure 1 uses the followig otatios: the first three compoets K, T, R simulate the rotatio of the vessel relative to three orthogoal axes, so that, ω K the agular velocity of rotatio of the dec at the agle of heel, ω T the agular velocity of rotatio of the dec at the corer of yaw, ω R the agular velocity of rotatio of the dec at the corer of pitch attitude. The poles of these compoets are the same ad are i the ceter of the ship rollig. The secod three compoets G, V, C simulate the optical device rotatio so that, ω G the agular velocity of the optical device rotatio by horizotal agle, ω V the agular velocity of the rotatio of the optical device by vertical agle, ω C the agular velocity of the optical device rotatio by agle of compesatio. The simplest ad most effective way to set up the equatios of such systems with geometric costraits o the locatio of its poits (ot free holoomic system) is based o the represetatio of their relatively geeralized forces determied by the eergy of the system, ad geeralized coordiates. The equatios of motio i this case as the Lagrage equatios of the 2d id ca be represeted i vector-matrix form: Aq Q 0,5Dq Q 1 2 aiq q i 1 q where A the matrix of derivatives of the ietic eergy by coordiates; Q geeralized forces matrix; q geeralized coordiates matrix (i=1, 2,,); D partial derivatives matrix; the umber of degrees of freedom; a i coefficiets of iertia. The advatage of this represetatio is its geerality ad uiform of costructig the equatios. However, due to the fact that the presetatio of (1) is complex i terms of the spatial orietatio of the system ad aalysis of its costituet parts, appears difficulty i simulatio modelig. I order to simplify the procedure of buildig up the equatios, we tae the cocept of quasi-velocities, i.e. we ll carry out the liearizatio of the geeralized velocities. Such a simplificatio allows a mathematical model of motio of the mechaical part of the system i a form suitable for its partial aalysis. I this case the stage of the withdrawal of differetial equatios for the dyamics of the body ca be omitted. I this way we ca preset the mechaical multitiered cotrol system of the optical device (see fig.1) i form of the graphic-vector descriptio as show i Figure 2, where the followig otatios are used : О the ceter of the earth rotatio; V K the liear speed of the ship; K, T, R, G, V, C poles of the relevat parts of the system; R K,T,R,G,V,C radius vectors of the correspodig poles; R C radius vector of the mass ceter of the correspodig part to the coordiate system associated with this part; M K,T,R,G,V,C momets of the forces actig o the relevat body; ω O,K,T,R,G,V,C agular velocity vector of the relevat parts; D T,R,G,V,C the poles of the drive uits of the correspodig parts; R D radius vectors of drive uit poles to the coordiate system associated with the respective parts; ω D agular velocities of the respective drive uits; R CD radius vectors of the momets of the forces correspodig to the appropriate drive uits; M D momets of forces of correspodig drive uits; γ T,R,G,V,C agles of bodies rotatio uder the actio of the correspodig drive uits. (1) Asia Busiess Cosortium ABC-JAR Page 45

Fig.2 Graphic-vector descriptio of the multitiered mechaical system of optical device cotrol A mathematical model of motio of the mechaical multitiered system ca be made o the basis of: the law of coservatio of eergy; virtual wor i couctio with the priciple d'alembert; Newto's ad Euler's equatios; Lagrage equatios; Hamilto's priciple. Accordig to this ratio as the iitial set-ow equatios of the free body motio i a vector form: m 0 W 0 + ω 0 r 0 + ω 0 ω 0 r 0 = J 0 0 + 0 J 0 0 + m 0 r 0 W 0 = M (2 ) W 0 = V 0 + 0 V 0, where m 0 и J 0 mass ad iertia tesor give to the top of the coordiate system fixed to the body X 0Y 0Z 0; F 0 и M 0 the mai vectors of exteral forces ad momets of exteral forces applied to the body accordig to the coordiate system X 0Y 0Z 0; 0 и V 0 agular ad liear velocity of top of the coordiate system X 0Y 0Z 0; r 0 radius vector of the body's mass ceter; W 0 acceleratio vector of the coordiate system X 0Y 0Z 0. The first equatio expresses the theorem quatity of motio (equatio of motio of the mass ceter). The secod equatio expresses the theorem quatity of motio relative to the coordiate system X 0Y 0Z 0 (equatio of "rotatio"). Ay mechaical system ca be decomposed ito a fiite set of solid bodies S i(i=,,n; =0,1, ). The positio of each body S i will be determied relative to the body S i-1 of this set, ad the positio ad movemet of the body S 0 will be determied relative to the iertial coordiate system XYZ. We will also assume that the body S i ca move relatively to the body S i-1 oly uder the ifluece of the force vectors F i from the body Si-1. I case of motio of the system of solid bodies, the positio of the ceter of iertia C i does ot chage correspodig to the pole O i, associated with the body S i coordiate system X iy iz i, but varies i relatio to the poles of the coordiate system with the idex other tha the «i». The positio of the ceter of iertia С О of the body S 0 relative to the iertial coordiate system XYZ is defied by the Asia Busiess Cosortium ABC-JAR Page 46

radius vector R 00 (3а), the positio of the ceter of iertia C 1 of the body S 1 relative to the XYZ ad X 0Y 0Z 0 is determied by the radius vector R 01 (3б), ad i geeral case (3c): R 00 = R 0 + R C0, a) R 01 = R 00 + R 1 + R C1 = R 0 + R C0 + R 1 + R C1 b) (3) R оп = i=0 R сп = R оп R i + R ci 1 i=0 ; R i + R ci + R п. where R 0 radius vector of the top of the coordiate system X 0Y 0Z 0 relative to the top of the coordiate system XYZ; R C0 radius vector of the ceter of the mass of the body S 0 i the X 0Y 0Z 0 coordiate system. To determie the velocity vectors of the ceters of iertia of the body relative to the appropriate coordiate system, you must perform the operatio of differetiatio with respect to time of the correspodig proectios of the radius vectors at the axis of coordiate systems, as show i 4a ad 4b i geeral case: V 00 = R 00 = V 0 + 0 R co ; а) (4) V 01 = R 01 = V 0 + V 1 + 0 R 1 + 0 + 1 R c1 ; V оп = R оп. V i = R i = V + ω l ω =i R + ω R c, where i=0, 1, 2,, N; =, +1,, N; V = R velocity vector of pole O, associated with the body S of the coordiate system X Y Z, relative to the coordiate system X -1Y - 1Z -1, attached to the body S -1. Determiatio of the acceleratio vector of ceters of iertia of bodies Si respect to the correspodig coordiate systems is similar. We ca show this i geeral case we get: W i = V i = R i = V + 2 ω l ω =i + ω l ω R l R + R c ω l ω V + ω l ω + ω R R + ω l ω ω R c, where V acceleratio vector of the pole O, of the coordiate system X Y Z attached to the body S relative to the coordiate system X -1Y -1Z -1, attached to the body S -1 proected o the axis X Y Z. The process of buildig a mathematical model of multitiered mechaical system is a sequece of operatios preseted i the form of the followig algorithm (see Fig.3). c) b) ω (5) Asia Busiess Cosortium ABC-JAR Page 47

Fig. 3 Algorithm of buildig a mathematical model of the multitiered mechaical system CONCLUSION The techique of buildig а mathematical model of arbitrary multitiered mechaical system has bee examied i this article. The relatioship betwee the compoets of the system is determied usig a graphic-vector scheme. Locatio of items i space is determied by the equatios of the motio of a free body i vector form. The algorithm of buildig a mathematical model of the multitiered mechaical system has bee also formed. The method is demostrated o the example. CONCLUSION Ahmad, M., Alam, M., & Taluder, M. (2013). No-Associativity of Loretz Trasformatio ad Associative Reciprocal Symmetric Trasformatio. Asia Joural Of Applied Sciece Ad Egieerig, 2(1), 48-56. Ozada, N., & Saeed Madai, S. (2012). A Novel Musculoseletal Joit Modelig for Orthopaedic Applicatios. ABC Joural Of Advaced Research, 1(1), 35-46. Taluder, M., & Ahmad, M. (2013). Relativistic Rule of Multiplicatio of Velocities Cosistet with Loretz Eistei Law of Additio ad Derivatio of the Missig Equatios of Special Relativity. Asia Joural Of Applied Sciece Ad Egieerig, 2(1), 26-41. Asia Busiess Cosortium ABC-JAR Page 48