A hybrid finite element formulation for flexible multibody dynamics

This work deals with the transient analysis of flexible multibody systems within a hybrid finite element framework. Hybrid finite elements are based on a two-field variational formulation in which the displacements and stresses are interpolated separately yielding very good coarse mesh accuracy. Most of the literature on flexible multibody systems uses beam-theory-based formulations. In contrast, the use of hybrid finite elements uses continuum-based elements, thus avoiding the problems associated with rotational degrees of freedom. In particular, any given three-dimensional constitutive relations can be directly used within the framework of this formulation. Since the coarse mesh accuracy as compared to a conventional displacement-based formulation is very high, the scheme is cost effective as well. A general formulation is developed for the constrained motion of a given point on a line manifold, using a total Lagrangian method. The multipoint constraint equations are implemented using Lagrange multipliers. Various kinds of joints such as cylindrical, prismatic, and screw joints are implemented within this general framework. Hinge joints such as spherical, universal, and revolute joints are obtained simply by using shared nodes between the bodies. In addition to joints, the formulation and implementation details for a DC motor actuator and for prescribed relative rotation are also presented. Several example problems illustrate the efficacy of the developed formulation.

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

A New Deformation Theory with Strain Gradient Effects

A new phenomenological deformation theory with strain gradient effects is proposed. This theory, which belongs to nonlinear elasticity, fits within the framework of general couple stress theory and involves a single material length scale l. In the present theory three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom u(i). omega(i) has no direct dependence upon ui and is called the micro-rotation, i.e. the material rotation theta(i) plus the particle relative rotation. The strain energy density is assumed to only be a function of the strain tensor and the overall curvature tensor, which results in symmetric Cauchy stresses. Minimum potential principle is developed for the strain gradient deformation theory version. In the limit of vanishing 1, it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in details. Comparisons between the present theory and the theory proposed by Shizawa and Zbib (Shizawa, K., Zbib, H.M., 1999. A thermodynamical theory gradient elastoplasticity with dislocation density Censor: fundamentals. Int. J. Plast. 15, 899) are given. With the same hardening law as Fleck et al. (Fleck, N.A., Muller, G.H., Ashby, M.F., Hutchinson, JW., 1994 Strain gradient plasticity: theory and experiment. Acta Metall. Mater 42, 475), the new strain gradient deformation theory is used to investigate two typical examples, i.e. thin metallic wire torsion and ultra-thin metallic beam bend. The results are compared with those given by Fleck et al, 1994 and Stolken and Evans (Stolken, J.S., Evans, A.G., 1998. A microbend test method for measuring the plasticity length scale. Acta Mater. 46, 5109). In addition, it is explained for a unit cell that the overall curvature tensor produced by the overall rotation vector is the work conjugate of the overall couple stress tensor. (C) 2002 Elsevier Science Ltd. All rights reserved.

Spectrophone measurements in sulfur-hexafluoride

The optoacoustic signal generated by pulsed 10.6 c infrared radiation incident upon a test cell filled with gaseous SF6 has been analyzed in detail. The effects ofm icroscopic energy transfer from the absorbing vibrational degrees of freedom, spontaneous emission, thermal conduction, and acoustic wave propagation are included. This complete treatment explains the experimental observations including a negative pressure response following irradiation at low gas pressure.

A rational fraction polynomials model to study vertical dynamic wheel-rail interaction

This paper presents a model designed to study vertical interactions between wheel and rail when the wheel moves over a rail welding. The model focuses on the spatial domain, and is drawn up in a simple fashion from track receptances. The paper obtains the receptances from a full track model in the frequency domain already developed by the authors, which includes deformation of the rail section and propagation of bending, elongation and torsional waves along an infinite track. Transformation between domains was secured by applying a modified rational fraction polynomials method. This obtains a track model with very few degrees of freedom, and thus with minimum time consumption for integration, with a good match to the original model over a sufficiently broad range of frequencies. Wheel-rail interaction is modelled on a non-linear Hertzian spring, and consideration is given to parametric excitation caused by the wheel moving over a sleeper, since this is a moving wheel model and not a moving irregularity model. The model is used to study the dynamic loads and displacements emerging at the wheel-rail contact passing over a welding defect at different speeds.

Long-lived oscillatory incoherent electron dynamics in molecules: trans-polyacetylene oligomers

