Semi-analytical study of the rotational motion stability of artificial satellites using quaternions

This study at aims performing the stability analysis of the rotational motion to artificial satellites using quaternions to describe the satellite attitude (orientation on the space). In the system of rotational motion equations, which is composed by four kinematic equations of the quaternions and by the three Euler equations in terms of the rotational spin components. The influence of the gravity gradient and the direct solar radiation pressure torques have been considered. Equilibrium points were obtained through numerical simulations using the softwares Matlab and Octave, which are then analyzed by the Routh-Hurwitz Stability Criterion.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

Spin-Stabilized Spacecrafts: Analytical Attitude Propagation Using Magnetic Torques

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Influência da suspensão na segurança e no conforto de um pulverizador autopropelido

Para avaliar o comportamento da suspensão do pulverizador autopropelido, foram desenvolvidos modelos físicos e matemáticos em função da excitação ocasionada pelas irregularidades do solo. Neste trabalho, estas irregularidades são representadas por obstáculos de uma pista normalizada segundo a norma ISO 5008. As equações do movimento são obtidas a partir dos modelos matemáticos de meio veículo. As simulações numéricas são executadas nos softwares Matlab® e Simulink®. A partir da entrada conhecida, podem-se determinar as características dos elementos da suspensão para obter níveis desejáveis de conforto e segurança. Foram analisadas quatro diferentes configurações do sistema, variando-se a relação de rigidez a partir de um modelo considerado padrão. Constatou-se que o aumento da relação de rigidez resulta na redução da aceleração vertical e no aumento do curso da suspensão, melhorando o conforto e diminuindo a segurança.

Analytical attitude propagation with non-singular variables and gravity gradient torque for spin stabilized satellite

An analytical approach for the spin stabilized satellite attitude propagation is presented using the non-singular canonical variables to describe the rotational motion. Two sets of variables were introduced for Fukushima in 1994 by a canonical transformation and they are useful when the angle between z-satellite axis of a coordinate system fixed in artificial satellite and the rotational angular momentum vector is zero or when the angle between Z-equatorial axis and rotation angular momentum vector is zero. Analytical solutions for rotational motion equations and torque-free motion are discussed in terms of the elliptic functions and by the application of some simplification to get an approximated solution. These solutions are compared with a numerical solution and the results show a good agreement for many rotation periods. When the mean Hamiltonian associated with the gravity gradient torque is included, an analytical solution is obtained by the application of the successive approximations' method for the satellite in an elliptical orbit. These solutions show that the magnitude of the rotation angular moment is not affected by the gravity gradient torque but this torque causes linear and periodic variations in the angular variables, long and short periodic variations in Z-equatorial component of the rotation angular moment and short periodic variations in x-satellite component of the rotation angular moment. The goal of this analysis is to emphasize the geometrical and physical meaning of the non-singular variables and to validate the approximated analytical solution for the rotational motion without elliptic functions for a non-symmetrical satellite. The analysis can be applied for spin stabilized satellite and in this case the general solution and the approximated solution are coincidence. Then the results can be used in analysis of the space mission of the Brazilian Satellites. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Swarm control designs applied to a non-ideal load transportation system

In this paper, a load transport system in platforms is considered. It is a transport device and is modelled as an inverted pendulum built on a car driven by a DC motor. The motion equations were obtained by Lagrange's equations. The mathematical model considers the interaction between the DC motor and the dynamic system. The dynamic system was analysed and a Swarm Control Design was developed to stabilize the model of this load transport system. ©2010 IEEE.

Nonlinear dynamic model and simulation of morphing wing profile actuated by shape memory alloys

Morphing aircraft have the ability to actively adapt and change their shape to achieve different missions efficiently. The development of morphing structures is deeply related with the ability to model precisely different designs in order to evaluate its characteristics. This paper addresses the dynamic modeling of a sectioned wing profile (morphing airfoil) connected by rotational joints (hinges). In this proposal, a pair of shape memory alloy (SMA) wires are connected to subsequent sections providing torque by reducing its length (changing airfoil camber). The dynamic model of the structure is presented for one pair of sections considering the system with one degree of freedom. The motion equations are solved using numerical techniques due the nonlinearities of the model. The numerical results are compared with experimental data and a discussion of how good this approach captures the physical phenomena associated with this problem. © The Society for Experimental Mechanics, Inc. 2012.

