随机结构弹性屈曲的递推求解方法

采用随机收敛的非正交的多项式展式表示未知的随机屈曲特征值和屈曲模态,利用摄动技巧,建立了随机结构弹性屈曲的递推求解方法。算例表明,和基于泰勒展开的摄动随机有限元方法相比,方法的结果能在较宽的随机涨落范围内更好地逼近蒙特卡洛模拟结果,即使只采用前四阶非正交多项式展式,逼近的结果仍然较好。     

随机结构重特征值分析的递推随机有限元法

利用递推随机有限元方法研究了具有随机参数结构的重特征值问题。采用随机收敛的非正交多项式展式表示未知的随机重特征值和随机特征向量，建立了和摄动法类似的一系列确定的递推方程，通过求解这些速推方程，得到了重特征值的统计值。算例表明，同基于二阶泰勒展开的摄动随机有限元法相比，递推随机有限元法的结果能在较宽的随机涨落范围内更好地逼近蒙特卡洛模拟结果。     

求解粘弹性问题的时域自适应等几何比例边界有限元法

提出一种基于分段时域自适应算法和等几何分析的求解粘弹性问题的数值方法。利用时域分段展开,建立了递推格式的比例边界元求解方程,环向比例边界采用等几何技术离散,在继承常规比例边界有限元半解析、便于处理应力奇异性/无限域问题等优点的同时,可更准确地描述几何边界,由此进一步提高了计算精度;在时域,通过分段时域自适应计算,保证不同时间步长下的计算精度。通过数值算例,从计算精度、收敛性等方面,对所提方法的有效性进行了验证。     

基于统计模型的结构损伤识别

提出了一种基于递推随机有限元方法(RSFEM)的随机结构损伤识别方法。在定义了随机损伤指数概念的基础上,考虑模型误差的不确定性和测量噪声的影响,建立了关于随机损伤指数的控制方程。然后,利用RSFEM得到了结构随机损伤指数的统计特性。数值算例的结果显示,新的方法能在考虑模型误差和测量噪声的情况下对结构损伤进行有效识别,且在结构随机参数有较大涨落情况下,该方法仍能有效识别出结构损伤,识别结果与蒙特卡洛模拟解非常吻合。     

A perturbation method for stochastic meshless analysis in elastostatics

A stochastic meshless method is presented for solving boundary‐value problems in linear elasticity that involves random material properties. The material property was modelled as a homogeneous random field. A meshless formulation was developed to predict stochastic structural response. Unlike the finite element method, the meshless method requires no structured mesh, since only a scattered set of nodal points is required in the domain of interest. There is no need for fixed connectivities between nodes. In conjunction with the meshless equations, classical perturbation expansions were derived to predict second‐moment characteristics of response. Numerical examples based on one‐ and two‐dimensional problems are presented to examine the accuracy and convergence of the stochastic meshless method. A good agreement is obtained between the results of the proposed method and Monte Carlo simulation. Since mesh generation of complex structures can be a far more time‐consuming and costly effort than the solution of a discrete set of equations, the meshless method provides an attractive alternative to finite element method for solving stochastic mechanics problems. Copyright © 2001 John Wiley & Sons, Ltd.     

基于同伦分析方法的随机结构静力响应求解

该文提出了一种基于同伦分析方法的求解含随机参数结构的静力响应的新方法。该方法将随机静力平衡方程重新进行同伦构造,利用含随机变量和趋近函数的同伦级数展式来表示结构的随机静力位移响应,该同伦级数的各阶确定性系数和趋近函数可通过对一系列的变形方程求解得到。由于趋近函数的引入,该同伦级数解相较于传统的摄动法有更大的收敛范围,对于含较大变异性随机参数的结构也能获得不错的求解精度。同时,该文提出了一种降维策略来提高该方法的计算效率。数值算例表明,与目前广泛应用的广义正交多项式展开法（GPC）相比,从计算精度上看,该文方法的3阶展开与GPC 2阶展开相当,该文方法的6阶展开与GPC 4阶展开相当,而计算时间上前者均明显少于后者。此外,该文方法也可以方便地应用到随机结构的几何非线性分析当中,并具有较好的计算精度和计算效率。     