We identify an intriguing feature of the electron-vibrational dynamics of molecular systems via a computational examination of trans-polyacetylene oligomers. Here, via the vibronic interactions, the decay of an electron in the conduction band resonantly excites an electron in the valence band, and vice versa, leading to oscillatory exchange of electronic population between two distinct electronic states that lives for up to tens of picoseconds. The oscillatory structure is reminiscent of beating patterns between quantum states and is strongly suggestive of the presence of long-lived molecular electronic coherence. Significantly, however, a detailed analysis of the electronic coherence properties shows that the oscillatory structure arises from a purely incoherent process. These results were obtained by propagating the coupled dynamics of electronic and vibrational degrees of freedom in a mixed quantum-classical study of the Su-Schrieffer-Heeger Hamiltonian for polyacetylene. The incoherent process is shown to occur between degenerate electronic states with distinct electronic configurations that are indirectly coupled via a third auxiliary state by vibronic interactions. A discussion of how to construct electronic superposition states in molecules that are truly robust to decoherence is also presented

From transistor to trapped-ion computers for quantum chemistry

Over the last few decades, quantum chemistry has progressed through the development of computational methods based on modern digital computers. However, these methods can hardly fulfill the exponentially-growing resource requirements when applied to large quantum systems. As pointed out by Feynman, this restriction is intrinsic to all computational models based on classical physics. Recently, the rapid advancement of trapped-ion technologies has opened new possibilities for quantum control and quantum simulations. Here, we present an efficient toolkit that exploits both the internal and motional degrees of freedom of trapped ions for solving problems in quantum chemistry, including molecular electronic structure, molecular dynamics, and vibronic coupling. We focus on applications that go beyond the capacity of classical computers, but may be realizable on state-of-the-art trapped-ion systems. These results allow us to envision a new paradigm of quantum chemistry that shifts from the current transistor to a near-future trapped-ion-based technology.

Trapped ion magnetic resonance: concepts and designs

A novel spectroscopy of trapped ions is proposed which will bring single-ion detection sensitivity to the observation of magnetic resonance spectra. The approaches developed here are aimed at resolving one of the fundamental problems of molecular spectroscopy, the apparent incompatibility in existing techniques between high information content (and therefore good species discrimination) and high sensitivity. Methods for studying both electron spin resonance (ESR) and nuclear magnetic resonance (NMR) are designed. They assume established methods for trapping ions in high magnetic field and observing the trapping frequencies with high resolution (<1 Hz) and sensitivity (single ion) by electrical means. The introduction of a magnetic bottle field gradient couples the spin and spatial motions together and leads to a small spin-dependent force on the ion, which has been exploited by Dehmelt to observe directly the perturbation of the ground-state electron's axial frequency by its spin magnetic moment.

A series of fundamental innovations is described m order to extend magnetic resonance to the higher masses of molecular ions (100 amu = 2x 10^5 electron masses) and smaller magnetic moments (nuclear moments = 10^(-3) of the electron moment). First, it is demonstrated how time-domain trapping frequency observations before and after magnetic resonance can be used to make cooling of the particle to its ground state unnecessary. Second, adiabatic cycling of the magnetic bottle off between detection periods is shown to be practical and to allow high-resolution magnetic resonance to be encoded pointwise as the presence or absence of trapping frequency shifts. Third, methods of inducing spindependent work on the ion orbits with magnetic field gradients and Larmor frequency irradiation are proposed which greatly amplify the attainable shifts in trapping frequency.

The dissertation explores the basic concepts behind ion trapping, adopting a variety of classical, semiclassical, numerical, and quantum mechanical approaches to derive spin-dependent effects, design experimental sequences, and corroborate results from one approach with those from another. The first proposal presented builds on Dehmelt's experiment by combining a "before and after" detection sequence with novel signal processing to reveal ESR spectra. A more powerful technique for ESR is then designed which uses axially synchronized spin transitions to perform spin-dependent work in the presence of a magnetic bottle, which also converts axial amplitude changes into cyclotron frequency shifts. A third use of the magnetic bottle is to selectively trap ions with small initial kinetic energy. A dechirping algorithm corrects for undesired frequency shifts associated with damping by the measurement process.

The most general approach presented is spin-locked internally resonant ion cyclotron excitation, a true continuous Stern-Gerlach effect. A magnetic field gradient modulated at both the Larmor and cyclotron frequencies is devised which leads to cyclotron acceleration proportional to the transverse magnetic moment of a coherent state of the particle and radiation field. A preferred method of using this to observe NMR as an axial frequency shift is described in detail. In the course of this derivation, a new quantum mechanical description of ion cyclotron resonance is presented which is easily combined with spin degrees of freedom to provide a full description of the proposals.

Practical, technical, and experimental issues surrounding the feasibility of the proposals are addressed throughout the dissertation. Numerical ion trajectory simulations and analytical models are used to predict the effectiveness of the new designs as well as their sensitivity and resolution. These checks on the methods proposed provide convincing evidence of their promise in extending the wealth of magnetic resonance information to the study of collisionless ions via single-ion spectroscopy.

