Global response sensitivity analysis using probability distance measures and generalization of Sobol's analysis

The study introduces two new alternatives for global response sensitivity analysis based on the application of the L-2-norm and Hellinger's metric for measuring distance between two probabilistic models. Both the procedures are shown to be capable of treating dependent non-Gaussian random variable models for the input variables. The sensitivity indices obtained based on the L2-norm involve second order moments of the response, and, when applied for the case of independent and identically distributed sequence of input random variables, it is shown to be related to the classical Sobol's response sensitivity indices. The analysis based on Hellinger's metric addresses variability across entire range or segments of the response probability density function. The measure is shown to be conceptually a more satisfying alternative to the Kullback-Leibler divergence based analysis which has been reported in the existing literature. Other issues addressed in the study cover Monte Carlo simulation based methods for computing the sensitivity indices and sensitivity analysis with respect to grouped variables. Illustrative examples consist of studies on global sensitivity analysis of natural frequencies of a random multi-degree of freedom system, response of a nonlinear frame, and safety margin associated with a nonlinear performance function. (C) 2015 Elsevier Ltd. All rights reserved.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Estimation of nonlinearities from pseudodynamic and dynamic responses of bridge structures using the Delay Vector Variance method

Analysis of the variability in the responses of large structural systems and quantification of their linearity or nonlinearity as a potential non-invasive means of structural system assessment from output-only condition remains a challenging problem. In this study, the Delay Vector Variance (DVV) method is used for full scale testing of both pseudo-dynamic and dynamic responses of two bridges, in order to study the degree of nonlinearity of their measured response signals. The DVV detects the presence of determinism and nonlinearity in a time series and is based upon the examination of local predictability of a signal. The pseudo-dynamic data is obtained from a concrete bridge during repair while the dynamic data is obtained from a steel railway bridge traversed by a train. We show that DVV is promising as a marker in establishing the degree to which a change in the signal nonlinearity reflects the change in the real behaviour of a structure. It is also useful in establishing the sensitivity of instruments or sensors deployed to monitor such changes. (C) 2015 Elsevier B.V. All rights reserved.

Response Dynamics of Recirculation Structures in Coaxial Nonpremixed Swirl-Stabilized Flames Subjected to Acoustic Forcing

This paper reports the time-mean and phase-locked response of nonreacting as well as reacting flow field in a coaxial swirling jet/flame (nonpremixed). Two distinct swirl intensities plus two different central pipe flow rates at each swirl setting are investigated. The maximum response is observed at the 105 Hz mode in the range of excitation frequencies (0-315 Hz). The flow/flame exhibited minimal response beyond 300 Hz. It is seen that the aspect ratio change of inner recirculation zone (IRZ) under nonreacting conditions (at responsive modes) manifests as a corresponding increase in the time-mean flame aspect ratio. This is corroborated by similar to 25% decrease in the IRZ transverse width in both flame and cold flow states. In addition, 105 Hz excited states are found to shed high energy regions (eddies) asymmetrically when compared to dormant 315 Hz pulsing frequency. The kinetic energy (KE) of the flow field is subsequently reduced due to acoustic excitation and a corresponding increase (similar to O (1)) in fluctuation intensity is witnessed. The lower swirl intensity case is found to be more responsive than the high swirl case as in the former flow state the resistance offered by IRZ to incoming acoustic perturbations is lower due to inherently low inertia. Next, the phase-locked analysis of flow and flame structure is employed to further investigate the phase dependence of flow/flame response. It is found that the asymmetric shifting of IRZ mainly results at 270 deg acoustic forcing. The 90 deg phase angle forcing is observed to convect the IRZ farther downstream in both swirl cases as compared to other phase angles. The present work aims primarily at providing a fluid dynamic view point to the observed nonpremixed flame response without considering the confinement effects.

Fabrication and Nonlinear Thermomechanical Analysis of SU8 Thermal Actuator

SU8-based micromechanical structures are widely used as thermal actuators in the development of compliant micromanipulation tools. This paper reports the design, nonlinear thermomechanical analysis, fabrication, and thermal actuation of SU8 actuators. The thermomechanical analysis of the actuator incorporates nonlinear temperature-dependent properties of SU8 polymer to accurately model its thermal response during actuation. The designed SU8 thermal actuators are fabricated using surface micromachining techniques and the electrical interconnects are made to them using flip-chip bonding. The issues due to thermal stress during fabrication are discussed and a novel strategy is proposed to release the thermal stress in the fabricated actuators. Subsequent characterization of the actuator using an optical profilometer reveals excellent thermal response, good repeatability, and low hysteresis. The average deflection is similar to 8.5 mu m for an actuation current of similar to 5 mA. The experimentally obtained deflection profile and the tip deflection at different currents are both shown to be in good agreement with the predictions of the nonlinear thermomechanical model. This underscores the need to consider nonlinearities when modeling the response of SU8 thermal actuators. 2015-0087]

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.

Material response and failure mechanism of unidirectional metal matrix composites under impulsive shear loading

