THE STOCHASTIC BOUNDARY ELEMENT METHOD IN STATISTICAL ANALYSIS OF MODERATELY THICK PLATES

In this paper a stochastic boundary element method(SBEM)isdeveloped to analyze moderately thick plates with random material parameters andrandom thickness.Based on the Taylor series expansion,the boundary integrationequations concerning the mean and deviation of the generalized displacements arederived,respectively.It is found that the randomness of material parameters isequivalent to a random load,so the mean and covariance matrices of unknowngeneralized boundary displacements and tractions can be obtained.Furthermore,themean and covariance of generalized displacements and forces at internal points can alsobe obtained.A numerical example has been worked out with the method proposed andnecessary analysis is made for the results.     

中厚板统计分析的随机边界元法

本文建立了分析含随机材料参数并具厚度不均匀性的中厚板问题的随机边界元法,基于Taylor级数展开技术,分析和到广义位移的均值和一阶偏差的积分方程,其中将材料参数的随机性和厚度的不均匀性作为等效荷载处理,从而得到广义边界位移或面力的均值和协方差,并进一步求出部点广义位移和内力的均值和协方差,最后用本文方法计算了两个数例,并对所得结果进行了分析,探讨。     

薄板统计分析的随机边界元法

本文建立了分析含随机材料参数的薄板弯曲问题的随机边界元法。基于Ｔａｙｌｏｒ级数展开技术，分别得到了位移的均值和一阶偏差的边界积分方程，发现材料参数的随机性可作为一个等效的随机荷载处理，从而得到边界位移或边界力的均值和协方差，并进一步求出内点位移和力矩的均值和协方差，最后用本文方法计算了两个算例，并对结果进行了必要的分析。     

概率断裂分析的随机边界元法

应用随机边界元法分析材料弹性常数的随机性和裂纹面随机性对应力强度因子的影响。文中首先简介了随机边界元法，给出了具有随机材料或几何参数的弹性体的边界位移或面力的协方差，进而给出了材料参数和裂纹面随机时应力强度因子均值和方差的计算公式。算例中详细讨论了杨氏模量、泊松比及裂纹面的随机性对应力强度因子的影响。     

On the local average of random field in stochastic finite element analysis

The distinguishing feature of stochastic finite element analysis is that it involves the discretization of the parameter space of random fields of material properties, the geometry of structure and / or the loads. It is shown in earlier investigations that a reasonable procedure of discretization is to take the local averages of the random fields on each element. In the present paper the formulae for the covariance of the local averages of a homogeneous random vector field on rectangular elements are generalized by relaxing the condition. For an inhomogeneous random field and /or non-rectangular elements, a procedure of using Gaussian quadrature to evaluate the means and covariances of the local averages is proposed. Thus, the stochastic finite element method (SFEM) based on the local averages of random fields is adapted to a structure with irregular shape and / or inhomogeneous random fields. The effects of the mesh geometry, the ratio of element size to the correlation scale as well as the number of Gaussian quadrature points on the convergence of SFEM are discussed. It is found that even better results could be obtained by utilizing appropriate Gaussian quadrature instead of exact local average.Project supported by National Natural Science Foundation of China.     

Stochastic finite element method based on local averages of random fields

The stochastic finite element method (SFEM) based on the local, averages of random fields, which was proposed in [5], is now generalized to analyze the structures with several correlated random parameters. The covariance matrix of the local averages of a random vector field is derived. The SFEM based on the local averages of random vector fields is formulated. The numerical examples show that the generalized SFEM preserves the advantages of the original one, i. e., rapid convergence, good accuracy and insensitivity to the correlation structures of random parameters.Project supported by National Natural Science Foundation of China.     

随机边界元法在薄板可靠性分析中的应用

本文利用Taylor级数展开，将结构的材料随机性化为一个等效随机外载荷来处理，从而利用相应同一结构确定性问题的基本解，分别建立关于响应的均值和偏差的边界积分方程；结合一阶二次矩法，应用随机边界元法分析了材料弹性常数的随机性和薄板厚度不均匀性对薄板结构可靠性的影响，对计算结果分析看出；（1）对于同一种单元划分和同一相关长度，可靠性指数保持了对具体相关模式的少敏感性；（2）对于同一种单元划分和同一相关模式，可靠性指数随相关长度的增大而减小。     

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation‐based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1‐D as well as the 3‐D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1‐D problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Stochastic finite element analysis of non-linear plane trusses

—This study considers the responses of geometrically and materially non-linear plane trusses under random excitations. The stress-strain law in the inelastic range is based on an explicit differential equation model. After a total Lagrangian finite element discretization, the nodal displacements satisfy a system of stochastic non-linear ordinary differential equations with right-hand-sides given by random functions of time. The exact solution of the above stochastic differential equation is generally difficult to obtain. To seek an approximate solution with good accuracy and reasonable computational effort, the stochastic linearization method is used to find the first and second statistical moments (i.e. the mean vector and the one-time covariance matrix) of the nodal displacements. Results of simple structures under Gaussian white-noise excitation indicate that the proposed method has good accuracy (generally underestimates the r.m.s. stationary response by 5–14%) and requires only a small fraction of the computation time of the time-history Monte-Carlo method.     

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

弯曲梁随机刚度有限元法

弯曲梁的随机材料参数和随机几何参数可结合为一个随机参数-弯曲刚度。以位移的均值和协方差函数的变分原理为基础,根据基本方程推导出具有随机刚度的静定弯曲梁的有限元法。这种方法不再采用摄动技术,可直接建立位移的均值和协方差函数。     

