Design of an Adaptive Neural Network Controller for Effective Position Control of Linear Pneumatic Actuators

Size: px
Start display at page:

Download "Design of an Adaptive Neural Network Controller for Effective Position Control of Linear Pneumatic Actuators"

Transcription

1 Researh Artile International Journal of Current Engineering and Tehnology E-ISSN , P-ISSN INPRESSCO, All Rights Reserved Available at Design of an Adaptive Neural Network Controller for Effetive Position Control of Linear Pneumati Atuators Osama. A. Montasser Ȧ* and B. A. El-Sayed Ḃ Ȧ Mehanial Power Engineering, Ain Shams University, Cairo, Egypt, on leave to join the British University in Egypt, BUE Ḃ Sakr Faroty, The Arab Organization for Industrialization (AOI, Cairo, Egypt Aepted 05 Ot 04, Available online Ot 04, Vol.4, No.5 (Ot 04 Abstrat The main target of this work is to design an appropriate position ontroller for a pneumati system using artifiial neural networks. The rod position of a double ating pneumati ylinder, ontrolled by proportional linear valves, was hosen as the present ontrol system. A mathematial dynami model for the pneumati system was derived. The model shows, as it is expeted, that the pneumati system is of highly nonlinear features. This is due to ylinder-piston mehanial frition, the ompressibility feature of air, and nonlinear harateristis of the flow through a valve orifie of variable area. The model shows, as well, that pneumati system is of time varying harateristis. A Proposed Neural Network Controller, PNNC, is designed and implemented. The PNNC is a rule-based ontroller, where both the slope and amplitude of the ativation funtion of eah neuron is adapted to enhane the ontrol system performane. A onsiderable improvement of the system response for different input onditions is ahieved by applying the PNNC on the present ontrol system. The robustness and effetiveness of the proposed ontroller were verified through omputer simulations using MATLAB pakage and SIMULINK toolbox. A omparison with the Conventional Neural Network Controller, CNNC and the typial PID ontroller, assured that the present PNNC is robust and more effiient in terms of both the system stability and speed of response. Keywords: Aurate Position Control, Adaptive Learning Algorithm, Neural Network Controllers, Pneumati Atuators, Sigmoid Ativation Funtion.. Introdution Hydrauli atuators are relatively easy to be positional ontrolled than the pneumati atuators. This is due to the high nonlinearity assoiated with pneumati media flow. A self-tuning fuzzy PID ontroller was suessfully developed and applied to an eletro-hydrauli atuator by (N. Ishak et al,. The proposed ontroller offers promising apabilities to guarantee the position traking auray of their hydrauli system. However, Pneumati atuators are now a reognized alternative tehnology to other hydrauli and eletri ounterparts in areas of high speed, medium power appliations (J. Moore and C. Wong, 000. Pneumati systems are apable of providing suitable power output at relatively low ost, in addition they are lean, of light weight and an be easily servied. Pneumati atuation has been mainly used in the form of logi ontrol systems. Many appliations in industry use pneumati power in open loop ontrol systems where the strokes of the moving parts are usually fixed using mehanial stops. However, many other appliations require the position ontrol of pneumati ylinders suh as, robots, gating systems and *Corresponding author: Osama. A. Montasser CNC mahines. (N. Tillet et al, 997 provided that pneumati systems have been limited to a simple end-stop appliation beause of the diffiulty of providing proportional ontrol of the rod position. This diffiulty stems from the nonlinear nature of air as power transmission medium and non-linearity assoiated with pneumati omponents (M. Shih and M. Ma, 998. (Verseveld and Bone, 997 proposed position ontrolled pneumati atuator. This pneumati system onsists of a double ating ylinder, with low frition option, onneted to a horizontal linear slide. A pneumati brake mehanism is inorporated within the linear slide to lok the atuator firmly in its position one the desired steady state auray is ahieved. Two standard three-way solenoid valves were used. The valves were atuated using a Pulse Width Modulation algorithm, PWM. A PID ontroller with added frition ompensation and position feed-forward is suessfully implemented. As it is well known, PID ontrollers are widely used for proess ontrol beause of their good performane. However, it is very diffiult to self tuning the ontrol parameters of these ontrollers. (K. Fujiwara and Y. Ishida, 995 desribed a method of using neural networks to self-tuning the PID ontrollers for pneumati ylinders. They show simulation and 3498 International Journal of Current Engineering and Tehnology, Vol.4, No.4 (Ot 04

2 experimental results to demonstrate the effetiveness of their model. (M. Shih and M. Ma, 998 have proposed the applying of a fuzzy ontrol tehnology to ontrol the position of a pneumati ylinder. In their work, a modified differential PWM ontrol method was applied. There are numerous works applied artifiial intelligene in the ontrol appliations but few of them are used in pneumati ontrol systems. (X. Cui and K. Shin, 99 used an Artifiial Neural Network Controller, ANNC, to ontrol the middle-point temperature in a thermal power plant. The ANNC was designed to deal with the nonlinearity due to negative effets of long response delays, saturation onditions and proess noise. The ontroller onsists of four-layer, with two hidden layers, pereptrons. Their work introdued a simple algorithm based on the Bak-Propagation, BP, for a lass of nonlinear systems, typified by proess ontrol appliations. The proposed NN ontroller is trained using the system s output errors diretly, with little a priori knowledge of the ontrolled plant. (M. Abdelhameed, 999 introdued an adaptive neural network ontroller for traking a robot of n-degree-offreedom. The learning algorithm is based on the adaptive updating of the weights of the network to minimizing the quadrant traking error and its derivative of eah robot arm. (G. Geng and G. Geary,997 proposed a method whih uses neural networks together with reursive least square algorithm to model nonlinear proesses. (G. Cembrano et al, 997 desribed the use of Cerebellar Model Artiulation Controller (CMAC networks for the adaptive dynami ontrol of an orange harvesting pneumati robot. A neural-network model was preferred sine it provided better approximation apabilities. Finally, it is notied that the appliation of onventional ontrol methods on position ontrol of pneumati systems have some restritions. Furthermore, there are relatively few works related to the appliation of artifiial neural network ontrollers in the field of pneumati atuation systems. It therefore, the main issue onsidered in the present work is to design an appropriate position ontrol system for a pneumati atuator using artifiial neural networks. The Artifiial Neural Network Controller, ANNC was hosen to position ontrol the piston movement of the present pneumati system. This is deided to get the benefits of its intelligene and learning apability needed to ope with the higher inflexibility and nonlinearity of the pneumati systems. A rule based NN ontroller is proposed, tested and found to suessfully ahieve a good traking ontrol performane and is robust to sever hanges in the applied load as well. The present paper is organized as follows. In the following setion, mathematial and state spae models of the present pneumati system are derived. An introdution to the artifiial neural network based ontrollers is reviewed afterward. In the next setion, the Conventional on-line self learning based Neural Network Controller, CNNC, is introdued. The omputational steps used to build a CNNC are listed in subsequent setion. A rule based neural network ontroller is proposed and designed in the following setion, PNNC. Results of the present CNN and PNN ontrollers are onsequently disussed. The onlusions of this study are summarized in the last setion.. Modeling of the present pneumati system In this setion a state spae model of the present pneumati system to be positional ontrolled is derived. Desription of the system is presented and the mathematial model of eah of its omponents is derived... Pneumati system desription Fig. Layout of the present pneumati system The present pneumati system is shown in Figure. It onsists of a double ating ylinder, linear pneumati atuator, ontrolled with two proportional valves, v and v, and arrying an external load M l. Air is disharged from a pressurized air supply to the right hand side of the ylinder through valve, v. The air is disharged from the ylinder to the atmosphere through valve, v. The pressure at the left hand side of the ylinder is held onstant by attahing it with a large aumulator, air pressure in whih, P l, is.9 bar. The funtion of the aumulator is to maintain a reasonable differene of air pressure between right and left sides of the ylinder. This redues the nonlinearity assoiated with large pressure differene between both sides of the ylinder. The value of the supply pressure, P s is fixed at six bar. When v is opened and v is losed, the air flows through v to the right side of the ylinder, whih is initially at fore balane state, to inrease its pressure, P. As P inreases, the piston of the ylinder moves to the left to the desired position, y ref. v is kept opened, as the position of the piston y(t is less than the desired piston position y ref. When the piston position y(t exeeds the desired piston position, y ref, v is losed and valve v is opened to reverse the diretion of the piston motion to the right. The air at the right side of the ylinder disharges to 3499 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

3 the atmosphere and P thus dereases. This lets the piston moves to the desired position at the right side and when reahed, the valve v is losed... Mathematial modeling In order to derive the governing equation for the motion of the ylinder piston, the momentum equation is onsidered as follows: M y P A P A f y ( l pr l pl Where: P, P l are the air pressure at the ylinder right and left side respetively and A pr and A pl are the ross setional areas of right and left sides of the ylinder respetively, f is the visous damping oeffiient, y and y are veloity and aeleration of the piston rod and M l is the equivalent mass of piston head and the load. The energy equation is applied to obtain a mathematial model for pressure dynamis of the air in the ylinder. Variation in air pressure at the ylinder left side, P l, is ignored sine it is onsidered onstant due to the aumulator effet. Considering the right side of the pneumati ylinder as a thermodynami ontrol volume, the governing equation for air pressure P is obtained by applying the energy equation as follows: d m u m C T m C T dv Q P p p s ( dt dt Where: m and m are the air mass flow rates through valves v and v respetively, Q is the rate of heat transfer aross the ylinder wall, P, V, m and u are pressure, volume, mass and speifi internal energy of air at the ylinder right side, P dv dt is the rate at whih work is done on the moving piston, T s is the supply air temperature and T is the air temperature at the ylinder right side, C m T m T is the enthalpy differene through the p s ontrol volume boundaries and dm u dt is the rate of hange of internal energy of the ontrol volume. Sine the mass transfer through the system is muh faster than the heat exhange through the system boundary, so, Q an be negleted. It was thus reasonably assumed that the system undergoes an adiabati proess. Applying thermodynamis relations for ideal gases, equation ( is redued to: dp dt P dv. V dt R V T m T m Where: is speifi heat ratio, R is the air gas onstant. While opening valve v, the status of air mass flow rate through it, m, hanges from moment to moment. It may be direted to inter the pneumati ylinder, inflow ( m, or may direted to leave it, outflow ( m. This is due to the high dynamis of air pressure at the right side, P, due to whih, it may take a lower or higher value than that of s (3 the supply pressure, P s. Outflow, m, should take a negative sign when it is substituted in energy equation (3. On ontrary, the air mass flow rate through valve v, m is always direted to the atmosphere, sine P never be lower the atmosphere pressure P a. Relations for m and m are onluded as follows: m (4 Av v v Cdv m (5 Av v v Cdv Where: and ρ is the air veloity and density at the valve orifie, A is the partially opened valve orifie area, A varies from 0% to 00% of the orifie area, and C d is the valve disharge oeffiient. Subsripts v or v refers to the valve number. Relations for the air density and veloity, ρ and, through both proportional valves are obtained by applying the ontinuity equation and the thermodynami relations for ideal gases as follows: Proportional valve Two possible onditions for the redued pressure ratio through the valve, P R, are onsidered, namely the subsoni ondition for P R > 0.58 and supersoni ondition for P R Relations for the two onditions are shown as follows: a Subsoni ondition a- Inflow ondition, P < P s, P R = P / P s > 0.58: P v (6 R T ( P s R R / v Ts ( PR ( ( vav / s As (7 a- Outflow ondition, P > P s, (P R = P s / P > 0.58 : Ps v (8 R T ( P R R / v T ( PR ( ( vav / A (9 where: ρ s is the density of supplied air and A s and A are the ross setional areas of the onnetions of air supply and the ylinder with the valve v respetively. b Supersoni ondition For inflow ondition, supersoni flow is possible and the air density relations are derived as follows: b- Inflow ondition, P < P s, P R = P / P s 0.58 : 0.58 Ps v (0 R Ts 3500 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

4 v R T s ( b- Outflow ondition, P > Ps: Outflow of air from the right side of the pneumati ylinder to the air supply is always subsoni and annot be supersoni sine P s / P never be possibly less than or equal Proportional valve Sine P > P a in all working onditions, flow of air from the ylinder to the atmosphere through valve is always outflow. a Subsoni ondition, P / P a > 0.58: v Pa ( R T ( P R R / v T ( PR ( ( vav / A (3 b Supersoni ondition P / P a 0.58: 0.58 P v (4 R T v R T Where: P a is the atmosphere pressure..3 State spae model (5 The state spae model of the pneumati system ould be derived as follows: Variables x, x, x 3 are defined as: x y, x x y, x3 P (6 From equations, 3 and 6, the following relations are obtained. x (7 y ax3 a a3x x x 3 x 3 P a4 a5 (8 x x Using the previous derived relations to 5, a to a 5 are alulated as follows: Apr a (9 M a a 3 l l Pl Apl (0 M f ( M l a ( 4 a 5. pr 350 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04. R ( Ts m T m (3 A The model, derived above, shows that dramati hanges always our in the density of air in pneumati systems. This is due to the high ompressibility essential feature of all gases. 3. Artifiial neural network based ontrollers A great attention has been devoted, in reent years, to the artifiial multilayer neural networks, ANNs, whih have proved to be extremely suessful in pattern reognition and ontrol problems. The main feature of the ANNs is their apability of dealing with nonlinear ontinuous funtions. That is why their appliations for ontrolling omplex nonlinear systems have beome an essential matter. In the present work, an ANN ontroller is proposed to be used in positional ontrol of the present pneumati system. Fig., below, shows the arhiteture of an artifiial neural network of one hidden layer. It onsists of three onseutive layers, namely, the input layer, of n dimensions, the neural network hidden layer, of N neurons, and the output layer, of m dimensions. Eah of the input and the output vetors is onneted with eah neuron of the hidden layer through a synapti weight, whih is a real fration number and is alulated aording to an ativation funtion. We, therefore, have two weight matries, n N dimensional matrix for input-neural network layers onnetions and N m dimensional matrix for neural network-output layers onnetions. The neurons are onneted aording to layered feedforward neural networks method, where a layer of neurons reeives input only from neurons of the previous layers. Various methods ould be used to update the weights of the neural network layer, what is so alled training the ANN, so that for a ertain input to the input layer, we an get the desired output from the output layer. One way is to set the weights expliitly using prior knowledge of the system under onsideration, open loop ontrol. The learning ategories are lassified into two distint types, namely the supervised learning and unsupervised learning. In supervised learning, assoiative learning type, the network, whih is alled self-supervised, is trained by providing it with inputs and the mathing outputs patterns. Fig. Arhiteture of an artifiial neural network of one hidden layer

5 Classi appliations of ANN ontrol systems are of supervised learning, open loop type. Reasonable auray of open loop ontrol using ANN ontrollers is attained for a ertain speified initial-desired output error e initial value. This is ahieved through pre-learning the ontroller for the speified value of e initial to get ontrol ations fairly fitting the inputs to the ontroller. Another way is to train the ANN by using a teahing pattern and letting the network sets its weights ON LINE, unsupervised learning. This is arried out aording to some predefined learning rules and an error performane index minimization algorithm. In the unsupervised learning, self-organizing type, the network is trained to respond to lusters of input patterns. Unlike the supervised learning pattern, there is no a priori ategories set, aording to whih the input-output patterns are to be lassified, rather, the system should develop its own representation of the input-output relationships. 4. The onventional on line self learning based NN ontroller, CNNC Reently, the ANN have been used to ontrol many industrial plants. The online learning by ontinuous mapping between the inputs and output fators of the ontroller enables the use of ANN ontrollers in feed-bak ontrol systems. A multilayer pereptron struture with bak propagation, BP, learning algorithm is one of the most effetive online learning algorithms. It an be used to approximate any ontinuous funtion mapping with any desired auray. For ontrol purposes, the ANN based ontroller utilizes the atual system input error to synthesize the ontrol ation. Sine, in ontrol systems, it is pratially to measure the output from the plant, the ontrolled variable y(t, rather than the output from the NN ontroller. It is therefore, the NN based ontrollers utilize the error in the ontrolled variable, the deviation between the desired and the atual output, e(t = y ref - y(t, to synthesize the ontrol ation. The error in plant output, e(t, is used also in training the NNC or in other words in updating the weights of the neural network layers. When designing the NN ontrol system, the major obstale is to train the NNC using the error in the systemoutput, e(t = y ref y(t, rather than the unknown networkoutput s output. (X. Cui and K. Shin, 99 presented a solution for this problem for onventional neural network based ontrollers CNNC. (X. Cui and K. Shin, 99, adopted the basi priniple of multilayer pereptrons with BP to update the weights of eah pereptron. Their formulas for updating the thresholds and the weights from the input to the output layers are utilized in this work to build a onventional neural network position ontroller CNNC. The piston rod position, y(t, of the present pneumati system, shown in Figure above, is ontrolled using a NN based ontroller. A blok diagram for the present position ontrol system is presented in Figure 3. A neural network ontroller, NNC, is asaded with the pneumati ontrolled system. The NNC output, ontrol ation (t, is used to modify the opening ratio, k, of the valves, v and v, to ahieve the desired ontrolled variable, y ref. Fig.3 A blok diagram for the present position ontrol system Figure 4 shows the basi struture of the present artifiial neural network. The hoie of the ANN s inputs should reflet the desired and atual status of the ontrolled variable. Therefore, the inputs of the ANN ontroller is usually the system s traking errors at the urrent time and at two preeding sampling time intervals, t, e(t, e(t - t, e(t - t. Inputs are seleted to onsider the history of the rate of hange of error signal, e(t, to simulate the integral ation in the onventional PID ontroller. Where t is the sampling time interval used in numerial alulations. t was seleted to be as small as possible to overome the numerial instability. Fig.4 Basi struture of the present artifiial neural network The number of the hidden layers depends on the ontrolled plant under onsideration, and is seleted to be two here in this work. The symbols used in Figure 4 are defined as; X: for inputs, I= N for input layer, N=3 in this work, J= N for the st hidden layer neurons, K= N for the nd hidden layer neurons, N and N are hosen to be 4 in this work, and subsripts, and 3 represent st and nd hidden layers and the output layer respetively. We have only one output from the present NNC, O =. 5. Computational steps of the CNNC Computations of the present NN ontroller, presented in Figure 4, are summarized in this setion. The on line self learned CNNC of (X. Cui and K. Shin, 99 are used here as a onventional position NN ontroller. Their 350 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

6 formulas for updating the thresholds and the weights from the input to the output layers are utilized in this work. The ativation funtion onsidered for eah neuron is of the sigmoid type. The sigmoid ativation funtion has the form: x t O e (4 Where: x(t is the output of a neuron, O(t is its input, the summation of the outputs of previous layer s neurons multiplied by the orresponding weight fators, (t is a threshold value, and is a slope parameter fixed at unity value for the present CNNC, =. Weight fators, W(t, and thresholds, (t, are initially defined and updated as desribed below. A satisfatory ontrol results ould not be attained before seleting the optimum values of the initial values of all the weight fators, W initials, thresholds, initials, learning rates, initials, and gain fators, initials. Learning rates and gain fators are used to update the weight fators and the thresholds. The seletion of these parameters is very tough and takes a very long time beause it depends on trial and error onept. To simplify the proess all members of eah of the above mentioned parameters are given same initial value. The optimum initial value of any of the parameters is obtained as follows; As the values of the other parameters are kept onstant, the value of the parameters under onsideration is hanged and at eah of its values ontrol performane measures are heked. Performane measures are onsidered as the settling time, the steady state error, and the root mean square of the error. The value of the parameters under onsideration whih gives the best ombination of the three performane measures is hosen as the optimum initial value of the parameters taken into onsideration. The same proedure is repeated for all the other parameters to get, at last, their optimum initial values. It is worth noting that there is a narrow range of initial values for eah of the NNC parameters out of whih the ontrol system is unstable. That is why it is very important to arry out the above mentioned seletion proess not only to get optimum ontrol performane but also to avoid system instability. The CNNC is trained to ahieve a ertain ontrol ondition, e initial = y ref - y initial. Computational steps are arried out aording to (X. Cui and K. Shin, 99 as listed in what follows: - For a ertain speified e initial and =, the initial values, W initials, initials, initials, and initials are seleted as mentioned above. - The outputs of the st HIDDEN layer are omputed at the urrent alulation step (t=0. 3- The outputs of the n HIDDEN layer are omputed at (t=0. 4- The output from output layer, ontrol ation, is omputed at (t=0. 5- Weights from the nd HIDDEN to OUTPUT layers are updated for the next alulation step (t = t + t. 6- Weights from the st HIDDEN to nd HIDDEN layers are updated for (t + t. 7- Weights from the INPUT to st HIDDEN layers are updated for (t + t. 8- The thresholds are updated for (t + t. 9- Steps from to 8 are repeated for the next alulation step (t = t + t. The onventional NN ontroller was found to be a single point ontroller and should be set up for eah different e initial value and for eah piston movement diretion. It is therefore it an t be used for traking position ontrol. A most smart ontroller was proposed by this study, so that a wide range of initial error values an be effiiently ontrolled. The Proposed NNC, PNNC, was heked for traking ontrol in both piston movement diretions and found to be suessful. The PNNC is presented in the next setion. 6. The proposed rule based neural network ontroller, PNNC The PNNC is proposed in this study to over a wide range of ontrol points in both piston movement diretions. The same struture, shown in Figure 4, for the CNNC is used. The slope parameter,, of the sigmoid ativation funtion, Equation 4, is fixed at unity value with the present CNNC, =. The slope value affets the degree of nonlinearity of the sigmoid funtion and thus affets the updating of the NNC parameters. It was found that when the value is hanged the optimum e initial of the CNNC is hanged than that at whih the optimum initial values of its parameters are seleted. So the proposed rule based NN ontroller, PNNC, was onstruted to ahieve a suessful traking ontrol as follows: - For forward piston movement and a unity slope parameter, =, a ertain initial error, e initial = y ref - y initial, is defined and for whih, the optimum initial values of the ontroller parameters are seleted as desribed in setion 5. - Keeping the optimum initial values unhanged, the value of the initial error, e initial, is hanged. At eah value of the initial error, the value of the slope parameter,, is hanged to meet the best ombination of the ontrol performane measures, namely, the settling time, the steady state error and the root mean square of the error value. This is arried out utilizing the trial and error tehnique. 3- A fairly good urve fitting, orrelation fator is almost equal to, is onluded for the optimum ( - e initial data for the forward piston movement. The urve fitting is resulted in a polynomial form as: s a i i0 e i initial of the polynomial terms., Where a i are onstants and s is the number 4- Steps from to 4 is repeated but for the piston bakward movement. So that two ( - e initial relations are obtained, one for eah piston movement s diretion. Having the relations the most suitable value of the slope parameter, is alulated at eah initial error value and for any piston movement diretion, and thus satisfatory ontrol results is ahieved International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

7 Displaement, m 7. Results of the present ontrol system The model of the pneumati ontrolled plant is built and all the simulation tests are arried out using MATLAB pakage and SIMULINK toolbox. forward and bakward position ontrol needs a different initial NN parameters seletion. A typial result for the piston bakward movement position ontrol is shown in Figure Results of the CNNC Typial results of the CNNC are shown in Figures 5 and 6. The figures present the ontrol system response with the onventional NN ontroller. Figure 5 shows the response when the initial error, e initial = y ref - y initial, has a value of 4 m for forward piston movement. This value is the value at whih the optimum values of the ontroller s parameters were determined and mentioned in setion 5 above. It is therefore, the system shows satisfatory results as Figure 5 illustrates. Fig.7 Control system response with CNNC set up for e initial = -.5 m 7.. Results of the PNNC 7... Single ontrol point omparison with the CNNC Fig.5 Control system response with CNNC set up for e initial = 4 m Figure 6 shows different results as the system failed to reah the referene position. The reason is that the initial error is hanged to 5.5 m, the value that is different from the value at whih the NN ontroller s parameters were determined. The PNNC showed a suessful ontrol results, as it is expeted, for single ontrol point in both movement diretions. A single universal set up is needed to over a wide range of ontrol point onditions. While the CNNC set up for e initial = 4 failed to over e initial = 5.5 m as shown in figure 6, the PNNC suesses to over multi ontrol onditions as figures 8 and 9, below, show Traking ontrol results for the PNNC Figures 8 and 9 show traking ontrol results of PNNC for two different traking trajetories. In both trajetories, reasonable settling time and steady state error are ahieved. The figures also show that, in both of movement diretions, the piston shows an osillatory movement forward and bakward on his way from its initial position to its referene position. This is beause of the high nonlinearity behavior of the pneumati systems. Y ref Fig.6 CNNC set up for e initial = 4 m failed to over e initial = 5.5 m The CNNC was tested for bakward piston movement and found to be also suessful for single ontrol point. Of ourse for the same e initial absolute value eah of the Time, s Fig.8 Traking ontrol results with PNNC, trak trajetory 3504 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

8 Displaement, m Load, % of 5.5kg Displaement, m Displaement, m Displaement, m Y ref Time, s Fig.9 Traking ontrol results with PNNC, trak trajetory 7..3Comparison with the results of the onventional PI ontroller The onventional PI ontroller is tested as a position ontroller of the present pneumati system. to be equal to 0.3, above whih the ontrol system beomes instable. Of ourse this leads to a large settling time. The integral effet was found to redue the settling time but it inreases the system instability, and thus it is exluded. Figure 0 shows a typial omparison between ontrol results of both the PNNC and the P ontroller of the largest possible gain, K = 0.3. The P ontroller shows stable results with large settling time, 5 seonds. The PNNC shows small settling time, seonds, with relatively allowable osillations whih are slowly onverged with time, as shown in the figure. The onventional P ontroller does not give satisfatory results when applied to pneumati systems. So, neural network ontrollers are proposed to do the job. 7..4Testing the PNNC against load disturbane The load disturbane effet on the PNNC performane is shown in Figure. The load value hanges in range of about ± 0.90% of a design load of 5.5 kg. The load disturbane is generated using a white noise generator built in faility of the MATLAB-SIMULINK ommerial pakage. The figure shows the robustness of the PNNC to severe flutuation in the applied load and still showing fairly good results. Y ref Time, s Time, s P ontroller, K =0.3 Y ref Time, s s Fig. Effet of external disturbane ± 90% around a design load of 5.5 kg Time, s PNN C Fig.0 Comparison between ontrol results for the PNNC and the P ontroller The proportional gain, K, has to be very small to avoid system instability. This is beause of the high nonlinearity of the ontrolled plant. The largest value of K was found Conlusions A mathematial nonlinear model is derived for a double ating pneumati ylinder plant. A rule based NN ontroller is proposed, PNNC, and applied to positional ontrol the piston movement of the present pneumati system. The robustness and effetiveness of the proposed ontroller were verified through omputer simulations using MATLAB pakage and SIMULINK toolbox. The performane of the present PNNC was tested against the PI and the onventional NN ontrollers. The main onlusions are extrated as follows: 3505 International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

9 The onventional NN ontroller, CNNC, gives satisfatory ontrol performane for a single ontrol ondition. The proposed NN ontroller, PNNC, suessfully ahieves a good traking ontrol performane for a wide range of ontrol point onditions. 3 The PNNC shows a satisfatory ontrol performane for load disturbane of about ±70% of the load design value. Appendix A The pneumati system initial onditions are set as follows: Initial pressure: 3.5 bar Air tank pressure:.94 bar Supply air pressure: 6.0 bar Supply air temperature: 98 K Air inside ylinder initial temperature: 98 K Air tank initial temperature: 98 K External load: 5.5 kg Cylinder diameter:.5 m Stem Diameter:.0 m Supply tube diameter:.0 m Disharge tube diameter:.0 m Piston stroke: 7.0 m Initial piston veloity: 0.0 m/s Gas onstant: 87 J/(kg K Speifi heat ratio:.4 ANN Controller Parameters The value of the ANN ontroller s parameters, W initials, initials, initials, and initials and the slope parameter were obtained as explained in the paper and listed below for both the forward and bakward piston s movement. Forward Piston Movement The seleted parameters are determined at = and for e initial =.50 m, y ref = 3 m and y initial = 0.5 m were: W initials = 0-6 initials = 0-3 initials = 0-3 initials = 0.. The slope parameter -initial error ( - e initial relation is obtained as 6 a i i0 e i initial Where: a 0 = 0.00, a = , a = 0.638, a 3 = , a 4 = -.466, a 5 = , a 6 = Bakward Piston Movement The seleted parameters are determined at = and for e initial = m, y ref = 3 m and y initial = 6.5 m were as follows; W initials =0-6 initials = 0. initials = 0. initials = The slope parameter -initial error ( - e initial relation is obtained as 6 a i i0 e i initial Where: a 0 = -0.6, a = , a = , a 3 =.58, a 4 = -9.93, a 5 = 6.643, a 6 = 0.83 List of symbols A Cross-setional area, m Subsripts C Cylinder inlet onnetion pl Piston left side pr Piston right side S Supply air onnetion v Valve orifie Valve orifie v ANN Artifiial neural network ANNC Artifiial neural network ontroller BP Bak Propagation C d Coeffiient of disharge, valve C d Coeffiient of disharge, valve C p Speifi heat at onstant pressure F Coeffiient of visous frition CNN Conventional Neural Network CNNC Conventional Neural Network Controller e Error signal y ref y(t m Mass flow rate through valve, kg/s m Mass flow rate through valve, kg/s M l External load, kg weight = 9.8 N NN Neural Network U Internal energy of system, J O Sum of node outputs = inputs to next layer nodes P Pressure, Pa/m Subsripts a Atmospheri Cylinder at right side l Cylinder at left side, Tank side R Redued pressure ratio through the valve s Supply air PNNC Proposed Neural Network Controller PWM Pulse width modulation Rate of heat transfer aross the wall of ylinder, Q W R Gas onstant, J/(kg o k t Time, se T Atmospheri Temperature, o K a T Temperature of air inside the ylinder, o K T s Temperature of supply air, o K V Volume of the right side of the ylinder W Weight of ANN node X Output of ANN node y Piston displaement, m y ref Piston referene position, m y initial Piston initial position, m y Piston veloity, m/s y Piston aeleration, m/s Greek letters Speifi heat ratio Slope parameter Air density at right side of ylinder, kg/m 3 s Supply air density, kg/m International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

10 v Air density after valve, kg/m 3 v Air density after valve, kg/m 3 Bias of ANN Learning rate of ANN Gain of ANN Referenes N. Ishak, M. Tajjudin, H. Ismail, M. Hezri, Y. Sam, R. Adnan, (0, PID Studies on Position Traking Control of an Eletro-Hydrauli Atuator, International Journal of Control Siene and Engineering, (5,pp Moore J. R., Wong C.B., (000, Smart Components-based servo pneumati atuation systems, Miroproessors and Mirosystems, Vol. 4, pp Tillet N.D., Vaughan N.D., Bowyer A., (997, A Non-linear model of a Rotary Pneumati Servo System, Pro. Instn Meh Engrs, Vol., Pert, pp Shih M.C., Ma M.A., (998, Position Control of a Pneumati Cylinder Using Fuzzy PWM Control Method, Int. Journal of Mehatronis, Vol. 8, pp Varseveld R.B., Bone G. M., (997, Aurate Position Control of a Pneumati Atuator Using On-Off Solenoid Valves, IEEE/ASME TRANS. On Mehatronis, Vol. No. 3, pp Fujiwara, K. K., Ishida Y., (995, Neural Network Based Adaptive I-PD Controller for Pneumati Cylinder, SICE, Sapporo, pp Cui, X., Shin, K.G., (99, Appliation of Neural Networks to Temperature Control in Thermal Power Plants, Eng. Appl. Artif. Intell., Vol. 5, No. 6, pp Abdelhameed, M.M., (999, Adaptive Neural Network Based Controller For Robots, International Journal of Mehatronis, Vol. 3, No., pp. -4. G. Geng, G. M. Geary, (997, Appliation of a Neural-Network- Based RLS Algorithm in the Generalized Preditive Control of a Nonlinear Air-Handling Plant, IEEE Trans. On Control Systems Tehnology, Vol. 5, No. 4, pp G. Cembrano, G. Welis, J. Surde, A. Ruggeri, (997, Dynami Control of a Robot Arm Using CMAC Neural Network, Control Eng. Pratie, Vol. 5, No. 4, pp International Journal of Current Engineering and Tehnology, Vol.4, No.5 (Ot 04

Control Theory association of mathematics and engineering

Control Theory association of mathematics and engineering Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology

More information

INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION

INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION Proeedings of ICFDP 8: 8 th International Congress of Fluid Dynamis & Propulsion Deember 14-17, 006, Sharm El-Shiekh, Sinai, Egypt ICFDP8-EG-154 INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE

More information

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.1 Tuning controllers, in the large view Quoting from Section 16.7: Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output

More information

Variation Based Online Travel Time Prediction Using Clustered Neural Networks

Variation Based Online Travel Time Prediction Using Clustered Neural Networks Variation Based Online Travel Time Predition Using lustered Neural Networks Jie Yu, Gang-Len hang, H.W. Ho and Yue Liu Abstrat-This paper proposes a variation-based online travel time predition approah

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are

More information

Subject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (

Subject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % ( 16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free

More information

Complexity of Regularization RBF Networks

Complexity of Regularization RBF Networks Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw

More information

Adaptive neuro-fuzzy inference system-based controllers for smart material actuator modelling

Adaptive neuro-fuzzy inference system-based controllers for smart material actuator modelling Adaptive neuro-fuzzy inferene system-based ontrollers for smart material atuator modelling T L Grigorie and R M Botez Éole de Tehnologie Supérieure, Montréal, Quebe, Canada The manusript was reeived on

More information

A simple expression for radial distribution functions of pure fluids and mixtures

A simple expression for radial distribution functions of pure fluids and mixtures A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.

More information

MultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates

MultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates Exerpt from the Proeedings of the COMSOL Conferene 9 Boston MultiPhysis Analysis of Trapped Field in Multi-Layer YBCO Plates Philippe. Masson Advaned Magnet Lab *7 Main Street, Bldg. #4, Palm Bay, Fl-95,

More information

The Nonlinear Accuracy Model of Electro-Hydrostatic Actuator

The Nonlinear Accuracy Model of Electro-Hydrostatic Actuator The Nonlinear Auray Model of Eletro-Hydrostati Atuator Rongjie Kang, Zongxia Jiao, Shuai Wu Yaoxing Shang Shool of Automation Siene and Eletrial Engineering Beihang University Beijing, China E-mail: angrongjie@vip.163.om

More information

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo

More information

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This

More information

A Queueing Model for Call Blending in Call Centers

A Queueing Model for Call Blending in Call Centers A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl

More information

Model-based mixture discriminant analysis an experimental study

Model-based mixture discriminant analysis an experimental study Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,

More information

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple

More information

Controller Design Based on Transient Response Criteria. Chapter 12 1

Controller Design Based on Transient Response Criteria. Chapter 12 1 Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors

More information

Simulation and Development of Trans-critical CO2 Rolling Piston Compressor

Simulation and Development of Trans-critical CO2 Rolling Piston Compressor Purdue University Purdue e-pubs International Compressor Engineering Conferene Shool of Mehanial Engineering 010 Simulation and Development of Trans-ritial CO Rolling Piston Compressor Yunfeng Chang Xi'an

More information

UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.

UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker. UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable

More information

3 Tidal systems modelling: ASMITA model

3 Tidal systems modelling: ASMITA model 3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

Assessing the Performance of a BCI: A Task-Oriented Approach

Assessing the Performance of a BCI: A Task-Oriented Approach Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,

More information

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012 INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability

More information

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3

More information

QCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines

QCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany

More information

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia

More information

Gluing Potential Energy Surfaces with Rare Event Simulations

Gluing Potential Energy Surfaces with Rare Event Simulations This is an open aess artile published under an ACS AuthorChoie Liense, whih permits opying and redistribution of the artile or any adaptations for non-ommerial purposes. pubs.as.org/jctc Gluing Potential

More information

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory

More information

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer Proeedings of the 4th IIAE International Conferene on Industrial Appliation Engineering 216 Stabilization of the Preision Positioning Stage Working in the Vauum Environment by Using the Disturbane Observer

More information

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential

More information

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION VIII International Conferene on Frature Mehanis of Conrete and Conrete Strutures FraMCoS-8 J.G.M. Van Mier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang Eds) MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP

More information

BINARY RANKINE CYCLE OPTIMIZATION Golub, M., Koscak-Kolin, S., Kurevija, T.

BINARY RANKINE CYCLE OPTIMIZATION Golub, M., Koscak-Kolin, S., Kurevija, T. BINARY RANKINE CYCLE OPTIMIZATION Golub, M., Kosak-Kolin, S., Kurevija, T. Faulty of Mining, Geology and Petroleum Engineering Department of Petroleum Engineering Pierottijeva 6, Zagreb 0 000, Croatia

More information

A model for measurement of the states in a coupled-dot qubit

A model for measurement of the states in a coupled-dot qubit A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:

More information

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2 Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque

More information

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that

More information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

STUDY OF INHERENT FREQUENCY OF HELMHOLTZ RESONATOR

STUDY OF INHERENT FREQUENCY OF HELMHOLTZ RESONATOR 005 WJTA Amerian Waterjet Conferene August -3, 005! Houston, Texas Paper 6B-4 STUDY OF INHERENT FREQUENCY OF HELMHOLT RESONATOR Gong Weili An Liqian Cui Longlian Xie Guixin Shool of Mehanis, Arhiteture

More information

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)

More information

Analysis of discretization in the direct simulation Monte Carlo

Analysis of discretization in the direct simulation Monte Carlo PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,

More information

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013

Wavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013 Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it

More information

Experiment 03: Work and Energy

Experiment 03: Work and Energy MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into

More information

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION 4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and

More information

13.Prandtl-Meyer Expansion Flow

13.Prandtl-Meyer Expansion Flow 3.Prandtl-eyer Expansion Flow This hapter will treat flow over a expansive orner, i.e., one that turns the flow outward. But before we onsider expansion flow, we will return to onsider the details of the

More information

Planning with Uncertainty in Position: an Optimal Planner

Planning with Uncertainty in Position: an Optimal Planner Planning with Unertainty in Position: an Optimal Planner Juan Pablo Gonzalez Anthony (Tony) Stentz CMU-RI -TR-04-63 The Robotis Institute Carnegie Mellon University Pittsburgh, Pennsylvania 15213 Otober

More information

The Effectiveness of the Linear Hull Effect

The Effectiveness of the Linear Hull Effect The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports

More information

arxiv: v1 [physics.gen-ph] 5 Jan 2018

arxiv: v1 [physics.gen-ph] 5 Jan 2018 The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the

More information

UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR

UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR Australia 006 Paper 5B UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR YAKOV S. KHODORKOVSKY & MICHAEL R. BELTRAN Beltran, In., U.S.A. ABSTRACT

More information

The Laws of Acceleration

The Laws of Acceleration The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the

More information

Active Magnetic Bearings for Frictionless Rotating Machineries

Active Magnetic Bearings for Frictionless Rotating Machineries Ative Magneti Bearings for Fritionless Rotating Mahineries Joga Dharma Setiawan Abstrat Ative magneti bearing (AMB systems an support a rotor without physial ontat and enable users to preisely ontrol rotor

More information

Nuclear Shell Structure Evolution Theory

Nuclear Shell Structure Evolution Theory Nulear Shell Struture Evolution Theory Zhengda Wang (1) Xiaobin Wang () Xiaodong Zhang () Xiaohun Wang () (1) Institute of Modern physis Chinese Aademy of SienesLan Zhou P. R. China 70000 () Seagate Tehnology

More information

Simplified Buckling Analysis of Skeletal Structures

Simplified Buckling Analysis of Skeletal Structures Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation

More information

SOME FUNDAMENTAL ASPECTS OF COMPRESSIBLE FLOW

SOME FUNDAMENTAL ASPECTS OF COMPRESSIBLE FLOW SOE FUNDAENAL ASECS OF CORESSIBLE FLOW ah number gas veloity mah number, speed of sound a a R < : subsoni : transoni > : supersoni >> : hypersoni art three : ah Number 7 Isentropi flow in a streamtube

More information

AC : A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN

AC : A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN AC 28-1986: A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN Minh Cao, Wihita State University Minh Cao ompleted his Bahelor s of Siene degree at Wihita State

More information

Sensitivity Analysis in Markov Networks

Sensitivity Analysis in Markov Networks Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores

More information

Theoretical and Experimental Research on the Yoke of the Axial Piston Pumps

Theoretical and Experimental Research on the Yoke of the Axial Piston Pumps Theoretial and Experimental Researh on the Yoke of the Axial Piston Pumps Vasiliu Niolae, P.E., Ph.D., Professor University "Politehnia" of Buharest, Fluid Power Department Rosu Cristian Adrian, M.E.,

More information

A Spatiotemporal Approach to Passive Sound Source Localization

A Spatiotemporal Approach to Passive Sound Source Localization A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,

More information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS

More information

Indirect Neural Control: A Robustness Analysis Against Perturbations

Indirect Neural Control: A Robustness Analysis Against Perturbations Indiret Neural Control: A Robustne Analysis Against Perturbations ANDRÉ LAURINDO MAITELLI OSCAR GABRIEL FILHO Federal Uversity of Rio Grande do Norte Potiguar Uversity Natal/RN, BRAZIL, ZIP Code 5907-970

More information

Time Domain Method of Moments

Time Domain Method of Moments Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations

More information

Part G-4: Sample Exams

Part G-4: Sample Exams Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

Bending resistance of high performance concrete elements

Bending resistance of high performance concrete elements High Performane Strutures and Materials IV 89 Bending resistane of high performane onrete elements D. Mestrovi 1 & L. Miulini 1 Faulty of Civil Engineering, University of Zagreb, Croatia Faulty of Civil

More information

Machining. Introduction

Machining. Introduction Mahining Introdution Mahining aims to generate the shape of the workpiee from a solid body, or to improve the toleranes and surfae finish of a previously formed workpiee, by removing exess material in

More information

The experimental plan of displacement- and frequency-noise free laser interferometer

The experimental plan of displacement- and frequency-noise free laser interferometer 7th Edoardo Amaldi Conferene on Gravitational Waves (Amaldi7) Journal of Physis: Conferene Series 122 (2008) 012022 The experimental plan of displaement- and frequeny-noise free laser interferometer K

More information

Neuro-Fuzzy Modeling of Heat Recovery Steam Generator

Neuro-Fuzzy Modeling of Heat Recovery Steam Generator International Journal of Mahine Learning and Computing, Vol. 2, No. 5, Otober 202 Neuro-Fuzzy Modeling of Heat Reovery Steam Generator A. Ghaffari, A. Chaibakhsh, and S. Shahhoseini represented in a network

More information

Dynamic Simulation and Composition Control in A 10 L Mixing Tank

Dynamic Simulation and Composition Control in A 10 L Mixing Tank roeeding of International Conferene on Chemial and Material Engineering ISBN : 978--97-8-7 SE.4 - Dynami Simulation and Composition Control in A L Mixing Tank Yulius Deddy Hermawan Chemial Engineering

More information

FREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia

FREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia FREQUENCY DOMAIN FEEDFORWARD COMPENSATION F.J. Pérez Castelo and R. Ferreiro Garia Dept. Ingeniería Industrial. Universidad de La Coruña javierp@ud.es, Phone: 98 7.Fax: -98-7 ferreiro@ud.es, Phone: 98

More information

Discrete Bessel functions and partial difference equations

Discrete Bessel functions and partial difference equations Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat

More information

Heat exchangers: Heat exchanger types:

Heat exchangers: Heat exchanger types: Heat exhangers: he proess of heat exhange between two fluids that are at different temperatures and separated by a solid wall ours in many engineering appliations. he devie used to implement this exhange

More information

FINITE WORD LENGTH EFFECTS IN DSP

FINITE WORD LENGTH EFFECTS IN DSP FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

Supplementary Materials

Supplementary Materials Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown

More information

Study on the leak test technology of spacecraft using ultrasonic

Study on the leak test technology of spacecraft using ultrasonic SINCE2013 Singapore International NDT Conferene & Exhibition 2013, 19-20 July 2013 Study on the test tehnology of spaeraft using ultrasoni Yan Rongxin, Li Weidan Beijing Institute of Spaeraft Environment

More information

A Novel Process for the Study of Breakage Energy versus Particle Size

A Novel Process for the Study of Breakage Energy versus Particle Size Geomaterials, 2013, 3, 102-110 http://dx.doi.org/10.4236/gm.2013.33013 Published Online July 2013 (http://www.sirp.org/journal/gm) A Novel Proess for the Study of Breakage Energy versus Partile Size Elias

More information

IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION

IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION 09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE

More information

Multicomponent analysis on polluted waters by means of an electronic tongue

Multicomponent analysis on polluted waters by means of an electronic tongue Sensors and Atuators B 44 (1997) 423 428 Multiomponent analysis on polluted waters by means of an eletroni tongue C. Di Natale a, *, A. Maagnano a, F. Davide a, A. D Amio a, A. Legin b, Y. Vlasov b, A.

More information

Directional Coupler. 4-port Network

Directional Coupler. 4-port Network Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding

More information

An iterative least-square method suitable for solving large sparse matrices

An iterative least-square method suitable for solving large sparse matrices An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method

More information

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept

More information

EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS

EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS -Tehnial Note- EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS H.-M. Chu * W.-L. Li ** Department of Mehanial Engineering Yung-Ta Institute of Tehnology & Commere Ping-Tung,

More information

Determination of the Aerodynamic Characteristics of Flying Vehicles Using Method Large Eddy Simulation with Software ANSYS

Determination of the Aerodynamic Characteristics of Flying Vehicles Using Method Large Eddy Simulation with Software ANSYS Automation, Control and Intelligent Systems 15; 3(6): 118-13 Published online Deember, 15 (http://www.sienepublishinggroup.om//ais) doi: 1.11648/.ais.1536.14 ISSN: 38-5583 (Print); ISSN: 38-5591 (Online)

More information

Developing Excel Macros for Solving Heat Diffusion Problems

Developing Excel Macros for Solving Heat Diffusion Problems Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper

More information

Sensor management for PRF selection in the track-before-detect context

Sensor management for PRF selection in the track-before-detect context Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,

More information

Maximum Likelihood Multipath Estimation in Comparison with Conventional Delay Lock Loops

Maximum Likelihood Multipath Estimation in Comparison with Conventional Delay Lock Loops Maximum Likelihood Multipath Estimation in Comparison with Conventional Delay Lok Loops Mihael Lentmaier and Bernhard Krah, German Aerospae Center (DLR) BIOGRAPY Mihael Lentmaier reeived the Dipl.-Ing.

More information

Relativity in Classical Physics

Relativity in Classical Physics Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of

More information

7 Max-Flow Problems. Business Computing and Operations Research 608

7 Max-Flow Problems. Business Computing and Operations Research 608 7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates

More information

ONLINE APPENDICES for Cost-Effective Quality Assurance in Crowd Labeling

ONLINE APPENDICES for Cost-Effective Quality Assurance in Crowd Labeling ONLINE APPENDICES for Cost-Effetive Quality Assurane in Crowd Labeling Jing Wang Shool of Business and Management Hong Kong University of Siene and Tehnology Clear Water Bay Kowloon Hong Kong jwang@usthk

More information

Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene

Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene Exerpt from the Proeedings of the OMSOL onferene 010 Paris Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene M. Bekmann-Kluge 1 *,. errero 1, V. Shröder 1, A. Aikalin and J. Steinbah

More information

Singular Event Detection

Singular Event Detection Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate

More information

Neuro-Fuzzy Control of Chemical Reactor with Disturbances

Neuro-Fuzzy Control of Chemical Reactor with Disturbances Neuro-Fuzzy Control of Chemial Reator with Disturbanes LENK BLHOÁ, JÁN DORN Department of Information Engineering and Proess Control, Institute of Information Engineering, utomation and Mathematis Faulty

More information

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures

Development of a user element in ABAQUS for modelling of cohesive laws in composite structures Downloaded from orbit.dtu.dk on: Jan 19, 2019 Development of a user element in ABAQUS for modelling of ohesive laws in omposite strutures Feih, Stefanie Publiation date: 2006 Doument Version Publisher's

More information

The simulation analysis of the bridge rectifier continuous operation in AC circuit

The simulation analysis of the bridge rectifier continuous operation in AC circuit Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of

More information

Seismic dip estimation based on the two-dimensional Hilbert transform and its application in random noise attenuation a

Seismic dip estimation based on the two-dimensional Hilbert transform and its application in random noise attenuation a Seismi dip estimation based on the two-dimensional Hilbert transform and its appliation in random noise attenuation a a Published in Applied Geophysis, 1, 55-63 (Marh 015) Cai Liu, Changle Chen, Dian Wang,

More information

General Equilibrium. What happens to cause a reaction to come to equilibrium?

General Equilibrium. What happens to cause a reaction to come to equilibrium? General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember

More information