结构动力分析的随机变分原理及随机有限元法

将结构动力系统的参数及激励的随机性直接引入结构的动力泛函变分表达式中，基于瞬时最小势能原理，应用小参数摄动法，建立了随机结构动力分析的随机变分列式及相应的确机有限元法。算例表明，应用此法分析随机结构动力响应，具有程序实施简便，计算效率高的优点。     

The random variational principle in finite deformation of elasticity and finite element method

ＴＨＥＲＡＮＤＯＭＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＩＮＦＩＮＩＴＥＤＥＦＯＲＭＡＴＩＯＮＯＦＥＬＡＳＴＩＣＩＴＹＡＮＤＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＧａｏＨａｎｇ－ｓｈａｎ（高行山）（ＮｏｒｔｈｗｅｓｔｅｎＰｏｌｙｔｅｃｈｎｉｃａｌＵ...     

断裂力学问题的杂交边界点方法

提出了一种求解断裂力学的新的边界类型无网格方法-杂交边界点法.以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性,求解时仅仅需要边界上离散点的信息.该文将杂交边界点方法应用到弹性断裂问题中,将移动最小二乘方法中的基函数扩充,能更好的模拟裂纹尖端应力场的奇异性,推导了求解断裂力学的杂交边界点法方程,与传统的元网格方法相比,文中方法具有后处理简单,计算精度高的优点.数值算例表明了该方法的稳定性和有效性.     

弹性理论中的临界变分及消除方法

临界变分现象是拉氏乘子法的固有特性，钱伟长应用高阶拉氏乘子消除了临界变分现象。本文将提出一种新的方法－凑合反推法，这种方法摒充了拉氏乘子法，把拉氏乘子所在的项目一个待定函数Ｆ代替。这样构成的泛函，作者称之为试泛函。而待定函数Ｆ的识别类似于拉氏乘子的识别。通过该法可以方便地构造出各种多变量广义变分原理，并且可以消除临界变分现象。     

New applications of variational iteration method using Adomian polynomials

In this paper, we discuss a new application of the variational iteration method considering Adomian’s polynomials on nonlinear physical equations. Two models of interest in physics are considered and solved by means of the variational iteration method. The behavior of the variational iteration method and the effects of different values of t are investigated. Comparisons are made among the standard Adomian decomposition method, exact solutions, and the proposed method. He’s variational iteration method is introduced to overcome the difficulty arising in calculating the Adomian polynomial in Adomian decomposition method. The results reveal that the proposed method is very effective and simple and can be applied to other nonlinear problems.     

第二类抛物型变分不等式中的MRM方法

本文讨论了第二类抛物型变分不等式中的MRM(多重互易方法)方法。首先采用时间项半离散和隐格式方法将抛物型变分不等式化解为一个椭圆变分不等式，然后利用MRM-边界积分方程，将其化解为MRM-边界混合变分不等式，并给出了MRM-边界混合变分不等式解的存在唯一性。说明了该MRM-边界混合变分不等式与常规边界积分方程得到的边界混合变分不等式是一致的，并且具有更容易编程实现。这为使用MRM边界元方法数值求解抛物型变分不等式提供了方法和理论依据。文末给出了数值算例。     

Some issues of turbulence statistics

The issue of dropping the random forcef _i and the arbitrariness of choosing the basic variable in the variational approach to turbulence closure problem, pointed out recently by the Russian scientists Bazdenkov and Kukharkin, are discussed. According to the mean-square estimation method, the random forcef _i should be dropped in the error expression of the LFP (Langevin-Fokker-Planck) model. However,f _i is not neglected, its effect has been taken into account by the variational approach. In order to optimize the perturbation solution of the Liouville equation, the LFP model requires that the basic variable is as near to Gaussian as possible. Hence, the velocity, instead of the vorticity, should be chosen as the basic variable in the three-dimensional turbulence. Although the LFP model and the zero-order Gaussian term of PDF (probability density function) imply whiteness assumption (zero correlation time off _i ), the higher-order non-Gaussian terms of PDF correspond to the nonwhiteness of turbulence dynamics, the variational approach does calculate the nonwhiteness effect properly. The work is supported by the National Basic Research Program “Non-linear Sciences” and the National Natural Science Foundation of China     

