Numerical modeling of the KdV random wave field

The evolution of the initially random wave field with a Gaussian spectrum shape is studied numerically within the Korteweg–de Vries (KdV) equation. The properties of the KdV random wave field are analyzed: transition to a steady state, equilibrium spectra, statistical moments of a random wave field, and the distribution functions of the wave amplitudes. Numerical simulations are performed for different Ursell parameters and spectrum width. It is shown that the wave field relaxes to the stationary state (in statistical sense) with the almost uniform energy distribution in low frequency range (Rayleigh–Jeans spectrum). The wave field statistics differs from the Gaussian one. The growing of the positive skewness and non-monotonic behavior of the kurtosis with increase of the Ursell parameter are obtained. The probability of a large amplitude wave formation differs from the Rayleigh distribution.     

Light beams and radiation patterns in random media

The purpose of this paper is to apply a numerical technique in solving problems involving light beams in a random medium. The technique starts by generating numerically sample media with prescribed statistical properties. Rays are then traced in these sample media and the ray statistics are compiled. These statistics are utilized to find the mean square displacement and distribution of a beam. Problems on beam broadening and distortion of radiation patterns are considered. Fluctations in wave amplitude and phase are also investigated. When possible, the numerical results are compared with the analytical results and experimental results.It seems that the numerical technique has a potential to solve a great variety of problems. This is because it does not have severe restrictive conditions as those imposed on the analytic formulation. For example the technique is equally applicable when the irregularities are anisotropic, or when the background medium is inhomogeneous, or when the background is anistropic, or when there exist background wind. Some of these are discussed.     

Reflection and transmission by a random medium

From ultrasonic inspections of heterogeneous materials the ultrasonic velocity in a test specimen and the ultrasonic reflection and transmission coefficients can be found over a certain range of frequencies. This paper presents a theoretical analysis of the frequency dependence of the reflection and transmission coefficients caused by random fluctuations of the material parameters.The procedure of using the effective frequency-dependent material parameters of an unbounded random medium to determine the wave propagation in a bounded random medium is not applicable at a distance of the boundary of the order of a correlation length. From this interface effect a contribution to the frequency dependence of the reflection and transmission coefficients arises.A previous consideration assuming the correlation length of the fluctuations to be small compared with the wavelength, is extended to shorter waves. Moreover, the wave propagation in a random half-space whose average material parameter differs from that of the exterior is considered. The consequences of the fluctuations of the elastic constants are of special interest, because they do not go with ω² in the governing wave equation as the density does.The interface effect is also discussed in the static limit of ‘homogeneous deformation’.     

Application of an Itô-based approximation method to random vibration of a pinching hysteretic system

In this paper, an extension of the Cumulant-Neglect closure scheme is utilized for the random vibration analysis of a single degree of freedom system with a general pinching hysteresis restoring force. The hysteresis element used in the system model can simulate commonly observed forms of stiffness, strength and pinching degradations. The second order statistics of the system response to a stationary Gaussian white noise input are derived using an Itô-based approximation technique. The validity of these response statistics are then verified by comparing them to Monte Carlo simulation results. The numerical studies performed for different combinations of degradation parameters and excitation levels show that the response estimates obtained by this solution method are in good agreement with Monte Carlo simulation. These studies also indicate the applicability of this technique for response analysis of complicated forms of non-linearities.     

Stochastic thermoelastic problem of a functionally graded plate under random temperature load

This study attempts to derive the statistics of temperature and thermal stress in functionally graded material (FGM) plates exposed to random external temperatures. The thermomechanical properties of the FGM plates are assumed to vary arbitrarily only in the plate thickness direction. The external temperatures are expressed as random functions with respect to time. The transient temperature field in the FGM plate is determined by solving a nonhomogeneous heat conduction problem for a multilayered plate with linear nonhomogeneous thermal conductivity and different homogeneous heat capacity in each layer. The autocorrelations and power spectrum densities (PSDs) of temperature and thermal stress are derived analytically. These statistics for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304) are computed under the condition that the fluctuation in the external temperature can be considered as white noise or a stationary Markov process.     

