Estimation of aircraft aerodynamic derivatives using Extended Kalman Filter

Design of flight control laws, verification of performance predictions, and the implementation of flight simulations are tasks that require a mathematical model of the aircraft dynamics. The dynamical models are characterized by coefficients (aerodynamic derivatives) whose values must be determined from flight tests. This work outlines the use of the Extended Kalman Filter (EKF) in obtaining the aerodynamic derivatives of an aircraft. The EKF shows several advantages over the more traditional least-square method (LS). Among these the most important are: there are no restrictions on linearity or in the form which the parameters appears in the mathematical model describing the system, and it is not required that these parameters be time invariant. The EKF uses the statistical properties of the process and the observation noise, to produce estimates based on the mean square error of the estimates themselves. Differently, the LS minimizes a cost function based on the plant output behavior. Results for the estimation of some longitudinal aerodynamic derivatives from simulated data are presented.

Prediction of transients and control reactions in a transonic wind tunnel

This work describes a lumped parameter mathematical model for the prediction of transients in an aerodynamic circuit of a transonic wind tunnel. Control actions to properly handle those perturbations are also assessed. The tunnel circuit technology is up to date and incorporates a novel feature: high-enthalpy air injection to extend the tunnel’s Reynolds number capability. The model solves the equations of continuity, energy and momentum and defines density, internal energy and mass flow as the basic parameters in the aerodynamic study as well as Mach number, stagnation pressure and stagnation temperature, all referred to test section conditions, as the main control variables. The tunnel circuit response to control actions and the stability of the flow are numerically investigated. Initially, for validation purposes, the code was applied to the AWT ("Altitude Wind Tunnel" of NASA-Lewis). In the sequel, the Brazilian transonic wind tunnel was investigated, with all the main control systems modeled, including injection.

A study of the facial aging - a multidisciplinary approach

This paper describes a mathematical and graphical model for face aging. It considers the possibility of predicting the aging process by offering an initial quantification of this process as it applies to the face. It is concerned with physical measurements and a general law of time dependence. After measuring and normalizing a photograph of a person, one could predict, with a known amount of error, the appearance of that person at a different age. The technique described has served its purpose successfully, with a representative amount of patient data behaving sufficiently near the general aging curve of each parameter. That model uses a warping technique to emulate the aging changes on the face of women. Frequently the warping methods are based on the interpolation between images or general mathematical functions to calculate the pixel attributes. The implemented process considers the age features of selected parts of a face such as the face outline and the shape of the lips. These age features were obtained by measuring the facial regions of women that have been photographed throughout their lives. The present work is first concerned with discussing a methodology to define the aging parameters that can be measured, and second with representing the age effects graphically.

On the numerical integration of rigid body nonlinear dynamics in presence of parameters singularities

One of the main complexities in the simulation of the nonlinear dynamics of rigid bodies consists in describing properly the finite rotations that they may undergo. It is well known that, to avoid singularities in the representation of the SO(3) rotation group, at least four parameters must be used. However, it is computationally expensive to use a four-parameters representation since, as only three of the parameters are independent, one needs to introduce constraint equations in the model, leading to differential-algebraic equations instead of ordinary differential ones. Three-parameter representations are numerically more efficient. Therefore, the objective of this paper is to evaluate numerically the influence of the parametrization and its singularities on the simulation of the dynamics of a rigid body. This is done through the analysis of a heavy top with a fixed point, using two three-parameter systems, Euler's angles and rotation vector. Theoretical results were used to guide the numerical simulation and to assure that all possible cases were analyzed. The two parametrizations were compared using several integrators. The results show that Euler's angles lead to faster integration compared to the rotation vector. An Euler's angles singular case, where representation approaches a theoretical singular point, was analyzed in detail. It is shown that on the contrary of what may be expected, 1) the numerical integration is very efficient, even more than for any other case, and 2) in spite of the uncertainty on the Euler's angles themselves, the body motion is well represented.

Gas-solid two-phase flow in the riser of circulating fluidized beds: mathematical modelling and numerical simulation