Numerical solution of plane and axially symmetric jet flow problems based on the variational inequality formulation

The numerical analysis of plane and axially symmetric jet flows of an incompressible inviscid fluid is treated. A new formulation of the variational inequality type is developed from the variational principle associated with jet problems. A successive approximation method is formulated by the combined use of variational inequality and the finite element method. Numerical examples based on the iterative method are presented. The results obtained agree well with those by other methods.     

分析力学初值问题的一种变分原理形式

明确了分析力学初值问题的控制方程,按照广义力和广义位移之间的对应关系,将 各控制方程卷乘上相应的虚量,代数相加,进而在 原空间中建立了分析力学初值问题的一种变分原理形式,即建立了分析力学初值问题的卷积 型变分原理和卷积型广义变分原理. 推导了分析力学初值问题卷积型变分原理和卷积型广义 变分原理的驻值条件. 在建立分析力学初值问题的一种变分原理形式的同时, 将变积方法推广为卷变积方法.     

Semi-inverse method and generalized variational principles with multi-variables in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎ 1 954,Ｈｕ[1,2 ]ｄｅｄｕｃｅｄＨｕ＿Ｗａｓｈｉｚｕｐｒｉｎｃｉｐｌｅｂｙｓｏ＿ｃａｌｌｅｄｔｒｉａｌ＿ａｎｄ＿ｅｒｒｏｒｍｅｔｈｏｄ ,ａｎｄｉｎ1 964 ,Ｃｈｉｅｎ[3]ｓｙｓｔｅｍａｔｉｃａｌｌｙｄｉｓｃｕｓｓｅｄｔｈｅＬａｇｒａｎｇｅｍｕｌｔｉｐｌｉｅｒｍｅｔｈｏｄ ,ｂｙｗｈｉｃｈｈｅｓｕｃｃｅｓｓｆｕｌｌｙｄｅｄｕｃｅｄＨｕ＿Ｗａｓｈｉｚｕｐｒｉｎｃｉｐｌｅ.Ａｆｔｅｒｔｈａｔｇｅｎｅｒａｌｉｚｅｄｖａｒｉａｔｉｏｎａｌｐｒｉｎｃｉｐｌｅｓｃａｎｂｅａｒｒｉｖｅｄａｔ…     

变分迭代算法：一种新的非线性分析方法

本文提出了求解非线性方程的一种新方法－变分迭代算法，这种方法的基本特点是：给定一个近似解（可以包含待定常数）然后用拉氏乘子来校正其近似解，拉氏乘子可用变分的概念最佳确定，这种方法不受小参数的限制，具有很大的应用前景，本文详细介绍了这种方法，并得到了Ｄｕｆｆｉｎｇ方程的两个新的近似解。     

HIERARCHICAL STOCHASTIC FINITE ELEMENT METHOD FOR STRUCTURAL ANALYSIS

In this paper, the hierarchical approach is adopted for series representation of the stochastic nodal displacement vector using the hierarchical basis vectors, while the Karhunen-Lòeve series expansion technique is employed to discretize the random field into a set of random variables. A set of hierarchical basis vectors are defined to approximate the stochastic response quantities. The stochastic variational principle instead of the projection scheme is adopted to develop a hierarchical stochastic finite element method (HSFEM) for stochastic structures under stochastic loads. Simplified expressions of coefficients of governing equations and the first two statistical moments of the response quantities in the schemes of the HSFEM are developed, so that the time consumed for computation can be greatly reduced. Investigation in this paper suggests that the HSFEM yields a series of stiffness equations with similar dimensionality as the perturbation stochastic finite element method (PSFEM). Two examples are presented for numerical study on the performance of the HSFEM in elastic structural problems with stochastic Young’s Modulus and external loads. Results show that the proposed method can achieve higher accuracy than the PSFEM for cases with large coefficients of variation, and yield results agreeing well with those obtained by the Monte Carlo simulation (MCS).     

A convolution-type semi-analytic DQ approach to transient response of rectangular plates