Effect of damping and wave parameters on offshore structure under random excitation

This paper describes the linearized and nonlinear dynamic response of a tension leg platform (TLP) to random waves and current forces. The forcing term of the equation of motion is inherently nonlinear due to the nonlinear drag force. Two analysis procedures are used: nonlinear time domain analysis and linear frequency domain analysis. For the nonlinear analysis, the random wave particle velocities and accelerations are simulated for a given wave spectrum. The nonlinear equation of motion is then integrated directly to obtain the system response statistics. For the linear frequency domain analysis, the nonlinear drag force is linearized through an introduction of linearization coefficients. The main objective of this paper is to investigate the effect of the structural damping and wave parameters on both nonlinear and linear dynamic response of the TLP by parametric studies. The results of stochastic nonlinear and linear dynamic response of the TLP, with and without the presence of current, are presented and compared.     

多铁性扇环截面柱体中的波传播

论文基于线性磁电弹性理论,研究了具有扇环形截面的多铁性柱形波导中的弹性波传播问题.利用波动势函数法,解析推导获得波动特征方程,进而得到弥散关系.通过算例研究了波传播的关键特性,深入分析了弥散曲线、相速度曲线和截止频率变化情况.结果显示,波的相速度和截止频率非常依赖于波导结构的扇环截面半顶角、内外径比和层合界面的弱界面系数,对于给定材料的波导结构,这些参数也是控制其弥散特性的重要影响因素.值得指出的是,在柱面应力自由的边界条件下相速度曲线中存在独特的频率带隙,而这通常是在周期结构中才会出现.     

Y波导多功能集成光学器件可靠性增长

为消除Y波导多功能集成光学器件潜在的薄弱环节提高其可靠性,通过＂试验-分析-改进-再试验＂对器件进行可靠性增长。分析了环境应力对Y波导多功能集成光学器件的影响,并对器件的可靠性进行摸底,暴露出潜在的薄弱环节为耦合胶耐高温高湿性能较差和尾纤结构不合理,在此基础上经多次试验对比选择了耐高温高湿的耦合胶与合理的尾纤结构,改进措施落实到器件上以后再次进行了高温高湿试验以验证改进措施的有效性。试验结果表明改进后的Y波导多功能集成光学器件在85℃、85%RH的条件下贮存500h其性能不会下降,这相当于40℃、85%RH下的20780h,通过材料和结构上的改进使得可靠性提高到改进前的39.6倍。     

Equivalent linearization for systems subjected to non-stationary random excitation

The method of eauivalent linearization is applied to the general problem of the response of non-linear discrete systems to non-stationary random excitation. Conditions for minimum equation difference are determined which do not depend explicitly on lime but only on the instantaneous statistics of the response process. Using the equivalent linear parameters, a deterministic non-linear ordinary differential equation for the covariance matrix is derived. An example is given of a damped Duffing oscillator subjected to modulated white noise.     

Parametric random excitation of nonlinear coupled oscillators

The stochastic bifurcation and response statistics of nonlinear modal interaction under parametric random excitation are studied analytically, numerically and experimentally. Two basic definitions of stochastic bifurcation are first introduced. These are bifurcation in distribution and bifurcation in moments. bifurcation in moments is examined for the case of a coupled oscillator subjected to parametric filtered white noise. The center frequency of the excitation is selected to be close to either twice the first mode or second mode natural frequencies or the sum of the two. The stochastic bifurcation in moments is predicted using the Fokker-Planck equation together with gaussian and non-Gaussian closures and numerically using the Monte Carlo simulation. When one mode is parametrically excited it transfers energy to the other mode due to nonlinear modal interaction. The Gaussian closure solution gives close results to those predicted numerically only in regions well remote from bifurcation points. However, bifurcation points predicted by the non-Gaussian closure are in good agreement with those estimated by numerical simulation. Depending on the excitation level, the probability density of the excited mode is strongly non-Gaussian and exhibits multi-maxima as predicted by Monte Carlo simulation. Experimental tests are carried out at relatively low excitation levels. In the neighborhood of stochastic bifurcation in mean square the measured results exhibit different regimes of response characteristics including zero motion and occasional small random motion regimes. These two regimes are characterized by the phenomenon of on-off intermittency. Both regimes overlap and thus it is difficult to locate experimentally the bifurcation point.     