A mathematical model is developed for gas-solids flows in circulating fluidized beds. An Eulerian formulation is followed based on the two-fluids model approach where both the fluid and the particulate phases are treated as a continuum. The physical modelling is discussed, including the formulation of boundary conditions and the description of the numerical methodology. Results of numerical simulation are presented and discussed. The model is validated through comparison to experiment, and simulation is performed to investigate the effects on the flow hydrodynamics of the solids viscosity.

Identification of flexural stiffness parameters of beams

A three degree of freedom model of the dynamic mass at the middle of a test sample, resembling a Stockbridge neutraliser, is introduced. This model is used to identify the hereby called equivalent complex cross section flexural stiffness (ECFS) of the beam element which is part of the whole test sample. This ECFS, once identified, gives the effective cross section flexural stiffness of the beam as well as its effective damping, measured as the loss factor of an equivalent viscoelastic beam. The beam element of the test sample may be of any complexity, such as a segment of stranded cable of the ACSR type. These data are important parameters for the design of overhead power transmission lines and other cable structures. A cost function is defined and used in the identification of the ECFS. An experiment, designed to measure the dynamic masses of two test samples, is described. Experimental and identified results are presented and discussed.

Identification of Flexural Stiffness Parameters of Beams

A three degree of freedom model of the dynamic mass at the middle of a test sample, resembling a Stockbridge neutraliser, is introduced. This model is used to identify the hereby called equivalent complex cross section flexural stiffness (ECFS) of the beam element which is part of the whole test sample. This ECFS, once identified, gives the effective cross section flexural stiffness of the beam as well as its effective damping, measured as the loss factor of an equivalent viscoelastic beam. The beam element of the test sample may be of any complexity, such as a segment of stranded cable of the ACSR type. These data are important parameters for the design of overhead power transmission lines and other cable structures. A cost function is defined and used in the identification of the ECFS. An experiment, designed to measure the dynamic masses of two test samples, is described. Experimental and identified results are presented and discussed.

A Design Methodology for the Compensation of Positioning Deviation in Gantry Manipulators

This work presents recent results concerning a design methodology used to estimate the positioning deviation for a gantry (Cartesian) manipulator, related mainly to structural elastic deformation of components during operational conditions. The case-study manipulator is classified as gantry type and its basic dimensions are 1,53m x 0,97m x 1,38m. The dimensions used for the calculation of effective workspace due to end-effector path displacement are: 1m x 0,5m x 0,5m. The manipulator is composed by four basic modules defined as module X, module Y, module Z and terminal arm, where is connected the end-effector. Each module controlled axis performs a linear-parabolic positioning movement. The planning path algorithm has the maximum velocity and the total distance as input parameters for a given task. The acceleration and deceleration times are the same. Denavit-Hartemberg parameterization method is used in the manipulator kinematics model. The gantry manipulator can be modeled as four rigid bodies with three degrees-of-freedom in translational movements, connected as an open kinematics chain. Dynamic analysis were performed considering inertial parameters specification such as component mass, inertia and center of gravity position of each module. These parameters are essential for a correct manipulator dynamic modelling, due to multiple possibilities of motion and manipulation of objects with different masses. The dynamic analysis consists of a mathematical modelling of the static and dynamic interactions among the modules. The computation of the structural deformations uses the finite element method (FEM).

Modeling of high pressure pretreatment process for gold leaching

Investigation of high pressure pretreatment process for gold leaching is the objective of the present master's thesis. The gold ores and concentrates which cannot be easily treated by leaching process are called "refractory". These types of ores or concentrates often have high content of sulfur and arsenic that renders the precious metal inaccessible to the leaching agents. Since the refractory ores in gold manufacturing industry take a considerable share, the pressure oxidation method (autoclave method) is considered as one of the possible ways to overcome the related problems. Mathematical modeling is the main approach in this thesis which was used for investigation of high pressure oxidation process. For this task, available information from literature concerning this phenomenon, including chemistry, mass transfer and kinetics, reaction conditions, applied apparatus and application, was collected and studied. The modeling part includes investigation of pyrite oxidation kinetics in order to create a descriptive mathematical model. The following major steps are completed: creation of process model by using the available knowledge; estimation of unknown parameters and determination of goodness of the fit; study of the reliability of the model and its parameters.

Optimization of parameters for fiber laser-MAG hybrid welding in shipbuilding applications

