Computer-aided generation of equations and structural diagrams for simulation of linear stationary mechanical dynamic systems/Lygciu ir strukturiniu schemu kompiuterizuotas sudarymas tiesinems stacionarioms mechaninems dinaminems sistemos modeliuoti.

Augustaitis, V.K. ; Gican, V. ; Sesok, N. 等

1. Introduction

Computer-aided methods of analysis and simulation are widely used for the investigation on dynamic properties of mechanic, mechatronic and other equipment [1-12].

In majority of cases, the application of computer-aided simulation in the environment of MATLAB / Simulink programs for such investigation is purposeful [13, 14]. For this purpose, the equations for description of the dynamic processes of the equipment under investigation, hereinafter referred to as the dynamic system or shortly system, and the structural diagram of the system formed on the base of the said equations according to the dynamic model of the system with discrete elements (widely used for examination of automatic control systems) should be available. Such a diagram formed according to the requirements of MATLAB/Simulink program package and "understandable" for a computer is referred to as Simulink-model [13, 14]. There are no substantial differences between it and the structural diagram: if any of them is available, the other is easily found on its base.

When a mechatronic system is examined, its electrical part, such as electric drive, components of sensors and so on, usually are provided as already known structural diagrams [15]. Cases of formation of structural diagrams of the mechanical part of mechanical or mechatronic systems appear to be more complicated. Because of the wide variety of structures of mechanical systems, their structural diagrams are not predictable in any specific case, so their formation in cases of complicated systems requires considerable attempts and errors are hardly avoidable if the process is not computerized.

The purpose of the paper is provision of a methodology of the application of MATLAB/Simulink program package to the available linearized mechanical system (or is a part of more complicated system) for computer-aided generation of equations describing movement of the system using Lagrange equations of the second type, transformation of the generated equations into the convenient structural diagram of the said system and Simulink-model, formation of literal and digital analytic expressions of transfer functions included in them.

Software package MATLAB Simulink program is not directly intended for creation of the said model. However, there are enough resources to solve this problem, which is devoted to this work.

The obtained structural diagram can be easily integrated in a structural diagram and Simulink-model of a more complicated system. In addition, the generated equations of the system are of independent value as well.

2. Generation of the equations

The equations describing movement of the linear mechanical dynamic system with lumped parameters required for the formation of structural diagram of the system are obtained from Lagrange equations of the second type

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where T, [PI] are kinetic and potential energy of the system, respectively; [PHI] dissipation function; [x.sub.i], [[??].sub.i] i-th generalized Lagrangian coordinate, shortly referred to as the coordinate, and its time (t) derivative; [F.sub.i] (t) the generalized outside force acting along the coordinate [x.sub.i].

It can be stated that each of coordinates [x.sub.i] corresponds to one equation of (1). In traditional case, the equation (1) include n generalized coordinates [x.sub.i]; their number n equals to the number m of degrees of freedom of the system (n = m) when a structural diagram for of self-contained mechanical system is generated of transformed Eq. (1). However, on generation of structural diagrams in more common case when the stationary linear mechanical system under discussion is a part (component) of a larger mechanical system with nonlinear or nonstationary components, a mechatronic system and so on, it can be m > n, i.e. we'll have n equations with m coordinates included in them and the structural diagram for such a system. Hereinafter, we'll mark such coordinates by [x.sub.s] (s = 1, ..., i, ..., m). In such a case, the Eq. (2) will include k = m - n redundant coordinates. For the simulation of a system described by equations with redundant coordinates, the values of such redundant coordinates are set or found from the "rejected" part of the system by connecting its structural diagram to the structural diagram of the mechanical part under discussion.

It is notable that in both cases, i.e. when m = n and when m > n, the methodology of formation of structural diagrams and Simulink-models remains the same. Incorporating of auxiliary coordinates in the model let us assume analyzed system as a part of another more complex system, which may have nonlinear and nonstationary elements, electrical elements and so on.

In the case under discussion, kinetic energy T is a function of derivatives of coordinates [[??].sub.s] and (more rarely) a function of the coordinates [x.sub.s] themselves; potential energy [PI] is a function of coordinates [x.sub.s], and dissipation function [PHI]--a function of derivatives [[??].sub.s]. In addition, the expressions of T, n, O can be direct functions of time t, when, for example, known kinematic excitations [u.sub.iz] (t) (z = 1, 2, ..., d) are included in the said expressions