Non-stationary random vibration of a linear structure

The non-stationary random vibration of a lightly damped linear structure subjected to white noise is considered. It is shown that the probability density function of the amplitude of the structural response can be approximated by a Rayleigh distribution. Analytical formulae for the time dependent statistics of the amplitude are presented. The analytical results are compared with data obtained by a numerical simulation     

Stochastic non-linear flutter of a panel subjected to random in-plane forces

—An analysis of non-linear flutter of a simply-supported panel exposed to supersonic gas flow and random in-plane forces is presented for two- and three-mode interactions. A first order quasi-steady state aerodynamic piston theory is used to model the aerodynamic loading. The Fokker-Planck equation is used to derive a general moment equation for two- and three-mode interactions. For stability analysis the moment equation is consistent and the mean square stability boundaries of the equilibrium are obtained in terms of the system parameters. The stability boundaries reveal common features to those predicted by the deterministic theory of panel nutter. For the non-linear response the moment equation is found inconsistent and a cumulant-neglect closure is used by setting cumulants of fifth and sixth orders to zero. This first order non-Gaussian closure is carried out to solve for the response statistics in terms of the air-to-plate mass ratio, aerodynamic pressure, modal damping, and in-plane random force spectral density. It is found that the non-Gaussian solution yields higher levels for the response statistics than those obtained by the Gaussian solution. The inclusion of more modes results in a reduction of the response levels and expands the stability region.     

Propagation and localization of Rossby waves over random topography (two-layer model)

We study propagation of Rossby waves over randomly stratified bottom topography in a two-layer system. The problem is reduced to coupled stochastic wave equations in latitudinal variable with suitable boundary conditions. When the two-dimensional system is broken into its basic barotropic and baroclinic modes, each one is scattered by the medium, and generates the like-wise and the cross-wise components. We study statistics of the reflection and transmission coefficients for both types of generated modes, and show localization for the like-wise components, and propagation for the cross-wise components. The localization lengths and the transmission rates are estimated in terms of the basic parameters of the system and the correlation-function of topographic fluctuations.     

非高斯随机激励下滞迟系统响应与首次穿越失效

针对由有界噪声、泊松白噪声和高斯白噪声共同构成的非高斯随机激励,通过Monte Carlo数值模拟方法研究了此激励作用下双线性滞迟系统和Bouc-Wen滞迟系统这两类经典滞迟系统的稳态响应与首次穿越失效时间。一方面,分析了有界噪声和泊松白噪声这两种分别具有连续样本函数和非连续样本函数的非高斯随机激励,在不同激励参数条件下对双线性滞迟系统和Bouc-Wen滞迟系统的稳态响应概率密度、首次穿越失效时间概率密度及其均值的不同影响;另一方面,揭示了在这类非高斯随机激励荷载作用下,双线性滞迟系统的首次穿越失效时间概率密度将出现与Bouc-Wen滞迟系统的单峰首次穿越失效时间概率密度截然不同的双峰形式。     

ON DOUBLE PEAK PROBABILITY DENSITY FUNCTIONS OF DUFFING OSCILLATOR TO COMBINED DETERMINISTIC AND RANDOM EXCITATIONS