The Arctic region becoming very active area of the industrial developments since it may contain approximately 15-25% of the hydrocarbon and other valuable natural resources which are in great demand nowadays. Harsh operation conditions make the Arctic region difficult to access due to low temperatures which can drop below -50 °C in winter and various additional loads. As a result, newer and modified metallic materials are implemented which can cause certain problems in welding them properly. Steel is still the most widely used material in the Arctic regions due to high mechanical properties, cheapness and manufacturability. Moreover, with recent steel manufacturing development it is possible to make up to 1100 MPa yield strength microalloyed high strength steel which can be operated at temperatures -60 °C possessing reasonable weldability, ductility and suitable impact toughness which is the most crucial property for the Arctic usability. For many years, the arc welding was the most dominant joining method of the metallic materials. Recently, other joining methods are successfully implemented into welding manufacturing due to growing industrial demands and one of them is the laser-arc hybrid welding. The laser-arc hybrid welding successfully combines the advantages and eliminates the disadvantages of the both joining methods therefore produce less distortions, reduce the need of edge preparation, generates narrower heat-affected zone, and increase welding speed or productivity significantly. Moreover, due to easy implementation of the filler wire, accordingly the mechanical properties of the joints can be manipulated in order to produce suitable quality. Moreover, with laser-arc hybrid welding it is possible to achieve matching weld metal compared to the base material even with the low alloying welding wires without excessive softening of the HAZ in the high strength steels. As a result, the laser-arc welding methods can be the most desired and dominating welding technology nowadays, and which is already operating in automotive and shipbuilding industries with a great success. However, in the future it can be extended to offshore, pipe-laying, and heavy equipment industries for arctic environment. CO2 and Nd:YAG laser sources in combination with gas metal arc source have been used widely in the past two decades. Recently, the fiber laser sources offered high power outputs with excellent beam quality, very high electrical efficiency, low maintenance expenses, and higher mobility due to fiber optics. As a result, fiber laser-arc hybrid process offers even more extended advantages and applications. However, the information about fiber or disk laser-arc hybrid welding is very limited. The objectives of the Master’s thesis are concentrated on the study of fiber laser-MAG hybrid welding parameters in order to understand resulting mechanical properties and quality of the welds. In this work only ferrous materials are reviewed. The qualitative methodological approach has been used to achieve the objectives. This study demonstrates that laser-arc hybrid welding is suitable for welding of many types, thicknesses and strength of steels with acceptable mechanical properties along very high productivity. New developments of the fiber laser-arc hybrid process offers extended capabilities over CO2 laser combined with the arc. This work can be used as guideline in hybrid welding technology with comprehensive study the effect of welding parameter on joint quality.

Glyphosate influence on the physiological parameters of Conyza bonariensis biotypes

The objective of this study was to determine changes in gas exchange and inhibition of EPSPs, based on the accumulation of shikimic acid in horseweed biotypes resistant and sensitive to glyphosate. Two experiments were conducted in a factorial model. The first one evaluated horseweed biotypes (one resistant and one susceptible to glyphosate), and herbicide rates (0 and 1,080 g a.e. ha ¹) applied on the weed. In the second experiment, the horseweed biotypes (susceptible and resistant to glyphosate) were evaluated in five periods as following: 0, 3, 7, 10, and 14 days after herbicide application (DAH). The photosynthetic rate, transpiration, carboxylation efficiency, and water efficiency were determined using an infrared gas analyzer (IRGA), and shikimic acid concentration by HPLC. The application of glyphosate damaged the photosynthetic parameters of the susceptible biotype, causing complete inhibition of the photosynthetic rate, transpiration rate, carboxylation efficiency and water use efficiency, starting from the 7 DAH. On the other hand, total inhibition of the photosynthetic parameters was not observed for the resistant biotype. Shikimic acid accumulation occurred in both biotypes after glyphosate application but the susceptible biotype had the highest concentrations, indicating greater sensitivity of the enzyme EPSPs. The accumulation of shikimic acid in the resistant biotype indicates that the mechanism of resistance is not related to the total insensitivity of the enzyme EPSPs to glyphosate and/or that other resistance mechanisms may be involved.

