重力场匹配导航的重力图生成技术

通过对当前用于海洋重力场格网插值的四种常用算法(距离倒数加权法、Kriging法、径向基函数法和改进的二次曲面Shepard方法)进行分析,以相对规则、不同密度的两组数据作为基准数据进行了插值比较。实验结果表明,基于改进的Shepard插值算法相对于其它三种算法具有速度快、精度高的优点,比较符合当前海洋重力数据获取的现状,是进行高精度重力图生成的有效方法。     

Advanced statistics to improve the physical interpretation of atomization processes

This paper reports an analysis of the physics of atomization processes using advanced statistical tools. Namely, finite mixtures of probability density functions, which best fitting is found using a Bayesian approach based on a Markov chain Monte Carlo (MCMC) algorithm. This approach takes into account eventual multimodality and heterogeneities in drop size distributions. Therefore, it provides information about the complete probability density function of multimodal drop size distributions and allows the identification of subgroups in the heterogeneous data. This allows improving the physical interpretation of atomization processes. Moreover, it also overcomes the limitations induced by analyzing the spray droplets characteristics through moments alone, particularly, the hindering of different natures of droplet formation. Finally, the method is applied to physically interpret a case-study based on multijet atomization processes.     

Exact stationary probability density functions for non-linear systems under Poisson white noise excitation

The paper presents exact stationary probability density functions for systems under Poisson white noise excitation. Two different solution methods are outlined. In the first one, a class of non-linear systems is determined whose state vector is a memoryless transformation of the state vector of a linear system. The second method considers the generalized Fokker-Planck (Kolmogorov-forward) equation. Non-linear system functions are identified such that the stationary solution of the system admits a prescribed stationary probability density function. Both methods make use of the stochastic integro-differential equations approach. This approach seems to have some computational advantages for the determination of exact stationary probability density functions when compared to the stochastic differential equations approach.     

采用重要抽样法的结构动力可靠度计算

首次对比分析了结构动力可靠度计算的三种重要抽样法,并对部分方法进行了补充修正.单元失效域法补充了依据随机教决定抽样区间的产生方法,根据单元失效域的条件概率和权重系数给出重要抽样密度函教.方差放大系数法直接通过激励过程的特性给出重要抽样密度函数的具体表达式.功率谱法的重要抽样密度函数仅为激励幅值的函数,根据结构反应的功率谱密度增大激励幅值的方差,建议幅值样本值的联合概率密度函数可表示为幅值样本值分量的概率密度函数的连乘形式.结果表明:对于线性体系三种方法的计算效率均比Monte-Carlo法有显著提高,而单元失效域法的计算效率又比另两种方法高.     

基于优化算法的串联体系可靠度分析

结构体系的失效概率数学上可以表示为结构体系失效域上联合概率密度函数的积分，一般情况下很难直接积分求解。近几十年来，结构体系可靠度分析一直是可靠度领域的一个研究热点，人们提出许多方法，如：Monte—Carlo法、重要性抽样法与界限法和概率网络估算技术等，这些算法在求解精度、计算效率、收敛性和易使用性等方面是不同的。本文采用优化算法(改进的可行方向法、序列线性规划和序列二次规划法)进行串联体系可靠度分析，并且与其他算法(HL—RF法、Monte—Carlo法和重要性抽样法)的结果以及一些精确解进行了比较。结果表明，相对于其他算法，基于优化算法的可靠度分析适用性广，在收敛性和健实性等方面具有明显的优势。     

A quadratic hierarchical Bayesian dynamic prediction model for infrastructure maintenance

The main objective of this study is the development of a correlation model in dynamic Bayesian belief networks (DBBNs) followed by an inverse economic analysis. This is based on a quadratic hierarchical Bayesian inference prediction method using Markov chain Monte Carlo simulations. The developed model is implemented to predict the future degradation and maintenance budget for a suspension bridge system. Bayesian inference is applied to find the posterior probability density function of the source parameters (damage indices and serviceability), given 10 years maintenance data. The simulated risk prediction under decreased serviceability conditions gives posterior distributions based on a prior distribution and likelihood data updated from annual maintenance tasks. Compared with a conventional linear prediction model, the proposed quadratic model provides highly improved convergence and closeness to the measured data. Finally, the developed inverse DBBN analysis method allows forecasts of future performance and the financial management of complex infrastructures by providing the sensitivity of serviceability and risky factors to the maintenance budgets of structural components and the overall system.     