The principal resonance of Duffing oscillator to combined deterministic and random external excitation was investigated. The random excitation was taken to be white noise or harmonic with separable random amplitude and phase. The method of multiple scales was used to determine the equations of modulation of amplitude and phase. The one peak probability density function of each of the two stable stationary solutions was calculated by the linearization method. These two one-peak-density functions were combined using the probability of realization of the two stable stationary solutions to obtain the double peak probability density function. The theoretical analysis are verified by numerical results.     

Radiation and attenuation of low- and high-frequency waves in a random medium

The wave field radiated by a point source in a weakly inhomogeneous, weakly dissipative, one-dimensional random medium is considered. Approximate expressions for the mean intensity and mean energy flux as a function of propagation range, based on general results obtained previously, are derived for the limiting cases of low and high frequencies. These expressions show that, in both cases, the mean intensity and mean energy flux decrease more rapidly with range than would be the case in the absence of randomness. This more rapid rate of decrease with range can be interpreted as an excess attenuation. These results are similar to those obtained previously for the intermediate-frequency case. The excess-attenuation results are found to be consistent with observations.     

Mapping closure approximation for reactive scalars in random flows

The Mapping Closure Approximation (MCA) approach is developed to describe the statistics of both conserved and reactive scalars in random flows. The statistics include Probability Density Function (PDF), Conditional Dissipation Rate (CDR) and Conditional Laplacian (CL). The statistical quantities are calculated using the MCA and compared with the results of the Direct Numerical Simulation (DNS). The results obtained from the MCA are in agreement with those from the DNS. It is shown that the MCA approach can predict the statistics of reactive scalars in random flows. The project supported by the National Committee of Science and Technology, China, under the Special Funds for Major Basic Research Project (G2000077305 and G1999032801), and the National Natural Science Foundation of China (10325211)     

Stochastic optimal active control of a half-car nonlinear suspension under random road excitation

The random response analysis and the stochastic optimal active control of a half-car model with nonlinear suspension stiffness and damping traversing a rough road are studied in this paper. The road roughness height is modeled as the output of a first-order linear filter to Gaussian white noise. Considering the hysteretic nonlinear stiffness and the square damping of the vehicle model, the response statistics of the nonlinear suspension with active control are obtained by using the equivalent linearization method. The performance indexes of the active suspension are evaluated and compared with those of the corresponding passive suspension. It is found that the nonlinear active suspension gives a better vehicle performances like ride comfort, suspension stroke and overall performance. Finally, the theoretical results are verified through Monte Carlo simulation.     

Generation of turbulent inflow and initial conditions based on multi‐correlated random fields

Direct or large eddy simulation of a turbulent flow field is strongly influenced by its initial or inflow boundary condition. This paper presents a new stochastic approach to generate an artificial turbulent velocity field for initial or inflow boundary condition based on digital filtering. Each velocity component of the artificial turbulent velocity field is generated by linear combination of individual uncorrelated random fields. These uncorrelated random fields are obtained by filtering random white‐noise fields. Using common elements in these linear combinations results in multi‐correlation among different velocity components. The generated velocity field reproduces locally desired Reynolds stress components and integral length scales including cross‐integral length scales. The method appears to be simple, flexible and more accurate in comparison with previously developed methods. The accuracy and performance of the method are demonstrated by numerical simulation of a homogeneous turbulent shear flow with high and low shear rates. To assess the accuracy and performance of the method, simulation results are compared with a reference simulation. Copyright © 2007 John Wiley & Sons, Ltd.     

Higher order linearization in non-linear random vibration

In this paper a higher order linearization method for analyzing non-linear random vibration problems is presented. The non-linear terms of the given equation are replaced by unknown linear terms. These are in turn described by extra non-linear differential equations. The combined system of equations is then linearized to arrive at a higher degree-of-freedom equation for the original system. The method is illustrated by considering the Duffing oscillator under white noise input. The equivalent two d.o.f linear system is derived by the present method. Numerical results on steady state variance and PSD functions are obtained. These are found to be better than the simple linearization results.