T = T(x, [??], t); [PI] = [PI] (x,t); [PHI] = [PHI] ([??], t) (2)

where x = [x.sub.1], [x.sub.2], ..., [x.sub.m]; [??] = [[??].sub.1], [[??].sub.2], ..., [[??].sub.m] are the totalities of coordinates [x.sub.s] and their derivatives [[??].sub.s].

In many cases [16, 17], striving to facilitate the generation of analytic expressions of the functions T, [PI], [PHI] for complicated systems, it is purposeful to use not only the generalized coordinates [x.sub.s] included in the Eq. (1), hereinafter referred to as the principal generalized coordinates, or shortly--to as the principal coordinates, but also auxiliary coordinates [[delta].sub.j]. They should be chosen in accordance with the below condition, i.e. they should be expressed unambiguously by the principal coordinates [x.sub.s] according to linear dependence in the following equations of connection

[[delta].sub.j] = [[alpha].sub.j,1] [x.sub.1] + [[alpha].sub.j,2] [x.sub.2] + ... + [[alpha].sub.j,m] [x.sub.m] + [[gamma].sub.j] (t) (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

where [[alpha].sub.j,1], ..., [[alpha].sub.j,m] are constant coefficients formed of the parameters of the system's dynamic model; [gamma](t) nonstationary members being predictable functions of time t (a part of the said coefficients or all nonstationary members or any part of them can be equal to zero).

When the auxiliary coordinates [delta] = [[delta].sub.1], [[delta].sub.2], ..., [[delta].sub.l] are applied, we obtain the following expression instead of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

Before differentiation of the functions T, [PI], [PHI] included in the Eq. (1), the auxiliary coordinates and their derivatives should be eliminated by using the Eqs. (3) and (4), i.e. T, [PI], [PHI] should be provided the expressions analogous to Eq. (2).

If the auxiliary coordinates are properly chosen according to the Eq. (3), the analytic Eq. (5) of the functions T, [PI], [PHI] developed by the investigator will be considerably simpler, as compared to Eq. (2). Then elimination of the auxiliary coordinates and their derivatives as well as differentiation of the functions Eq. (2) shall be carried out in a computer-aided way.

When the final Eq. (2) of the functions T, [PI], [PHI] F after differentiation and other relevant procedures are available, linear differential equations with constant coefficients of a degree not higher than second (in some cases, linear algebraic equations) are obtained for description of the dynamic model of the system under discussion. As it was mentioned, each i-th generalized coordinate [x.sub.i] corresponds to i-th equation from the Eq. (1). After introducing the differentiation operator p = d/dt, the equation will be transformed as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

where

[d.sub.i,j](p) = [a.sub.i,j] [p.sup.2] + [b.sub.i,j] p + [c.sub.i,j] (i = 1, 2, ..., n ; j = 1, 2, ..., m) (7)

[d.sub.i,i] (p) = [a.sub.i,i] [p.sup.2] + [b.sub.i,i] p + [a.sub.i,i] 8)

polynomials of the second degree in respect of p; a ,b, c with relevant indexes--constant coefficients (any of them can be equal to zero); [F.sub.i](t)--the generalized outside force acting along the [x.sub.i]--th generalized coordinate; [H.sub.i] (t) is the component of the generalized force acting along the [x.sub.i]--th generalized coordinate obtained on differentiation of the functions T, [PI], [PHI] in the cases of their direct dependence on time t (for example, when the equations (3) and (4) include nonstationary members [[gamma].sub.j] and [[??].sub.j]); [x.sub.1], ..., [x.sub.n] are principal generalized coordinates (each of them corresponds to one of the equations (1)), and; [x.sub.n+1], ..., [x.sub.m] are redundant generalized coordinates when they are used. It is accepted that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

[g.sub.i,e] (p) = [[alpha].sub.i,e] [p.sup.2] + [[eta].sub.i,e] p + [[xi].sub.i,e] (10)