Correction of wandering smoothing effects on static measurements of a wing-tip vortex

Wandering is a typical feature of wing-tip vortices and it consists in random fluctuations of the vortex core. Consequently, vortices measured by static measuring techniques appear to be more diffuse than in reality, so that a correction method is needed. In the present paper statistical simulations of the wandering of a Lamb-Oseen vortex are first performed by representing the vortex core locations through bi-variate normal probability density functions. It is found that wandering amplitudes smaller than 60% of the core radius are well predicted by using the ratio between the RMS value of the mean cross-velocity and its slope measured at the mean vortex center. Furthermore, the principal axes of wandering can be accurately evaluated from the opposite of the cross-correlation coefficient between the spanwise and the normal velocities measured at the mean vortex center. The correction of the wandering smoothing effects is then carried out through four different algorithms that perform the deconvolution of the mean velocity field with the probability density function that represents the wandering. The corrections performed are very accurate for the simulations with wandering amplitudes smaller than 60% of the core radius, whereas errors become larger with increasing wandering amplitudes. Subsequently, the whole procedure to evaluate wandering and to correct the mean velocity field is applied to static measurements, carried out with a fast-response five-hole pressure probe, of a tip vortex generated from a NACA 0012 half-wing model. It is found that the wandering is predominantly in the upward-outboard to downward-inboard direction. Furthermore, the wandering amplitude grows with increasing streamwise distance from the wing, whereas it decreases with increasing angle of attack and free-stream velocity.     

Nonlinear dynamics and performance analysis of modified snap-through vibration energy harvester with time-varying potential function

Vibration energy harvesting has emerged as a promising method to harvest energy for small-scale applications. Enhancing the performance of a vibration energy harvester(VEH) incorporating nonlinear techniques, for example, the snap-through VEH with geometric non-linearity, has gained attention in recent years. A conventional snap-through VEH is a bi-stable system with a time-invariant potential function, which was investigated extensively in the past. In this work, a modified snap-through VEH wit...     

New approximate solution for N-point correlation functions for heterogeneous materials

Statistical N-point correlation functions are used for calculating properties of heterogeneous systems. The strength and the main advantage of the statistical continuum approach is the direct link to statistical information of microstructure. Two-point correlation functions are the lowest order of correlation functions that can describe the morphology and the microstructure-properties relationship. Experimentally, statistical pair correlation functions are obtained using SEM or small x-ray scattering techniques. Higher order correlation functions must be calculated or measured to increase the precision of the statistical continuum approach. To achieve this aim a new approximation methodology is utilized to obtain N-point correlation functions for non-FGM (functional graded materials) heterogeneous microstructures. Conditional probability functions are used to formulate the proposed theoretical approximation. In this approximation, weight functions are used to connect subsets of (N?1)-point correlation functions to estimate the full set of N-point correlation function. For the approximation of three and four point correlation functions, simple weight functions have been introduced. The results from this new approximation, for three-point probability functions, are compared to the real probability functions calculated from a computer generated three-phase reconstructed microstructure in three-dimensional space. This three-dimensional reconstruction was based on an experimental two-dimensional microstructure (SEM image) of a three-phase material. This comparison proves that our new comprehensive approximation is capable of describing higher order statistical correlation functions with the needed accuracy.     

A finite element solution of the time-dependent incompressible Navier–Stokes equations using a modified velocity correction method

A finite element solution of the two-dimensional incompressible Navier–Stokes equations has been developed. The present method is a modified velocity correction approach. First an intermediate velocity is calculated, and then this is corrected by the pressure gradient which is the solution of a Poisson equation derived from the continuity equation. The novelty, in this paper, is that a second-order Runge–Kutta method for time integration has been used. Discretization in space is carried out by the Galerkin weighted residual method. The solution is in terms of primitive variables, which are approximated by polynomial basis functions defined on three-noded, isoparametric triangular elements. To demonstrate the present method, two examples are provided. Results from the first example, the driven cavity flow problem, are compared with previous works. Results from the second example, uniform flow past a cylinder, are compared with experimental data.     