Hybridization of stochastic reduced basis methods with polynomial chaos expansions

We propose a hybrid formulation combining stochastic reduced basis methods with polynomial chaos expansions for solving linear random algebraic equations arising from discretization of stochastic partial differential equations. Our objective is to generalize stochastic reduced basis projection schemes to non-Gaussian uncertainty models and simplify the implementation of higher-order approximations. We employ basis vectors spanning the preconditioned stochastic Krylov subspace to represent the solution process. In the present formulation, the polynomial chaos decomposition technique is used to represent the stochastic basis vectors in terms of multidimensional Hermite polynomials. The Galerkin projection scheme is then employed to compute the undetermined coefficients in the reduced basis approximation. We present numerical studies on a linear structural problem where the Youngs modulus is represented using Gaussian as well as lognormal models to illustrate the performance of the hybrid stochastic reduced basis projection scheme. Comparison studies with the spectral stochastic finite element method suggest that the proposed hybrid formulation gives results of comparable accuracy at a lower computational cost.     

Interval spectral stochastic finite element analysis of structures with aggregation of random field and bounded parameters

This paper presents the study on non‐deterministic problems of structures with a mixture of random field and interval material properties under uncertain‐but‐bounded forces. Probabilistic framework is extended to handle the mixed uncertainties from structural parameters and loads by incorporating interval algorithms into spectral stochastic finite element method. Random interval formulations are developed based on K–L expansion and polynomial chaos accommodating the random field Young's modulus, interval Poisson's ratios and bounded applied forces. Numerical characteristics including mean value and standard deviation of the interval random structural responses are consequently obtained as intervals rather than deterministic values. The randomised low‐discrepancy sequences initialized particles and high‐order nonlinear inertia weight with multi‐dimensional parameters are employed to determine the change ranges of statistical moments of the random interval structural responses. The bounded probability density and cumulative distribution of the interval random response are then visualised. The feasibility, efficiency and usefulness of the proposed interval spectral stochastic finite element method are illustrated by three numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

非线性拟双曲方程的有限元配置法数值分析

基于有限元配置法,采用分片双三次Hermite插值多项式空间作为逼近函数空间,本文对粘性振动及神经传播过程中涉及的一类非线性拟双曲方程的初边值问题建立了二维半离散和全离散格式.并对两种格式证明了数值解的存在唯一性,应用微分方程先验估计的理论和技巧得到了L2模最佳阶误差估计.数值实验结果表明:所提方法在保证整体误差估计要求且不增加计算量的前提下,比传统有限元方法有更高的逼近精度,并扩展了配置法的应用范围.     

The scaled boundary finite element method in structural dynamics

The scaled boundary finite element method is extended to solve problems of structural dynamics. The dynamic stiffness matrix of a bounded (finite) domain is obtained as a continued fraction solution for the scaled boundary finite element equation. The inertial effect at high frequencies is modeled by high‐order terms of the continued fraction without introducing an internal mesh. By using this solution and introducing auxiliary variables, the equation of motion of the bounded domain is expressed in high‐order static stiffness and mass matrices. Standard procedures in structural dynamics can be applied to perform modal analyses and transient response analyses directly in the time domain. Numerical examples for modal and direct time‐domain analyses are presented. Rapid convergence is observed as the order of continued fraction increases. A guideline for selecting the order of continued fraction is proposed and validated. High computational efficiency is demonstrated for problems with stress singularity. Copyright © 2008 John Wiley & Sons, Ltd.     

Mixed Approximation Scheme of the Finite-Element Method for the Solution of Two-Dimensional Problems of the Elasticity Theory

For the solution of two-dimensional boundary-value problems of the elasticity theory, a triangular finite element, ensuring stability and convergence of mixed approximation, is proposed. The system of resolving equations of the mixed method is derived with account of strict satisfaction of static boundary conditions at the surface. To solve matrix equations of the mixed method, various algorithms of the conjugate-gradient method with the pre-conditional matrix have been considered. Numerical data on convergence and accuracy of the solution for a number of test problems of the elasticity theory and fracture mechanics are given. The results obtained by the conventional and mixed finite-element method approaches are compared.     

Use of higher‐order shape functions in the scaled boundary finite element method

The scaled boundary finite element method is a novel semi‐analytical technique, whose versatility, accuracy and efficiency are not only equal to, but potentially better than the finite element method and the boundary element method for certain problems. This paper investigates the possibility of using higher‐order polynomial functions for the shape functions. Two techniques for generating the higher‐order shape functions are investigated. In the first, the spectral element approach is used with Lagrange interpolation functions. In the second, hierarchical polynomial shape functions are employed to add new degrees of freedom into the domain without changing the existing ones, as in the p‐version of the finite element method. To check the accuracy of the proposed procedures, a plane strain problem for which an exact solution is available is employed. A more complex example involving three scaled boundary subdomains is also addressed. The rates of convergence of these examples under p‐refinement are compared with the corresponding rates of convergence achieved when uniform h‐refinement is used, allowing direct comparison of the computational cost of the two approaches. The results show that it is advantageous to use higher‐order elements, and that higher rates of convergence can be obtained using p‐refinement instead of h‐refinement. Copyright © 2005 John Wiley & Sons, Ltd.     

Gaussian process emulators for the stochastic finite element method

This paper explores a method to reduce the computational cost of stochastic finite element codes. The method, known as Gaussian process emulation, consists of building a statistical approximation to the output of such codes based on few training runs. The incorporation of emulation is explored for two aspects of the stochastic finite element problem. First, it is applied to approximating realizations of random fields discretized via the Karhunen–Loève expansion. Numerical results of emulating realizations of Gaussian and lognormal homogeneous two‐dimensional random fields are presented. Second, it is coupled with the polynomial chaos expansion and the partitioned Cholesky decomposition in order to compute the response of the typical sparse linear system that arises due to the discretization of the partial differential equations that govern the response of a stochastic finite element problem. The advantages and challenges of adopting the proposed coupling are discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

On the dual iterative stochastic perturbation‐based finite element method in solid mechanics with Gaussian uncertainties

The main idea is a dual mathematical formulation and computational implementation of the iterative stochastic perturbation‐based finite element method for both linear and nonlinear problems in solid mechanics. A general‐order Taylor expansion with random coefficients serves here for the iterative determination of the basic probabilistic characteristics, where linearization procedure widely applicable in stochastic perturbation technique is replaced with the iterative one. The expected values and, in turn, the variances are derived first, and then, they are substituted into the equations for higher central probabilistic moments and additional probabilistic characteristics. The additional formulas for up to the fourth‐order probabilistic characteristics are derived thanks to the 10th‐order Taylor expansion. Computational implementation of this idea in the stochastic finite element method is provided by using the direct differentiation method and, independently, the response function method with polynomial basis. Numerical experiments include the simple tension of the elastic bar, nonlinear elasto‐plastic analysis of the aluminum 2D truss, and solution to the homogenization problem of periodic fiber‐reinforced composite with random elastic properties. The expected values, coefficients of variation, skewness, and kurtosis of the structural response determined via this iterative scheme are contrasted with these estimated by the Monte Carlo simulation as well as with the results of the semi‐analytical probabilistic technique following the response function method itself. Although the entire methodology is illustrated here by using the Gaussian variables where all odd‐order terms simply vanish, it can be extended towards non‐Gaussian processes as well and completed with all the perturbation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

旋翼桨叶结构载荷计算方法比较研究

