热荷载作用下Timoshenko功能梯度夹层梁的静态响应

在精确考虑轴线伸长和一阶横向剪切变形的基础上建立了Timoshenko功能梯度夹层梁在热载荷作用下的几何非线性控制方程.采用打靶法数值求解所得强非线性边值问题,获得了两端固支功能梯度夹层梁在横向非均匀升温作用下的静态热过屈曲和热弯曲变形数值解.分析了功能梯度材料参数变化、不同表层厚度和升温参数对夹层梁弯曲变形、拉-弯耦...     

功能梯度截顶圆锥壳的热弹性弯曲精确解

研究了功能梯度材料截顶圆锥壳在横向机械载荷与非均匀热载荷同时作用下的变形问题. 基于经典线性壳体理论推导出了以横向剪力和中面转角为基本未知量的功能梯度薄圆锥壳轴对称变形的混合型控制方程. 假设功能梯度圆锥壳的材料性质为沿厚度方向按照幂函数连续变化的形式. 然后采用解析方法求解,得到了问题的精确解. 分别就两端简支和两端固支边界条件,给出了圆锥壳的变形随其载荷、材料参数等变化的特征关系曲线,重点分析和讨论了载荷参数与材料梯度变化参数对变形的影响.      

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

S型功能梯度扁球壳非线性屈曲的渐近分析

本文利用渐近迭代法获得了边界弹性支撑的功能梯度扁球壳的非线性屈曲问题的理论解．假设材料组分体积分数沿壳体厚度方向呈sigmoid幂函数变化,边界上考虑一般的弹性支撑约束．基于经典的薄壳理论和几何非线性关系,导出了S型功能梯度扁球壳的非线性屈曲问题的控制方程．采用渐近迭代法通过两次迭代得到了无量纲挠度和均布荷载之间的非线性特征关系．通过与已有有限元方法和其他数值方法的结果对比,验证了本文解的有效性和高精度．同时,通过计算阐述了缺陷因子、材料参数、边界约束系数及特征几何参数对壳体临界屈曲荷载的影响．     

矩形脉冲作用下复合材料层合扁球壳的冲击屈曲分析

研究了计及横向剪切的复合材料层合扁球壳在矩形脉冲载荷作用下的非线性动力屈曲问题;采用Galerkin方法得到以顶点挠度表达的动力响应方程,并用Runge-Kutta方法进行数值求解,应用Budiansky-Roth准则(简称B-R准则)确定冲击屈曲的临界荷载;讨论了壳体几何尺寸和物理参数对复合材料层合扁球壳冲击屈曲的影响;数值算例表明,该方法是可行的.     

均布压力作用下圆底扁锥壳的大挠度问题

1.引言圆底扁锥壳是工程中的常见壳型结构,与扁球壳相比,对圆底扁锥壳的大挠度理论和实验研究都不多,仅有少数文献进行过研究,文献[1]和[2]用的是伽辽金方法,文献[3]用摄动法,文献[4]用修正迭代法,它们分别在壳体几何参数λ不大的范围内确定了失稳临界载荷。由于这些方法本身存在的困难,所得解答当λ稍大时就不准确了。本文利用牛顿-样条函数方法对均布压力作用下圆底扁锥壳(图1)的非线性弯曲和稳定性进行研究,获得一些有意义的结果。2.基本方程及其求解考虑均布压力作用下圆底扁锥壳的轴对称非线性方程:     

爆炸冲击下复合材料层合扁球壳的动力屈曲

研究计及横向剪切的复合材料层合扁球壳在爆炸冲击载荷作用下的非线性轴对称动力屈曲问题。通过在复合材料层合扁球壳非线性稳定性的基本方程中增加横向转动惯量项并引入R.H.Cole理论的爆炸冲击力,得到爆炸冲击下复合材料层合扁球壳的动力控制方程,应用Galerkin方法得到用顶点挠度表达的爆炸冲击动力响应方程,并采用Runge-Kutta方法进行数值求解,采用Budiansky-Roth准则确定冲击屈曲的临界载荷,讨论了壳体几何尺寸对复合材料层合扁球壳冲击屈曲的影响;数值算例表明,此方法是可行的。     

U型波纹管大挠度变形的积分方程描述及其迭代算法

采用格林函数法,导出了U型波纹管圆环壳部分和截头扁锥壳部分的非线性的积分方程,其中的四个未知参数由圆环壳和截头扁锥壳的连接条件确定,联合应用梯度法和积分方程迭代法建立了U型波纹管大挠度分析的迭代算法,开发了相应的程序系统,数值结果表明,本文方法具有较高的精度,压缩角对峰值应力和刚度影响十分显著,应用作为波纹管设计的重要参数。     

夹层扁球壳的非线性稳定性

基于Reissner假设和变分原理,给出夹层扁球壳在均布压力作用下的大挠度方程,采用修正迭代法求得了夹层扁球壳非线性稳定问题的解析解,得到两类边界条件下临界屈曲载荷的表达式,讨论了几何参数和物理参数对临界屈曲载荷的影响     

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

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

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

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