Stability of large flexible damped spacecraft modeled as elastic continua

Vibrational stability of a large flexible, structurally damped spacecraft subject to large rigid body rotations is analysed modelling the system as an elastic continuum. Using solution of rigid body attitude motion under torque free conditions and modal analysis, the vibrational equations are reduced to ordinary differential equations with time-varying coefficients. Stability analysis is carried out using Floquet theory and Sonin-Polya theorem. The cases of spinning and non-spinning spacecraft idealized as a flexible beam plate undergoing simple structural vibration are analysed in detail. The critical damping required for stabilization is shown to be a function of the spacecraft's inertia ratio and the level of disturbance.

NUVIEW:software for display and interactive manipulation of nucleic acid models

The NUVIEW software package allows skeletal models of any double helical nucleic acid molecule to be displayed out a graphics monitor and to apply various rotations, translations and scaling transformations interactively, through the keyboard. The skeletal model is generated by connecting any pair of representative points, one from each of the bases in the basepair. In addition to the above mentioned manipulations, the base residues can be identified by using a locator and the distance between any pair of residues can be obtained. A sequence based color coded display allows easy identification of sequence repeats, such as runs of Adenines. The real time interactive manipulation of such skeletal models for large DNA/RNA double helices, can be used to trace the path of the nucleic acid chain in three dimensions and hence get a better idea of its topology, location of linear or curved regions, distances between far off regions in the sequence etc. A physical picture of these features will assist in understanding the relationship between base sequence, structure and biological function in nucleic acids.

A real-time algorithm for simulation of flexible objects and hyper-redundant manipulators

Flexible objects such as a rope or snake move in a way such that their axial length remains almost constant. To simulate the motion of such an object, one strategy is to discretize the object into large number of small rigid links connected by joints. However, the resulting discretised system is highly redundant and the joint rotations for a desired Cartesian motion of any point on the object cannot be solved uniquely. In this paper, we revisit an algorithm, based on the classical tractrix curve, to resolve the redundancy in such hyper-redundant systems. For a desired motion of the `head' of a link, the `tail' is moved along a tractrix, and recursively all links of the discretised objects are moved along different tractrix curves. The algorithm is illustrated by simulations of a moving snake, tying of knots with a rope and a solution of the inverse kinematics of a planar hyper-redundant manipulator. The simulations show that the tractrix based algorithm leads to a more `natural' motion since the motion is distributed uniformly along the entire object with the displacements diminishing from the `head' to the `tail'.

Hamilton's theory of turns revisited

We present a new approach to Hamilton's theory of turns for the groups SO(3) and SU(2) which renders their properties, in particular their composition law, nearly trivial and immediately evident upon inspection. We show that the entire construction can be based on binary rotations rather than mirror reflections.

Anomalous dynamical charges, phonons, and the origin of negative thermal expansion in Y2W3O12

We use a combination of classical model and first-principles density functional theory calculations to study lattice dynamics of Y2W3O12 and identify phonons responsible for its negative thermal expansion (NTE). Born dynamical charges of various atoms are found to deviate anomalously from their nominal values. We find that the phonons with energy from 4 to 10 meV are the primary contributors to its NTE. These phonons involve rotations of the YO6 octahedra and WO4 tetrahedra in mutually opposite sense and collective translational atomic displacements, reflecting a strong mixing between acoustic and optic modes.

Consistent quaternion interpolation for objective finite element approximation of geometrically exact beam

We explore an isoparametric interpolation of total quaternion for geometrically consistent, strain-objective and path-independent finite element solutions of the geometrically exact beam. This interpolation is a variant of the broader class known as slerp. The equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations. We show that, for a two-noded beam element, the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient. The interpolation of total quaternion, which is computationally more efficient than the one based on local rotations, converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold. This interpolation, unlike the additive interpolation of total rotation, corresponds to a geodesic on the rotation manifold. For beam elements with more than two nodes, however, a consistent extension of the proposed quaternion interpolation is difficult. Alternatively, a quaternion-based procedure involving interpolation of relative rotations is proposed for such higher order elements. We also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions, proposed in [I. Romero, The interpolation of rotations and its application to finite element models of geometrically exact rods, Comput. Mech. 34 (2004) 121–133]. The strain-objectivity and path-independence of solutions are justified theoretically and then demonstrated through numerical experiments. This study, being focused only on the interpolation of rotations, uses a standard finite element discretization, as adopted by Simo and Vu-Quoc [J.C. Simo, L. Vu-Quoc, A three-dimensional finite rod model part II: computational aspects, Comput. Methods Appl. Mech. Engrg. 58 (1986) 79–116]. The rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector. Nodal rotations are stored in terms of rotation vectors and no secondary storages are required.

