热/机械载荷下功能梯度材料矩形厚板的弯曲行为

采用Reddy高阶剪切板理论，考虑材料物性参数随坐标和温度变化的特性，研究在均匀变化的温度场内功能梯度材料矩形板在面内与横向载荷共同作用下的横向弯曲问题，基于一维DQ法和Galerkin技术，给出了一对边固支，另对边任意约束时板弯曲问题的半解析解，以Si3N4/SUS304板为例考察了材料组份，温度场，面内载荷及边界约束条件等对功能梯度材料板弯曲行为的影响。     

弹性地基中含孔隙的功能梯度圆锥壳的振动分析

为研究弹性地基中多孔功能梯度材料圆锥壳的振动特性,基于经典薄壳理论建立了弹性地基中含均匀和非均匀分布孔隙的功能梯度材料圆锥薄壳的振动方程,并用伽辽金法求得了自由振动和动力响应的解。通过参数分析讨论了孔隙、弹性地基参数、半锥角等因素对功能梯度圆锥壳自由振动和动力响应的影响。结果表明,弹性地基的压缩和剪切刚度的增大提高了圆锥壳的振动频率而显著减小了动力响应;当半锥角增大时,圆锥壳的动力响应显著增大。与非均匀分布孔隙壳体相比,均匀分布孔隙壳体的自振频率和动力响应随孔隙率的变化更为敏感。     

轴对称自由振动功能梯度材料微圆板中的热弹性阻尼

基于经典弹性薄板理论和单向耦合热传导理论,研究了材料性质沿厚度连续变化的功能梯度微圆板的热弹性阻尼特性．首先,考虑热力耦合效应,建立了功能梯度微圆板轴对称横向自由振动微分方程．然后,忽略温度梯度在面内的变化,建立了单向耦合变系数一维热传导方程．采用分层均匀化近似方法,将变系数热传导方程转化为一系列常系数的微分方程,利用上下表面的热边界条件和层间连续性条件获得了微圆板温度场解析解．将所得温度场代入微圆板的自由振动微分方程,得到了包含热弹性阻尼的复频率,从而获得了反映热弹性阻尼水平的逆品质因子．最后,针对材料性质沿板厚按幂函数变化的陶瓷-金属功能梯度微圆板,定量地分析材料梯度指数、几何尺寸、边界条件、温度环境等对微圆板热弹性阻尼的影响．     

功能梯度材料扁锥壳的热弹性大变形分析

研究了功能梯度材料扁薄锥壳在横向非均匀升温场中的几何非线性大变形问题.基于von Kármán几何非线性理论推导出了以中面位移为基本未知量的功能梯度扁薄锥壳在横向非均匀热载荷作用下的轴对称大挠度控制方程.采用打靶法数值求解所得非线性常微分方程边值问题,得到了锥壳的大挠度弯曲变形数值解.给出了锥壳的变形与其形状参数、载荷和材料参数等变化的特征关系曲线,分析和讨论了温度参数和材料梯度变化参数对变形的影响.     

界面特性对功能梯度智能梁静动态响应的影响研究

采用状态空间法分析了两边简支的含压电夹层的功能梯度梁的静力弯曲和自由振动问题．为了考虑中间压电层与上、下功能梯度层之间的粘结效果,采用线性弹簧模型以模拟界面性能．假设上下功能梯度层的材料参数沿厚度连续变化,而压电层则是均匀材料,并且它们都是正交各向异性的．由于功能梯度梁的不均匀性使得直接求解比较困难,文中用层合模型来进行近似．数值算例中,分别考虑了压电层用于传感器或作动器的情形,分析了粘结界面完美程度对组合梁静力弯曲和自由振动频率的影响．     

功能梯度压电、压磁柱对称问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下, 假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布时, 研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题. 首先在柱坐标系下, 由功能梯度材料的参数、本构、梯度和平衡方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解, 进而对空间柱对称的压电、压磁功能梯度材料动力控制问题进行了理论分析. 可以看出, 当梯度参数时, 即完全退化为横观各项同性压电、压磁柱对称的振动问题, 与文献[16]的基本方程为柱坐标下得出的结果一致. 最后给出数值算例, 结果表明材料不均匀性对沿径向振动各物理量的显著影响, 且用一个特定不均匀性参数值可以优化电磁力耦合的性能, 这在现代工程设计中尤为重要.     

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

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

内含弹性介质功能梯度压电球壳的径向振动调控

针对内含弹性介质的功能梯度压电球壳的径向振动问题,利用变量替换技术,求得了解析解。该解适用于材料参数沿厚度按幂律变化且密度的梯度指标可与其它物理量的梯度指标不同的功能梯度压电球壳,克服了以往分析中所有材料参数梯度指标均假设为相等的局限。数值结果表明,功能梯度压电球壳的径向振动特性可利用材料参数梯度指标和内部弹性介质的刚度进行有效调控。     

功能梯度压电、压磁柱对称问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下, 假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布时, 研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题. 首先在柱坐标系下, 由功能梯度材料的参数、本构、梯度和平衡方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解, 进而对空间柱对称的压电、压磁功能梯度材料动力控制问题进行了理论分析. 可以看出, 当梯度参数时, 即完全退化为横观各项同性压电、压磁柱对称的振动问题, 与文献[16]的基本方程为柱坐标下得出的结果一致. 最后给出数值算例, 结果表明材料不均匀性对沿径向振动各物理量的显著影响, 且用一个特定不均匀性参数值可以优化电磁力耦合的性能, 这在现代工程设计中尤为重要.     

功能梯度材料明德林矩形微板的热弹性阻尼

功能梯度材料微板谐振器热弹性阻尼的建模和预测是此类新型谐振器热?弹耦合振动响应的新课题. 本文采用数学分析方法研究了四边简支功能梯度材料中厚度矩形微板的热弹性阻尼. 基于明德林中厚板理论和单向耦合热传导理论建立了材料性质沿着厚度连续变化的功能梯度微板热弹性自由振动控制微分方程. 在上下表面绝热边界条件下采用分层均匀化方法求解变系数热传导方程, 获得了用变形几何量表示的变温场的解析解. 从而将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的偏微分方程. 然后,利用特征值问题在数学上的相似性,求得了四边简支条件下功能梯度材料明德林矩形微板的复频率解析解, 进而利用复频率法获得了反映谐振器热弹性阻尼水平的逆品质因子. 最后, 给出了材料性质沿板厚按幂函数变化的陶瓷?金属组分功能梯度矩形微板的热弹性阻尼数值结果. 定量地分析了横向剪切变形、材料梯度变化以及几何参数对热弹性阻尼的影响规律. 结果表明, 采用明德林板理论预测的热弹性阻尼值小于基尔霍夫板理论的预测结果, 而且两者的差别随着相对厚度的增大而变得显著.      

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

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

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

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