Optimum design of a helicopter rotor for low vibration using aeroelastic analysis and response surface methods

An aeroelastic analysis based on finite elements in space and time is used to model the helicopter rotor in forward flight. The rotor blade is represented as an elastic cantilever beam undergoing flap and lag bending, elastic torsion and axial deformations. The objective of the improved design is to reduce vibratory loads at the rotor hub that are the main source of helicopter vibration. Constraints are imposed on aeroelastic stability, and move limits are imposed on the blade elastic stiffness design variables. Using the aeroelastic analysis, response surface approximations are constructed for the objective function (vibratory hub loads). It is found that second order polynomial response surfaces constructed using the central composite design of the theory of design of experiments adequately represents the aeroelastic model in the vicinity of the baseline design. Optimization results show a reduction in the objective function of about 30 per cent. A key accomplishment of this paper is the decoupling of the analysis problem and the optimization problems using response surface methods, which should encourage the use of optimization methods by the helicopter industry. (C) 2002 Elsevier Science Ltd. All rights reserved.

Symmetry-breaking transitions in equilibrium shapes of coherent precipitates: Effect of elastic anisotropy and inhomogeneity

We examine the symmetry-breaking transitions in equilibrium shapes of coherent precipitates in two-dimensional (2-D) systems under a plane-strain condition with the principal misfit strain components epsilon(xx)*. and epsilon(yy)*. For systems with cubic elastic moduli, we first show all the shape transitions associated with different values of t = epsilon(yy)*/epsilon(xx)*. We also characterize each of these transitions, by studying its dependence on elastic anisotropy and inhomogeneity. For systems with dilatational misfit (t = 1) and those with pure shear misfit (t = -1), the transition is from an equiaxed shape to an elongated shape, resulting in a break in rotational symmetry. For systems with nondilatational misfit (-1 < t < 1; t not equal 0), the transition involves a break in mirror symmetries normal to the x- and y-axes. The transition is continuous in all cases, except when 0 < t < 1. For systems which allow an invariant line (-1 less than or equal to t < 0), the critical size increases with an increase in the particle stiffness. However, for systems which do not allow an invariant line (0 < t less than or equal to 1), the critical size first decreases, reaches a minimum, and then starts increasing with increasing particle stiffness; moreover, the transition is also forbidden when the particle stiffness is greater than a critical value.

Nonlinear Dynamics of Beams with Stochastic Parameter Variations

The aim of this paper is to investigate the steady state response of beams under the action of random support motions. The study is of relevance in the context of earthquake response of extended land based structures such as pipelines and long span bridges, and, secondary systems such as piping networks in nuclear power plant installations. The following complicating features are accounted for in the response analysis: (a) differential support motions: this is characterized in terms of cross power spectral density functions associated with distinct support motions, (b) nonlinear support conditions, and (c) stochastically inhomogeneous stiffness and mass variations of the beam structure; questions on non-Gaussian models for these variations are considered. The method of stochastic finite elements is combined with equivalent linearization technique and Monte Carlo simulations to obtain response moments.

Nonlinear Finite-Element Modeling of Batter Piles under Lateral Load

In this paper, a finite-element model is developed in which the nonlinear soil behavior is represented by a hyperbolic relation for static load condition and modified hyperbolic relation, which includes both degradation and gap for a cyclic load condition. Although batter piles are subjected to lateral load, the soil resistance is also governed by axial load, which is incorporated by considering the P-Δ moment and geometric stiffness matrix. By adopting the developed numerical model, static and cyclic load analyses are performed adopting an incremental-iterative procedure where the pile is idealized as beam elements and the soil as elastoplastic spring elements. The proposed numerical model is validated with published laboratory and field pile test results under both static and cyclic load conditions. This paper highlights the importance of the degradation factor and its influence on the soil resistance-displacement (p-y) curve, number of cycles of loading, and cyclic load response.

Temperature profiles of accretion discs around rapidly rotating strange stars in general relativity: A comparison with neutron stars

We compute the temperature profiles of accretion discs around rapidly rotating strange stars, using constant gravitational mass equilibrium sequences of these objects, considering the full effect of general relativity. Beyond a certain critical value of stellar angular momentum (J), we observe the radius ( $r_{\rm orb}$) of the innermost stable circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes nor in rotating neutron stars). The reason for this is traced to the crucial dependence of ${\rm d}r_{\rm orb}/{\rm d}J$ on the rate of change of the radial gradient of the Keplerian angular velocity at $r_{\rm orb}$ with respect to J. The structure parameters and temperature profiles obtained are compared with those of neutron stars, as an attempt to provide signatures for distinguishing between the two. We show that when the full gamut of strange star equation of state models, with varying degrees of stiffness are considered, there exists a substantial overlap in properties of both neutron stars and strange stars. However, applying accretion disc model constraints to rule out stiff strange star equation of state models, we notice that neutron stars and strange stars exclusively occupy certain parameter spaces. This result implies the possibility of distinguishing these objects from each other by sensitive observations through future X-ray detectors.