Crack Tip Fields in a Single Edge Notched Aluminum Single Crystal Specimen

We report a combined experimental and computational study of a low constraint aluminum single crystal fracture geometry and investigate the near-tip stress and strain fields. To this end, a single edge notched tensile (SENT) specimen is considered. A notch, with a radius of 50 µm, is taken to lie in the (010) plane and its front is aligned along the [101] direction. Experiments are conducted by subjecting the specimen to tensile loading using a special fixture inside a scanning electron microscope chamber. Both SEM micrographs and electron back-scattered diffraction (EBSD) maps are obtained from the near-tip region. The experiments are complemented by performing 3D and 2D plane strain finite element simulations within a continuum crystal plasticity framework assuming an isotropic hardening response characterized by the Pierce–Asaro–Needleman model. The simulations show a distinct slip band forming at about 55 deg with respect to the notch line corresponding to slip on (11-bar 1)[011] system, which corroborates well with experimental data. Furthermore, two kink bands occur at about 45 deg and 90 deg with respect to the notch line within which large rotations in the crystal orientation take place. These predictions are in good agreement with the EBSD observations. Finally, the near-tip angular variations of the 3D stress and plastic strain fields in the low constraint SENT fracture geometry are examined in detail.

Secondary flow induced by the rotation of two concentric spheres about their common diameter in micropolar fluids

In this paper, we have studied the secondary flow induced in a micropolar fluid by the rotation of two concentric spheres about a fixed diameter. The secondary flow exhibits behaviour commonly observed in visco-elastic fluids. In particular we have obtained the expressions for microrotation vector. Numerical results have been obtained for a number of values of relative rotations of the two spheres for a chosen set of values of fluid parameters. The results are presented graphically and compared with the previous investigations.

Secondary flows induced by the rotation of a sphere or coaxial cones in micropolar fluids

The micropolar fluids like Newtonian and Non-Newtonian fluids cannot sustain a simple shearing motion, wherein only one component of velocity is present. They exhibit both primary and secondary motions when the boundaries are subject to slow rotations. The primary motion, as in Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., resembles that of Newtonian fluid for slow steady rotation. We further notice that the micro-rotation becomes identically equal to the vorticity present in the fluid and the condition b) of "Wall vorticity" can alone be satisfied at the boundaries. As regards, the secondary motion, we notice that it can be determined by the above procedure for a special class of fluids, namely that for which j0(n2-n3)=4 n3/l2. Moreover for this class of fluids, the micro-rotation is identical with the vorticity of the fluid everywhere. Also the stream function for the secondary flow is identical with that for the Newtonian fluid with a suitable definition of the Reynolds number. In contrast with the Non-Newtonian fluids, characterized by the equation due to Rivlin-Ericksen, Oldroyd, Walters etc., this class of micropolar fluids does not show separation. This is in conformity with the statement of Condiff and Dahler (3) that in any steady flow, internal spin matches the vorticity everywhere provided that (i) spin boundary conditions are satisfied, (ii) body torques and non-conservative body forces are absent, and (iii) inertial and spin-inertial terms are either negligible or vanish identically.

3D shape extraction of human face in presence of facial hair: A profilometric approach

The increasing use of 3D modeling of Human Face in Face Recognition systems, User Interfaces, Graphics, Gaming and the like has made it an area of active study. Majority of the 3D sensors rely on color coded light projection for 3D estimation. Such systems fail to generate any response in regions covered by Facial Hair (like beard, mustache), and hence generate holes in the model which have to be filled manually later on. We propose the use of wavelet transform based analysis to extract the 3D model of Human Faces from a sinusoidal white light fringe projected image. Our method requires only a single image as input. The method is robust to texture variations on the face due to space-frequency localization property of the wavelet transform. It can generate models to pixel level refinement as the phase is estimated for each pixel by a continuous wavelet transform. In cases of sparse Facial Hair, the shape distortions due to hairs can be filtered out, yielding an estimate for the underlying face. We use a low-pass filtering approach to estimate the face texture from the same image. We demonstrate the method on several Human Faces both with and without Facial Hairs. Unseen views of the face are generated by texture mapping on different rotations of the obtained 3D structure. To the best of our knowledge, this is the first attempt to estimate 3D for Human Faces in presence of Facial hair structures like beard and mustache without generating holes in those areas.

