Gauss白噪声外激下Rayleigh振子的平稳响应与首次穿越

研究了Rayleigh振子在Gauss白噪声外激下的平稳响应和首次穿越。首先利用随机平均法给出了系统随机平均It^o微分方程,对平均方程的稳态概率密度做了数值分析;然后建立了条件可靠性函数的后向Kolmogorov方程及首次穿越时间条件矩的Pontragin方程;最后对三组不同的参数值分析了首次穿越的概率统计特性。     

截尾Gauss概率密度函数的性质及待定参数的求解

截尾Gauss 分布是湍流燃烧中常用的一种瞬时标量概率密度函数形式,但其中的待定参数较难以计算. 本文对截尾Gauss 分布进行了理论分析,获得了截尾Gauss 分布的若干性质. 表明当标量脉动均方值(g) 较大时,待定参数μ和σ之间呈线性关系,这使得求解待定参数的方程组退化为包含单一变量的方程. 在一般的g 取值条件下,通过提取标量平均值(f) 的等值线并沿其进行插值,可快速地求解出不同f与g 值下的待定参数μ和σ,并建立了相应的表格.     

非Gauss分布随机波浪力对海洋平台的非线性作用

为避免求解决定Maikov过程转移概率密度的Fokker—Planck方程，基于尺度分离的假设导出了一组描述非线性海洋平台受非Gauss分布随机波浪载荷作用所产生响应的矩量的常微分方程组。矩量方程清楚地反映出分别对应随机载荷和结构响应的两种不同统计特性的相互关系。由于矩量方程不依赖载荷的概率分布的具体细节，以它来模拟随机激励作用下的非线性系统将免于Monte Carlo方法所面临的正确模拟载荷概率分布的困难任务。将摄动法用于矩量方程可使线性化不再需要，这样就不会因为线性化而产生不可预料的误差。     

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想,进一步采用分部积分将δ移出奇异积分式中积分核之外,转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和,可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

非连续边界元积分的精确表达式及相关问题

以二维位势问题边界元分析为例,给出了利用线性非连续边界元离散边界积分方程时系数矩阵积分计算的精确表达式,通过和利用Gauss积分方法计算系数矩阵所得数值结果的比较表明：配位点选择不同对数值计算结果精度影响的主要原因是积分计算的精度,尤其当配位因子选择较大时,存在的准奇异积分（Nearly Singular Integrals)很难利用常规Gauss积分方法准确求得。     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

二维位势边界元法高阶单元几乎奇异积分半解析算法

准确计算几乎奇异积分是边界元法难题之一。目前,对于一般的高阶单元的几乎奇异积分尚缺乏通用高效的计算方法。本文在单元局部坐标系中表征了二维高阶单元的几何特征,提出了源点相对高阶单元的接近度概念。针对二维位势边界元法的3节点二次等参单元,构造出与单元积分核具有相同几乎奇异性的近似奇异核函数。从二维位势几乎奇异积分单元积分核中扣除近似奇异核函数,把几乎奇异积分项转换为规则积分和奇异积分两部分之和,规则积分部分用常规Gauss数值积分计算,奇异积分部分由导出的解析公式计算,从而建立了二维位势问题高阶单元几乎强奇异和超奇异积分的半解析算法。算例结果表明了本文半解析算法的有效性和计算精度。     

三维问题边界元法中几乎奇异积分的正则化算法

当源点靠近边界单元时,边界积分方程通常存在几乎奇异积分的计算难题.基于三角形单元,将源点到单元的距离与单元特征长度比值定义为接近度,用于度量边界单元中积分奇异性的程度.将单元上的面积分在局部的极坐标系ρθ下表示,利用一些初等函数的积分公式,获得对变量ρ作单层积分的解析表达式.几乎强奇异和超奇异面积分被转化为沿单元围道上一系列线积分,而Gauss数值积分能够有效计算这些线积分.应用该算法分析三维弹性薄壁结构获得了成功.     

多点非均匀调制演变随机激励下结构地震的响应

针对大跨度结构在非均匀调制演变随机激励作用下,考虑行波效应时的非平衡随机地震响应问题,应用虚拟激励法进行了分析,由于虚拟激励法自动计及了参振振型的互相关项以及激励之间的互相关项,理论上是精确解,时变功率谱的计算采用精细逐步积分格式,使计算效率进一步得到提高。     