Modeling of anisotropic polydispersed composites by progressive micromechanical approach

A progressive micromechanical method is presented in order to predict the elastic constants of polydispersed composites including multi-directional or randomly oriented reinforcement particles. Heterogeneities of various types are introduced into the matrices in a gradual manner. At each step, the Mori-Tanaka method is used to obtain the stiffness tensor of the intermediate medium used as a matrix of the following step. The proposed method is capable of introducing any kind of heterogeneities based on their dimensions, orientations, mechanical properties, and volume fractions to the matrix. Furthermore, suitable probability density functions can be defined for physical and structural parameters of the composite, including the level of the filler-matrix interfacial bonding, the aspect ratio, and the orientation of reinforcement particles. The efficiency of the iterative approach and the convergence of the solution are studied by computing the stiffness tensors of unidirectional and bidirectional particulate composites. The results of the present study are also compared with the literature data for a randomly oriented particulate composite.     

轿车钢板表面形貌特征参数的提取

通过对典型轿车钢板表面形貌的分析,以形貌与其摩擦学性能的关系为出发点,对粗糙表面、粗糙峰和粗糙谷的分布特性进行了统计分析。结果表明：粗糙表面、粗糙峰和粗糙谷数的统计分布特性与润滑及接触等性能密切相关;针对轿车钢板表面、粗糙峰和粗糙谷数在高度方向概率分布曲线的非Gauss特性,提出采用Weibull概率密度函数的非线性曲线拟合对分布特征进行分析处理。所得拟合参数具有非常明确的物理意义,所此可以采用少量参数来较全面地描述其对应的轿车钢板表面形貌特征分布。     

Enhancing water level prediction through model residual correction based on Chaos theory and Kriging

Hydrodynamic models based on the physical processes are indispensable tools for predicting water levels in ocean environment. Nonetheless, their accuracies are limited by various factors such as simplifying assumptions, complex ocean bathymetry, and so on. Residual correction, as one of the data assimilation techniques, can extract information from observation and assimilate it into a numerical model to correct the model output directly. Such correction is often performed in two steps: prediction of the model residuals at measured stations followed by spatial distribution at non‐measured locations. For long‐term residual forecast, the accuracy of prediction usually deteriorates with the forecast horizon. In addition to the residual correction at measurement locations, in this paper, we address the critical question as to how to effectively update outputs for computational points without measurements. We develop a hybrid data assimilation procedure, which combines a modified local model (MLM) and an approximated ordinary kriging (AOK). This technique improves the forecasts over a long horizon over the entire computational domain. Using the proposed residual correction technique, the hybrid procedure is examined on a case study of Singapore Regional Model for correcting the water level outputs at locations with and without measurements. In order to provide a comparison, the analysis is carried out throughout prediction horizon of model residuals, tidal residuals, and sea level anomaly, respectively. The comparisons show that the proposed method can successfully assimilate and forecast all variables. Results indicate that resulting prediction accuracy can be significantly improved for all locations of interest independently of the forecast horizon. Copyright © 2014 John Wiley & Sons, Ltd.     

模态参数不确定性分析的贝叶斯方法研究

在结构损伤诊断和参数识别中,实测结构模态参数不可避免地存在误差。本文将模态参数视为随机变量,采用贝叶斯方法对模态参数的不确定性进行分析。分析中选用高斯联合概率密度函数作为先验密度函数,通过多次独立的模态参数测试,得到传递函数的条件概率密度函数和模态参数的后验估计表达式,再利用拉普拉斯渐近方法求解边缘概率密度函数,得到模态参数的最大后验估计。在钢筋混凝土框架结构的模态试验中,利用本文方法给出了结构模态参数的估计值,结果表明,本文方法具有良好的收敛性。     

基于显示连通贝叶斯网络的机场塔台地震易损性分析