The convolution-type Gurtin variational principle is known as the only variational principle that is, from the mathematics point of view, totally equivalent to the initial value problem system. In this paper, the equation of motion of rectangular thin plates is first transformed to a new governing equation containing initial conditions by using a convolution method. A convolution-type semi-analytical DQ approach, which involves differential quadrature （DQ） approximation in the space domain and an analytical series expansion in the time domain, is proposed to obtain the transient response solution. This approach offers the same advantages as the Gurtin variational principle and, at the same time, is much simpler in calculation. Numerical results show that it is very accurate yet computationally efficient for the dynamic response of plates.     

A generalized variational principle in micromorphic thermoelasticity

Recently Nappa obtains a Gurtin-type variational principle for micromorphic thermoelasticity. However, the use of convolutions is of course restricted to the linear case, which sets a limit to the rang of applicability. This paper establishes a classic variational principles for the discussed problem by the semi-inverse method.     

饱和多孔介质耦合系统的几类变分原理

采用变积方法，建立了一组等温准静态下饱和多孔介质的六类变量的广义变分原理，在此基础上，引入约束条件得到五类变量，四类变量，三类变量和二类变量的变分原理，为建立饱和多孔介质的有限元模型提供了基础。     

Variational analysis of thermomechanically coupled steady-state rolling problem

A steady-state, rigid-plastic rolling problem for temperature and strain-rate dependent materials with nonlocal friction is considered. A variational formulation is derived, coupling a nonlinear variational inequality for the velocity and a nonlinear vari- ational equation for the temperature. The existence and uniqueness results are obtained by a proposed fixed point method.     

Fuzzy variational principle and its applications

Linear and non-linear peaky fuzzy numbers and their arithmetic operations are constructed for the analysis of engineering structures with fuzzy characteristic quantities. Fuzziness of the corresponding quantities is consistently incorporated into the functional of the total potential energy. A set of deterministic recursive equations is obtained as the alternative expressions of the fuzzy variational principle by means of the second-order perturbation technique. The fuzzy Ritz method and the fuzzy finite element method are presented as the applications of the fuzzy variational principle. Accordingly, the roundabout procedures frequently used in the formulations of the fuzzy finite element method are avoided. A benchmark problem of a bending beam with fuzzy Young's modulus under fuzzy external loading is solved by the developed fuzzy numerical methods. Numerical examples show that results determined by these two fuzzy methods are both little conservative, and are in good agreement with those obtained by the analytical method. Moreover, the fuzzy Ritz method or the fuzzy finite element method can provide more valuable information than the conventional deterministic methods.     

考虑畸变影响的箱梁横向弯矩计算方法

为研究偏心垂向荷载作用下梯形截面单室箱梁的横向弯矩,对框架分析法计算箱梁横向弯矩的方程进行优化,并在刚性支承法的基础上提出一种更加简单的横向弯矩计算方法;与框架分析法不同,横向弯矩可采用能量变分法求解,建立以箱梁顶板剪力差为未知量的四阶控制微分方程,采用比拟的弹性地基梁解法解出剪力差,得出梯形截面单室箱梁横向弯矩的能量变分法解。对几种箱梁横向弯矩计算方法用两个算例进行验证,结果表明,能量变分法解将箱梁底板上的弯矩误差绝对值由15.41%降到了9.68%;本文方法得出的横向弯矩结果和有限元结果吻合较好,弯矩误差绝对值最大不超过6.01%;本文方法和能量变分法可有效降低箱梁底板上的弯矩误差,计算精度得到提高。     

求解三维裂纹前缘SIF分布时间历程的通用权函数法和有限变分法

将三维热权函数法扩展为适用于表面力、体积力和温度载荷的通用权函数法(UWF).推导出以变分型积分方程表达的UWF法基本方程,从变分的角度,将求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE)改造为可以适用于一般的变分型积分方程的一类新型数值方法--有限变分法(FVM).在FVM中可以引入无穷多种线性无关的局部变分模式,可以根据计算要求在求解域中插入任意多个计算节点,单一型裂纹问题FVM所得到的最终方程组的系数矩阵总是一个对称的窄带矩阵,而且对角元总是大数,具有良好的数值计算性能.FVM对于SIF沿裂纹前缘急剧变化的复杂情况具有较好的数值模拟能力和较高的计算精度,利用自身一致性,可以求得三维裂纹前缘SIF的高精度解.