The material response and failure mechanism of unidirectional metal matrix composite under impulsive shear loading are investigated in this paper. Both experimental and analytical studies were performed. The shear strength of unidirectional C-f/A356.0 composite and A356.0 aluminum alloy at high strain rate were measured with a modified split Hopkinson torsional bar technique. The results indicated that the carbon fibers did not improve the shear strength of aluminum matrix if the fiber orientation aligned with the shear loading axis. The microscopic inspection of the fractured surface showed a multi-scale zigzag feature which implied a complicated shear failure mechanism in the composite. In addition to testing, the micromechanical stress field in the composite was analyzed by the generalized Eshelby equivalent method (GEEM). The influence of cracking in matrix on the micromechanical stress field was investigated as well. The results showed that the stress distribution in the composite is quite nonhomogeneous and very high shear stress concentrations are found in some regions in the matrix. The high shear stress concentration in the matrix induces tensile cracking at 45 degrees to the shear direction. This in turn aggravates the stress concentration at the fiber/matrix interface and finally leads to a catastrophic failure in the composite. From the correlation between the analysis and experimental results, the shear failure mechanism of unidirectional C-f/A356.0 composite can be elucidated qualitatively.

Identification of Nonlinear Systems Through Time-Frequency Filtering Technique

In the previous paper, a class of nonlinear system is mapped to a so-called skeleton linear model (SLM) based on the joint time-frequency analysis method. Behavior of the nonlinear system may be indicated quantitatively by the variance of the coefficients of SLM versus its response. Using this model we propose an identification method for nonlinear systems based on nonstationary vibration data in this paper. The key technique in the identification procedure is a time-frequency filtering method by which solution of the SLM is extracted from the response data of the corresponding nonlinear system. Two time-frequency filtering methods are discussed here. One is based on the quadratic time-frequency distribution and its inverse transform, the other is based on the quadratic time-frequency distribution and the wavelet transform. Both numerical examples and an experimental application are given to illustrate the validity of the technique.

Numerical Analyses of Discrete Gust Response for an Aircraft

Based on Navier-Stokes equations and structural and flight dynamic equations of motion, dynamic responses in vertical discrete gust flow perturbation are investigated for a supersonic transport model. A tightly coupled method was developed by subiterations between aerodynamic equations and dynamic equations of motion. First, under the assumption of rigid-body and single freedom of motion in the vertical plunging, the results of a direct-coupling method are compared with the results of quasi-steady model method. Then, gust responses for the one-minus-cosine gust profile arc analyzed with two freedoms of motion in plunging and pitching for the airplane configurations with and without the consideration of structural deformation.

Impact response and failure mechanisms of ultrahigh modulus polyethylene fiber composites and polyethylene fiber-carbon fiber hybrid composites

The impact response and failure mechanisms of ultrahigh modulus polyethylene (UHMPE) fiber composites and UHMPE fiber-carbon fiber hybrid composites have been investigated. Charpy impact, drop weight impact and high strain rate impact experiments have been performed in order to study the impact resistance, notch sensitivity, strain rate sensitivity and hybrid effects. Results obtained from dynamic and quasi-static measurements have been compared. Because of the ductility of UHMPE fibers, the impact energy absorption of UHMPE fiber composites is very high, thereby leading to excellent damage tolerance. By hybridizing with UHMPE fibers, the impact properties of carbon fiber composites can be greatly improved. The impact and shock failure mechanisms of these composites are discussed.

A dynamic view of the spread and intracellular distribution of Salmonella enterica.

The events that determine the dynamics of proliferation, spread and distribution of microbial pathogens within their hosts are surprisingly heterogeneous and poorly understood. We contend that understanding these phenomena at a sophisticated level with the help of mathematical models is a prerequisite for the development of truly novel, targeted preventative measures and drug regimes. We describe here recent studies of Salmonella enterica infections in mice which suggest that bacteria resist the antimicrobial environment inside host cells and spread to new sites, where infection foci develop, and thus avoid local escalation of the adaptive immune response. We further describe implications for our understanding of the pathogenic mechanism inside the host.

Time-Frequency Analysis of Nonlinear Systems: The Skeleton Linear Model and the Skeleton Curves

The joint time-frequency analysis method is adopted to study the nonlinear behavior varying with the instantaneous response for a class of S.D.O.F nonlinear system. A time-frequency masking operator, together with the conception of effective time-frequency region of the asymptotic signal are defined here. Based on these mathematical foundations, a so-called skeleton linear model (SLM) is constructed which has similar nonlinear characteristics with the nonlinear system. Two skeleton curves are deduced which can indicate the stiffness and damping in the nonlinear system. The relationship between the SLM and the nonlinear system, both parameters and solutions, is clarified. Based on this work a new identification technique of nonlinear systems using the nonstationary vibration data will be proposed through time-frequency filtering technique and wavelet transform in the following paper.

Modelling nonlinear vehicle dynamics with neural networks

The nonlinear modelling ability of neural networks has been widely recognised as an effective tool to identify and control dynamic systems, with applications including nonlinear vehicle dynamics which this paper focuses on using multi-layer perceptron networks. Existing neural network literature does not detail some of the factors which effect neural network nonlinear modelling ability. This paper investigates into and concludes on required network size, structure and initial weights, considering results for networks of converged weights. The paper also presents an online training method and an error measure representing the network's parallel modelling ability over a range of operating conditions. Copyright © 2010 Inderscience Enterprises Ltd.

A Nonlinear Threshold Model Applied to Spallation Analysis

Spallation in heterogeneous media is a complex, dynamic process. Generally speaking, the spallation process is relevant to multiple scales and the diversity and coupling of physics at different scales present two fundamental difficulties for spallation modeling and simulation. More importantly, these difficulties can be greatly enhanced by the disordered heterogeneity on multi-scales. In this paper, a driven nonlinear threshold model for damage evolution in heterogeneous materials is presented and a trans-scale formulation of damage evolution is obtained. The damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that some dimensionless numbers govern the whole process of deformation and damage evolution. The effects of heterogeneity in terms of Weibull modulus on damage evolution in spallation process are also investigated.