将基于性能的多维易损性分析方法,结合显示连通贝叶斯网络,应用于机场塔台的多维易损性分析。考虑地震激励的不确定性,通过非线性时程分析获得结构响应数据;将塔台结构分为三个层次,每个层次按包含的层数分为相应的子层次。根据功能特性确定子层次的评价指标和极限状态,建立服从多元对数正态分布的概率地震需求模型;考虑各种极限状态之间的相关性,建立极限状态方程,确定失效域,通过蒙特卡洛法求得构件的超越概率;建立塔台结构的显示连通贝叶斯网络模型,利用层次分析法获得中间节点的条件概率表,利用MATLAB进行贝叶斯网络的推理计算,实现从单一层次的易损性到整体易损性的推理。     

Multi-parameter accelerated modified Newton-Raphson methods for rigid-plastic FE analysis

Applications of certain multi-parameter acceleration techniques used with the modified Newton-Raphson (mN-R) methods to solve the nonlinear equations arising from rigid-plastic finite element analysis are investigated. New modified multi-parameter techniques, developed from Crisfield's multi-parameter methods, are utilized to solve these nonlinear equations. The numerical performance of these techniques is compared with the standard Newton-Raphson method (sN-R), Crisfield's single parameter method (C1), Crisfield's two parameter method (C2) and Crisfield's three parameter method (C3). The new techniques do not involve additional residual force calculation and require little extra computational effort. In addition, they are more robust and efficient than other existing acceleration techniques.     

Investigation of bed shear stresses subject to external turbulence

The fluctuating bed shear stress has largely been investigated only for uniform channel flows and boundary layers. In practical engineering, the flow conditions are often modified due to the presence of various hydraulic structures. To what extent the modification affects the known characteristics of the bed shear stress is not clear. This study presents experimental results of the bed shear stress fluctuations, which are obviously subjected to external turbulence generated by superimposing artificial structures in the open channel flows. The statistical analysis of the measurements shows that the probability density function of the bed shear stress can be described by the lognormal function. Some associated relations concerning higher-order moments, skewness and kurtosis, which are derived from the lognormal function, are further compared with the experimental data. Physical implication of the skewed probability density distribution is finally discussed.     

A method for reconstruction of residual stress fields from measurements made in an incompatible region

A method is introduced by which the complete state of residual stress in an elastic body may be inferred from a limited set of experimental measurements. Two techniques for carrying out this reconstruction using finite element analysis are compared and it is shown that for exact reconstruction of the stress field via this method, the stress field must be measured over all eigenstrain-containing regions of the object. The effects of error and incompleteness in the measured part of the stress field on the subsequent analysis are investigated in a series of numerical experiments using synthetic measurement data based on the NeT TG1 round-robin weld specimen. It is hence shown that accurate residual stress field reconstruction is possible using measurement data of a quality achievable using current experimental techniques.     

Transitions in a Duffing oscillator excited by random noise

We investigate a Duffing oscillator driven by random noise which is assumed to be a harmonic function of the Wiener process. We show that the correlation time of the noise has a strong effect on the form of the response stationary probability density functions. It represents the so-called reentrance transitions, i.e. for the same noise intensity the probability density function has an identical modality for both the small and the large correlation time but a different modality for the moderate correlation time. The transitions are observed for both the single-well and twin-well potential case. A new approach is used to study the response probability density function. It is based on analysis of hyperbolic systems.     

Gas-kinetic description of shock wave structures by solving Boltzmann model equation

On the basis of the mesoscopic theory of Boltzmann-type velocity distribution function, the modified Boltzmann model equation describing the one-dimensional gas flows from various flow regimes is presented by incorporating the molecular interaction models relating to the viscosity and diffusion cross-sections, density, temperature and the dependent exponent of viscosity into the molecular collision frequency. The gas-kinetic numerical method for directly solving the molecular velocity distribution function is studied by introducing the reduced distribution functions and the discrete velocity ordinate method, in which the unsteady time-splitting method and the NND finite difference scheme are applied. To study the inner flows of non-equilibrium shock wave structures, the one-dimensional unsteady shock-tube problems with various Knudsen numbers and the steady shock wave problems at different Mach numbers are numerically simulated. The computed results are found to give good agreement with the theoretical, DSMC and experimental results. The computing practice has confirmed the good precision and reliability of the gas-kinetic numerical algorithm in solving the highly nonequilibrium shock wave disturbances from various flow regimes.