圆薄板材料性能随温度变化的大挠度弯曲分析

记入材料性能参数随温度的变化，导出了横向力作用下圆薄板轴对称大挠度弯曲的位移型控 制方程，并利用伽辽金方法求出了圆薄板大挠度弯曲的近似位移解；针对一个简例进行了理 论计算和有限元数值模拟，二者结果一致. 在此基础上分析了圆薄板大挠度弯曲的规律及材 料参数随温度变化的热软化效应，给出了相关结果.     

轴对称圆薄板大挠度微分方程的数值分析方法

根据轴对称问题的特点,利用级数展开和求极限法则,证明了轴对称大挠度圆薄板在圆心处应满足的边界条件,并以圆薄板轴对称大挠度弯曲变形微分方程为基础,建立了圆心处非奇异的轴对称大挠度圆板弯曲微分方程,从而可以方便地利用现有的常微分方程数值求解方法(如变步长龙格-库塔法)对实心圆板的轴对称问题进行数值求解,又不必像摄动法那样推导复杂的公式。在数值求解轴对称圆板大挠度弯曲变形微分方程时,将非线性微分方程的求解主要归结为迭代求解圆心处三个未知边界条件的问题,即圆心处的径向膜力、圆心处的挠度、圆心处挠度的二阶导数,并提出了相应的求解方法。实例中,对于圆薄板受均布横向荷载的问题,分析了周边固支边界条件下的非线性弯曲问题,给出了中心挠度参数大范围变化时的荷载和部分边界值变化曲线,并与经典摄动解进行了对比。对比结果可见,本文方法和摄动法的解非常接近,在量纲归一化中心挠度不超过4.0时,两种方法解的相对误差均小于5.0%。另外,本文还分析了与挠度有关的液体压力作用下和集中荷载作用下周边固支圆板的非线性弯曲问题。通过算例可见:本文方法可以灵活处理不同的荷载问题;对于不同的问题,计算过程相似,不必推导复杂的计算公式,计算精度容易控制。     

轴对称变形的弹性圆板和圆扁壳问题的联合解法

复杂载荷作用下圆板和圆扁壳的大挠度和稳定问题,目前还缺少工程上的有效解法。下面我们给出,在定常温度和轴对称径向力、横向力的作用下,弹性圆板和圆扁壳问题的差分法与摄动法(或逐次逼近法)的联合解法。     

复合载荷作用下圆板大挠度弯曲问题的伽辽金法

1.引言 关于复合载荷作用下圆板的大挠度弯曲问题,文献[1,2]曾分别选取不同的参数用摄动方法给予求解。本文用一个简便的方法来分析此问题。即先假设一个挠度试函数,使相容方程完全满足,求出薄膜力;然后再用伽辽金加权残数法求解平衡微分方程。 已知均布荷载及中心集中力联合作用下圆板的大挠度方程为     

功能梯度与均匀圆板弯曲解的线性转换关系

基于一阶剪切理论, 研究了功能梯度材料圆板与均匀圆板轴对称弯曲解之间的线性转换关系. 通过理论分析和比较 功能梯度材料圆板和均匀圆板在一阶剪切理论下的位移形式的轴对称弯曲控制方程, 发现了功能梯度材料圆板的转角与均匀圆板的转角之间的相似转换关系. 在假设材料性质沿板厚连续变化的情况下, 给出了相似转换系数的解析表达式. 在此基础上, 进一步导出了一阶剪切理论下功能梯度圆板的挠度与经典理论下, 均匀圆板的挠度之间的线性关系. 从而, 可将功能梯度材料圆板在一阶剪切理论下的弯曲问题求解, 转化为相应均匀薄圆板在经典理论下的弯曲问题求解, 以及转换系数的计算问题. 这一方法为功能梯度非均匀中厚度圆板的求解提供了简捷途径, 而且更便于工程应用. 作为例子, 采用上述方法分别求得了周边简支和夹紧条件下, 梯度圆板在均布横向载荷作用下的弯曲解析解, 该解答与Reddy得到的结果完全吻合.      