Solução analítica da equação de movimento clássica para o potencial de Morse generalizado e potenciais de Morse simétricos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Suspensões de veículos - Um estudo de caso

Pós-graduação em Engenharia Mecânica - FEG

Analytical prediction of the spin stabilized satellite's attitude using the solar radiation torque

The aim of this paper is to present an analytical solution for the spin motion equations of spin-stabilized satellite considering only the influence of solar radiation torque. The theory uses a cylindrical satellite on a circular orbit and considers that the satellite is always illuminated. The average components of this torque were determined over an orbital period. These components are substituted in the spin motion equations in order to get an analytical solution for the right ascension and declination of the satellite spin axis. The time evolution for the pointing deviation of the spin axis was also analyzed. These solutions were numerically implemented and compared with real data of the Brazilian Satellite of Data Collection - SCD1 an SCD2. The results show that the theory has consistency and can be applied to predict the spin motion of spin-stabilized artificial satellites.

Controle de dissociação molecular com ferramentas de dinâmica não linear

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

O Método de Kane aplicado a um modelo de corpo humano para obtenção da pressão intradiscal na coluna lombar

Nos últimos anos tem surgido vários métodos para se estudar sistemas multicorpos. A vantagem destes métodos é a obtenção das equações do movimento de forma mais fácil e precisa. Robôs, manipuladores, mecanismos, satélites, máquinas em geral e sistemas biomcânicos, pelas suas características, podem ser vantajosamente tratados por formalismo multicorpo, mas precisamente, o método de Kane a um modelo do corpo humano, com a finalidade de se determinar a pressão intradiscal entre as vértebras L4/L5 da coluna lombar, em algumas atividades normais de vida diária, como: sentar, abaixar, pegar e elevar cargas, realizar movimentos de flexão e extensão. Pretendeu-se, também obter os esforços internos nas diversas articulações deste modelo. A pressão intradiscal representa um parâmetro importante para caracterizar a sobrecarga na coluna, podendo ter relação com a degeneração do disco intervertebral, o que produz muita dor em indivíduos. Inicialmente, o modelo foi proposto e o método de Kane foi aplicado para se obter as equações do movimento. Para resolvê-las, foi necessário realiza a filmagem de um sujeito realizando as atividades citadas, a fim de se obter as equações das trajetórias dos diversos segmentos do corpo humano. Após a simulação, os parâmetros desejados foram obtidos e comparados com valores experimentais citados na literatura.

Satélites estabilizados por rotação e torque de radiação solar direta

The aims of this work are to analyze the direct solar radiation pressure torque (TPRS) in the rotational motion of spin-stabilized artificial satellites, to numerically implement these solutions and to compare the results with real data of the Brazilian Satellite Data Collection – SCD1 and SCD2, supplied by INPE. The mathematical model for this torque is determined for a cylindrical satellite, and the components of this torque are determined in a fixed system in the satellite. An analytical solution for the spin motion equations is proposed, in which TPRSD does not affect the spin velocity of the satellite. Two approaches are adopted in the numerical implementation of the developed theory: the first one considers the proposed theory and the second introduces a variation in the spin velocity based on its real variation. The results obtained indicate that the solar radiation pressure torque has little influence in the right ascension and declination axis of rotation due to the small dimension of the satellite and altitude in which it is found. To better validate the application of the presented theory, the angular deviation of the spin axis and solar aspect angle were also analyzed. The comparison of the results of the approaches conducted with real data show good precision in the theory, which can be applied in the prediction of the rotational motion of the spin-stabilized artificial satellites, when others external torques are considered besides the direct solar radiation pressure torque

Modelagem matemática e comportamento dinâmico da suspensão passiva de um pulverizador agrícola autopropelido

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos