Neumann stochastic finite element method for calculating temperature field of frozen soil based on random field theory

To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory （MATLAB） sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.

High-efciency improved symmetric successive over-relaxation preconditioned conjugate gradient method for solving large-scale finite element linear equations

Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.

High-Order Spectral Stochastic Finite Element Analysis of Stochastic Elliptical Partial Differential Equations

This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.

Finite Difference Preconditioners for Legendre Based Spectral Element Methods on Elliptic Boundary Value Problems

Finite difference type preconditioners for spectral element discretizations based on Legendre-Gauss-Lobatto points are analyzed. The latter is employed for the approximation of uniformly elliptic partial differential problems. In this work, it is shown that the condition number of the resulting preconditioned system is bounded independently of both of the polynomial degrees used in the spectral element method and the element sizes. Several numerical tests verify the h-p independence of the proposed preconditioning.

Interval Analysis of the Finite Element Method for Stochastic Structures

A random parameter can be transformed into an interval number in the structural analysis with the concept of the confidence interval. Hence, analyses of uncertain structural systems can be used in the traditional FEM software. In some cases, the amount of solutions in stochastic structures is nearly as many as that in the traditional structural problems. In addition, a new method to evaluate the failure probability of structures is presented for the needs of the modern engineering design.

A MPI PARALLEL PRECONDITIONED SPECTRAL ELEMENT METHOD FOR THE HELMHOLTZ EQUATION

Spectral element method is well known as high-order method, and has potential better parallel feature as compared with low order methods. In this paper, a parallel preconditioned conjugate gradient iterative method is proposed to solving the spectral element approximation of the Helmholtz equation. The parallel algorithm is shown to have good performance as compared to non parallel cases, especially when the stiffness matrix is not memorized. A series of numerical experiments in one dimensional case is carried out to demonstrate the efficiency of the proposed method.

Stochastic Finite Element Method for Mechanical Vibration Based on Conjugate Gradient(CG)

When material properties, geometry parameters and applied loads are assumed to be stochastic, the vibration equation of a system is transformed to static problem by using Newmark method. In order to improve the computational efficiency and to save storage, the Conjugate Gradient （CG） method is presented. The CG is an effective method for solving a large system of linear equations and belongs to the method of iteration with rapid convergence and high precision. An example is given and calculated results are compared to validate the proposed methods.

A Sub-element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods

In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.

A METHOD OF SOLVING THE FUZZY FINITE ELEMENT EQUATIONS IN MONOSOURCE FUZZY NUMBERS

For same cases the rules of monosource fuzzy numbers con be used into the solution of fuzzy stochastic finite element equations in engineering. This method can reduce the computing quantity of the solution. It can be proved that the amount of the solution is nearly as much as that with the general stochastic finite element method (SFEM). In addition, a new method to appreciate the structural fuzzy failure probability is presented for the needs of the modem engineering design.

Multiple linear system techniques for 3D finite element method modeling of direct current resistivity

The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.

RELIABILITY OF ELASTO-PLASTIC STRUCTURE USING FINITE ELEMENT METHOD

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.

Stochastic Finite Element Analysis of Plate and Shell Construction

The response of random plate and shell construction is analyzed with the stochastic finite element method (SFEM). Random material properties and geometric dimensions of construction are involved in this paper. A simplified isoparametric local average model is used to describe the random field. Numerical results of the examples indicate that the approach presented herein is an economical and efficient solution for such an analysis compared with Monte Carlo simulation (MCS).

Well-Posedness and Finite Element Approximations for Elliptic SPDEs with Gaussian Noises

The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises.The noise term is approximated through the spectral projection of the covariance operator,which is not required to be commutative with the Laplacian operator.Through the convergence analysis of SPDEs with the noise terms replaced by the projected noises,the well-posedness of the SPDE is established under certain covariance operator-dependent conditions.These SPDEs with projected noises are then numerically approximated with the finite element method.A general error estimate framework is established for the finite element approximations.Based on this framework,optimal error estimates of finite element approximations for elliptic SPDEs driven by power-law noises are obtained.It is shown that with the proposed approach,convergence order of white noise driven SPDEs is improved by half for one-dimensional problems,and by an infinitesimal factor for higher-dimensional problems.