where [[alpha].sub.i,e], [[eta].sub.i,e], [[xi].sub.i,e] are constant coefficients; [N.sub.ie] (t) is the known functions of time t, i.e. the nonstationary members [[gamma].sub.j] (t) included in the Eq. (3) or kinematic excitations [u.sub.i,z] (t) included in the T, [PI], [PHI] Eqs. (2) and (5) but not assessed in the Eq. (3).

The expressions of the generalized forces [F.sub.i](t) are developed by the investigator, so they are not discussed in details herein.

For the formation of structural diagram, each Eq. (6) is solved in respect of the coordinate [x.sub.i] included in it

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

According to the equations (11), the structural diagrams corresponding to them are formed; they are connected into the complete typical structural diagram of the whole system shown in Fig. 1 and on the base of the latter, Simulink-model for simulation of the system is developed.

[FIGURE 1 OMITTED]

3. The example

Let's suppose that low frequency vibrations of a vehicle moving on an uneven road are examined (Fig. 2). The vehicle is considered as absolutely solid body standing on elastic supports with damping; the supports simulate units of its wheels. The masses of the wheels and the suspension are neglected. For generation of equations for describing movement of the vehicle, two systems of coordinates are chosen. One of them, i.e. the system [O.sub.1], [X.sub.1], [Y.sub.1], [Z.sub.1] is fixedly connected to the vehicle and moves together with it. The point of the origin of coordinates of the said system of coordinates coincides with the center of stiffness of the vehicle located in the plane of fixing four elastic supports to the vehicle. In the plane, the points of fixation of the supports 1, 2, 3, and 4 are situated. All these points are situated in the same distance from the point of origin of coordinates [O.sub.1] and are symmetrical in respect of it (their positions are defined by the distances [l.sub.1] and [l.sub.2], respectively). The other (immovable) system of coordinates O, X, Y, Z, upon the rest of the vehicle when it is affected by the force of gravity only coincides with the system [O.sub.1], [X.sub.1], [Y.sub.1], [Z.sub.1], (in such a state, both systems of coordinates are shown in Fig. 2). It is conditionally supposed that the velocity of the vehicle's movement equals to zero and the uneven road is "moving" with its velocity u. It is considered that vibrations of the vehicle in respect of the system of coordinates O, X, Y, Z are small. During the vibrations, it rotates by small angles [[phi].sub.x], [[phi].sub.y] about the axes passing the point O1 and parallel to the axes X and Y ; in addition, it moves in vertical direction to the distance z along the axis of coordinates Z. For simplifications of the example, the shifts along the axes X, Y and rotation about the axis Z are neglected. It is considered that coordinates of the center S of the mass of the vehicle in the system of coordinates [O.sub.1], [X.sub.1], [Y.sub.1], [Z.sub.1] are [x.sub.s], [y.sub.s], [z.sub.s]. On simulation of the wheels and units of their suspensions, all four elastic supports of the vehicle are considered alike (their coefficients of stiffness and resistance are k and h, respectively). It was supposed that vibrations of the vehicle are excited in a kinematic way by the road's inequalities [z.sub.1] (t), [z.sub.2] (t), [z.sub.3] (t) and [z.sub.4] (t) that are known functions of time t (kinematic excitation). In addition, excitation of vibrations of the vehicle by the forces and moments of forces that impact the engine is assessed as well; after reduction to the origin of coordinates of the system [O.sub.1], [X.sub.1], [Y.sub.1], [Z.sub.1], the said forces and moments of forces are expressed by the vertical force [P.sub.z] (t) and the moments [M.sub.x] (t) and [M.sub.y] (t) of rotational forces. The mass of the vehicle is m, its moments of inertia about the axes [X.sub.1], [Y.sub.1] are [J.sub.x], [J.sub.y], respectively, and its combined moment of inertia is [J.sub.x,y] (no other moments of inertia exist in the case under discussion).

Thus, the dynamic model discussed upon in the example has three degrees of freedom and its movement is defined by the principal generalized coordinates [[phi].sub.x], [[phi].sub.y], z. The following specific values of the parameters of the dynamic model presented in the example were accepted: m = 5000 kg; [J.sub.x] = 3600kg[m.sup.2]; [J.sub.y] - 6000kg [m.sup.2]; [J.sub.xy] = 600kg [m.sup.2]; [l.sub.1] = 1.35m; [l.sub.2] = 0.85m; k = 80000 N/m; h = 1920Ns/m; [x.sub.s] = 0.2m; [y.sub.s] = 0.05m.