三维常数势边界元中的精确积分

对三维问题边界元方法中应用最广泛的常数边界元的积分提出一种精确积分方法。借助于一个假想的闭合曲面,将特定的势场应用于边界积分方程,发现对于三维问题,常数势项的积分可以化作球面三角型的面积计算,而导数项的积分则可在平面域用极坐标进行。本文方法结果精确,公式简单,同一计算公式可以用来计算非奇异、几乎奇异和奇异积分,统一了积分算法。     

二维热弹性力学边界元法中几乎奇异积分的正则化

针对二维热弹性力学边界元法中近边界点的几乎强奇异和超奇异积分,采用一种通用算法,将其实施正则化．该方法适用于线性单元,与近边界点邻近的单元上的积分采用正则化积分公式计算,远处单元的积分仍保持常规高斯积分．算例证明了该法的有效性和精确性．     

Characterization of Hamiltonian symmetries and their first integrals

We provide explicit criteria when the Hamiltonian symmetries for a finite dimensional canonical Hamiltonian system correspond to their first integrals. There are two approaches used for the construction of the first integrals once the symmetry is known. In the standard classical approach the first integrals are obtained up to a distinguished function of time t. In the second, which is recent, the integrals are given by a formula which involves the determination of the divergence terms. In both methods utilized, the first integrals are not determined uniquely. Firstly we show what conditions need to be imposed on the Hamiltonian symmetry in order that it constructively and uniquely yields a first integral. Secondly we provide the extra condition on the first integral for the first approach and the integrability conditions on the divergence term for the second. As a consequence, we show that both methods are in fact equivalent. Furthermore, it is shown that when the Hamiltonian symmetries provide first integrals they form a Lie algebra. Moreover, we prove that the Hamilton first integral is invariant under the Hamilton action symmetry. Several examples taken from the literature are given to illustrate our results and conditions.     

Analytical removal of singularity and direct evaluation of singular integrals in 3-D elastoplastic finite deformation analysis by BEM

As a further development of the present authors' research work [1,2], in this paper a method of the so-called quadratic pentahedron polar co-ordinate transformation and analytical removal of singularity of Cauchy principal value singular integrals is proposed to evaluate the strongly singular integrals in the sense of Cauchy principal values and the weakly singular integrals over quadratic internal cells in 3-D elastoplastic finite deformation analysis by BEM. First, a quadratic pentahedron polar co-ordinate transformation technique is used to reduce the order of singularity of the singular integrals. Then, a form of Gauss' theorem is introduced to remove the singularity in the Cauchy principal value singular integrals analytically. Therefore, the evaluation of all those strongly and weakly singular integrals can be carried out by standard Gaussian quadrature accurately and efficiently. Numerical examples of the 3-D elastoplastic problem and 3-D finite deformation problem are given to demonstrate that the method possesses good accuracy and numerical stability, and is convenient to implement. The method in this paper can be applied extensively to evaluating the singular integrals over cubic and higher order elements.     

Complementary energy release rates and dual conservation integrals in micropolar elasticity

The complementary energy momentum tensor, expressed in terms of the spatial gradients of stress and couple-stress, is used to construct the and conservation integrals of infinitesimal micropolar elasticity. The derived integrals are related to the release rates of the complementary potential energy associated with a defect translation or rotation. A nonconserved integral is also derived and related to the energy release rate that is associated with a self-similar cavity expansion. The results are compared to those obtained on the basis of the classical energy momentum tensor, expressed in terms of the spatial gradients of displacement and rotation, and the release rates of the potential energy. It is shown that the evaluation of the complementary conservation integrals is of similar complexity to that of the classical conservation integrals, so that either can be effectively used in the energetic analysis of the mechanics of defects. The two-dimensional versions of the dual conservation integrals are then derived and applied to an out-of-plane shearing of a long cracked slab.     

The tallest column problem: New first integrals and estimates

We analyze the problem of finding the shape of the tallest column. For the system of equations that determine the optimal shape we construct a variational principle and two new first integrals. From the first integrals we are able to determine, analytically, the size of the cross-sectional area of the optimal column at the bottom, as well as the corresponding bending moment and curvature of the elastic line. Our result for critical load is compared with the results obtained by other methods.     

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。