A Theory and Method for Modeling of Structures with Stochastic Parameters

In order to reflect the stochastic characteristics of structures more comprehensively and accurately, a theory and method for modeling of structures with stochastic parameters is presented by using probability finite element method and stochastic experiment data of structures based on the modeling of structures with deterministic parameters. Double-decker space frame is taken as an example to validate this theory and method, good results are gained.

IMPROVED PRECONDITIONED CONJUGATE GRADIENT METHOD AND ITS APPLICATION IN F.E.A.FOR ENGINEERING

In this paper two theorems with theoretical and practical significance are given in respect to the preconditioned conjugate gradient method (PCCG). The theorems discuss respectively the qualitative property of the iterative solution and the construction principle of the iterative matrix. The authors put forward a new incompletely LU factorizing technique for non-M-matrix and the method of constructing the iterative matrix. This improved PCCG is used to calculate the ill-conditioned problems and large-scale three-dimensional finite element problems, and simultaneously contrasted with other methods. The abnormal phenomenon is analyzed when PCCG is used to solve the system of ill-conditioned equations, ft is shown that the method proposed in this paper is quite effective in solving the system of large-scale finite element equations and the system of ill-conditioned equations.

重塑黄土直剪试验误差分析及对边坡可靠性影响评价

为通过室内快剪试验误差来源分析强度参数的随机特性,比较了击实法和压样法对抗剪强度离散性的影响,并在两种制样方式的基础上,进行再现性试验,通过对抗剪强度及正应力修正,采用Grubbs准则检验可疑值,求得了能反映强度参数随机特性的概率分布,并使用K-L级数展开法离散强度参数正态随机场,分析了两种制样方式对边坡可靠性的影响.研究表明:压样法制样比击实法获得数据离散性更小,使用再现性法能够体现强度参数的随机特性;两种制样方法在边坡稳定性分析时安全系数接近,由于击实法统计数据的离散性更大,可靠性分析时比压样法产生更大的失效概率.

随机振动下管径壁厚对航空液压管路的疲劳寿命的影响

根据目前存在由流量设计的不同管径液压管路的疲劳寿命达不到要求的情况,不同管径和壁厚的航空液压导管在随机振动工况下进行疲劳寿命的影响分析。将管路模型简化处理并进行模态分析和随机振动响应分析,最后结合疲劳损伤累积理论和材料的S-N曲线数据得出不同管径壁厚下的疲劳寿命时间。结果表明:当管径和壁厚越大时,对疲劳寿命时间的影响越大。

双洞公路隧道中夹岩柱可靠度研究

目的为了探究双洞隧道开挖后中夹岩柱的变形规律及破坏特征,改进岩柱稳定性定量评估方法。方法以某市公路双洞隧道为例,建立了非侵入式随机有限元-蒙特卡罗分析框架,提出基于中夹岩柱中部轴线变形的应变失稳准则。通过改变岩土体的物理力学参数变异系数、岩柱宽度以及隧道埋深研究双洞隧道开挖后中夹岩柱的失效概率,并分析其变化规律,进而采用系统稳定性分析法对影响岩柱的参数进行敏感性分析。结果在该工程背景下的中夹岩柱失效概率为0.2204,需定期进行稳定性监测。比较各影响因素归一化敏感度因子,岩柱失效概率对各因素敏感性由高到低依次为弹性模量变异系数、岩柱宽度、黏聚力变异系数、内摩擦角变异系数与隧道埋深,后续设计相似岩体参数分布的双洞隧道时,应关注弹性模量变异系数对岩柱稳定性的影响。结论非侵入式随机有限元-蒙特卡罗分析框架能够对中夹岩柱可靠度进行计算,实现结构的稳定性定量评估,为双洞隧道设计及稳定性控制提供一定的理论依据。