Response of nonlinear systems to narrow-band excitation

The hardening cubic spring oscillator is studied under narrow-band gaussian excitation. Equivalent linearization leads to multiple steady states. The realizability of the solution is discussed through stochastic stability analysis. Theoretical results are supported by numerical simulation.

Superharmonic vibrations of order 3 in stretched strings

Superharmonic vibrations of order 3 in stretched strings driven by a single-mode planar simple-harmonic force are investigated. It is shown that planar as well as nonplanar superharmonic resonances can occur. After giving a few analytical results, the problem is thoroughly investigated numerically. The stability analysis shows that the region of stable nonplanar oscillations is much smaller than that of stable planar oscillation. It is observed in the case of nonplanar oscillations that the amplitude of the superharmonic in the plane normal to the plane of excitation is larger than that in the plane of excitation.

Narrowband random excitation of a limit cycle system

The response of the Van der Pol oscillator to stationary narrowband Gaussian excitation is considered. The central frequency of excitation is taken to be in the neighborhood of the system limit cycle frequency. The solution is obtained using a non-Gaussian closure approximation on the probability density function of the response. The validity of the solution is examined with the help of a stochastic stability analysis. Solution based on Stratonovich''s quasistatic averaging technique is also obtained. The comparison of the theoretical solutions with the digital simulations shows that the theoretical estimates are reasonably good.

Instability and dewetting of ultrathin solid viscoelastic films on homogeneous and heterogeneous substrates

Instability and dewetting engendered by the van der Waals force in soft thin (<100 nm) linear viscoelastic solid (e. g., elastomeric gel) films on uniform and patterned surfaces are explored. Linear stability analysis shows that, although the elasticity of the film controls the onset of instability and the corresponding critical wavelength, the dominant length-scale remains invariant with the elastic modulus of the film. The unstable modes are found to be long-wave, for which a nonlinear long-wave analysis and simulations are performed to uncover the dynamics and morphology of dewetting. The stored elastic energy slows down the temporal growth of instability significantly. The simulations also show that a thermodynamically stable film with zero-frequency elasticity can be made unstable in the presence of physico-chemical defects on the substrate and can follow an entirely different pathway with far fewer holes as compared to the viscous films. Further, the elastic restoring force can retard the growth of a depression adjacent to the hole-rim and thus suppress the formation of satellite holes bordering the primary holes. These findings are in contrast to the dewetting of viscoelastic liquid films where nonzero frequency elasticity accelerates the film rupture and promotes the secondary instabilities. Thus, the zero-frequency elasticity can play a major role in imposing a better-defined long-range order to the dewetted structures by arresting the secondary instabilities. (C) 2011 American Institute of Physics. doi: 10.1063/1.3554748]

3-D Warping in Four-Bar Laminated Linkages

This paper deals with the evaluation of the component-laminate load-carrying capacity, i.e., to calculate the loads that cause the failure of the individual layers and the component-laminate as a whole in four-bar mechanism. The component-laminate load-carrying capacity is evaluated using the Tsai-Wu-Hahn failure criterion for various layups. The reserve factor of each ply in the component-laminate is calculated by using the maximum resultant force and the maximum resultant moment occurring at different time steps at the joints of the mechanism. Here, all component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically nonlinear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and nonlinear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional nonlinearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the nonlinear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the nonlinear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict more quickly and accurately than would otherwise be possible. Local 3-D stress, strain and displacement fields for representative sections in the component-bars are recovered, based on the stress resultants from the 1-D global beam analysis. A numerical example is presented which illustrates the failure of each component-laminate and the mechanism as a whole.

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.

Instabilities of soft elastic microtubes filled with viscous fluids: Pearls, wrinkles, and sausage strings

A linear stability analysis is presented to study the self-organized instabilities of a highly compliant elastic cylindrical shell filled with a viscous liquid and submerged in another viscous medium. The prototype closely mimics many components of micro-or nanofluidic devices and biological processes such as the budding of a string of pearls inside cells and sausage-string formation of blood vessels. The cylindrical shell is considered to be a soft linear elastic solid with small storage modulus. When the destabilizing capillary force derived from the cross-sectional curvature overcomes the stabilizing elastic and in-plane capillary forces, the microtube can spontaneously self-organize into one of several possible configurations; namely, pearling, in which the viscous fluid in the core of the elastic shell breaks up into droplets; sausage strings, in which the outer interface of the mircrotube deforms more than the inner interface; and wrinkles, in which both interfaces of the thin-walled mircrotube deform in phase with small amplitudes. This study identifies the conditions for the existence of these modes and demonstrates that the ratios of the interfacial tensions at the interfaces, the viscosities, and the thickness of the microtube play crucial roles in the mode selection and the relative amplitudes of deformations at the two interfaces. The analysis also shows asymptotically that an elastic fiber submerged in a viscous liquid is unstable for Y = gamma/(G(e)R) > 6 and an elastic microchannel filled with a viscous liquid should rupture to form spherical cavities (pearling) for Y > 2, where gamma, G(e), and R are the surface tension, elastic shear modulus, and radius, respectively, of the fiber or microchannel.