Winkler地基上变厚度圆板的轴对称弯曲

本文提出了Winkler地基上变厚度圆板轴对称弯曲的传递矩阵算法。首先,根据贝塞尔函数理论获得了等厚度圆板和环板单元在任意荷载作用下轴对称弯曲的解析解,这些解均由通解和特解两部分组成。基于这些解析解,导出了等厚度圆板和环板单元的传递矩阵。然后沿径向将变厚度圆板划分成一个等厚度圆板单元和一系列等厚度环板单元,应用传递矩阵算法原理获得了变厚度圆板的整体传递矩阵。引入圆板的边界条件,给出了该板每条节线上的挠度、径向转角、径向弯矩和径向剪力。最后,讨论了受均布荷载作用的简支线性变厚度圆板的弯曲,将本文数值解与解析解进行比较,证实了本文方法的有效性,并简要地讨论了地基参数对板挠度和径向弯矩的影响。     

THERMAL BUCKLING ANALYSIS OF CIRCULAR PLATE EMBEDDED IN AN INFINITE ELASTIC PLATE 1)

分析了嵌入无限大弹性板中的圆板在变温时的热屈曲问题。由于圆板的热膨胀系数与无限大弹性板的热膨胀系数不同，温度变化时圆板中会产生压应力。当压应力达到其临界值时，圆板会发生热屈曲。首先，基于弹性力学平面应力问题的基本理论，得到圆板和无限大弹性板的应力和位移；然后建立圆板热屈曲的控制微分方程，求得临界屈曲温度的解析解和数值解，着重讨论圆板和无限大弹性板的材料物性参数的关系对圆板临界屈曲温度的影响。     

非均匀圆柱型正交各向异性圆板的弯曲

研究了非均匀圆柱型正交各向异性圆板在均布横向载荷作用下的弯曲问题，求得了折算刚度随半径按指数规律变化的非均匀圆柱型正交各向异性圆板弯曲问题的渐近解，给出了周边固支和简支条件下的渐近解．通过算例可以看出，这种非均匀性对圆板中心挠度的影响是显著的     

粘弹性大挠度圆板的轴对称弯曲

本文探讨粘弹性大挠度圆板的轴对称弯曲的基本方程和求解方法.用半逆解和摄动法分析挠度与膜力,对标准线性固体进行数例计算,并与小挠度理论相比较.全部方程与解答可退化得相应的弹性大挠度板的结果.     

考虑地基水平摩阻的Winkler地基上板模型解析

在Winkler地基模型基础上,采用板与地基之间的水平摩阻应力正比于板与地基接触点间水平位移假设,推演了包含板弯曲和轴向拉压效应的非线性地基板微分方程;应用摄动法分析竖向集中力作用下,板与地基水平摩阻应力对板内力（弯曲力、轴向力和剪应力）、位移（挠度、水平位移）的影响规律,实证了板与地基水平摩阻应力对板剪应力的非线性影响可以忽略不计,从而得到了轴对称条件下计入地基水平摩阻的Winkler地基上薄板的线性微分方程．随后通过汉克尔变换得到了轴对称无限大板的解析解,以及包含四项复宗量贝塞尔函数的轴对称圆形板的解析解,给出了板内力、位移的计算式;最后,通过算例分析了板与地基间水平摩阻状况对板挠度、截面弯矩的影响规律．结果表明：地基板挠度与弯矩随着板与地基间水平摩阻的增大而减少;当地基板水平摩阻参数大于0.01时,地基水平摩阻力对板挠度和弯矩影响有可能超过2%,应予考虑;无限大板作用圆形均布荷载时,水平摩阻的存在最大可使板最大挠度（即荷载圆中心点处挠度）下降约50%,板最大截面弯矩（即荷载圆中心点处弯矩）下降约70%.     

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

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

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

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