Kinetic energy of the machine according to [18] without auxiliary coordinates

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

For the simplification of the initial expression of potential energy [PI] and dissipation function [PHI], auxiliary coordinates are introduced. Let's suppose that on vibrations of the vehicle, the elastic elements used for simulation of the wheels and their suspensions are deformed by the values [[delta].sub.1], [[delta].sub.2], [[delta].sub.3], [[delta].sub.4] that are considered auxiliary generalized coordinates. So, potential energy

[PI] = k/2 ([[delta].sup.2.sub.1] + [[delta].sup.2.sub.2] + [[delta].sup.2.sub.3] + [[delta].sup.2.sub.4] (13)

The auxiliary coordinates are expressed by the principal ones using the Eq. (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

After insertion of the values of the auxiliary coordinates [[delta].sub.j] in the Eq. (13), potential energy expressed by principal coordinates only is found

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

The expression of dissipative function is found in an analogous way

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

[FIGURE 3 OMITTED]

The generalized outside forces

[F.sub.1] = [P.sub.z](t); [F.sub.2] = [M.sub.x](t); [F.sub.3] = [M.sub.y] (t) (17)

Then, using Lagrange Eq. (1), the equations of the vehicle's movement in principal generalized coordinates [[phi].sub.x], [[phi].sub.y],z (there are no auxiliary coordinates in the example under discussion) are generated. The expressions of T, [PI], [PHI] Eqs. (12), (15) and (16) included in the Eq. (1) are differentiated and if the generalized outside forces Eq. (17) are known, the differential equations of the vehicle's vibrations are found; from the latter, Eq. (11) type are obtained

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19)

The structural diagram shown in the Fig. 3 is formed on the base of the solutions of (18) and the typical structural diagram (Fig. 1).

4. Computerized realization of the proposed methodology

The application of the computerized methodology of computation of Eq. (11) type equations describing mechanical dynamic systems and the parameters of Simulink--model formed on their base in literal and digital form can be divided to the following phases.

1. Familiarization with the system, formation of its dynamic model, computation of the values of its parameters, choosing the principal generalized coordinates, choosing the auxiliary generalized coordinates (if they are used) generation of their Eqs. (3) and (4) in the principal generalized coordinates, generation of expressions of kinetic energy T, potential energy [PI] and dissipative function [PHI] included in the Eq. 1), formation of expressions of the generalized outside forces [F.sub.i], kinematic excitations [N.sub.i,e] and [u.sub.i,z] that affect the system. All materials are developed by the investigator.

2. Entering the initial data in to a computer. The data to be entered: data names; control matrix K; names of auxiliary coordinates Eqs. (3) and, (4), if they are used, and their time t derivatives; analytic expressions of the functions T, [PI], [PHI], numerical values of the parameters of the dynamic model.

First of all, data names are entered. They are written in lines; each line starts from the word "syms". The order of entering is chosen freely. The entered data names are distributed in a line, for example, as follows.

In the first line, all principal generalized coordinates [x.sub.1], [x.sub.2], ..., [x.sub.m], the kinematic excitations [u.sub.iz] (if they are used) in the expressions of T, [PI], [PHI] and nonstationary members [[gamma].sub.j] included in the expressions of the auxiliary coordinates Eq. (3) are listed.

In the second line, the derivatives of all values listed in the first line are provided; they are marked with the letter D before any of such value, for example, Dz, Dx and so on.

In the third line, the names of all parameters included in the dynamic model of the system under discussion (all constants included in the expressions of the functions T, [PI], [PHI], the connection Eq. (3), kinematic and outside excitations and generalized forces) are provided.

In the fourth line, the names of generalized outside forces [F.sub.i] (t), the differentiation operator p, and time t included in the equation (1) are provided.

For the example under discussion, the data names are written as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20)

The control matrix K is entered in four lines. Each line should include the same number of elements.

In the first line, the names of the principal generalized coordinates [x.sub.1], [x.sub.2], ..., [x.sub.m], kinematic excitations and nonstationary members are entered in any order.

