线性随机结构在随机激励下动力响应分析

利用虚拟激励法对随机结构正交展开理论进行扩展,并在Ritz向量子空间中对扩阶系统方程进行动力聚缩,提出了一类可以快速高效地进行线性随机结构复合随机振动分析的计算方法．算例分析表明,该法可以方便地分析随机结构在平稳或非平稳随机激励下的复合随机振动问题,且分析结果与 Monte Carlo模拟分析结果符合良好;与均值参数确定性结构传统随机振动分析计算结果相比,随机结构在相同随机激励下响应自谱密度曲线具有峰值降低、谱宽增大的特点．     

不确定结构动力区间分析方法研究

将结构中的不确定性参数用区间数来表示,建立了系统的区间动力特征方程和动力控制方程.通过将Monte Carlo方法引入到区间分析中来,并改进了一般Monte Carlo方法在区间模拟中存在的缺陷,对区间动力特征方程和动力控制方程进行了求解,得到了结构的固有频率、位移响应和应力响应的区间计算结果.算例表明改进的Monte Carlo方法求解简单,具有普适性,能够给出较高精度的响应区间.     

精细积分时域平均法和随机扩阶系统法

讨论含随机参数结构的动力响应的计算问题,发展了精细积分时域平均法（ＴＡＰＩＭ）,它可以用来计算确定性系统受到随机激励时的动力响应;结合随机扩阶系统方法与随机有限元法,将ＴＡＰＩＭ方法应用于计算随机参数结构的动力响应,取得了较好的结果。结出了数值算例,结果表明随机扩阶系统法,随机有限元法与精细积分时域平均法的结合是计算 随机参数结构动力响应的有效方法。     

随机结构动力分析的递归聚缩算法

首先介绍了一种新型的随机结构动力分析方法：扩阶系统方法，然后建议了用于扩阶系统动力分析的速归聚缩算法．在不降低计算精度的条件下，递归聚缩算法可以大幅度地提高随机结构动力反应分析的速度，使关于扩阶系统动力分析的计算速度接近于相应的确定性结构系统的计算速度．从而，为随机结构分析的扩阶系统方法进入到实用阶段铺平了道路．     

随机结构动力分析的递归聚食缩算法

首先介绍一种新型的随机结构动力分析方法，扩阶系统方法，然后建议了用于扩阶系统动力分析的递归聚缩算法，在不降低计算精度的条件下，递归聚缩算法可以大幅度地提高随机结构动力反应分析的速度，使关于扩阶系统动力分析的计算速度接近于相应的确定性结构系统的计算速度，从而，为随机结构分析的扩阶系统方法进入到实用了阶段铺平了道路。     

CFRP加固RC梁的模糊随机可靠度数值模拟

CFRP加固RC结构可靠度的不确定性,不仅要考虑随机性,也要考虑模糊性.本文根据模糊数学与可靠性理论,研编了引入模糊示性函数的Monte Carlo计算程序,通过试验讨论模糊示性函数的确定方法.利用Matlab软件,直接产生随机变量数组,解决Monte Carlo 抽样困难.CFRP加固梁破坏试验获得全部随机变量的统计参量.编制基于Matlab的Monte Carlo计算程序,完成了对CFRP加固RC梁抗弯模糊随机可靠度数值模拟,模拟计算结果得到传统FOSM法验证.     

复合材料拉伸断裂过程的Monte-Carlo模拟

本文假定材料的断裂韧性Г是一个按概率密度分布的随机参数，其概率密度服从Weibull分布。应用Monte—Carlo方法模拟了层合板在拉伸条件下的断裂过程。计算结果表明：理论预测值与实验结果吻合较好。并得出结论：用Monte—Carlo模型能够较好模拟层合板的拉伸断裂过程，对层合板中基体开裂的演化进行相关的预测是可行并且有效的。     

管线随机有限元的抗震可靠性研究

基于动力分析的摄动有限元法及管线可靠度分析模型 ,建立管线的动态随机有限元方程 ,分析了管线质量随机性和几何特征随机性 ,得出了在Elcentro地震波的激励下 ,管线结构的动力响应及具有非穿透性裂纹的抗震可靠度 ;并与Monte Carlo有限元法进行计算比较 ,表明了方法的正确性和可行性     

考虑颗粒碰撞的多重Monte Carlo算法

从减少计算代价和改进碰撞算法出发, 提出了考虑颗粒碰撞的多重Monte Carlo算法, 它采用直接模拟Monte Carlo算法来考虑颗粒碰撞, 并与求解颗粒拉氏Langevin方程的Monte Carlo算法耦合起来, 跟踪比实际颗粒数目小得多的虚拟颗粒. 提出了时间步长选定标准、虚拟碰撞伙伴所在控制容积的判断准则、颗粒碰撞发生的判 断准则、虚拟碰撞伙伴的选择、基于随机碰撞角度的碰撞动力学, 构成了考虑颗粒碰撞的完整多重Monte Carlo算法. 对理想工况的细微颗粒流和粗重颗粒流进行了数值模拟, 颗粒碰撞率的模拟结果与理论分析解和DNS结果均符合很好, 颗粒场演变的细节信息, 如时间平均和特定时刻的颗粒数密度, 速度和颗粒湍动能等, 均与DNS结果符合很好. 数值模拟结果证明该算法不仅具有较低的计算代价, 而且能够达到足够的计算精度.     

DYNAMICS ANALYSIS OF STOCHASTIC SPATIAL FLEXIBLE MULTIBODY SYSTEM 1)

轻质、高精度的柔性多体系统被广泛应用于实际工程领域中.由于实际设计公差、制造误差及环境温度等多种不确定因素的存在,使得柔性多体系统的结构参数(物理参数和几何参数)表现出随机性.具有随机结构参数的动力学模型能够客观地反映出真实系统的动力学行为,且结构参数的不确定性对空间柔性多体系统动力学响应的影响是不容忽视的.针对具有多个随机参数的空间柔性多体系统,提出了一种基于广义alpha算法的非侵入式随机柔性多体系统动力学计算方法.采用绝对节点坐标公式(absolute node coordinate formulation, ANCF)来描述柔性体, 推导建立多体系统动力学模型.利用混沌多项式展开(polynomial chaos expansion, PCE)法构建系统随机动力学方程的代理模型,然后将随机响应面法(stochastic response surface method, SRSM)嵌入广义-alpha方法中,分别采用改进抽样的回归方法(regression method of improved sampling, RMIS)和单项求容积法则(Monte Carlo simulation, MCR)来确定样本点.将数值计算结果与蒙特卡洛模拟(Monte Carlo simulation, MCS)结果进行对比, 验证了所提算法的有效性.在相同的定积分精度的条件下,根据单项求容积法则确定的样本点的计算结果稳定性更强, 且其计算效率更高.     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。