Quantum error correction with biased noise

Quantum computing offers powerful new techniques for speeding up the calculation of many classically intractable problems. Quantum algorithms can allow for the efficient simulation of physical systems, with applications to basic research, chemical modeling, and drug discovery; other algorithms have important implications for cryptography and internet security.

At the same time, building a quantum computer is a daunting task, requiring the coherent manipulation of systems with many quantum degrees of freedom while preventing environmental noise from interacting too strongly with the system. Fortunately, we know that, under reasonable assumptions, we can use the techniques of quantum error correction and fault tolerance to achieve an arbitrary reduction in the noise level.

In this thesis, we look at how additional information about the structure of noise, or "noise bias," can improve or alter the performance of techniques in quantum error correction and fault tolerance. In Chapter 2, we explore the possibility of designing certain quantum gates to be extremely robust with respect to errors in their operation. This naturally leads to structured noise where certain gates can be implemented in a protected manner, allowing the user to focus their protection on the noisier unprotected operations.

In Chapter 3, we examine how to tailor error-correcting codes and fault-tolerant quantum circuits in the presence of dephasing biased noise, where dephasing errors are far more common than bit-flip errors. By using an appropriately asymmetric code, we demonstrate the ability to improve the amount of error reduction and decrease the physical resources required for error correction.

In Chapter 4, we analyze a variety of protocols for distilling magic states, which enable universal quantum computation, in the presence of faulty Clifford operations. Here again there is a hierarchy of noise levels, with a fixed error rate for faulty gates, and a second rate for errors in the distilled states which decreases as the states are distilled to better quality. The interplay of of these different rates sets limits on the achievable distillation and how quickly states converge to that limit.

Quantum optomechanics with silicon nanostructures

Mechanical resonators are the most basic and ubiquitous physical systems known. In on-chip form, they are used to process high frequency signals in every cell phone, television, and laptop. They have also been in the last few decades in different shapes and forms, a critical part of progress in quantum information sciences with kilogram-scale mirrors for gravitational wave detection measuring motion at its quantum limits, and the motion of single ions being used to link qubits for quantum computation.

Optomechanics is a field primarily concerned with coupling light to the motion of mechanical structures. This thesis contains descriptions of recent work with mechanical systems in the megahertz to gigahertz frequency range, formed by nanofabricating novel photonic/phononic structures on a silicon chip. These structures are designed to have both optical and mechanical resonances, and laser light is used to address and manipulate their motional degrees of freedom through radiation pressure forces. We laser cool these mechanical resonators to their ground states, and observe for the first time the quantum zero-point motion of a nanomechanical resonator. Conversely, we show that engineered mechanical resonances drastically modify the optical response of our structures, creating large effective optical nonlinearities not present in bulk silicon. We experimentally demonstrate aspects of these nonlinearities by proposing and observing ``electromagnetically induced transparency'' and light slowed down to 6 m/s, as well as wavelength conversion, and generation of nonclassical optical radiation. Finally, the application of optomechanics to longstanding problems in quantum and classical communications are proposed and investigated.

Universal quantum computation with trapped ions in thermal motion by adiabatic passage

We propose a universal quantum computation scheme for trapped ions in thermal motion via the technique of adiabatic passage, which incorporates the advantages of both the adiabatic passage and the model of trapped ions in thermal motion. Our scheme is immune from the decoherence due to spontaneous emission from excited states as the system in our scheme evolves along a dark state. In our scheme the vibrational degrees of freedom are not required to be cooled to their ground states because they are only virtually excited. It is shown that the fidelity of the resultant gate operation is still high even when the magnitude of the effective Rabi frequency moderately deviates from the desired value.

A bayesian probabilistic approach to structural health monitoring

A Bayesian probabilistic methodology for on-line structural health monitoring which addresses the issue of parameter uncertainty inherent in problem is presented. The method uses modal parameters for a limited number of modes identified from measurements taken at a restricted number of degrees of freedom of a structure as the measured structural data. The application presented uses a linear structural model whose stiffness matrix is parameterized to develop a class of possible models. Within the Bayesian framework, a joint probability density function (PDF) for the model stiffness parameters given the measured modal data is determined. Using this PDF, the marginal PDF of the stiffness parameter for each substructure given the data can be calculated.

Monitoring the health of a structure using these marginal PDFs involves two steps. First, the marginal PDF for each model parameter given modal data from the undamaged structure is found. The structure is then periodically monitored and updated marginal PDFs are determined. A measure of the difference between the calibrated and current marginal PDFs is used as a means to characterize the health of the structure. A procedure for interpreting the measure for use by an expert system in on-line monitoring is also introduced.