基于有限转角假设,建立了刚柔耦合旋翼动力学模型。该模型考虑了刚体转动与弹性变形之间的耦合效应,相较于基于小转角假设的传统有限元模型具有明显的优势。气动力以广义力形式与桨叶刚体转动及弹性变形耦合组建方程。在方程求解的单步上,分别采用力积分法、反力法以及曲率法计算桨叶剖面结构振动载荷。以BO105模型桨叶及SA349/2小铃羊直升机为仿真对象,比较研究了这三种载荷计算方法的预测精度与适用范围。对于不考虑气动力的纯结构振动载荷,三种计算方法具有相同的精度。在气弹瞬态计算中,力积分法对桨根载荷的预测精度不足。曲率法与反力法在桨叶有限元节点处得到了相近的结果。反力法预测精度取决于有限元建模精度,且只对节点处载荷有效。由于曲率法只计入弹性桨叶的弯曲曲率,该方法需要更高阶次的形函数以满足自由度二阶导数的连续性。此外,为加速收敛及减少累积误差,本文开发了基于外推法的数值积分算法。     

Application of the generalized perturbation-based stochastic boundary element method to the elastostatics

This paper is entirely devoted to the demonstration of a solution for some boundary value problems of isotropic linear elastostatics with random parameters using the boundary element method. The stochastic perturbation technique in its general nth-order Taylor series expansion version is used to express all the random parameters and the state functions of the problem. These expansions inserted in the classical deterministic equilibrium statement return up to the nth-order (both PDEs and matrix) equations. Contrary to the previous implementations of the stochastic perturbation technique, any order partial derivatives with respect to the random input are derived from the deterministic structural response function (SRF) at a given point. This function is approximated using polynomials by the least-squares method from the multiple solution of the initial deterministic problem solved for the expectations of random structural parameters. First two probabilistic moments have been computed symbolically here using the computational MAPLE environment, also as the polynomial expressions including perturbation parameter ε. It should be mentioned that such a generalized perturbation approach makes it possible to analyze all types of random variables (not only Gaussian) and to compute even higher probabilistic moments with a priori given accuracy. The entire methodology can be implemented after minor modifications to analyze nonlinear phenomena for both statics and dynamics of even heterogeneous domains.     

各向异性中厚度板壳的弹塑性大变形分析

本文提出了一个解各向异性中厚度板壳弹塑性大变形问题的一般有限元方法。几何非线性的描述采用Ｖ．Ｋａｒｍａｎ的位移应变关系；材料的弹塑性分析采用Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应。特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了计算塑性大变形和后屈曲路径的效率。几个数值算例表明本文给出的有限元方法对于各种各向异性壳体的弹塑性大变形分析有较好的精度。     

Design sensitivity analysis of nonlinear structures using endochronic constitutive model

The continuum formulation for design sensitivity analysis incorporating the unified endochronic constitutive model described in Part 1 is discretized and implemented using the standard isoparametric finite element procedure. Only the total Lagrangian approach is investigated since it has been concluded to be surperior to the updated Lagrangian approach. Solution of the sensitivity equation and its convergence criterion are reported. Several cases of a ten-member truss under static and dynamic loads are investigated for sensitivity verification for both nonshape and shape design problems. It is concluded that the incremental solution procedure for design sensitivity analysis of history dependent nonlinear problems can give quite accurate sensitivities. However, the equation of motion, constitutive equations and the sensitivity equations must be integrated very accurately. Thus small increments are generally needed and the computational effort for such problems can be substantial.     

Interval Finite Element Analysis using Interval Factor Method

A new method called the interval factor method for the finite element analysis of truss structures with interval parameters is presented in this paper. The structural parameters and applied forces can be considered as interval variables by using the interval factor method, the structural stiffness matrix can then be divided into the product of two parts corresponding to the interval factors and the deterministic value. From the static governing equations of interval finite element method of structures, the structural displacement and stress responses are expressed as the functions of the interval factors. The computational expressions for lower and upper bounds, mean value and interval change ratio of structural static responses are derived by means of the interval operations. The effect of the uncertainty of the structural parameters and applied forces on the structural displacement and stress responses is demonstrated by truss structures.