In the second line, time t derivatives of all the values listed in the first line are provided; they are marked with the letter D before any of such value. Under the name of each value listed in the first line, the name of its derivative should be specified in the second line, i.e. the name of a value mentioned in the first line and the name of its derivative should be in the same column of the matrix K.

In the third line, names of the generalized outside forces [F.sub.i](t) acting along the generalized coordinates [x.sub.i] shall be written. The said names shall be written in those columns with coordinate's [x.sub.i] of the matrix K that correspond to the directions of acting of the forces. Other elements of the third line are equal to zero.

In the fourth line, it is specified that the functions T, [[PI], [PHI] included in Lagrange equations (1) should be differentiated according to the coordinates [x.sub.1], [x.sub.2], [x.sub.n] and their derivatives. In the elements of the line situated in the columns of the names of the mentioned coordinates and their derivatives, any positive whole numbers, such as 1, is written. Other elements of the line are equal to zero. For the example under discussion:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

After the matrix K, the analytic (literal) expressions of auxiliary coordinates [[delta].sub.j] and their derivatives [[delta].sub.j] are entered.

They are entered in lines according to the order specified in the Eqs. (3) and (4).

For the example under discussion

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]; (22)

Then the analytic expressions of the functions T, [PI], [PHI] are entered (on entering, they are marked as TM, PM, FM, respectively). They consist of the sum of the summonds of the line [1]. The summonds of each of said functions with the auxiliary coordinates (when they are used) are entered as summarized elements of the column of the vector.

For the example under discussion

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

Then the original program "Functions" is called.

Functions

Then the numeral values of the parameters of the dynamic model are entered. For the example under discussion

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

The Program Functions carries out the following actions: develops the analytic (literal) expressions of the functions T, [PI], [PHI] (2) without auxiliary coordinates; differentiates them according to coordinates [x.sub.i], according their derivatives [[??].sub.i] and time t (Eq. (1)); transforms them into the operator form and solves in respect of the coordinates [x.sub.i], i.e. generates the Eq. (11); forms the matrices Mark, Den and Denc (Table) with all the data required for formation of the structural diagram and Simulink-model of the system.

In the matrix Den, the specific analytic (literal) expressions of the coefficients [d.sub.i,j], [d.sub.i,i], [g.sub.i,e] in the parameters of the dynamic model of the system under discussion are provided, and in the matrix Denc their numeral values are provided. The matrix Mark shows the order of dislocations of the values in the matrices Den and Denc (each line of the matrix Mark corresponds to the lines of the matrices with the same series number). Let's analyze structure of these matrices.

The four-column matrix Mark consists of n groups of lines. Each i-th group of lines (i = 1, 2, ..., n) conforms to the [x.sub.i] -th generalized coordinate in respect of which the i-th Eq. (11) is solved. The number of lines of a group is [r.sub.i] = 1 + [k.sub.i] + [u.sub.i]; where [k.sub.i]--the number of the coordinates [x.sub.s] in the right part of the i - th Eq. (11), (1 [less than or equal to] [k.sub.i] [less than or equal to] m - 1), [u.sub.i]--the total number of generalized output forces [F.sub.i] and kinematic excitations [N.sub.i,e] in the same Eq. (11) (0 [less than or equal to] [u.sub.i] [less than or equal to] [v.sub.i] + 1). The total number of all lines of the matrix Mark is: [psi] = [r.sub.1] + ... + [r.sub.n]. Groups of the lines are situated in this matrix in the same order as the coordinates [x.sub.i] in the first line of the matrix K, i.e. the group of lines that corresponds to the coordinate [x.sub.1] is in the beginning of the matrix, then the group of lines that corresponds to the coordinate [x.sub.2] follows and so on.

The first element of a line of any group of lines is the number of the line in the matrix (the lines are numbered consecutively starting from the number 1 and ending by the number [psi]). Other three elements of the line depend on the group the line belongs to and on its place in the group.

Let's suppose that we have the i-th group of lines. All third elements of this group of lines (the third column of the group) are the serial number of the coordinate [x.sub.i] in the matrix K, i.e. the number "i" (i = 1, 2, ..., n). The second element of the first line of the group is the name of the coordinate [x.sub.i] and the fourth element is symbolic inscription "den " showing that the line with such inscription is to be used as a denominator of the coordinate [x.sub.i].

