脉动流条件下血管壁的应力分布

为了分析血液-血管耦合运动所产生血液脉动压力载荷对血管壁应力分布的影响,利用线性化的血液-血管耦合运动方程的Womersley解,导得血液脉动压力载荷下的血管壁Green应变,同时利用Fung的血管壁应变能密度函数,导得相应血管壁应力分布的一般表达式.数值结果表明,在脉动流条件下,当考虑血液-血管耦合运动时,血管壁中周向应力最大,轴向应力居中,径向应力最小;血管壁的残余应力将明显减小血管内壁的应力集中;脉动压力载荷将导致血管壁周向应力在一个心动周期中随时间的脉动,而且随着Womersley数α和血管轴向约束参数K*的增大,血管壁周向应力的脉动将明显加剧,提示在分析动脉重建时必须计及血液-血管耦合运动对血管壁应力分布的影响.     

用p-V指数关系确定血管壁的周向应力

通过检测，得到在体轴向伸长比条件下血管段的压力－容积（p-V)数据，利用指数函数关系，求得在一定内压作用下血管壁周向应力沿壁厚的分布。并以正常大鼠腹主动脉为例进行计算，结果表明，在生理压力作用下，残余应力对血管壁周向应力的分布有很大影响，不考虑残余应力将使血管内壁处出现明显的周向应力集中，当考虑存在于血管壁内的残余应力后，周向应力沿壁厚的分布将变得较均匀，例如当内压p=10kPa时，血管内、外壁周向应力值之比，对前者将达10以上，而对后者差不多只有2左右。所得结果与用Fung提出的应变能密度函数的计算结果作比较，发现二者符合得相当好。     

脉动血压对在体轴向极限强约束血管壁应力分布的影响

本文详细分析正弦脉动压力载荷对动脉管壁应力分布的影响。结果表明,对某确定时刻,包含在运动方程中的惯性项对血管壁应力的影响很小,可忽略不计;在体轴向极限强约束下,血管壁周向应力随时间的变化曲线与压力曲线同相位和同周期,而且振幅不随频率的快慢而改变。一个心动周期内,血管内、外壁处周向应力随时间变化的平均值只与脉动压力载荷的平均值有关,而振幅主要与脉压有关,在内壁处随平均压的改变有一定变化,但在外壁处受平均压的影响很小。     

在体轴向伸长比条件下血管壁的周向应力分布

本文通过检测在体轴向伸长比条件下血管段的压力－容积数据,利用应变能密度函数,求得血管壁周向应力沿壁厚的分布。文章的正常大鼠腹主动脉为例进行计算,结果发现,在生理压力作用下（例如ｐ＝１３．３ｋＰａ）,不考虑残余应力确实出现在血管内壁处的周向应力集中,血管内、外壁处的应力值之比可达１５左右,而存在于实际动脉管中的残余应力会使血管壁中的周向应力沿壁厚的分布明显趋于均匀,血管处的周向应力之比差不多只有２左     

血管残余应力的一种确定方法

动脉无载荷状态下存在残余应力,分析血管壁在载荷状态下的应力分布必须计及血管壁的残余应力。本文通过检测动脉不同轴向伸长比条件下的压力—容积(p—V)数据,根据其p—V实验曲线的“S”型特征,提出了一种用反正切函数关系对试验数据点进行拟合,进而导得动脉管壁周向残余应力沿壁厚分布的方法。文章对10只正常大鼠颈动脉的分析结果表明,大鼠颈动脉周向残余应力沿血管壁厚呈单调上升趋势。残余应力在内壁处为负呈压缩状态,在外壁处为正呈拉伸状态,并且内外壁处应力绝对值大小基本相等,大约在2．5kPa左右。文章结果对进一步分析动脉管壁在载荷作用下的应力分布提供必要的信息。     

动脉压对静脉血管壁周向应力分布的影响

自体静脉是病变动脉管段常用的替代物。移植后因承受压力急剧升高引起的静脉管壁应力改变是影响移植手术的主要因素之一。为了比较移植前(静脉压作用下)和移植后(动脉压作用下)静脉管壁的周向应力分布，本文通过检测一定轴向伸长比条件下静脉管的p(压力)——V(容积)试验数据，利用3参数的应变能密度函数对实验数据进行拟合，进而求得静脉管壁的残余应力和沿血管壁的周向应力分布。对狗的股静脉和颈静脉的分析结果表明，在动脉压作用下静脉管壁周向应力将急剧增大。与处于静脉压环境相比，处于动脉压环境中的颈静脉管周向应力将增大差不多2个数量级。计算结果还显示，静脉管壁残余应力的数值虽然比动脉管壁的相应值小很多，但是与动脉管相同，血管壁残余应力依然对静脉管壁上的周向应力分布影响显著，残余应力的存在将大大削弱在静脉管内壁处的周向应力集中，使周向应力沿静脉管壁厚的变化梯度明显减小。     

局部扩张动脉管血液振荡流的切应力分布

本文求解局部缓慢扩张动脉管中血液振荡流的基本方程，得到血管内血液的流速与压力梯度的关系。通过导出压力梯度沿局部扩张管轴向的变化特性。建立利用扩张段上游血管均匀段中心流速波形确定局部扩张管中血液流的速度和切应力分布的方法，文章以人体颈动脉余弦扩张为例进行分析。详细讨论了局部扩张对血管壁切应力及其梯度分布的影响。数值结果表明，在与刚性均匀管中管壁切应力沿轴向保持不变不同，在局部扩张段，管壁切应力将随着血管半径的增大而减小，因而管壁切应力梯度一般不为零，甚至在某些位置达到相当大的数值。另外，随着血管扩张程度的增加，管壁切应力还将进一步减小，而且管壁切应力梯度也将进一步增大，血管扩张导致管壁切应力的这些变化将直接影响血管壁的结构和功能，使其产生适应性的变化。     

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

分级加卸载硬岩短时蠕变特性实验研究

由单向抗压强度实验及8级加卸载短时蠕变实验,得到蠕变下限为27MPa的细砂岩试样典型的坚硬岩石脆性断裂特征.应力-应变等时曲线线性回归函数的相关系数均高于0.92,长期强度与瞬时强度之比达94.39%,证明细砂岩试样的整体蠕变特性不强.应力-轴向应变等时曲线线性回归函数的平均相关系数高出应力-径向应变等时曲线线性回归函数的平均相关系数3.92%,因此其轴向非线性蠕变特性相对于径向非线性蠕变特性更弱.随着时间延续,细砂岩试样的非线性蠕变特性总体服从负Gauss分布规律,具有明显的时间效应.随加载应力水平的提高,应力-轴向应变等时曲线和应力-径向应变等时曲线线性拟合函数的平均相关系数下降幅度分别为0.97%和0.67%,线性相关性普遍降低.所以细砂岩试样的非线性蠕变程度随加载应力水平的提高而提高,载荷效应明显.     

考虑轴向力和渗透力时软化围岩隧道解析

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明:在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.      

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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