Dielectric behaviour and ferroelectricity in alums

Introduction The alum family of double salts with the general formula M1+M3+(RO4)2.12H2O where M1+ is a monovalent ion (M1+ = K, Rb, Cs, Tl, NH4, CH3NH3, NH3OH3 NH3NH2, etc.) and M3+ is a trivalent metal such as Al, Fe, Cr, V, In, Ga, etc. and R is S or Se, form an isomorphous series and their general features indicate a common cubic space group Pa3. Lipson1 showed subsequently that there exist three different structure types agr, β and γ and the structure of a particular alum is dependent on the radius of the monovalent atom. The agr structure is typical of medium sized ions, the β of the larger ones and the γ of the small Na atom.2 Ferroelectricity has been reported only in alums containing NH4, CH3NH3, NH3NH3 and NH3OH. Their hindered rotations as well as the influence of sulphate group disorder on the dielectric behaviour of alums is still not clear.3 No study of the temperature dependence of the low frequency dielectric constant of some of the alums, particularly those of Cs, Rb and Tl, have been made so The present investigation was undertaken to correlate their dielectric behaviour with their composition and structural differences. Under the same experimental conditions, methyl ammonium and ammonium alums also were studied and compared with the known data.

An unsteady flow due to eccentrically rotating porous disk and a fluid at infinity

An exact solution of the unsteady Navier-Stokes equations is obtained for the flow due to non-coaxial rotations of a porous disk, executing non-torsional oscillations in its own plane, and a fluid at infinity. It is shown that the infinite number of solutions existing for a flow confined between two disks reduce to a single unique solution in the case of a single disk. The adjustment of the unsteady flow near the rotating disk to the flow at infinity rotating about a different axis is explained.

Characterization of hot deformation behavior of brasses using processing maps: Part II. β Brass and α-β brass

The hot deformation behaviors of β brass in the temperature range of 550°C to 800°C and α-β brass in the temperature range of 450°C to 800°C have been characterized in the strain rate range of 0.001 to 100 s−1 using processing maps developed on the basis of the Dynamic Materials Model. The map for β brass revealed a domain of superplasticity in the entire temperature range and at strain rates lower than 1 s−1, with a maximum efficiency of power dissipation of about 68 pct. The temperature variation of the efficiency of power dissipation in the domain is similar to that of the diffusion coefficient for zinc in β brass, confirming that the diffusion-accommodated flow controls the superplasticity. The material undergoes microstructural instability in the form of adiabatic shear bands and strain markings at temperatures lower than 700°C and at strain rates higher than 10 s−1. The map for α-β brass revealed a wide domain for processing in the temperature range of 550°C to 800°C and at strain rates lower than 1 s−1, with a maximum efficiency of 54 pct occurring at about 750°C and 0.001 s−1. In the domain, the α phase undergoes dynamic recrystallization and controls the hot deformation of the alloy, while the β phase deforms superplastically. At strain rates greater than 1 s−1, α-β brass exhibits microstructural instabilities manifested as flow rotations at lower temperatures and localized shear bands at higher temperatures.

Vibrational Spectra and Normal Coordinate Analysis of N-Methyl Thiourea

The infrared spectra of N-methylthiourea (NMTU) and its N-deuterated and S-methylated species were measured. Assignment of the infrared and Raman spectra of NMTU has been accomplished by correlation with thiourea and by use of infrared band shifts on N-deuteration as well as S-methylation. Normal coordinate analysis was performed for all the fundamentals of NMTU and NMTU-d3, the assignments obtained from the force field calculations being discussed in relation to those in other related thioureas and thioamides. The potential barriers to the internal rotations for the �NH2, �CH3, and �CN groups were estimated from the force constants.

On the relation between rotation increments in different tangent spaces

In computational mechanics, finite rotations are often represented by rotation vectors. Rotation vector increments corresponding to different tangent: spaces are generally related by a linear operator, known as the tangential transformation T. In this note, we derive the higher order terms that are usually left out in linear relation. The exact nonlinear relation is also presented. Errors via the linearized T are numerically estimated. While the concept of T arises out of the nonlinear characteristics of the rotation manifold, it has been derived via tensor analysis in the context of computational mechanics (Cardona and Geradin, 1988). We investigate the operator T from a Lie group perspective, which provides a better insight and a 1-1 correspondence between approaches based on tensor analysis and the standard matrix Lie group theory. (C) 2010 Elsevier Ltd. All rights reserved.