The other lines of the group are usable for marking the names of the values included in the numerator of the right part of the Eq. (11). The fourth element of all these lines is the symbol "num" showing that the value mentioned in the line is in the numerator of the right part of the Eq. (11). The second element of the second line of the group is the name of the generalized force [F.sub.i]. Starting from the third line, total ki lines are used for the names of the coordinates included in the right part of the i-th Eq. (11). The second elements of all said lines are names of coordinates; the lines that correspond to the said coordinates are situated in the same order as the names of the said coordinates in the matrix K. The lines of the last, i.e. i-th group (the number of them is [v.sub.i] - Eq. (9)) are used for listing the kinematic excitations [N.sub.i,e] that impact the system. Their names are the second elements of the said lines.

In matrix Den, the analytic (literal) expressions of the coefficients [d.sub.i,j], [d.sub.i,i], [g.sub.i,e] are provided; their dislocation in this matrix is coordinated with dislocation of the elements of the second column of the matrix Mark. The said coefficients are polynomials of the second degree in respect of the operator p according to the Eqs. (7), (8) and (10) and are described by the following matrices lines

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

Three-column and [psi]-line matrix Den consists of the lines (25); they are supplemented by the lines

[f.sub.i] = [0 0 1] (26)

for description the individual coefficients at forces [F.sub.i] in the Eq. (11); their analytic expressions are developed by the investigator. The elements [a.sub.i,j], ..., [[xi].sub.i,e] of the lines Eq. (25) in the matrix Den are provided in the analytic (literal) form expressed in the parameters of dynamical model of the system.

The name of the coordinate [x.sub.s], force [F.sub.i] or kinematic excitation [N.sub.i,e] written in the second elements (the elements of the second column) of lines of the matrix Mark corresponds to the one of the lines Eqs. (25) or (26) of the matrix Den (that includes one of the said names) dislocated according to the same order as the lines of the matrix Mark.

Thus, both the matrix Den and the matrix Mark consist of n groups of lines. Each group includes [r.sub.i] lines and corresponds to one of the coordinates [x.sub.i] (i = 1, ..., n). The first line of the i-th group of the matrix Mark (its second element is the name of the coordinate [x.sub.i]) corresponds to the line [d.sub.i,i] of the same group of the matrix Den. The second line includes the force [F.sub.i] and it corresponds to the line [f.sub.i] of the matrix Den. Analogously, one of lines Eq. (25) of the matrix Den corresponds to relevant other lines of the i-th group of the matrix Mark.

Matrix Denc differs from matrix Den only by numeral values of the coefficients [a.sub.ii], ..., [g.sub.i,e] provided in the Eq. (25) instead of their analytic (literal) expressions.

When the data of the matrices Mark, Den, Denc are available, structural diagram of the system under discussion is developed in accordance with the instructions of MATLAB / Simulink set of programs and its simulation is carried out. In addition, the data provided in the matrices are sufficient for generation of the equations of the type Eq. (11) for the systems under discussion that is required for the application of various other methods of computation and research.

In addition to the data provided in the Table, pictograms of blocks of transfer functions that present a basis of the structure of Simulink-model of the object under discussion are developed.

For convenience of simulation, the Simulink-model of the whole object is divided to subsystems that correspond to relevant coordinates [x.sub.i]. In Fig. 4, a, the subsystem that corresponds to the coordinate [[phi].sub.x] of the example (Subsystem 2) is shown.

The name Trans is automatically classified to the blocks of transfer functions; after it, a fraction follows; the numerator and the denominator of the fraction specify the numbers of the lines of the matrices Mark, Den, Denc used in the block (Table).

Blocks of transfer functions are provided in a compact form.

The pictogram of block Trans 2/1 in Fig. 4, a is enlarged to show more clearly the mark Num(s)/Den(s) (inside its contour) showing that the numerator and the denominator of the transfer functions are lines of the matrix.

The fraction following the word Trans, for example, 2/1 (Fig. 4, a), shows the lines of the matrix Denc where the values of the coefficients of the numerator and the denominator of the transfer function are specified.

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