Contact instabilities of anisotropic and inhomogeneous soft elastic films

Anisotropy plays important roles in various biological phenomena such as adhesion of geckos and grasshoppers enabled by the attachment pods having hierarchical structures like thin longitudinal setae connected with threads mimicked by anisotropic films. We study the contact instability of a transversely isotropic thin elastic film when it comes in contact proximity of another surface. In the present study we investigate the contact stability of a thin incompressible transversely isotropic film by performing linear stability analysis. Based on the linear stability analysis, we show that an approaching contactor renders the film unstable. The critical wavelength of the instability is a function of the total film thickness and the ratio of the Young's modulus in the longitudinal direction and the shear modulus in the plane containing the longitudinal axis. We also analyze the stability of a thin gradient film that is elastically inhomogeneous across its thickness. Compared to a homogeneous elastic film, it becomes unstable with a longer wavelength when the film becomes softer in going from the surface to the substrate.

Real-time monitoring of critical nodes with minimal number of Phasor Measurement Units

This paper presents a method for placement of Phasor Measurement Units, ensuring the monitoring of vulnerable buses which are obtained based on transient stability analysis of the overall system. Real-time monitoring of phase angles across different nodes, which indicates the proximity to instability, the very purpose will be well defined if the PMUs are placed at buses which are more vulnerable. The issue is to identify the key buses where the PMUs should be placed when the transient stability prediction is taken into account considering various disturbances. Integer Linear Programming technique with equality and inequality constraints is used to find out the optimal placement set with key buses identified from transient stability analysis. Results on IEEE-14 bus system are presented to illustrate the proposed approach.

Amplitude death in complex networks induced by environment

We present a mechanism for amplitude death in coupled nonlinear dynamical systems on a complex network having interactions with a common environment like external system. We develop a general stability analysis that is valid for any network topology and obtain the threshold values of coupling constants for the onset of amplitude death. An important outcome of our study is a universal relation between the critical coupling strength and the largest nonzero eigenvalue of the coupling matrix. Our results are fully supported by the detailed numerical analysis for different network topologies.

Safety Assessment of a Masonry Arch Bridge: Field Testing and Simulations

The safety of an in-service brick arch railway bridge is assessed through field testing and finite-element analysis. Different loading test train configurations have been used in the field testing. The response of the bridge in terms of displacements, strains, and accelerations is measured under the ambient and design train traffic loading conditions. Nonlinear fracture mechanics-based finite-element analyses are performed to assess the margin of safety. A parametric study is done to study the effects of tensile strength on the progress of cracking in the arch. Furthermore, a stability analysis to assess collapse of the arch caused by lateral movement at the springing of one of the abutments that is elastically supported is carried out. The margin of safety with respect to cracking and stability failure is computed. Conclusions are drawn with some remarks on the state of the bridge within the framework of the information available and inferred information. DOI: 10.1061/(ASCE)BE.1943-5592.0000338. (C) 2013 American Society of Civil Engineers.

Evaluation of strength of component-laminates in strip-based mechanisms

This paper deals with the evaluation of the component-laminate load-carrying capacity, i.e., to calculate the loads that cause the failure of the individual layers and the component-laminate as a whole in four-bar mechanism. The component-laminate load-carrying capacity is evaluated using the Tsai-Wu-Hahn failure criterion for various lay-ups. The reserve factor of each ply in the component-laminate is calculated by using the maximum resultant force and the maximum resultant moment occurring at different time steps at the joints of the mechanism. Here, all component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (strip-like beam). Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional nonlinearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the non-linear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the non-linear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict more quickly and accurately than would otherwise be possible. Local 3-D stress, strain and displacement fields for representative sections in the component-bars are recovered, based on the stress resultants from the 1-D global beam analysis. A numerical example is presented which illustrates the failure of each component-laminate and the mechanism as a whole.

A multifold reduction in the transition Reynolds number, and ultra-fast mixing, in a micro-channel due to a dynamical instability induced by a soft wall

