MATHEMATICAL MODELING OF THE TRACTOR-GRADER AGRICULTURAL SYSTEM CINEMATIC DURING LAND IMPROVING WORKS

Similar documents
1. VELOCITY AND ACCELERATION

Suggested Practice Problems (set #2) for the Physics Placement Test

SPH3U: Projectiles. Recorder: Manager: Speaker:

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

PHYSICS 149: Lecture 9

The motions of the celt on a horizontal plane with viscous friction

Chapter 3 (Lectures 12, 13 and 14) Longitudinal stick free static stability and control

Structural Dynamics and Earthquake Engineering

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams

Q.1 Define work and its unit?

Electrical and current self-induction

Curling Stress Equation for Transverse Joint Edge of a Concrete Pavement Slab Based on Finite-Element Method Analysis

Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture

Lab #2: Kinematics in 1-Dimension

CLASS XI SET A PHYSICS. 1. If and Let. The correct order of % error in. (a) (b) x = y > z (c) x < z < y (d) x > z < y

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

Some Basic Information about M-S-D Systems

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

At the end of this lesson, the students should be able to understand

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

KINEMATICS IN ONE DIMENSION

L1, L2, N1 N2. + Vout. C out. Figure 2.1.1: Flyback converter

Comparative study between two models of a linear oscillating tubular motor

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Sliding Mode Controller for Unstable Systems

Non-uniform circular motion *

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

2002 November 14 Exam III Physics 191

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Module 4: Time Response of discrete time systems Lecture Note 2

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification

Transverse and Longitudinal Damped Vibrations of Hydraulic Cylinder in a Mining Prop

STATE-SPACE MODELLING. A mass balance across the tank gives:

MECHANICS OF MATERIALS Poisson s Ratio

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Advanced Organic Chemistry

Exponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives

Reasonable compensation coefficient of maximum gradient in long railway tunnels

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration

Elements of Computer Graphics

Chapter 10 INDUCTANCE Recommended Problems:

v A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

Stability and Bifurcation in a Neural Network Model with Two Delays

Phys 221 Fall Chapter 2. Motion in One Dimension. 2014, 2005 A. Dzyubenko Brooks/Cole

- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

SUPPLEMENTARY INFORMATION

Final Spring 2007

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

IB Physics Kinematics Worksheet

LabQuest 24. Capacitors

and v y . The changes occur, respectively, because of the acceleration components a x and a y

Sub Module 2.6. Measurement of transient temperature

Section 3.8, Mechanical and Electrical Vibrations

Best test practice: Take the past test on the class website

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

A car following model for traffic flow simulation

ψ(t) = V x (0)V x (t)

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Physics 180A Fall 2008 Test points. Provide the best answer to the following questions and problems. Watch your sig figs.

Two Coupled Oscillators / Normal Modes

0 time. 2 Which graph represents the motion of a car that is travelling along a straight road with a uniformly increasing speed?

A FAMILY OF THREE-LEVEL DC-DC CONVERTERS

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

Unit 1 Test Review Physics Basics, Movement, and Vectors Chapters 1-3

Estimation of Kinetic Friction Coefficient for Sliding Rigid Block Nonstructural Components

NEWTON S SECOND LAW OF MOTION

Flow-Induced Vibration Analysis of Supported Pipes with a Crack

Solutions from Chapter 9.1 and 9.2

Today: Graphing. Note: I hope this joke will be funnier (or at least make you roll your eyes and say ugh ) after class. v (miles per hour ) Time

Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member

Constant Acceleration

Lab 10: RC, RL, and RLC Circuits

A Robust Exponentially Weighted Moving Average Control Chart for the Process Mean

Chapter 3 Kinematics in Two Dimensions

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES

Displacement ( x) x x x

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Development of a new metrological model for measuring of the water surface evaporation Tovmach L. Tovmach Yr. Abstract Introduction

This exam is formed of 4 obligatory exercises in four pages numbered from 1 to 4 The use of non-programmable calculators is allowed

th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3256

Shells with membrane behavior

Qualitative behavior of mixing phenomena - the case of axisymmetric extensional flows

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

Integration of the equation of motion with respect to time rather than displacement leads to the equations of impulse and momentum.

We may write the basic equation of motion for the particle, as

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections

Traveling Waves. Chapter Introduction

Chapter Q1. We need to understand Classical wave first. 3/28/2004 H133 Spring

Physics 5A Review 1. Eric Reichwein Department of Physics University of California, Santa Cruz. October 31, 2012

Transcription:

Bullein of he Transilvania Universiy of Braşov Series II: Foresry Wood Indusry Agriculural Food Engineering Vol. 5 (54) No. 1-2012 MATHEMATICA MODEING OF THE TRACTOR-GRADER AGRICUTURA SYSTEM CINEMATIC DURING AND IMPROVING WORKS Gheorghe BRĂTUCU 1 Absrac: This paper sudies he influence of he pulled leveling machine on he general dynamics of he wheel racor. The obained daa are employed in he design of he ransmission and of oher racor sub-assemblies, represening he disurbance facors of he mahemaical and dynamic models se up of esablishing he srain on various componens. I is shown ha he rear axle of he racor is loaded up o 53%, while he fron one is released o values below he admissible minimum (20% of he saic load). The loading coefficien of he rear wheels can reach A 2 = 0.76, while he longiudinal sabiliy angle ils during climbing and worsens by abou 57% in relaion o he solo racor. Key words: agriculural racors on wheels, general dynamic, graders. 1. Symbols is Z 01, Z 02 - saic reacions on he fron and rear racor axels; Z 1, Z 2 - dynamic reacions on he fron and rear racor axels; G, G m - racor weigh, grader weigh respecively; F a, F j, F m, F f, F - forces of: air resisance; racor ineria; driving; rolling resisance; resisance on he working movemen of he echnical sysem; r m - dynamic radius of he driving wheel; M r1, M r2 - rolling resisance momens on he fron and rear racor wheels;, a, b, c, h c, h, h a - geomerical dimensions, according o Figure 1; f - rolling resisance coefficien; α, β - longiudinal and ransversal angles of he slopes; γ - angle made by he force F wih he horizonal; k u - soil leveling resisance coefficien; A - working surface (in soil) of he grader blade. 2. Inroducion Graders are working machines which cu he high irregulariies of he errain using blades or boomless buckes which ranspor rough creeping he released maerial o lower places [2]. From he echnical sysem formaion poin of view of he graders are owed machines, wih rigid or aricled frame (Figure 1). In Figure 1 we presened he echnical sysem scheme, made from U-650 M racor and he owed grader NT-2.25, a ramp ascension (a) and moving on ransversal slopes (b) [3]. 1 Dep. of Engineering and Managemen in Food and Tourism, Transilvania Universiy of Braşov.

130 Bullein of he Transilvania Universiy of Braşov Series II Vol. 5 (54) No. 1-2012 Fig. 1. Movemen on longiudinal ramps (a) and ransversal slopes (b) of he racor-grader agriculural sysem The racion resisance opposed by he grader is calculaed wih he relaion: F = k A + G f. (1) u m Under he dynamic aspec, he grader influences he echnical sysem rough F forces magniude and he γ angle, on which forces are made wih an axel parallel wih he soil [1]. For his echnical sysem here will be considered in all he cases, he following consan values: racor weigh (G = 36200 N), racor wheel base ( = 2430 mm), graviy cener posiion (b = 850 mm, b = a = 1580 mm, h c = 900 mm) and racor gauge (B = 1920 mm). Variable values will be considered he ones which are relaed o he rolling resisance condiions and laeral slip, grader weigh, grader graviy cener posiion from he rear racor axel, racion resisance force and is posiion in relaion o he horizonal surface. The simulaion of he echnical sysem dynamics is made in MATHCAD language, using mahemaical models, corresponding o differen funcional siuaions and characerisic elemens of he general dynamics. 3. Mahemaical Modeling and Analysis of Tracor Axels Reacions As a reference elemen in analyzing he axels reacions of he racors which form he echnical sysems of land melioraion machines, are he saic reacions resuled from he racor weigh repariion on he wo axels, on horizonal errain: G b Z01 = = 12700 [ N], (2) G ( b) Z02 = = 23600 [ N]. (3) The dynamic reacions on he fron and rear wheels Z 1 and Z 2 are calculaed wih relaions: G b cos α ( G sin α + Fj ) hc Fa h a + F + F c g γ + M r1 + M r2 Z1 =, (4)

Brăucu, Gh.: Mahemaical Modeling of he Tracor-Grader Agriculural Sysem 131 Z 2 G ( b) cosα + ( G = M r1 + M r2 +. sin α + F j ) h c Fa ha + F + F ( + c) g γ + (5) The variaions of he dynamic reacions on hose wo racor axels are calculaed wih he relaions: Z. (6) 1 = Z1 Z01; Z 2 = Z 2 Z02 A he sudy of dynamic reacions i has been considered he case of a racor moving on a horizonal errain wih consan speed, wih he racion force parallel o he soil and neglecing air resisance. Resuls are presened in Table 1, from which we can find he influences of he main facors (F, γ and f) on he variaion of racor axels loads. Table 1 Reacion forces on he racor wheels a consan speed on a horizonal surface, wih he racion force parallel wih he soil surface Variables Z 01, Z 02, Z 1, Z 1 Z 2, Z 2 Symbols Values [%] [%] [N] [N] [%] [N] [N] [%] 5000 9860 2800 23 26340 +2800 12 F, [N] 15000 6480 6220 49 29470 +6370 27 35000 590 1260 96 36140 +12600 53 γ, [ ] 0 8920 3740 29 27280 +3710 15.9 12700 23600 14 6230 6470 51 29970 +6370 27 100 100 30 3320 9930 78 32870 +9930 42 0.06 3200 9460 74 33000 +9460 40 f 0.1 2800 9900 78 33510 +9910 42 0.2 1560 1110 87 34640 +11000 47 Observaions γ = 30 f = 0.l F = 25000 f = 0.1 F = 25000 γ = 30 In he condiions of mainaining all he echnical sysem parameers consan, he resisan racion force from he grader loads on he rear axel of he racor wih 12% for F = 5000 N, wih 27% for F = 15000 N and wih 53% if F would have he value F = 35000 N. If racion resisance remains consan, bu he angle γ is modified, an angle which his force makes wih he horizonal, he load on he rear axel increases by 15.9% for γ = 0 and F = 25000 N, from F = 0, wih 27% for γ = 14 and by 42% for γ = 30. For he rolling resisance coefficien f = 0.06, F = 25000 N and γ = 30, he variaion of he loads on he racor axels owards he unloaded racor siuaion, shows a rear axle load increasing by 40% and a download of he fron axel by 74%. Increasing he value of f o f = 0.1 makes he rear axle load increase by 42% and he fron axel o discharge by 78%. In he hypoheical case where F = 35000 N, he rear axle load is of 47%, and he fron axel download of 87%. 4. Mahemaical Modeling and Analyses of oading Coefficiens for Tracor Wheels Tracor axels loading a he leveling work need o consider also he ire lif, expressed hrough wheel loading coefficiens λ. For he saic posiion of he racor hese coefficiens have he values given by he relaion: b λ fs = λ01 = = 0.35, (7)

132 Bullein of he Transilvania Universiy of Braşov Series II Vol. 5 (54) No. 1-2012 b λ ss = λ02 = = 0.65. (8) For he siuaion when he racor is moving wih racion force on a horizonal errain, he coefficiens λ f and λ s have he values: b F + f rm G λ f =, (9) G b F + f rm G λ s = +. (10) G The changing of hese coefficiens in he condiions of echnical sysem movemen on a plain errain, wih he variable force F of he grader, bu arranged in an angle of 30 from he horizonal and a rolling resisance coefficien f = 0.01, is represened in Figure 2. We found a decrease of he dynamic loading coefficien of he seering wheels λ f = λ 1 (curve 3) compared wih he coefficien λ fs = λ 01 (curve 1) of he same wheels, from he wheels λ f = λ 1 = 0.305 from F = 5000 N, o he value λ f = λ 1 = 0.231 for F = 35000 N. Fig. 2. The variaion of racor wheels loading coefficiens during leveling work A he same ime, an increase of he driving wheel loading coefficien (curve 4) from he value λ s = λ 2 = 0.695 for F = 5000 N, o he value λ s = λ 2 = 0.769 for F = 35000 N (curve 2 represens λ ss ) is noiced. 5. Mahemaical Modeling and Analyses of he ongiudinal Sabiliy of he Tracor For he deerminaion of he criical angles of saic sabiliy a boarding ramp (α 02max ) and ramp descending (α 01max ), he racor is considered slowed in boh siuaions and he following relaions are applied: b α02 max = arcg = h c b α01 max = arcg = h c o 43.363[ ], (11) o 60.33 [ ]. (12) The criical angle of racor dynamic sabiliy a boarding ramp is calculaed wih he relaion: b ( ϕ f ) ( h + c g γ) gαcr =, (13) h ( h + c g γ) c if he angle γ has a significan value wih he relaion: b ( ϕ f ) h g αcr =, (14) hc h if angle γ = 0. The case of his agriculural sysem is no considered relevan in deermining he criical angles of longiudinal dynamic sabiliy a ramp descending, as a resul of his significanly improving wih he owed echnical sysem. In Table 2 he angle variaion of he owed echnical sysem longiudinal sabiliy is presened according o he resisance force F of he grader movemen. In he case of he racor-grader agriculural sysem s movemen wih a racion resisance F parallel wih he soil surface (γ = 0), he influence of he F force on back overhrow ( α 2 ) is in he limi of α 2 = 18.93%, if F = 20000 N and of α 2 = 9.2%, if F = 5000 N. However, if γ = 30, α 2 = 57.27%,

Brăucu, Gh.: Mahemaical Modeling of he Tracor-Grader Agriculural Sysem 133 The angle variaion of he owed echnical sysem longiudinal sabiliy according o he resisance forces F of he grader movemen Table 2 Variables α 01, α 02, α 1, α 1 α 2, α 2 Symbols Values [%] [%] [ ] [ ] [%] [ ] [ ] [%] 5000 63.7 +3.4 +5.6 37.77 5.53 12.7 F, [N] 10000 65.6 +5.3 +8.7 34.70 8.6 19.8 15000 67.3 +7.0 +11.6 31.40 11.91 27.5 20000 68.75 +8.45 +14.0 27.83 15.47 35.7 5000 65.46 +5.16 +8.5 36.0 7.3 16.8 60.3 43.3 10000 68.47 +8.17 +13.5 30.78 12.52 28.9 F, [N] 100 100 15000 70.85 +10.55 +17.5 24.94 18.36 42.4 20000 72.78 +12.48 +20.9 18.5 24.8 52.3 5000 62.0 + 1.7 +2.8 39.3 4.0 9.2 F, [N] 10000 62.5 +2.2 +3.6 38.0 5.3 12.3 15000 63.0 +2.7 +4.4 36.6 6.7 15.5 20000 63.5 +3.2 +5.3 35.1 8.2 18.9 Observaions G = 36200 N h = 0.25 f = 0.1 γ = 15 G = 36200 N h = 0,25 f = 0.l γ = 30 G = 36200 N h = 0.25 f = 0.1 γ = 0 for F = 20000 N and α 2 = 16.8%, for F = 5000 N. More suggesive is he represenaion from Figure 3 of he γ angle influence over he longiudinal racor sliding sabiliy angles, in which 1 represens he posiion of longiudinal saic sabiliy criical angle a slopes descending (α 01max ); 2 - he posiion of he saic sabiliy criical angle a ramps ascending (α 02max ); 4 - he posiion of dynamic sabiliy angle a ramps ascending, if γ = 0, and 3 - he posiion of longiudinal dynamic sabiliy criical angle o sliding if γ = 0...30. I is noed ha if he angle γ increases, he rear axel is geing loaded and he dynamic sliding sabiliy is improved. 6. The Mahemaical Modeling and Analyze of Tracor Transversal Sabiliy The saic sabiliy a agriculural sysem overhrow on ransversal slopes is calculaed wih he relaion: 0.5 B βmax = arcg. (15) h c For he maximum value of he racor U-650M rack widh (B = l.92 m) he saic sabiliy overhrow angle has he variaion represened in Figure 4, dependen only upon he racor graviy cener posiion h c. I is noiced ha for h c = 0.6 m, β max = 58, and for h c = l m, β max = 43, reaching a β max = 35, for h c = 1.4 m. The racor dynamic sabiliy a laeral overhrow shall be ensured if: G 2 0.5B cosβ hc sin β v R G, (16) g h c Fig. 3. The influence of angle γ over he saic and dynamic sabiliy angles a movemen on ramps and slopes in which he firs par (he cenrifugal force) is dependen on he speed v (for G and R consan on a racor), and he second par is dependen upon angle β of slope inclinaion.

134 Bullein of he Transilvania Universiy of Braşov Series II Vol. 5 (54) No. 1-2012 7. Conclusions Fig. 4. The variaion of racor overhrow saic sabiliy angle depending on he posiion of he racor graviy cener h c Fig. 5. The variaion of dynamic sabiliy a laeral overhrow The condiion required by he relaion (13) is graphically represened in Figure 5, where on he abscissa, he movemen speed v and he slope angle β are represened, and on he ordinae, he cenrifugal force and he oher par of he inequaion are rendered. I is noiced ha he agriculural sysem s sabiliy is assured only if he movemen speed v < 2 m/s (curve 1), and he ransversal slope angle β < 30 (curve 2), a differen siuaion from he one sudied in he case of sabiliy a saic overhrow. 1. The moving resisance force F of he operaing racor-grader agriculural sysem modifies he loads on he racor axels by almos +27% on rear axel and by 49% an fron axel (a F = 15000N), he angle γ loads he rear axel by almos 26% and unloads he fron axel by 60%, and he rolling resisance coefficien causes he rear axel loading by a maximum of 7% and a fron axel unloading of 13%. 2. The fron axel unloading can no be acceped wih more han 80%, because he agriculural sysem s movemen direcions can no be safely conrolled. 3. The dynamic loading coefficiens of he racor wheels have similar variaions wih he axels loads, in he sense ha a he fron wheels (seering wheels) hose coefficiens decrease from 0.35 o 0.231 ( 34%) and heir increase on he rear wheels from 0.65 o 0.769 (+18.3%) because of force F and in smaller proporions if only he influences of γ and f are aken ino consideraion. References 1. Ionescu, E., e al.: Tracoare şi auomobile (Tracors and Auomobiles). Braşov. Transilvania Universiy Publishing House, 1979. 2. Mihăilescu, S., e al.: Maşini de consrucţii (Building Machines). Vol. II. Bucureşi. Technical Publishing House, 1985. 3. Năsăsoiu, S., Popescu, S., e al.: Tracoare (Tracors). Bucureşi. Didacic and Pedagogic Publishing House, 1983.