Variance-Reduced Particle Filters for Structural System Identification Problems

A few variance reduction schemes are proposed within the broad framework of a particle filter as applied to the problem of structural system identification. Whereas the first scheme uses a directional descent step, possibly of the Newton or quasi-Newton type, within the prediction stage of the filter, the second relies on replacing the more conventional Monte Carlo simulation involving pseudorandom sequence with one using quasi-random sequences along with a Brownian bridge discretization while representing the process noise terms. As evidenced through the derivations and subsequent numerical work on the identification of a shear frame, the combined effect of the proposed approaches in yielding variance-reduced estimates of the model parameters appears to be quite noticeable. DOI: 10.1061/(ASCE)EM.1943-7889.0000480. (C) 2013 American Society of Civil Engineers.

Biomechanics of Substrate Boring by Fig Wasps: Role of Zinc in Insect Cuticle

There are many biomechanical challenges that a female insect must meet to successfully oviposit and ensure her evolutionary success. These begin with selection of a suitable substrate through which the ovipositor must penetrate without itself buckling or fracturing. The second phase corresponds to steering and manipulating the ovipositor to deliver eggs at desired locations. Finally, the insect must retract her ovipositor fast to avoid possible predation and repeat this process multiple times during her lifetime. From a materials perspective, insect oviposition is a fascinating problem and poses many questions. Specifically, are there diverse mechanisms that insects use to drill through hard substrates without itself buckling or fracturing? What are the structure-property relationships in the ovipositor material? These are some of the questions we address with a model system consisting of a parasitoid fig wasp - fig substrate system. To characterize the structure of ovipositors, we use scanning electron microscopy with a detector to quantify the presence of transition elements. Our results show that parasitoid ovipositors have teeth like structures on their tips and contain high amounts of zinc as compared to remote regions. Sensillae are present along the ovipositor to aid detection of chemical species and mechanical deformations. To quantify the material properties of parasitoid ovipositors, we use an atomic force microscope and show that tip regions have higher modulus as compared to remote regions. Finally, we use videography to show that ovipositors buckle during oviposition and estimate the forces needed to cause substrate boring based on Euler buckling analysis. Such methods may be useful for the design of functionally graded surgical tools.

Rotating beams isospectral to axially loaded nonrotating uniform beams

In this paper, the stiffness and mass per unit length distributions of a rotating beam, which is isospectral to a given uniform axially loaded nonrotating beam, are determined analytically. The Barcilon-Gottlieb transformation is extended so that it transforms the governing equation of a rotating beam into the governing equation of a uniform, axially loaded nonrotating beam. Analysis is limited to a certain class of Euler-Bernoulli cantilever beams, where the product between the stiffness and the cube of mass per unit length is a constant. The derived mass and stiffness distributions of the rotating beam are used in a finite element analysis to confirm the frequency equivalence of the given and derived beams. Examples of physically realizable beams that have a rectangular cross section are shown as a practical application of the analysis.

On the Chern number of I-admissible filtrations of ideals

Let I be an m-primary ideal of a Noetherian local ring (R, m) of positive dimension. The coefficient e(1)(I) of the Hilbert polynomial of an I-admissible filtration I is called the Chern number of I. A formula for the Chern number has been derived involving the Euler characteristic of subcomplexes of a Koszul complex. Specific formulas for the Chern number have been given in local rings of dimension at most two. These have been used to provide new and unified proofs of several results about e(1)(I).

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

We undertake a systematic, direct numerical simulation of the twodimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. Firstly, there are transients that depend on the initial conditions. In the second regime, powerlaw scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power laws and the extents of the scaling regions change with time and depend on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers k > kc and partial thermalization takes place for modes with k < kc; the self-truncation wave number kc(t) depends on the initial conditions and it grows either as a power of t or as log t. Finally, in the fourth regime, complete thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms. Our work is a natural generalization of recent studies of thermalization in the Euler and other hydrodynamical equations; it combines ideas from fluid dynamics and turbulence, on the one hand, and equilibrium and nonequilibrium statistical mechanics on the other.

Generalized model predictive static programming and its application to 3D impact angle constrained guidance of air-to-surface missiles

