双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

连续体结构屈曲约束的ICM方法拓扑优化

基于ICM(独立、连续、映射)方法解决具有屈曲约束的连续体拓扑优化问题。建立以结构重量为目标,以屈曲临界力为约束的拓扑优化模型;采用独立的连续拓扑变量,借助泰勒展式将目标函数作二阶近似展开;借助瑞利商、泰勒展式、过滤函数将约束化为近似显函数,避免了灵敏度的计算;将优化模型转化为对偶规划,并利用序列二次规划求解,减少了设计变量的数目,缩小了模型的求解规模。给出三个算例,结果表明:该方法可有效地解决屈曲约束的连续体拓扑优化问题,能够得到合理的拓扑结构,并有较高的计算效率。     

基于响应面方法的桁架截面敏度分析和优化

把响应面方法引入桁架截面优化中,将应力和位移约束近似表达为桁架截面倒变量的线性函数。为拟合响应面,基于中心复合和单纯形试验设计方法开发了中心对称和拟单纯形试验设计两种方法,既可保证约束近似精度,又降低了结构分析计算量。对于桁架结构重量目标函数,直接推出倒变量的二阶形式,以桁架重量最小为目标的优化问题构造为标准的二次规划模型。将响应面方法计算的位移对设计变量的敏度与莫尔积分方法的近似显式进行了对比。以MSC.Patran为平台的数值算例表明:结合响应面方法,应用序列二次规划对问题进行寻优,其收敛精度及稳定性都可获得保障。     

基于快速Fourier变换的可靠性分析方法

提出了一种用于非线性隐式功能函数可靠性分析的快速Fourier变换方法.首先利用高效多维数值积分方法或基于响应面法的卷积积分求解功能函数的分布函数的特征函数,也就是进行Fourier正变换;然后利用快速Fourier逆变换数值求解功能函数的概率密度函数,进而求得失效概率.给出的算例表明此法具有精度和计算效率较高的优点.     

改进的基于雅可比椭圆函数的随机平均法

改进了基于雅可比椭圆函数的随机平均法,用于预测高斯白噪声激励下硬弹簧及软弹簧系统的随机响应.引入包含雅可比椭圆正弦函数、余弦函数及delta函数的雅可比椭圆函数变换,导出关于响应幅值和相位的随机微分方程.应用随机平均原理,将响应幅值近似为Markov扩散过程,建立其平均的It随机微分方程.响应幅值的稳态概率密度由相应的简化Fokker-Planck-Kolmogorov方程解出,进而得到系统位移和速度的稳态概率密度.以Duffing-Van der Pol振子为例,研究了硬弹簧及软弹簧情形下的随机响应,通过与Monte-Carlo数值模拟结果比较证实了本文方法的可行性及精度.     

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法，首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身成更好的线性关系。因此以广义中间变量做一介泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算，近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明本文的方法。是十分有效的。     

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法。首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身（即节点坐标和截面积）成更好的线性关系。因此，以广义中间变量做一阶泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算。近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明了本文的方法是十分有效的。     

桁架拓扑优化的多点逼近遗传算法

提出一种基于多点逼近函数和遗传算法的桁架拓扑优化方法。该方法建立了包含连续尺寸和离散拓扑两类变量的优化模型，并通过构造多点逼近函数建立了结构优化问题的第一级序列显式近似，然后采用分层优化方法：在外层对拓扑变量采用遗传算法进行优化；在内层对尺寸变量通过可由对偶法求解的第二级序列近似问题进行优化。几个经典的桁架拓扑优化算例表明该方法能以较少的结构分析次数获得比较理想的概率意义上的最优解。     

含模糊随机变量的多模式系统广义失效概率计算的内插迭代线抽样方法

通过正则化基本变量的度量空间,定义了单个基本变量同时具有模糊和随机双重不确定性时的广义失效概率.在广义失效概率的计算中,模糊随机变量被等价变换为随机变量.从而使得广义失效概率的计算变换为随机失效概率的计算.当模糊随机变量的密度函数和隶属函数均为正态型时,推导了其等价概率密度函数的形式和参数.采用自适应线抽样方法对基本变量同时具有模糊和随机不确定性时的多模式广义失效概率进行了计算,并采用数值算例对自适应线抽样广义失效概率计算方法的效率和精度进行了验证.算例分析表明该方法的计算结果是合理的,并且由于自适应线抽样法具有较高的效率和精度,因而所提方法具有一定的工程意义.     

结构可靠性分析高效自适应重要抽样方法

提出了一种基于自适应Metropolis算法和快速高斯变换技术的结构可靠性分析高效自适应重要抽样方法. 该方法首先利用自适应Metropolis算法高效生成结构失效域样本, 然后运用自适应宽核密度估计方法构造重要抽样密度函数, 最后采用快速高斯变换加速重要抽样过程中核函数的计算. 与传统方法相比, 自适应Metropolis算法能够在相同计算量下提供更多结构失效域信息从而改善计算精度, 即为求得给定精度问题的解, 可有效减少样本生成过程中的结构分析次数, 提高方法的计算效率; 快速高斯(Gauss)变换大幅降低核密度估计的计算复杂度从而大幅缩减重要抽样的计算耗时. 通过数值算例可以看出该方法具有较高的计算精度和效率.      

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

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

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


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).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).