开口厚壁截面短构件的约束扭转

为了分析开口厚壁截面短构件的约束扭转问题,采用统一分析梁模型与有限节线法,对T形和L形厚壁截面短构件约束扭转时横截面的翘曲和应力分布情况等问题进行了分析研究.算例计算结果表明:开口厚壁截面短构件存在与其横截面形心位置不一致的扭转(弯曲)中心,构件在不过扭转中心的外力作用下会产生弯扭耦合变形,其横截面将产生不均匀翘曲,横截面上的翘曲正应力和扭转剪应力均呈非线性分布.     

开口厚壁截面短构件的约束扭转1）

为了分析开口厚壁截面短构件的约束扭转问题,采用统一分析梁模型与有限节线法,对T形和L形厚壁截面短构件约束扭转时横截面的翘曲和应力分布情况等问题进行了分析研究.算例计算结果表明：开口厚壁截面短构件存在与其横截面形心位置不一致的扭转（弯曲）中心,构件在不过扭转中心的外力作用下会产生弯扭耦合变形,其横截面将产生不均匀翘曲,横截面上的翘曲正应力和扭转剪应力均呈非线性分布.     

带裂纹角钢形截面杆抗扭附度和第三型应力强度因子的计算

本文是作者前一工作的继续,文中讨论了开裂角钢形截面抗扭刚度和第三型应力强度因子的计算方法.对于图1所示截面,若令扭转问题中的应力函数为φ(x,y)=-x~2 u(x,y) (1)不难导出对于函数u(x,y)的定解问题为((?)~2u/(?)x~2) ((?)~2u/(?)y~2)=0,u|L=x~2 (2)其中L 表示截面的周边.同时抗扭刚度为D=μJ J=2(?)(-x~2 u(x,y))dxdy (3)     

扇形截面杆扭转问题的差分线法

导出了扇形截面杆扭转问题偏微分方程的差分线法常微分方程组, 并解析求解了该方程组, 得到了扭转应力函数的半解析解, 计算了扭转应力及扭转刚度. 计算过程中, 用追赶法计算 常微分方程组的特解, 用公式计算三对角矩阵的特征值与特征向量, 利用实对阵矩阵的特征 向量相互正交的特性避免矩阵求逆计算, 利用复化梯形公式计算扭转刚度. 整个求解过程在 角度方向离散微分方程和用复化梯形公式进行面积积分时引入了误差, 其他求解过程是精确 的. 计算结果与已有结果进行了对比, 显示了算法的正确性. 该算法对工程中扇形截面扭 转杆的设计有一定的实用价值.     

曲面拟合在平面等达因法数据处理中的应用

1 前言平面等达因法(Insodyne method)是散光法的一种发展,能得到全场的清晰条纹.与一般由片光得到的散光条纹图相比,物体的形状和条纹都不发生畸变,两幅正交方向的等达因条纹图,就是分解受力模型内各点的三个应力分量,但必须通过求导的数据处理.通常的作法,对从实验中得到的离散的等达因条纹数据,沿某一截面,用曲线拟合后再求导,以得到某一截面上的应力分量,如果需要另一方向的导数,需另外选择截面,再用曲线拟合,再求导.事实上等达因条纹级次的分布是(x,y)坐标的函数,对应的是一个曲面,在进行应力分析时,同时需要x 和y 方向的导数,因此在对实验中测得的离散数据,用     

TORSIONAL CENTER OF BEAMS OF ARBITRARY COMPLEX CROSS-SECTION 1)

为了计算任意复杂非圆截面梁横截面扭转中心的位置，用节线法将其约束受扭后所有横截面面外变形的形状用一族包含节线未知函数的曲面表示，建立梁约束受扭时的控制方程后，再用常微分方程求解器分别求出单纯扭矩与横向载荷单独作用时节线未知函数的数值解，最后用刚度等效原理导出复杂截面梁横截面扭转中心的位置。算例计算结果表明：该方法是合理的、有效的，是计算任意复杂非圆截面梁横截面扭转中心位置的可靠方法。     

任意复杂截面梁的扭转中心

为了计算任意复杂非圆截面梁横截面扭转中心的位置,用节线法将其约束受扭后所有横截面面外变形的形状用一族包含节线未知函数的曲面表示,建立梁约束受扭时的控制方程后,再用常微分方程求解器分别求出单纯扭矩与横向载荷单独作用时节线未知函数的数值解,最后用刚度等效原理导出复杂截面梁横截面扭转中心的位置。算例计算结果表明:该方法是合理的、有效的,是计算任意复杂非圆截面梁横截面扭转中心位置的可靠方法。     

薄壁弹性梁在截面突变处的连接条件

本文对开口薄壁弹性直梁在截面突变处两侧截面,利用扇性坐标转换公式,提出了在一个截面内部分截面的主扇性坐标的概念,建立了轴向位移连续条件;据此,得到了在截面突变处轴向位移、弯曲和扭转耦合的一般的连接条件(位移连续条件和力平衡条件)。     

一种快速收敛的简单矩形板弯曲元素

1.引言文献[1]中对具有12个节点位移参数的简单矩形板弯曲元素,采用三个特定的变形假设,构造了两个分别沿元素局部坐标x、y方向的位移插值函数.这两个位移插值函数都完整到二次多项式,并包含有坐标x、y高达三——六次的项,这在同     

变截面圆轴扭转问题用非正交曲线坐标的新解法

本文是在飞机发动机涡轮轴应力分析的基础上提高为变截面圆轴扭转问题的一个新解法.利用向量的散度和旋度对不同坐标系是不变量的特点,通过张量分析推导出在任意非正交曲线坐标系中变截面圆轴扭转问题的平衡和协调方程,包括用应力函数表达的协调方程和应力函数与应力分量的关系式.用任意非正交曲线坐标和差分法求解应力函数.本文计算得到了全轴的等应力函数线和剪应力分布,并得到沿小凹槽边任意点的应力,计算结果和光弹试验结果接近.本文还计算了有解析解的空心锥轴,误差小于百分之一.通过计算说明本文提出的新解法收敛性很好,并且所需计算机容量少(可在容量32k的TQ-16机上同时计算800多个节点),应用方便,便于编排通用程序,计算量较有限元法少;另一方面,由于采用了任意非正交曲线坐标,因此,适用于解决复杂曲线边界的问题,提高了通常用的差分法的适应性和灵活性;此外,本方法用应力函数作为未知量,从所得的等应力函数线和等位移函数线图可以看出全轴应力分布的概况,对于改进设计很有帮助.     

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

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

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

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