A new beam finite element for the analysis of functionally graded materials

A new beam element is developed to study the thermoelastic behavior of functionally graded beam structures. The element is based on the first-order shear deformation theory and it accounts for varying elastic and thermal properties along its thickness. The exact solution of static part of the governing differential equations is used to construct interpolating polynomials for the element formulation. Consequently, the stiffness matrix has super-convergent property and the element is free of shear locking. Both exponential and power-law variations of material property distribution are used to examine different stress variations. Static, free vibration and wave propagation problems are considered to highlight the behavioral difference of functionally graded material beam with pure metal or pure ceramic beams. (C) 2003 Elsevier Science Ltd. All rights reserved.

A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

An Improved Compact Compliant Mechanism for an External Pipe-Crawler

An external pipe-crawling device presented in this paper aids the inspection of pipes in hazardous environments and areas inaccessible to humans. The principal component of our design, which uses inchworm type motion, is a compliant ring mechanism actuated using shape memory alloy (SMA) wire. It was fabricated and tested and was reported in our earlier work. But this device had a drawback of low crawling speed (about 1 mm/min) owing to the delay in heating and cooling of the SMA strips in the linear actuation. Additionally, that design also had the difficulties of mounting on pipes with closed ends, large radial span, and the need for housing for electrical insulation and guiding of the SMA wire. In this paper we present a compact design that overcomes the difficulties of the earlier design. In particular, we present a compact compliant mechanism with two halves so as to enable mounting and un-mounting on any closed or open pipe. Another feature is the presence of insulation and guiding of the SMA wire without housing. This design results in a reduction of the radial span of the ring from 22 mm to 12 mm, and the stiffness of the mechanism and the SMA wire are matched. An SMA helical spring is to used in the place of an SMA strip to increase the crawling speed of the device. A microcontroller-based circuitry is also fitted to cyclically.activate the SMA wires and springs.

A Compliant Mechanism Kit with Flexible Beams and Connectors

We present the concept, prototypes, and an optimal design method for a compliant mechanism kit as a parallel to the kits available for rigid-body mechanisms. The kit consists of flexible beams and connectors that can be easily hand-assembled using snap fits. It enables users, using their creativity and mechanics intuition, to quickly realize a compliant mechanism. The mechanisms assembled in this manner accurately capture the essential behavior of the topology, shape, size and material aspects and thereby can lead the way for a real compliant mechanism for practical use. Also described in this paper are the design of the connector to which flexible beams can be added in eight different directions; and prototyping of the spring steel connectors as well as beams using wire-cut electro discharge machining. It is noted in this paper that the concept of the kit also resolves a discrepancy in the finite element (FE) modeling of beam-based compliant mechanisms. The discrepancy arises when two or more beams are joining at one point and thus leading to increased stiffness. After resolving this discrepancy, this work extends the topology optimization to automatically generate designs that can be assembled with the kit. Thus, the kit and the accompanying analysis and optimal synthesis procedures comprise a self-contained educational as well as a research and pragmatic toolset for compliant mechanisms. The paper also illustrates how human creativity finds new ways of using the kit beyond the original intended use and how it is useful even for a novice to design compliant mechanisms.

Doubly spectral stochastic finite element method (DSSFEM) for random field problems

Uncertainties in complex dynamic systems play an important role in the prediction of a dynamic response in the mid- and high-frequency ranges. For distributed parameter systems, parametric uncertainties can be represented by random fields leading to stochastic partial differential equations. Over the past two decades, the spectral stochastic finite-element method has been developed to discretize the random fields and solve such problems. On the other hand, for deterministic distributed parameter linear dynamic systems, the spectral finite-element method has been developed to efficiently solve the problem in the frequency domain. In spite of the fact that both approaches use spectral decomposition (one for the random fields and the other for the dynamic displacement fields), very little overlap between them has been reported in literature. In this paper, these two spectral techniques are unified with the aim that the unified approach would outperform any of the spectral methods considered on their own. An exponential autocorrelation function for the random fields, a frequency-dependent stochastic element stiffness, and mass matrices are derived for the axial and bending vibration of rods. Closed-form exact expressions are derived by using the Karhunen-Loève expansion. Numerical examples are given to illustrate the unified spectral approach.

