Inverse dynamics based predictive control for unified power flow controllers without DC bus

This paper presents a new predictive digital control method applied to Matrix Converters (MC) operating as Unified Power Flow Controllers (UPFC). This control method, based on the inverse dynamics model equations of the MC operating as UPFC, just needs to compute the optimal control vector once in each control cycle, in contrast to direct dynamics predictive methods that needs 27 vector calculations. The theoretical principles of the inverse dynamics power flow predictive control of the MC based UPFC with input filter are established. The proposed inverse dynamics predictive power control method is tested using Matlab/Simulink Power Systems toolbox and the obtained results show that the designed power controllers guarantees decoupled active and reactive power control, zero error tracking, fast response times and an overall good dynamic and steady-state response.

A simple model of aggregate pension expenditure

This paper develops a simple model that can be used to analyze the long-term sustainability of the contributive pension system and the steady-state response of pension expenditure to changes in some key demographic and economic variables, in the characteristics of the average pensioner and in the parameters that describe how pensions are calculated in Spain as a function of workers' Social Security contribution histories.

Sliding Mode Control for Switched Mode Power Supplies

Control applications of switched mode power supplies have been widely investigated. The main objective ofresearch and development (R&D) in this field is always to find the most suitable control method to be implemented in various DC/DC converter topologies. Inother words, the goal is to select a control method capable of improving the efficiency of the converter, reducing the effect of disturbances (line and load variation), lessening the effect of EMI (electro magnetic interference), and beingless effected by component variation. The main objective of this research work is to study different control methods implemented in switched mode power supplies namely (PID control, hysteresis control, adaptive control, current programmed control, variable structure control (VSC), and sliding mode control (SMC). The advantages and drawbacks of each control method are given. Two control methods, the PID and the SMC are selected and their effects on DC/DC (Buck, Boost, and Buck-Boost) converters are examined. Matlab/SimulinkTM is used to implement PID control method in DC/DC Buck converter and SMC in DC/DC (Buck, and Buck Boost) converters. For the prototype, operational amplifiers (op-amps) are used to implement PID control in DC/DC Buck converter. For SMC op-amps are implemented in DC/DC Buck converter and dSPACETM is used to control the DC/DC Buck-Boost converter. The SMC can be applied to the DC/DC (Buck, Boost, and Buck-Boost) converters. A comparison of the effects of the PID control and the SMC on the DC/DC Buck converter response in steady state, under line variations, load variations, and different component variations is performed. Also the Conducted RF-Emissions between the PID and SMC DC/DC Buck Converter are compared. The thesis shows that, in comparison with the PID control, the SMC provides better steady-state response, better dynamic response, less EMI, inherent order reduction, robustness against system uncertainty disturbances, and an implicit stability proof. Giving a better steady-state and dynamic response, the SMC is implemented in a DC/DC resonant converter. The half-wave zero current switching (HWZCS) DC/DC Buck converter is selected as a converter topology. A general guideline to select the tank component values, needed for the designing of a HWZCS DC/DC Buck, is obtained. The implementation of the SMC to a HWZCS DC/DC Buck converter is analysed. The converter response is investigated in the steady-state region and in the dynamic region.

Päällystyskoneiden runkojen värähtelyiden FE-analysoinnin kehittäminen

Työssä kehitettiin päällystyskoneen runkojen FE-analysointia erityisesti vastelaskennan kannalta, joskin myös ominaistaajuuslaskennan tarkentamiseksi esitettiin parannuksia. Työssä rajoitutaan harmonisten telaherätteiden aiheuttamien vasteiden analysointiin. . Työssä käsitellään värähtelyjen teoriaa ja tarkastellaan Abaqus-ohjelmiston (versio 5.8) laskentamenetelmiä vakiotilan värähtelyvasteen laskemiseksi. FEM-mallin rakenteeseen liittyen käsitellään perustuksien ja maaperän, konepalkin ja telaherätteiden mallintamista ja mallinnuksen laajuutta. Telaherätteitä käsitellään yksittäisen telan ja telojen samanvaiheisuuden kannalta. Samanvaiheisuutta tutkitaan työssä kehitetyllä summafunktiolla. Värähtelyjen mittaukseen FEM-mallien verifioimiseksi esitetään parannuksia. Nykyinen mallinnustapa käsitellään lyhyesti. Parannusehdotuksia kokeiltiin mallintamalla Rauma 400-päällystyskone ja vertaamalla tuloksia mitattuihin. Laskennan tulokset vastasivat vaihtelevasti mittaustuloksia, mittaustulosten puutteellisuus vaikeutti vertailua. Tulosten perusteella herätetiedon parantaminen on perusteltua ja mallin laajennus lisää todenmukaisuutta. Maaperän huomiointi vaikuttaa ennen kaikkea ominaistaajuuksiin ja muotoihin ja suora ratkaisutapa on käyttökelpoinen vasteen laskentamenetelmä otettaessa maaperä huomioon.

Chaos and Unpredictability in the Vibration of an Elasto-Plastic Beam

This contribution discusses the nonlinear dynamics of a pin-ended elasto-plastic beam with both kinematic and isotropic hardening. An iterative numerical procedure based on the operator split technique is developed in order to deal with the nonlinearities in the equations of motion. Free and forced responses for harmonic sinusoidal and square wave excitations are investigated. Numerical simulations present many interesting behaviors such as jump phenomena, sensitivity to initial conditions, chaos and transient chaos. These results indicate that there are practical problems in predicting the response of the beam even when periodic steady state response is expected.

Hopping-tunneling model to describe electric charge injection at metal/organic semiconductor heterojunctions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Current collection by an active spherical electrode in an unmagnetized plasma

A theoretical model for the steady-state response of anodic contactors that emit a plasma current Ii and collect electrons from a collisionless, unmagnetized plasma is presented. The use of a (kinetic) monoenergetic population for the attracted species, well known in passive probe theory, gives both accuracy and tractability to the theory. The monoenergetic population is proved to behave like an isentropic fluid with radial plus centripetal motion, allowing direct comparisons with ad hoc fluid models. Also, a modification of the original monoenergetic equations permits analysis of contactors operating in orbit-limited conditions. Besides that, the theory predicts that, only for plasma emissions above certain threshold current a presheath/double layer/core structure for the potential is formed (the core mode), while for emissions below that threshold, a plasma contactor behaves exactly as a positive-ion emitter with a presheath/sheath structure (the no-core mode). Ion emitters are studied as a particular case. Emphasis is placed on obtaining dimensionless charts and approximate asymptotic laws of the current-voltage characteristic.

Método computacional para la aceleración de análisis paramétricos de modificaciones locales en estructuras complejas sometidas a cargas dinámicas

Dentro del análisis y diseño estructural surgen frecuentemente problemas de ingeniería donde se requiere el análisis dinámico de grandes modelos de elementos finitos que llegan a millones de grados de libertad y emplean volúmenes de datos de gran tamaño. La complejidad y dimensión de los análisis se dispara cuando se requiere realizar análisis paramétricos. Este problema se ha abordado tradicionalmente desde diversas perspectivas: en primer lugar, aumentando la capacidad tanto de cálculo como de memoria de los sistemas informáticos empleados en los análisis. En segundo lugar, se pueden simplificar los análisis paramétricos reduciendo su número o detalle y por último se puede recurrir a métodos complementarios a los elementos .nitos para la reducción de sus variables y la simplificación de su ejecución manteniendo los resultados obtenidos próximos al comportamiento real de la estructura. Se propone el empleo de un método de reducción que encaja en la tercera de las opciones y consiste en un análisis simplificado que proporciona una solución para la respuesta dinámica de una estructura en el subespacio modal complejo empleando un volumen de datos muy reducido. De este modo se pueden realizar análisis paramétricos variando múltiples parámetros, para obtener una solución muy aproximada al objetivo buscado. Se propone no solo la variación de propiedades locales de masa, rigidez y amortiguamiento sino la adición de grados de libertad a la estructura original para el cálculo de la respuesta tanto permanente como transitoria. Adicionalmente, su facilidad de implementación permite un control exhaustivo sobre las variables del problema y la implementación de mejoras como diferentes formas de obtención de los autovalores o la eliminación de las limitaciones de amortiguamiento en la estructura original. El objetivo del método se puede considerar similar a los que se obtienen al aplicar el método de Guyan u otras técnicas de reducción de modelos empleados en dinámica estructural. Sin embargo, aunque el método permite ser empleado en conjunción con otros para obtener las ventajas de ambos, el presente procedimiento no realiza la condensación del sistema de ecuaciones, sino que emplea la información del sistema de ecuaciones completa estudiando tan solo la respuesta en las variables apropiadas de los puntos de interés para el analista. Dicho interés puede surgir de la necesidad de obtener la respuesta de las grandes estructuras en unos puntos determinados o de la necesidad de modificar la estructura en zonas determinadas para cambiar su comportamiento (respuesta en aceleraciones, velocidades o desplazamientos) ante cargas dinámicas. Por lo tanto, el procedimiento está particularmente indicado para la selección del valor óptimo de varios parámetros en grandes estructuras (del orden de cientos de miles de modos) como pueden ser la localización de elementos introducidos, rigideces, masas o valores de amortiguamientos viscosos en estudios previos en los que diversas soluciones son planteadas y optimizadas, y que en el caso de grandes estructuras, pueden conllevar un número de simulaciones extremadamente elevado para alcanzar la solución óptima. Tras plantear las herramientas necesarias y desarrollar el procedimiento, se propone un caso de estudio para su aplicación al modelo de elementos .nitos del UAV MILANO desarrollado por el Instituto Nacional de Técnica Aeroespacial. A dicha estructura se le imponen ciertos requisitos al incorporar un equipo en aceleraciones en punta de ala izquierda y desplazamientos en punta de ala derecha en presencia de la sustentación producida por una ráfaga continua de viento de forma sinusoidal. La modificación propuesta consiste en la adición de un equipo en la punta de ala izquierda, bien mediante un anclaje rígido, bien unido mediante un sistema de reducción de la respuesta dinámica con propiedades de masa, rigidez y amortiguamiento variables. El estudio de los resultados obtenidos permite determinar la optimización de los parámetros del sistema de atenuación por medio de múltiples análisis dinámicos de forma que se cumplan de la mejor forma posible los requisitos impuestos con la modificación. Se comparan los resultados con los obtenidos mediante el uso de un programa comercial de análisis por el método de los elementos .nitos lográndose soluciones muy aproximadas entre el modelo completo y el reducido. La influencia de diversos factores como son el amortiguamiento modal de la estructura original, el número de modos retenidos en la truncatura o la precisión proporcionada por el barrido en frecuencia se analiza en detalle para, por último, señalar la eficiencia en términos de tiempo y volumen de datos de computación que ofrece el método propuesto en comparación con otras aproximaciones. Por lo tanto, puede concluirse que el método propuesto se considera una opción útil y eficiente para el análisis paramétrico de modificaciones locales en grandes estructuras. ABSTRACT When developing structural design and analysis some projects require dynamic analysis of large finite element models with millions of degrees of freedom which use large size data .les. The analysis complexity and size grow if a parametric analysis is required. This problem has been approached traditionally in several ways: one way is increasing the power and the storage capacity of computer systems involved in the analysis. Other obvious way is reducing the total amount of analyses and their details. Finally, complementary methods to finite element analysis can also be employed in order to limit the number of variables and to reduce the execution time keeping the results as close as possible to the actual behaviour of the structure. Following this third option, we propose a model reduction method that is based in a simplified analysis that supplies a solution for the dynamic response of the structure in the complex modal space using few data. Thereby, parametric analysis can be done varying multiple parameters so as to obtain a solution which complies with the desired objetive. We propose not only mass, stiffness and damping variations, but also addition of degrees of freedom to the original structure in order to calculate the transient and steady-state response. Additionally, the simple implementation of the procedure allows an in-depth control of the problem variables. Furthermore, improvements such as different ways to obtain eigenvectors or to remove damping limitations of the original structure are also possible. The purpose of the procedure is similar to that of using the Guyan or similar model order reduction techniques. However, in our method we do not perform a true model order reduction in the traditional sense. Furthermore, additional gains, which we do not explore herein, can be obtained through the combination of this method with traditional model-order reduction procedures. In our procedure we use the information of the whole system of equations is used but only those nodes of interest to the analyst are processed. That interest comes from the need to obtain the response of the structure at specific locations or from the need to modify the structure at some suitable positions in order to change its behaviour (acceleration, velocity or displacement response) under dynamic loads. Therefore, the procedure is particularly suitable for parametric optimization in large structures with >100000 normal modes such as position of new elements, stiffness, mass and viscous dampings in previous studies where different solutions are devised and optimized, and in the case of large structures, can carry an extremely high number of simulations to get the optimum solution. After the introduction of the required tools and the development of the procedure, a study case is proposed with use the finite element model (FEM) of the MILANO UAV developed by Instituto Nacional de Técnica Aeroespacial. Due to an equipment addition, certain acceleration and displacement requirements on left wing tip and right wing tip, respectively, are imposed. The structure is under a continuous sinusoidal wind gust which produces lift. The proposed modification consists of the addition of an equipment in left wing tip clamped through a rigid attachment or through a dynamic response reduction system with variable properties of mass, stiffness and damping. The analysis of the obtained results allows us to determine the optimized parametric by means of multiple dynamic analyses in a way such that the imposed requirements have been accomplished in the best possible way. The results achieved are compared with results from a commercial finite element analysis software, showing a good correlation. Influence of several factors such as the modal damping of the original structure, the number of modes kept in the modal truncation or the precission given by the frequency sweep is analyzed. Finally, the efficiency of the proposed method is addressed in tems of computational time and data size compared with other approaches. From the analyses performed, we can conclude that the proposed method is a useful and efficient option to perform parametric analysis of possible local modifications in large structures.

The relationship between CA/C ratio and individual differences in dynamic accommodative responses while viewing stereoscopic images

The oculomotor synergy as expressed by the CA/C and AC/A ratios was investigated to examine its influence on our previous observation that whereas convergence responses to stereoscopic images are generally stable, some individuals exhibit significant accommodative overshoot. Using a modified video refraction unit while viewing a stereoscopic LCD, accommodative and convergence responses to balanced and unbalanced vergence and focal stimuli (BVFS and UBVFS) were measured. Accommodative overshoot of at least 0.3 D was found in 3 out of 8 subjects for UBVFS. The accommodative response differential (RD) was taken to be the difference between the initial response and the subsequent mean static steady-state response. Without overshoot, RD was quantified by finding the initial response component. A mean RD of 0.11 +/- 0.27 D was found for the 1.0 D step UBVFS condition. The mean RD for the BVFS was 0.00 +/- 0.17 D. There was a significant positive correlation between CA/C ratio and RD (r = +0.75, n = 8, p <0.05) for only UBVFS. We propose that inter-subject variation in RD is influenced by the CA/C ratio as follows: an initial convergence response, induced by disparity of the image, generates convergence-driven accommodation commensurate with the CA/C ratio; the associated transient defocus subsequently decays to a balanced position between defocus-induced and convergence-induced accommodations.

Fopid Controller Design for Robust Performance Using Particle Swarm Optimization

Mathematics Subject Classification: 26A33; 93C15, 93C55, 93B36, 93B35, 93B51; 03B42; 70Q05; 49N05

Nuevo método de identificación usando la respuesta escalón para sistemas de fase no mínima de segundo orden

Este artículo presenta un nuevo método de identificación para sistemas de fase no mínima basado en la respuesta escalón. El enfoque propuesto provee un modelo aproximado de segundo orden evitando diseños experimentales complejos. El método propuesto es un algoritmo de identificación cerrado basado en puntos característicos de la respuesta escalón de sistemas de fase no mínima de segundo orden. Él es validado usando diferentes modelos lineales. Ellos tienen respuesta inversa entre 3,5% y 38% de la respuesta en régimen permanente. En simulaciones, ha sido demostrado que resultados satisfactorios pueden ser obtenidos usando el procedimiento de identificación propuesto, donde los parámetros identificados presentan errores relativos medios, menores que los obtenidos mediante el método de Balaguer.

Blind separation of auditory event-related brain responses into independent components

Averaged event-related potential (ERP) data recorded from the human scalp reveal electroencephalographic (EEG) activity that is reliably time-locked and phase-locked to experimental events. We report here the application of a method based on information theory that decomposes one or more ERPs recorded at multiple scalp sensors into a sum of components with fixed scalp distributions and sparsely activated, maximally independent time courses. Independent component analysis (ICA) decomposes ERP data into a number of components equal to the number of sensors. The derived components have distinct but not necessarily orthogonal scalp projections. Unlike dipole-fitting methods, the algorithm does not model the locations of their generators in the head. Unlike methods that remove second-order correlations, such as principal component analysis (PCA), ICA also minimizes higher-order dependencies. Applied to detected—and undetected—target ERPs from an auditory vigilance experiment, the algorithm derived ten components that decomposed each of the major response peaks into one or more ICA components with relatively simple scalp distributions. Three of these components were active only when the subject detected the targets, three other components only when the target went undetected, and one in both cases. Three additional components accounted for the steady-state brain response to a 39-Hz background click train. Major features of the decomposition proved robust across sessions and changes in sensor number and placement. This method of ERP analysis can be used to compare responses from multiple stimuli, task conditions, and subject states.

Differential regulation of PGC-1 alpha expression in rat liver and skeletal muscle in response to voluntary running

Background: The beneficial actions of exercise training on lipid, glucose and energy metabolism and insulin sensitivity appear to be in part mediated by PGC-1 alpha. Previous studies have shown that spontaneously exercised rats show at rest enhanced responsiveness to exogenous insulin, lower plasma insulin levels and increased skeletal muscle insulin sensitivity. This study was initiated to examine the functional interaction between exercise-induced modulation of skeletal muscle and liver PGC-1 alpha protein expression, whole body insulin sensitivity, and circulating FFA levels as a measure of whole body fatty acid (lipid) metabolism. Methods: Two groups of male Wistar rats (2 Mo of age, 188.82 +/- 2.77 g BW) were used in this study. One group consisted of control rats placed in standard laboratory cages. Exercising rats were housed individually in cages equipped with running wheels and allowed to run at their own pace for 5 weeks. At the end of exercise training, insulin sensitivity was evaluated by comparing steady-state plasma glucose (SSPG) concentrations at constant plasma insulin levels attained during the continuous infusion of glucose and insulin to each experimental group. Subsequently, soleus and plantaris muscle and liver samples were collected and quantified for PGC-1 alpha protein expression by Western blotting. Collected blood samples were analyzed for glucose, insulin and FFA concentrations. Results: Rats housed in the exercise wheel cages demonstrated almost linear increases in running activity with advancing time reaching to maximum value around 4 weeks. On an average, the rats ran a mean (Mean +/- SE) of 4.102 +/- 0.747 km/day and consumed significantly more food as compared to sedentary controls (P < 0.001) in order to meet their increased caloric requirement. Mean plasma insulin (P < 0.001) and FFA (P < 0.006) concentrations were lower in the exercise-trained rats as compared to sedentary controls. Mean steady state plasma insulin (SSPI) and glucose (SSPG) concentrations were not significantly different in sedentary control rats as compared to exercise-trained animals. Plantaris PGC-1 alpha protein expression increased significantly from a 1.11 +/- 0.12 in the sedentary rats to 1.74 +/- 0.09 in exercising rats (P < 0.001). However, exercise had no effect on PGC-1 alpha protein content in either soleus muscle or liver tissue. These results indicate that exercise training selectively up regulates the PGC-1 alpha protein expression in high-oxidative fast skeletal muscle type such as plantaris muscle. Conclusion: These data suggest that PGC-1 alpha most likely plays a restricted role in exercise-mediated improvements in insulin resistance (sensitivity) and lowering of circulating FFA levels.

The first non-turnover voltammetric response from a molybdenum enzyme: direct electroscopy of dimethylsulfoxide reductase from Rhodobacter capsulatus

The first direct voltammetric response from a molybdenum enzyme under non-turnover conditions is reported. Cyclic voltammetry of dimethylsulfoxide reductase from Rhodobacter capsulatus reveals a reversible Mo-VI/V response at + 161 mV followed by a reversible Mo-V/IV response at -102 mV versus NHE at pH 8. The higher potential couple exhibits a pH dependence consistent with protonation upon reduction to the Mo-V state and we have determined the pK(a) for this semi-reduced species to be 9.0. The lower potential couple is pH independent within the range 5 < pH < 10. The optical spectrum of the Mo chromophore has been investigated with spectroelectrochemistry. At high potential, in its resting state, the enzyme exhibits a spectrum characteristic of the Mo-VI form. This changes significantly following bulk electrolysis (-400 mV versus NHE) at an optically transparent, indium-doped tin oxide working electrode, where a single visible electronic maximum at 632 nm is observed, which is comparable with spectra reported previously for the dithionite-reduced enzyme. This two-electron process is chemically reversible by reoxidizing the enzyme at the electrode in the absence of mediators or promoters. The activity of the enzyme has been established by observation of a catalytic current in the presence of DMSO at pH 8, where a sigmoidal (steady state) voltammogram is seen. Electronic supplementary material to this paper (Fig. S 1) can be obtained by using the Springer Link server located at http://dx.doi.org/10.1007/s00775-002-0374-y.

Response characteristics of arsenic-sensitive bioreporters expressing the gfp reporter gene

This paper describes the development of an analytical technique for arsenic analyses that is based on genetically-modified bioreporter bacteria bearing a gene encoding for the production of a green fluorescent protein (gfp). Upon exposure to arsenic (in the aqueous form of arsenite), the bioreporter production of the fluorescent reporter molecule is monitored spectroscopically. We compared the response measured as a function of time and concentration by steady-state fluorimetry (SSF) to that measured by epi-fluorescent microscopy (EFM). SSF is a bulk technique; as such it inherently yields less information, whereas EFM monitors the response of many individual cells simultaneously and data can be processed in terms of population averages or subpopulations. For the bioreporter strain used here, as well as for the literature we cite, the two techniques exhibit similar performance characteristics. The results presented here show that the EFM technique can compete with SSF and shows substantially more promise for future improvement; it is a matter of research interest to develop optimized methods of EFM image analysis and statistical data treatment. EFM is a conduit for understanding the dynamics of individual cell response vs. population response, which is not only a matter of research interest, but is also promising in the practical terms of developing micro-scale analysis.