均布压力作用下圆弧双铰拱对称和反对称失稳的折迭突变和尖点突变模型

建立了径向均布压力作用下圆弧双铰拱正对称失稳和反对称失稳的折迭和尖点突变模型，模型可以很好地描述这两个失稳过程的各种性态，求解拱失稳时的临界荷载，给出拱屈曲后路径与荷载的初步关系。     

薄壳失稳机理浅析

对薄壳失稳问题研究的理论与实验成果进行了总结和讨论，对薄壳后屈曲理论研究结果提出了不同的看法，同时应用动力学原理对薄壳失稳问题进行了探讨，并建立了计算模型。文中应用动力学原理描述了从加载初期的一个呈现静力学特征的薄壳随荷载的增加而逐渐成为一个呈现动力学特征的薄壳的过程，从薄壳受扰振动乃至共振的角度解释了失稳临界荷载实验数据值及其离散并低于失稳临界荷载理论值的原因。     

基于数值模拟的混凝土系列有效断裂判据

基于ANSYS平台,对某系列混凝土断裂破坏试验的加载过程进行有限元数值模拟,获得了荷载-加力 点位移曲线.利用获得的曲线,在线性渐进叠加假定的基础上,确定了试件在失稳断裂时的临界有效裂缝长度.将模拟峰值荷载和临界有效裂缝长度代入由LEFM推导出的公式计算了考虑断裂过程区影响的有效断裂韧度.最后,通过对系列试验有效断裂韧...     

横向冲击作用下高桥墩的几何非线性动力稳定性

考虑大变形的影响,建立了高桥墩在横向冲击荷载作用下的非线性动力学基本方程式;通过位移形函数假设,采用伽辽金积分方法得到了时间变率的动力学控制方程;对时间变率的非线性微分方程进行数值求解,给出了不同冲击荷载作用下不同柔度高桥墩的位移响应曲线以及其从产生横向振动到失稳的全过程,得到了横向冲击时高桥墩失稳的临界冲击荷载和失稳时刻;通过数值算例比较了三角形和矩形冲击荷载作用下高桥墩的荷载-位移响应曲线;分析了横向冲击力幅值、冲击区域大小、冲击持续时间、高桥墩柔度、桥面质量大小对位移幅值响应曲线、临界冲击荷载、失稳时刻的影响。结果表明:无论是在矩形还是三角形形式横向冲击作用下,随着冲击区域的增大高桥墩的振动幅值变大;桥面简化质量越大,高桥墩失稳的临界冲击载荷越小;柔度越大时,高桥墩失稳的临界冲击荷载幅值越小;对于矩形冲击,当冲击持续时间大于1s后,冲击持续时间的增加对高桥墩的稳定性无明显影响;对于三角形冲击荷载,随着冲击持续时间的增大高桥墩的振动幅值变大。     

斜靠式拱桥结构侧倾失稳分析的Ritz法

采用Ritz法,首次推导了主拱肋拱脚铰接、稳定拱肋拱脚固接边界条件下,斜靠式拱桥的侧倾失稳临界荷载系数的计算公式,并通过有限元法验证了该计算公式的正确性.进一步分析了稳定拱肋弯曲变形能和扭转变形能、主拱肋与稳定拱肋间横撑切向和径向弯曲变形能及吊杆非保向力位势和桥面系侧向弯曲变形能对斜靠式拱桥侧倾失稳临界荷载系数λcr的影响.研究结果表明:稳定拱肋的侧向弯曲变形能对斜靠式拱桥侧倾失稳临界荷载系数λcr有一定的影响,而扭转变形对侧倾失稳临界荷载系数λcr影响甚微;若忽略横撑的所有变形能,侧倾失稳临界荷载系数λcr降低65％～83％,即横撑对斜靠式拱桥侧倾失稳临界荷载系数λcr的影响非常显著;随着稳定拱肋倾角的增加,侧倾失稳临界荷载系数λcr呈增大趋势,且增设稳定拱肋对提高斜靠式拱桥的侧向稳定性大有帮助.     

纯压钢管拱稳定临界荷载计算的等效柱法

以均布荷载下的抛物线钢管拱为研究对象,在考虑双重非线性的有限元分析基础上,讨论了完善拱和有初始几何缺陷的拱的弹性失稳和弹塑性失稳的特性,提出纯压钢管拱稳定临界荷载计算的等效柱法.分析结果表明,矢跨比是计算拱临界荷载的重要影响因素,而现有等效柱法中没有考虑这一因素的影响,为此,提出等效柱的稳定系数中考虑矢跨比影响的计算方法.有初始几何缺陷的拱将发生极值点失稳,且极值点荷载要小于分支屈曲临界荷载,为此提出缺陷拱等效柱法考虑缺陷影响的计算方法.给出了钢管拱失稳临界荷载等效柱法计算的相应公式和实用表格.与双重非线性有限元计算结果对比表明,提出的等效柱法能方便且较精确地估算钢管拱的非线性临界荷载.     

弹性屈曲临界力测试方法的理论及实验研究

本文研讨了测定弹性屈曲临界力值的一个静力法和一个动力法.对于前者,论证了Southwell 图解法的广泛适用性;后者的要点是分别测定荷载作用前后结构的自振频率并据之推算弹性屈曲临界力.     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。     

圆柱壳纯弯曲时塑性失稳临界曲率半径模型

薄壁管材在等曲率矫直生产中,塑性失稳临界曲率半径作为重要的工艺参数,直接决定了设备结构和产品质量。而目前现场仍沿用经验图表结合人工经验和反复试矫对其进行估定,亟待建立针对性的临界曲率半径数学模型以指导生产。在力学建模和分析时,就是确定具有初始曲率的圆柱壳体在纯弯曲条件下塑性失稳的临界曲率半径,为此从旋转壳体一般几何方程出发,基于J2形变理论和能量理论,运用里茨法建立了圆柱壳体在纯弯曲条件下塑性失稳时的临界弯矩,以此确定了临界曲率半径模型,并给出了数值解法。应用ANSYS/LS-DYNA进行了有限元动态仿真试验,证明了模型是近似正确的,并通过仿真对比分析证明了轴向起皱先于截面畸变是圆柱壳体在纯弯曲条件下塑性失稳的主要模态。     

钢网壳模型的动力稳定性振动台试验研究

对两个钢扁网壳结构模型在地震荷载作用下的动力稳定性进行了振动台试验研究,观察了模型的局部动力失稳和整体动力失稳现象,通过对试验测量数据结果的处理,给出了模型中杆件的内力时程曲线和结点位移时程曲线,对这些时程曲线进行分析,总结了结构产生动力失稳时的特征和规律,同时也给出了部分理论分析的结果,并结合试验结果对理论研究中给出的动力失稳定义和判定准则以及其它结论分别进行了验证和说明．     

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

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

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

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