Viscoelastic phenomenology based structure assignment for closed-loop vibration control of a beam with sensors and actuators

In this paper we incorporate a novel approach to synthesize a class of closed-loop feedback control, based on the variational structure assignment. Properties of a viscoelastic system are used to design an active feedback controller for an undamped structural system with distributed sensor, actuator and controller. Wave dispersion properties of onedimensional beam system have been studied. Efficiency of the chosen viscoelastic model in enhancing damping and stability properties of one-dimensional viscoelastic bar have been analyzed. The variational structure is projected on a solution space of a closed-loop system involving a weakly damped structure with distributed sensor and actuator with controller. These assign the phenomenology based internal strain rate damping parameter of a viscoelastic system to the usual elastic structure but with active control. In the formulation a model of cantilever beam with non-collocated actuator and sensor has been considered. The formulation leads to the matrix identification problem of two dynamic stiffness matrices. The method has been simplified to obtain control system gains for the free vibration control of a cantilever beam system with collocated actuator-sensor, using quadratic optimal control and pole-placement methods.

Axial wave propagation in coupled nanorod system with nonlocal small scale effects

This article deals with the axial wave propagation properties of a coupled nanorod system with consideration of small scale effects. The nonlocal elasticity theory has been incorporated into classical rod/bar model to capture unique features of the coupled nanorods under the umbrella of continuum mechanics theory. Nonlocal rod model is developed for coupled nanorods. The strong effect of the nonlocal scale has been obtained which leads to substantially different wave behavior of nanorods from those of macroscopic rods. Explicit expressions are derived for wavenumber, cut-off frequency and escape frequency of nanorods. The analysis shows that the wave characteristics of nanorods are highly over estimated by the classical rod model, which ignores the effect of small-length scale. The studies also shows that the nonlocal scale parameter introduces certain band gap region in axial or longitudinal wave mode, where no wave propagation occurs. This is manifested in the spectrum cures as the region, where the wavenumber tends to infinite or wave speed tends to zero. The effect of the coupled spring stiffness is also capture in the present analysis. It has been also shown that the cut-off frequency increases as the stiffness of the coupled spring increases and also the coupled spring stiffness has no effect on escape frequency of the axial wave mode in the nanorod. This cut-off frequency is also independent of the nonlocal small scale parameter. The present study may bring in helpful insights while investigating multiple-nanorod-system-models for future nano-optomechanical systems applications. The results can also provide useful guidance for the study and design of the next generation of nanodevices that make use of the wave propagation properties of coupled single-walled carbon nanotubes or coupled nanorods. (C) 2011 Elsevier Ltd. All rights reserved.

Design of Geocell Reinforcement for Supporting Embankments on Soft Ground

The methods of design available for geocell-supported embankments are very few. Two of the earlier methods are considered in this paper and a third method is proposed and compared with them. The first method is the slip line method proposed by earlier researchers. The second method is based on slope stability analysis proposed by this author earlier and the new method proposed is based on the finite element analyses. In the first method, plastic bearing failure of the soil was assumed and the additional resistance due to geocell layer is calculated using a non-symmetric slip line field in the soft foundation soil. In the second method, generalpurpose slope stability program was used to design the geocell mattress of required strength for embankment using a composite model to represent the shear strength of geocell layer. In the third method proposed in this paper, geocell reinforcement is designed based on the plane strain finite element analysis of embankments. The geocell layer is modelled as an equivalent composite layer with modified strength and stiffness values. The strength and dimensions of geocell layer is estimated for the required bearing capacity or permissible deformations. These three design methods are compared through a design example.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

Wave dispersion in a cellular composite with modulated microstructure

In this paper we report a modeling technique and analysis of wave dispersion in a cellular composite laminate with spatially modulated microstructure, which can be modeled by parameterization and homogenization in an appropriate length scale. Higher order beam theory is applied and the system of wave equations are derived. Homogenization of these equations are carried out in the scale of wavelength and frequency of the individual wave modes. Smaller scale scattering below the order of cell size are filtered out in the present approach. The longitudinal dispersion relations for different values of a modulation parameter are analyzed which indicates the existence of stop and pass band patterns. Dispersion relations for flexural-shear case are also analyzed which indicates a tendency toward forming the stop and pass bands for increasing values of a shear stiffness modulation parameter. The effect the phase angle (θ) of the incident wave indicates the existence more number of alternative stop bands and pass bands for θ = 45°.