Chlorophyll fluorescence parameters in populations of two legume trees: Stryphnodendron adstringens (Mart.) Coville (Mimosoideae) and Cassia ferruginea (Schrad.) Schrad. ex DC. (Caesalpinoideae)

The aim of this study was to investigate the photosynthetic performance in populations of two legume tree species, Stryphnodendron adstringens (Mimosoideae), typical from Cerrado, and Cassia ferruginea (Caesalpinoideae) from the Atlantic Rain Forest. The photosynthetic traits were assessed by measures of chlorophyll fluorescence in progenies of naturally pollinated plants from three populations of S. adstringens and a population of C. ferruginea. Plants of S. adstringens growing under similar conditions of C. ferruginea plants demanded higher light values for photosynthesis saturation, 600 µmol.m-2.s-1 and 350 µmol.m-2.s-1 respectively, and showed higher intrinsic photosynthetic efficiency of photosystem II, Fv/Fm of 0.814 versus 0.783 in C. ferruginea. The highest values of Fv/Fm observed in S. adstringens can explain the highest electron transport rates (ETR) obtained for this species. No significant differences were found among progenies from different C. ferruginea trees nor among populations of S. adstringens, and only in few cases, variation among progenies within populations were found for S. adstringens plants. The fact that fluorescence parameters distinguished species but not populations or most of progenies may be related to low intraspecific genetic variation of these chlorophyll fluorescence traits or due to lack of expression on genetic differences in plants under no stressful conditions.

Biomass and mineralmass estimates in a "cerrado" ecosystem

This work aimed to develop allometric equations for tree biomass estimation, and to determine the site biomass in different "cerrado" ecosystems. Destructive sampling in a "campo cerrado" (open savanna) was carried out at the Biological Reserve of Moji-Guaçu, State of São Paulo, southeastern Brazil. This "campo cerrado" (open savanna) grows under a tropical climate and on acid, low nutrient soils. Sixty wood plants were cut to ground level and measurements of diameter, height and weight of leaves and stems were taken. We selected the best equations among the most commonly used mathematical relations according to R² values, significance, and standard error. Both diameter (D) and height (H) showed good relationship with plant biomass, but the use of these two parameters together (DH and D²H) provided the best predictor variables. The best equations were linear, but power and exponential equations also showed high R² and significance. The applicability of these equations is discussed and biomass estimates are compared with other types of tropical savannas. Mineralmass was also estimated. "Cerrados" proved to have very important carbon reservoirs due to their great extent. In addition, high land-use change that takes place nowadays in the "cerrado" biome may significantly affect the global carbon cycle.

Warranty Prediction Using Weibull Simulation

The main objective of this master’s thesis is to examine if Weibull analysis is suitable method for warranty forecasting in the Case Company. The Case Company has used Reliasoft’s Weibull++ software, which is basing on the Weibull method, but the Company has noticed that the analysis has not given right results. This study was conducted making Weibull simulations in different profit centers of the Case Company and then comparing actual cost and forecasted cost. Simula-tions were made using different time frames and two methods for determining future deliveries. The first sub objective is to examine, which parameters of simulations will give the best result to each profit center. The second sub objective of this study is to create a simple control model for following forecasted costs and actual realized costs. The third sub objective is to document all Qlikview-parameters of profit centers. This study is a constructive research, and solutions for company’s problems are figured out in this master’s thesis. In the theory parts were introduced quality issues, for example; what is quality, quality costing and cost of poor quality. Quality is one of the major aspects in the Case Company, so understand-ing the link between quality and warranty forecasting is important. Warranty management was also introduced and other different tools for warranty forecasting. The Weibull method and its mathematical properties and reliability engineering were introduced. The main results of this master’s thesis are that the Weibull analysis forecasted too high costs, when calculating provision. Although, some forecasted values of profit centers were lower than actual values, the method works better for planning purposes. One of the reasons is that quality improving or alternatively quality decreasing is not showing in the results of the analysis in the short run. The other reason for too high values is that the products of the Case Company are com-plex and analyses were made in the profit center-level. The Weibull method was developed for standard products, but products of the Case Company consists of many complex components. According to the theory, this method was developed for homogeneous-data. So the most im-portant notification is that the analysis should be made in the product level, not the profit center level, when the data is more homogeneous.