A dynamical instability is observed in experimental studies on micro-channels of rectangular cross-section with smallest dimension 100 and 160 mu m in which one of the walls is made of soft gel. There is a spontaneous transition from an ordered, laminar flow to a chaotic and highly mixed flow state when the Reynolds number increases beyond a critical value. The critical Reynolds number, which decreases as the elasticity modulus of the soft wall is reduced, is as low as 200 for the softest wall used here (in contrast to 1200 for a rigid-walled channel) The instability onset is observed by the breakup of a dye-stream introduced in the centre of the micro-channel, as well as the onset of wall oscillations due to laser scattering from fluorescent beads embedded in the wall of the channel. The mixing time across a channel of width 1.5 mm, measured by dye-stream and outlet conductance experiments, is smaller by a factor of 10(5) than that for a laminar flow. The increased mixing rate comes at very little cost, because the pressure drop (energy requirement to drive the flow) increases continuously and modestly at transition. The deformed shape is reconstructed numerically, and computational fluid dynamics (CFD) simulations are carried out to obtain the pressure gradient and the velocity fields for different flow rates. The pressure difference across the channel predicted by simulations is in agreement with the experiments (within experimental errors) for flow rates where the dye stream is laminar, but the experimental pressure difference is higher than the simulation prediction after dye-stream breakup. A linear stability analysis is carried out using the parallel-flow approximation, in which the wall is modelled as a neo-Hookean elastic solid, and the simulation results for the mean velocity and pressure gradient from the CFD simulations are used as inputs. The stability analysis accurately predicts the Reynolds number (based on flow rate) at which an instability is observed in the dye stream, and it also predicts that the instability first takes place at the downstream converging section of the channel, and not at the upstream diverging section. The stability analysis also indicates that the destabilization is due to the modification of the flow and the local pressure gradient due to the wall deformation; if we assume a parabolic velocity profile with the pressure gradient given by the plane Poiseuille law, the flow is always found to be stable.

Suppression of gravitational instabilities by dominant dark matter halo in low-surface-brightness galaxies

The low-surface-brightness galaxies are gas rich and yet have a low star formation rate; this is a well-known puzzle. The spiral features in these galaxies are weak and difficult to trace, although this aspect has not been studied much. These galaxies are known to be dominated by the dark matter halo from the innermost regions. Here, we do a stability analysis for the galactic disc of UGC 7321, a low-surface-brightness, superthin galaxy, for which the various observational input parameters are available. We show that the disc is stable against local, linear axisymmetric and non-axisymmetric perturbations. The Toomre Q parameter values are found to be large (>> 1) mainly due to the low disc surface density, and the high rotation velocity resulting due to the dominant dark matter halo, which could explain the observed low star formation rate. For the stars-alone case, the disc shows finite swing amplification but the addition of dark matter halo suppresses that amplification almost completely. Even the inclusion of the low-dispersion gas which constitutes a high disc mass fraction does not help in causing swing amplification. This can explain why these galaxies do not show strong spiral features. Thus, the dynamical effect of a halo that is dominant from inner regions can naturally explain why star formation and spiral features are largely suppressed in low-surface-brightness galaxies, making these different from the high-surface-brightness galaxies.

Homeostasis of functional maps in active dendrites emerges in the absence of individual channelostasis

The maintenance of ion channel homeostasis, or channelostasis, is a complex puzzle in neurons with extensive dendritic arborization, encompassing a combinatorial diversity of proteins that encode these channels and their auxiliary subunits, their localization profiles, and associated signaling machinery. Despite this, neurons exhibit amazingly stereotypic, topographically continuous maps of several functional properties along their active dendritic arbor. Here, we asked whether the membrane composition of neurons, at the level of individual ion channels, is constrained by this structural requirement of sustaining several functional maps along the same topograph. We performed global sensitivity analysis on morphologically realistic conductance-based models of hippocampal pyramidal neurons that coexpressed six well-characterized functional maps along their trunk. We generated randomized models by varying 32 underlying parameters and constrained these models with quantitative experimental measurements from the soma and dendrites of hippocampal pyramidal neurons. Analyzing valid models that satisfied experimental constraints on all six functional maps, we found topographically analogous functional maps to emerge from disparate model parameters with weak pairwise correlations between parameters. Finally, we derived a methodology to assess the contribution of individual channel conductances to the various functional measurements, using virtual knockout simulations on the valid model population. We found that the virtual knockout of individual channels resulted in variable, measurement and location-specific impacts across the population. Our results suggest collective channelostasis as a mechanism behind the robust emergence of analogous functional maps and have significant ramifications for the localization and targeting of ion channels and enzymes that regulate neural coding and homeostasis.