A new `generalized model predictive static programming (G-MPSP)' technique is presented in this paper in the continuous time framework for rapidly solving a class of finite-horizon nonlinear optimal control problems with hard terminal constraints. A key feature of the technique is backward propagation of a small-dimensional weight matrix dynamics, using which the control history gets updated. This feature, as well as the fact that it leads to a static optimization problem, are the reasons for its high computational efficiency. It has been shown that under Euler integration, it is equivalent to the existing model predictive static programming technique, which operates on a discrete-time approximation of the problem. Performance of the proposed technique is demonstrated by solving a challenging three-dimensional impact angle constrained missile guidance problem. The problem demands that the missile must meet constraints on both azimuth and elevation angles in addition to achieving near zero miss distance, while minimizing the lateral acceleration demand throughout its flight path. Both stationary and maneuvering ground targets are considered in the simulation studies. Effectiveness of the proposed guidance has been verified by considering first order autopilot lag as well as various target maneuvers.

FUZZY CONIFOLD Y-F(6) AND MONOPOLES ON S-F(2) x S-F(2)

In this paper, we construct the fuzzy (finite-dimensional) analogs of the conifold Y-6 and its base X-5. We show that fuzzy X-5 is (the analog of) a principal U(1) bundle over fuzzy spheres S-F(2) x S-F(2) and explicitly construct the associated monopole bundles. In particular, our construction provides an explicit discretization of the spaces T-k,T-k and T-k,T-0.

Fractal dimension method for delamination detection in composite beams

A new delaminated composite beam element is formulated for Timoshenko as well as Euler-Bernoulli beam models. Shape functions are derived from Timoshenko functions; this provides a unified formulation for slender to moderately deep beam analyses. The element is simple and easy to implement, results are on par with those from free mode delamination models. Katz fractal dimension method is applied on the mode shapes obtained from finite element models, to detect the delamination in the beam. The effect of finite element size on fractal dimension method of delamination detection is quantified.

Shock-boundary layer interaction and transonic flutter

The transonic flutter dip of an aeroelastic system is primarily caused by compressibility of the flowing fluid. Viscous effects are not dominant in the pre-transonic dip region. In fact, an Euler solver can predict this flutter boundary with considerable accuracy. However with an increase in Mach number the shock moves towards the trailing edge causing shock induced separation. This shock-boundary layer interaction changes the flutter boundary in the transonic and post-transonic dip region significantly. We discuss the effect of viscosity in changing the flutter boundary in the post-transonic dip region using a RANS solver coupled to a two-degree of freedom model of the structural dynamics of a wing.

Transonic flutter, shocks, and parameter effects

There is a drop in the flutter boundary of an aeroelastic system placed in a transonic flow due to compressibility effects and is known as the transonic dip. Viscous effects can shift the lo-cation of the shock and depending on the shock strength the boundary layer may separate leading to changes in the flutter speed. An unsteady Euler flow solver coupled with the structural dynamic equations is used to understand the effect of shock on the transonic dip. The effect of various system parameters such as mass ratio, location of the center of mass, position of the elastic axis, ratio of uncoupled natural frequencies in heave and pitch are also studied. Steady turbulent flow results are presented to demonstrate the effect of viscosity on the location and strength of the shock.

A parametric study on the buckling of functionally graded material plates with internal discontinuities using the partition of unity method

In this paper, the effect of local defects, viz., cracks and cutouts on the buckling behaviour of functionally graded material plates subjected to mechanical and thermal load is numerically studied. The internal discontinuities, viz., cracks and cutouts are represented independent of the mesh within the framework of the extended finite element method and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction and the effective properties are estimated using the Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation theory. The influence of various parameters, viz., the crack length and its location, the cutout radius and its position, the plate aspect ratio and the plate thickness on the critical buckling load is studied. The effect of various boundary conditions is also studied. The numerical results obtained reveal that the critical buckling load decreases with increase in the crack length, the cutout radius and the material gradient index. This is attributed to the degradation in the stiffness either due to the presence of local defects or due to the change in the material composition. (C) 2013 Elsevier Masson SAS. All rights reserved.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Freeform Skeletal Shape Optimization of Compliant Mechanisms

Compliant mechanisms are elastic continua used to transmit or transform force and motion mechanically. The topology optimization methods developed for compliant mechanisms also give the shape for a chosen parameterization of the design domain with a fixed mesh. However, in these methods, the shapes of the flexible segments in the resulting optimal solutions are restricted either by the type or the resolution of the design parameterization. This limitation is overcome in this paper by focusing on optimizing the skeletal shape of the compliant segments in a given topology. It is accomplished by identifying such segments in the topology and representing them using Bezier curves. The vertices of the Bezier control polygon are used to parameterize the shape-design space. Uniform parameter steps of the Bezier curves naturally enable adaptive finite element discretization of the segments as their shapes change. Practical constraints such as avoiding intersections with other segments, self-intersections, and restrictions on the available space and material, are incorporated into the formulation. A multi-criteria function from our prior work is used as the objective. Analytical sensitivity analysis for the objective and constraints is presented and is used in the numerical optimization. Examples are included to illustrate the shape optimization method.

Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Edge states of a three-dimensional topological insulator

We use the bulk Hamiltonian for a three-dimensional topological insulator such as Bi-2 Se-3 to study the states which appear on its various surfaces and along the edge between two surfaces. We use both analytical methods based on the surface Hamiltonians (which are derived from the bulk Hamiltonian) and numerical methods based on a lattice discretization of the bulk Hamiltonian. We find that the application of a potential barrier along an edge can give rise to states localized at that edge. These states have an unusual energy-momentum dispersion which can be controlled by applying a potential along the edge; in particular, the velocity of these states can be tuned to zero. The scattering and conductance across the edge is studied as a function of the edge potential. We show that a magnetic field in a particular direction can also give rise to zero energy states on certain edges. We point out possible experimental ways of looking for the various edge states.