Parameters inversion of fluid-saturated porous media based on neural networks

The multi-layers feedforward neural network is used for inversion of material constants of fluid-saturated porous media. The direct analysis of fluid-saturated porous media is carried out with the boundary element method. The dynamic displacement responses obtained from direct analysis for prescribed material parameters constitute the sample sets training neural network. By virtue of the effective L-M training algorithm and the Tikhonov regularization method as well as the GCV method for an appropriate selection of regularization parameter, the inverse mapping from dynamic displacement responses to material constants is performed. Numerical examples demonstrate the validity of the neural network method. Project supported by the National Natural Science Foundation of China (Nos. 19872002 and 10272003) and Climbing Foundation of Northern Jiaotong University.     

弹性、热弹性物性值反问题的边界元分析

本文采用边界元法和卡尔曼滤波,对弹性、热弹性问题物性值进行反分析,由有限个观测点的位移值,同时反算出材料的拉压弹性模量Ε、泊桑比ν和线膨胀系数α。     

An extended stochastic response surface method for random field problems

An efficient and accurate uncertainty propagation methodology for mechanics problems with random fields is developed in this paper. This methodology is based on the stochastic response surface method (SRSM) which has been previously proposed for problems dealing with random variables only. This paper extends SRSM to problems involving random fields or random processes fields. The favorable property of SRSM lies in that the deterministic computational model can be treated as a black box, as in the case of commercial finite element codes. Numerical examples are used to highlight the features of this technique and to demonstrate the accuracy and efficiency of the proposed method. A comparison with Monte Carlo simulation shows that the proposed method can achieve numerical results close to those from Monte Carlo simulation while dramatically reducing the number of deterministic finite element runs. The project supported by the National Natural Science Foundation of China (10602036). The English text was polished by Yunming Chen.     

Computing elastic moduli on 3-D X-ray computed tomography image stacks

A numerical task of current interest is to compute the effective elastic properties of a random composite material by operating on a 3D digital image of its microstructure obtained via X-ray computed tomography (CT). The 3-D image is usually sub-sampled since an X-ray CT image is typically of order 1000³ voxels or larger, which is considered to be a very large finite element problem. Two main questions for the validity of any such study are then: can the sub-sample size be made sufficiently large to capture enough of the important details of the random microstructure so that the computed moduli can be thought of as accurate, and what boundary conditions should be chosen for these sub-samples? This paper contributes to the answer of both questions by studying a simulated X-ray CT cylindrical microstructure with three phases, cut from a random model system with known elastic properties. A new hybrid numerical method is introduced, which makes use of finite element solutions coupled with exact solutions for elastic moduli of square arrays of parallel cylindrical fibers. The new method allows, in principle, all of the microstructural data to be used when the X-ray CT image is in the form of a cylinder, which is often the case. Appendix A describes a similar algorithm for spherical sub-samples, which may be of use when examining the mechanical properties of particles. Cubic sub-samples are also taken from this simulated X-ray CT structure to investigate the effect of two different kinds of boundary conditions: forced periodic and fixed displacements. It is found that using forced periodic displacements on the non-geometrically periodic cubic sub-samples always gave more accurate results than using fixed displacements, although with about the same precision. The larger the cubic sub-sample, the more accurate and precise was the elastic computation, and using the complete cylindrical sample with the new method gave still more accurate and precise results. Fortran 90 programs for the analytical solutions are made available on-line, along with the parallel finite element codes used.     

随机结构有限元分析的递推求解方法

提出了一种新的谱随机有限元分析方法——递推求解方法。该方法将随机结构的随机响应表示成非正交多项式展式,建立了和摄动法类似的一系列确定的递推方程,并通过确定性有限元方法对这些递推方程进行静力问题求解。算例表明,当随机量出现较大涨落时,计算结果相对于传统摄动法有不小的改进。     

基于摄动法的不确定性有限元模型修正方法研究

开展了考虑不确定性的有限元模型修正方法的研究。基于摄动法推导了待修正参数均值和协方差矩阵的迭代格式,其中协方差的迭代格式包括是否考虑试验数据与修正参数之间相关性的两种形式。在理论研究基础上开展数值仿真研究,实现了不确定性有限元模型修正的摄动法,并研究了试验数据样本数量对修正误差的影响。仿真结果表明,该方法适用于解决系统参数与试验数据存在不确定性的模型修正问题,试验样本数量对待修正参数标准差的修正精度影响较大;忽略试验模态参数与待修正参数不确定性之间的相关性,能够避免计算二阶灵敏度矩阵,在保证修正结果准确性的前提下减少计算量。     

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead of the finite element method is used to perform lower bound limit and shakedown analysis of structures. The self-equilibrium stress fields are constructed by a linear combination of several basic self-equilibrium stress fields with parameters to be determined. These basic self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains and obtained through elastic-plastic incremental analysis. The complex method is used to solve nonlinear programming and determine the maximal load amplifier. The limit analysis is treated as a special case of shakedown analysis in which only the proportional loading is considered. The numerical results show that SGBEM is efficient and accurate for solving limit and shakedown analysis problems. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctorial Dissertation of China (No. 200025) and the Basic Research Foundation of Tsinghua University.     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.     

悬臂轴拉杆一阶随机分析的解析解

针对悬臂轴拉杆在自由端受集中力问题进行了随机分析，其中考虑了轴拉杆杨氏模量的变异，在小变异情况下获得了不同相关结构与相关长度下问题的解析解，并就结果对相关结构的敏感性进行了讨论．