In Fig. 4, b, the window for entering parameters of block Trans 2/1 is shown (such entering considerably reduces the probability of entering wrong coefficients of the transfer function). In the column "Numerator coefficients", the program shows: [Denc(2,:)]. This means that all elements of the second line of matrix Denc are used in the numerator of the transfer function. In the column "Denominator coefficients", the program shows: [Denc(1,:)]. This means that all elements of the first line of matrix Denc are used in the denominator.

The connecting lines of the whole model and its subsystems in Simulink-models are connected by the investigator using the graphical editor of Simulink set. The fully connected system of the example is shown in Fig. 5, a. The simulated transitional process obtained on an abrupt change of the road unevenness [z.sub.1] = 0.1m is shown in Fig. 5, b.

5. Conclusions

1. The methodology of investigation of linear stationary dynamic models of the mechanical part of mechanical and mechatronic equipment in the environment MATLAB/Simulink programs required for computer-aided generation of structural diagrams (Simulink-models) is provided.

2. On the base of the analytic (literal) expressions of kinetic energy, potential energy, dissipation function and generalized forces included in Lagrange equations of second type, the literal or numerical expressions of transfer functions included in the structural diagram are generated in the computer-aided way.

3. In course of generation of the structural diagrams, data of independent value required for the generation of differential equations of the second degree for the system under discussion and development of the normal form of the said equations are obtained.

Received February 11, 2011

Accepted June 15, 2011

References

[1.] Augustaitis, V.; Gican, V.; Sesok, N.; Iljin, I. 2006. Application of MATLAB programmes set for individual drive component investigation of a section of the printing press, Mechanika 6(62): 17-24.

[2.] Jonusas, R.; Juzenas, E.; Juzenas, K. 2010. Analysis of some extreme situations in exploitation of complex rotary systems, Mechanika 1(81): 53-57.

[3.] Tumonis, L.; Schneider, M.; Kacianauskas, R.; Kaceniauskas, A. 2009. Comparison of dynamic behaviour of EMA-3 railgun under differently induced loadings, Mechanika 4(78): 31-37.

[4.] Zavadskas, E.K.; Turskis, Z.; Tamosaitiene, J. 2010. Risk RISK assessment of construction projects, Journal of Civil Engineering and Management 16(1): 33-46.

[5.] Rimkeviciene, J.; Ostasevicius, V.; Jurenas, V.; Gaidys, R. 2009. Experiments and simulations of ultrasonically assisted turning tool, Mechanika 1(75): 42-46

[6.] Daunys, M.; Putnaite D.; Bazaras, Z. 2009. Principles for modeling technological processes investigation into the strength and durability of automatic coupler SA-3 in railway carriages, Transport 24(2): 83-92.

[7.] Vansauskas, V.; Bogdevicius, M. 2009. Investigation into the stability of driving an automobile on the road pavement with ruts, Transport 24(2): 170-179.

[8.] Mamcic, Stanislav; Bogdevicius, Marijonas. 2010. Simuliation of dinamic processes in hydraulic accumulators, Transport 25(2): 215-221.

[9.] Sesok, N.; Iljin, I. 2010. The investigation on the impact of the contact rings of cylinders and on transversal vibrations of blanket and plate cylinders in offset printing press, Journal of Vibroengineering 12(1): 95-101.

[10.] Bilyk, V.; Korobko, E.; Binshtok, A.; Bubulis, A. 2010. Simulation of oscillation dynamics of mechanical system with the electrorheological shock-absorber, Journal of Vibroengineering 12(4): 419-429.

[11.] Volkovas, V.; Perednis, A. 2010. Adaptable vibration monitoring in rotor system, Journal of Vibroengineering 12(4): 396-406.

[12.] Dorosevas, V. The analytical method of acoustic field estimation in the cylindrical shape, Journal of Vibroengineering 12(4): 589-593.

[13.] Chernykh, I.V. 2003. Simulink: Creation Environment Engineering Applications. Moscow, Dialog--MIFI. 496p. (in Russian).

[14.] Dabney, J.; Harman, T.L. 2004. Mastering Simulink. Pearson/Prentice Hall. 376p.

[15.] Dorf, R.C.; Bishop, R.H. 1998. Modern Control Systems. Addision-Wesley. 832p.

