Mechatronic Approach To The Machine Tool Design
Mechatronic Approach To The Machine Tool Design
Abstract: In the paper we consider the possibilities of complex approach to the design of ma-
chine tools in the light of the newest means which the designer can use on the earlier stages of the
design procedure. Same of components is delivered with the controllers in the form of the micro-
processors. In that case it is useful to change nominal software to adjust the control law to our
purposes. So more we can add additional control loops to improve dynamics of the machine
control methods are applied to the mechanical sys- more flexible and easier for implementation of “intel-
tems. They are called robust methods since such ligence”. At the same time their stiffness may be
system is robust to the parameter changes and to different than in the case of pure mechanical connec-
external excitations. Robust controller has constant tions.
parameters and assures sufficient stability and per- An important possibility of “mechatronic ma-
formance of the closed-loop system in the broad chine tools” is reconfiguration. Modules bought
range of parameters of the operation points. from the catalogued items and used in modern ma-
In the paper we consider the possibilities of chine tools must have unified (and frequently stan-
complex approach to the design of machine tools in dardized) interfaces. It makes possible changing
the light of the newest means which the designer can configuration of the machine tool after some time of
use on the earlier stages of the design procedure. its use. The reconfiguration may change the machin-
Same of components is delivered with the control- ing parameters and allow the user of the machine
lers in the form of the microprocessors. In that case tool to be much more flexible and respond to the
it is useful to change nominal software to adjust the changing demands of the market quicker and with
control law to our purposes. So more we can add smaller expenses.
additional control loops to improve dynamics of the But the best way to reconfigure system is the
machine. Such problems will be considered in our introduction of proper software to the controllers. It
paper. allows the low cost modernization of the machine
tools. If we add the control systems of structure
2. CONFIGURATION OF MACHINE parameters we can design the smart machine tools.
TOOLS
3. MECHANICAL STRUCTURE
A designer of a new machine tool has much more
MODELING
freedom than a designer of a new car. In the car
design there are just a few configurations of the Classic motion equations of flexible or mass-
main modules and a designer has a limited choice at lumped structure with a finite number of degrees of
the stage of a conceptual design. In the field of ma- freedom are as follows:
chine tools there are many possible configurations of
modules which can fulfill: the required parameters Mx
&& + Cx& + kx = f (1)
of workpieces to be machined, NC movements and
additional functions. where x and f are the vectors of generalized dis-
The freedom of choice requires responsibility of placements (translations and displacements) and
a machine tool designer. It is very important for the forces (point forces and torques) and M, K and C are
final result what configuration of the machine tool has respectively the mass, stiffness, and damping matri-
been chosen at the conceptual stage of design. The ces, they are symmetric and semi positive definite.
know-how of the design team and a heuristic ap- M and K arise from the discretization of the struc-
proach may be not enough. A certain method of multi ture, usually with finite elements. A lumped mass
body simulation, BMS, is propose as a tool for the system has diagonal mass matrix while the finite
configuration choice based of estimation of the work- element method usually leads to non-diagonal mass
ing stiffness of different mechanical structures [3]. matrix. In the complex machine tools the lamped
In modern manufacturing monitoring is more mass system also have non-diagonal mass matrix
and more popular, frequently as a part of automatic [3].
supervision [15]. Automatic supervision attempts to The damping matrix C represents the various
eliminate the influence of disturbances, to guard dissipation mechanisms in the structure, which are
machine tool or/and workpiece and to optimize pro- usually poorly known. Therefore it is customary to
duction. The term monitoring was adopted from assume hypothesis about Rayleigh damping:
medicine where it is as a kind of automatic, simpli- C = α M + β K , where coefficients α, β are selected
fied diagnostics. Diagnostics, which is directed only to fit the structure under consideration.
at a chosen features of monitored object, but is ac- Let us introduce modal transformation:
tive all the time. In hospitals the special sensors are x = Φz , where z is the vector of modal coordinates.
mounted on seriously ill patients and the signals The transformation leads to decoupled modal equa-
from the sensors are continuously and automatically tions:
analyzed from the point of view of selected features
z + 2ξΩz& + Ω 2 = μ -1Φ T f
&& (2)
(e.g. pulse, blood pressure... ). When the value of
selected feature is identified as critical the hospital where: Ω = diag (ωi ) is the matrix of natural fre-
personnel is automatically called to intervene. quencies, μ = diag ( μi ) is the matrix of modal
Mechatronic modules with their own drive and 1 α
sensor units are connected by information signals masses, and ξ = diag (ξi = [ + βωi ]) is the matrix
with other mechatronic units and with the control 2 ωi
system of the machine tool. Such structure is much of model damping.
EDIProD’2004 167
The transfer function between the force f as an As we show later the bandwidth can be
input and z as an output is in the following form: wider in the case of proper choice of control law. A
robust control system seems to be a perfect choice in
G (ω ) = ⎡⎣ −ω 2 M + jωC + K ⎤⎦ = ∑
−1 n φφ i i
T
the case of loaded machines. To design the robust
i =1 μ (ω − ω + 2 jξ ωω )
i i
2 2
i i
controller we should consider all interactions of the
(3) control plant with surroundings and take into acount
Typical Bode characteristics of the transfer function all spectrum of nonlinear plant parameters (which
(3) are shown in the Fig.1. are conected for example with operation points)
when we are designing of linear model of the control
plant [4].
The interactions and changes of parameters
are represented by weigth functions in th form of
transfer functions as it is shown in the Fig. 2. The
weight functions are chosen in such a way to create
the systems parameters which are able to fulfill the
assumptions. So, the weight disturbances Wo(s) i
Wi(s) describe the contents of the disturbances di i do
or they may be used to model disturbance power
spectrum which depends on the nature of signals
involved in the practical systems. The weight Wn(s)
describes the frequency model of the sensor noise.
From the opposite side, the weight We(s) reflects the
requirements put on the shape of the closed-loop
Fig.1. Typical frequency characteristics of a me- transfer functions, i.e. on the shape of the output
chanical structure. sensitivity function. Similarly, the weight Wu(s)
describes some restrictions put on the control or
For control purposes we usually reduce the model to actuator signals. The last weight Wr(s) is an optional
m lowest modes: element used to achieve desired shape of command
φ (k ) 2
or to represent a nonunity feedback system in an
G (ω ) = ∑
m
+R (4)
i
One of the aims of vibration control system of ma- Fig.2. Functional diagram of system for robust
chine tools is to make the static flexible matrix ele- analysis purposes.
ments as small as possible.
The way of choice of the weights is shown in Fig. 3
4. DRIVES and is described by equations (7). The control error
weight We(s) is associated with the sensitivity
Electric and hydraulic drives are used in machine function S(s). The sensitivity function |S| should
tools. Usully we assume that drives work perfectly. have the low value for the low frequencies. In this
It means that they have amplifications equal one (0 range the control signal and disturbance signal play
[dB]) in the sufficient bandwith.. It is true in the case a meaningful role. So, the sensitivity function must
of slowly working machine tools. In this case we can fulfill the requirement |We(s)⋅S(s)| ≤ 1, where the
omit the dynamics of the drives. We should take into weight We(s)⋅is described by equation (7). The form
account the dynamcs of servo-drives when the ma- of the weight is realized in systems with steady-state
chines work with acceleration reaching 3[g]. The error. In considered systems we assume the error
bandwidth of servo-drives depends on many causes: lower than ε. This is ensured when |S(0)| ≤ ε. So, for
loads, operation point, control law, external distur- |We(0)|≥1 the norm is ||We(s)S(s)||∞ ≤ 1. The similar
bances, and so on.
168 EDIProD’2004
analyses can be used for other weights and functions troller were optimised in Matlab by using of NCD
of the system. Blockset Toolbox [9]. As a result we obtain the
following PID controller parameters: Kp=1 – pro-
k
⎛ ωbc ⎞
k portional gain; Ki=0.25 – gain of integration part;
⎛ 1 ⎞
⎜ k M s + ωb ⎟ ⎜ s+ k ⎟ Kd=0.01 – gain of differential part.
⎜ Mu ⎟ We assume in our case that the following
We ( s ) = ⎜ s ⎟ Wu ( s ) = (7)
⎜ s + ωb k ε ⎟ ⎜ k ε1 s + ωbc ⎟ parameters change during the electric drive opera-
⎜ ⎟ ⎜⎜ ⎟⎟ tion: k E - voltage constant of drive; Tm - mechani-
⎝ ⎠ ⎝ ⎠
where: k ≥ 1, ωb – the cut-off frequency for function cal time constant; Te - electrical time constant; in
S(jω), ε – the lower restriction of S(jω), Ms – the range ±10% of nominal value. For example,
the upper restriction of S(jω), ωbc – the cut-off electrical time constant can have values from the set
frequency for function R(jω), {0.252, 0.28, 0.308}. Transfer function of closed-
ε1 – the lower restriction of R(jω), Mu – the loop system with any nonnominal parameters will be
upper restriction of R(jω). denoted by GP(s).
Let us introduce relative error in the follow-
ing form: e( jω ) = GP ( jω ) − G ( jω ) . The choice of the
G ( jω )
weight is determined by condition that for all fre-
quencies the Bode diagram of weight should be
above the plot of relative error. To do this we chose
the weight transfer function in the
−3
form: ω u (s ) =
1 .5 ⋅ 10 ⋅ s
Fig.3. Limitations put on the weights functions. 1.5 ⋅ 10− 7 ⋅ s + 1
The next considered weight is connected with
The nominal model of the plant together with the scaling of the controller output signal u. For
weight functions forms an extended model of the simplicity we assume a scalar weight where its mag-
open-loop system. For such model we can obtain the nitude is as follows: ω I = 7 .
optimal controllers using of commercial software, In our case after calculations and simplifications
for example, Matlab-Simulink. we have obtained the controller described by the
following transfer function:
4.1. Electric drive
49.86 ⋅ s 2 + 2600.11 ⋅ s + 9096.58
Modern control theory has recently been widely K (s ) = .
applied to design the controllers for electric drives s 3 + 773.35 ⋅ s 2 + 4141.89 ⋅ s + 4.14
[7,8]. The reasons for this are following: 1) the ne-
cessity of meeting increasingly stringent require- The results of computer simulations of closed-loop
ments on the performance of drive control systems, system for both controllers are shown in the Fig. 4.
ROZRUCH SILNIKA
2) easy access to modern power semiconductor 5000
3000
practically applicable to microprocessor-based drive 2500
control by taking into account physical constraints 2000
such as input delay times and input and state vari- 1500
able constraints. 1000 PID
The object of consideration is DC motor D – 500
Hinf
trollers the start phase during which electric drive piston and servo-valve displacements ΔXV=ΔX =0 ,
reaches its operation speed was recorded on the lab and for three flow rates.
stand and transient processes are shown in Fig.5. 5600
ZMIANA nzad Z 4400 NA 5400 [obr/min]
OBROTY [obr/min]
4000
5000
3500
4800
3000
OBROTY [obr/min]
2500 4600
nzad
2000 PID
4400 Hinf
1500
nzad=4400 [obr/min]
1000 4200
PID -0.5 0 0.5 1 1.5
Hinf
500
CZAS [s]
OBROTY [obr/min]
3800
3200
Parameters
3000
Fig. 17. Root locus for the system with inertial con-
troller.
s 2 2 2
s + 0.5s + 5.29 s + 14s + 400 s + 10s + 22500
150s + 1 2
s + 3s + 100
2
s + 15s + 3600 2
s + 4s + 32400
1
s +1
2
1
s + 2s + 1
5.2. Experimental results
Fig.15. Model of system with second-order filter as
the controller. The inertial controller was verified in experimental
way by using of it in the vibration control system of
The closed-loop system with inertial controller beam from Fig. 9. Durig transient response the
is awans stable since branches of plots are in the left control system was lock-in and the results are seen
half-plane for all controller gains. System with sec- in Fig. 19. It was calculated that logarithmic decre-
ond-order filter becmes unstable for some controller ment of damping increased from δ=0.2 to δ=0.51.
gains. For uncorrect model the system designed as The amplitude-frequency characteristics is pre-
stable one can appear unstable. sented in Fig. 20. We see the high level of the damp-
ing in the whole range of frequency. Characteristics
is similar to respective characteristics of the simu-
lated closed-loop system from Fig. 13.
172 EDIProD’2004