The probabilistic framework is developed in order to address the model parameter uncertainty issue inherent in the health monitoring problem. To illustrate this issue, consider a very simplified deterministic structural health monitoring method. In such an approach, the model parameters which minimize an error measure between the measured and model modal values would be used as the "best" model of the structure. Changes between the model parameters identified using modal data from the undamaged structure and subsequent modal data would be used to find the existence, location and degree of damage. Due to measurement noise, limited modal information, and model error, the "best" model parameters might vary from one modal dataset to the next without any damage present in the structure. Thus, difficulties would arise in separating normal variations in the identified model parameters based on limitations of the identification method and variations due to true change in the structure. The Bayesian framework described in this work provides a means to handle this parametric uncertainty.

The probabilistic health monitoring method is applied to simulated data and laboratory data. The results of these tests are presented.

Efficient methods for stochastic optimal control

The Hamilton Jacobi Bellman (HJB) equation is central to stochastic optimal control (SOC) theory, yielding the optimal solution to general problems specified by known dynamics and a specified cost functional. Given the assumption of quadratic cost on the control input, it is well known that the HJB reduces to a particular partial differential equation (PDE). While powerful, this reduction is not commonly used as the PDE is of second order, is nonlinear, and examples exist where the problem may not have a solution in a classical sense. Furthermore, each state of the system appears as another dimension of the PDE, giving rise to the curse of dimensionality. Since the number of degrees of freedom required to solve the optimal control problem grows exponentially with dimension, the problem becomes intractable for systems with all but modest dimension.

In the last decade researchers have found that under certain, fairly non-restrictive structural assumptions, the HJB may be transformed into a linear PDE, with an interesting analogue in the discretized domain of Markov Decision Processes (MDP). The work presented in this thesis uses the linearity of this particular form of the HJB PDE to push the computational boundaries of stochastic optimal control.

This is done by crafting together previously disjoint lines of research in computation. The first of these is the use of Sum of Squares (SOS) techniques for synthesis of control policies. A candidate polynomial with variable coefficients is proposed as the solution to the stochastic optimal control problem. An SOS relaxation is then taken to the partial differential constraints, leading to a hierarchy of semidefinite relaxations with improving sub-optimality gap. The resulting approximate solutions are shown to be guaranteed over- and under-approximations for the optimal value function. It is shown that these results extend to arbitrary parabolic and elliptic PDEs, yielding a novel method for Uncertainty Quantification (UQ) of systems governed by partial differential constraints. Domain decomposition techniques are also made available, allowing for such problems to be solved via parallelization and low-order polynomials.

The optimization-based SOS technique is then contrasted with the Separated Representation (SR) approach from the applied mathematics community. The technique allows for systems of equations to be solved through a low-rank decomposition that results in algorithms that scale linearly with dimensionality. Its application in stochastic optimal control allows for previously uncomputable problems to be solved quickly, scaling to such complex systems as the Quadcopter and VTOL aircraft. This technique may be combined with the SOS approach, yielding not only a numerical technique, but also an analytical one that allows for entirely new classes of systems to be studied and for stability properties to be guaranteed.

The analysis of the linear HJB is completed by the study of its implications in application. It is shown that the HJB and a popular technique in robotics, the use of navigation functions, sit on opposite ends of a spectrum of optimization problems, upon which tradeoffs may be made in problem complexity. Analytical solutions to the HJB in these settings are available in simplified domains, yielding guidance towards optimality for approximation schemes. Finally, the use of HJB equations in temporal multi-task planning problems is investigated. It is demonstrated that such problems are reducible to a sequence of SOC problems linked via boundary conditions. The linearity of the PDE allows us to pre-compute control policy primitives and then compose them, at essentially zero cost, to satisfy a complex temporal logic specification.

Diseño y evaluación de un sistema de medida de aceleración en el elemento terminal de un manipulador de cinemática paralela

[ES]Este proyecto tiene como objetivo desarrollar una línea de investigación de opciones de sensorización de un mecanismo mediante acelerómetros. Se construirá para ello un sistema de adquisición y tratamiento de señales destinado a la sensorización de un mecanismo de cinemática paralela en base a los conocimientos adquiridos durante el curso. Se trabajará además con otros alumnos para llevar a cabo el diseño y montaje de un robot prototipo de cinemática paralela de dos grados de libertad sobre el que se experimentará y llevará a cabo el proyecto. Se plantean de este modo dos líneas de trabajo que se desarrollarán en este proyecto: Elaboración de un sistema de adquisición y tratamiento de señales adaptable a distintos sensores. Utilización de señales de múltiples acelerómetros para conocer en primer lugar aceleración, y de ser posible, posición de puntos de interés del mecanismo.