[16.] Augustaitis, V.K. 1994. Analysis and Synthesis of Vibration within the Machines As Applied to Drive Units of Machine Tools. Vilnius, Technika. 308p. (in Russian).

[17.] Augustaitis, V.K. 2000. Fundamentals of Mechanical Vibrations. Vilnius, Ziburis. 319p. (in Lithuanian)

[18.] Lojcanskij, L.G.; Lurie, L.I. 1983. Course of Theoretical Mechanics. Vol. II. Dynamics. 2nd ed. Moscow, Nauka. 640p. (in Russian).

Table
The elements of the matrices Mark, Den
and Denc found on solving the example

       Mark                     Den                  Denc

[1, phix, 1, den]    [Jx, 4*h*l2^2, 4*k*l2^2]      3.6000e+003
[2, Mx, 1, num]                     [0, 0, 1]                0
[3, phiy, 1, num]                 [Jxy, 0, 0]      6.0000e+002
[4, Z, 1, num]                   [m*ys, 0, 0]      2.5000e+002
[5, z1, 1, num]             [0, -h*l2, -k*l2]                0
[6, z2, 1, num]             [0, -h*l2, -k*l2]                0
[7, z3, 1, num]               [0, h*l2, k*l2]                0
[8, z4, 1, num]               [0, h*l2, k*l2]                0
[9, phiy, 2, den]    [Jy, 4*h*l1^2, 4*k*l1^2]      6.0000e+003
[10, My, 2, num]                    [0, 0, 1]                0
[11, phix, 2, num]                [Jxy, 0, 0]              600
[12, Z, 2, num]                 [-m*xs, 0, 0]     -1.0000e+003
[13, z1, 2, num]            [0, -h*l1, -k*l1]                0
[14, z2, 2, num]              [0, h*l1, k*l1]                0
[15, z3, 2, num]            [0, -h*l1, -k*l1]                0
[16, z4, 2, num]              [0, h*l1, k*l1]                0
[17, Z, 3, den]                 [m, 4*h, 4*k]      5.0000e+003
[18, P3, 3, num]                    [0, 0, 1]                0
[19, phix, 3, num]               [m*ys, 0, 0]      2.5000e+002
[20, phiy, 3, num]              [-m*xs, 0, 0]     -1.0000e+003
[21, z1, 3, num]                    [0, h, k]                0
[22, z2, 3, num]                    [0, h, k]                0
[23, z3, 3, num]                    [0, h, k]                0
[24, z4, 3, num]                    [0, h, k]                0

       Mark                       Denc

[1, phix, 1, den]      5.5488e+007     2.3120e+005
[2, Mx, 1, num]                  0               1
[3, phiy, 1, num]                0               0
[4, Z, 1, num]                   0               0
[5, z1, 1, num]       -1.6320e+003    -6.8000e+004
[6, z2, 1, num]       -1.6320e+003    -6.8000e+004
[7, z3, 1, num]        1.6320e+003     6.8000e+004
[8, z4, 1, num]        1.6320e+003     6.8000e+004
[9, phiy, 2, den]      1.3997e+004     5.8320e+005
[10, My, 2, num]                 0               1
[11, phix, 2, num]               0               0
[12, Z, 2, num]                  0               0
[13, z1, 2, num]      -2.5920e+003    -1.0800e+005
[14, z2, 2, num]       2.5920e+003     1.0800e+005
[15, z3, 2, num]      -2.5920e+003    -1.0800e+005
[16, z4, 2, num]       2.5920e+003     1.0800e+005
[17, Z, 3, den]        7.6800e+003     3.2000e+005
[18, P3, 3, num]                 0               1
[19, phix, 3, num]               0               0
[20, phiy, 3, num]               0               0
[21, z1, 3, num]       1.9200e+003     8.0000e+004
[22, z2, 3, num]       1.9200e+003     8.0000e+004
[23, z3, 3, num]       1.9200e+003     8.0000e+004
[24, z4, 3, num]       1.9200e+003     8.0000e+004

V. K. Augustaitis, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E- mail: vytautas.augustaitis@vgtu.lt

V. Gican, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E-mail: vladimir.gican@vgtu.lt

N. Sesok, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E-mail: pgses@vgtu.lt

I. Iljin, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E-mail: pgilj@vgtu.lt