一种适用于简支多边形薄板问题的有限元单元法

简支多边形薄板双调和定解问题分解为二个互不耦合的Ｐｏｉｓｓｏｎ方程的定问题，通过比拟，可用二维平面应力单元进行有限元分析。数值算例表明，本文采用的方法在较粗的网格下即可得到较高的精度，而且具有速度快，占用计算机内存少等优点。     

平面Poisenille流动的二维Hopf分叉

本文用直接数值模拟的方法计算了二维Ｐｏｉｓｅｕｉｌｌｅ流动中扰动波的演化问题。得到了二维平衡态，在一定的波数下，Ｒｅ３９５０时，这种平衡态，将变得不稳定，模拟发现出现第二周期解，即二次分叉。     

平面Poisenille流动的二维Hopf分叉

本文用直接数值模拟的方法计算了二维Ｐｏｉｓｅｕｉｌｌｅ流动中扰动波的演化问题。得到了二维平衡态，在一定的波数下，Ｒｅ３９５０时，这种平衡态，将变得不稳定，模拟发现出现第二周期解，即二次分叉。     

位移模式对Cross－PLY矩形板热弹性响应的影响

本文建立了以Ｒｅｉｓｓｎｅｒ混合变分原理为基础的迭层板剪切变形理论来求解Ｃｒｏｓｓ－ＰＬＹ矩形板的热弯曲问题，用二种方式来模拟平面位移沿厚度的变化：将一阶剪切变形理论的平面位移叠加一项交错线性函数；设平面住移是分段线性连续函数。为保证层间应力的连续性，将横向剪切应力处理成厚度坐标的二次函数，温度沿厚度线性变化。通过对称及反对称Ｃｒｏｓｓ－ＰＬＹ矩形板的热应力及挠度分析比较，指出了位移模式对热弹性响应的影响。     

位移模式对Cross—PLY矩形板热弹性响应的影响

本文建立了以Ｒｅｉｓｓｎｅｒ混合变分原理为基础的迭层板剪切变形理论来求解Ｃｒｏｓｓ－ＰＬＹ矩形板的热弯曲问题，用二种方式来模拟平面位移沿厚度的变化：将一阶剪切变形理论的平面位移叠加一项交错线性函数；设平面位移是分线性连续函数。为保证层间应力的连续性，将横向剪切应力处理成厚度坐标的二次函数，温度沿厚度线性变化。通过对称及反对称Ｃｒｏｓｓ－ＰＬＹ矩形板的热应力及挠度分析比较，指出了位移模式对热弹性响应     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

一类新变量的Reissner板弯曲边界元法

文［４］导出了二维弹性力学平面问题的一类新型边界积分方程，本文将该理论和方法推广到三变量的Ｒｅｉｓｓｎｅｒ板弯曲中，给出边界场变量含广义位移和新型广义力的边界积分方程。从而边界弯矩应力张量可直接由离散边界积分方程求出。     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献〔２，３〕的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几     

天然气管道裂纹稳定扩展问题的分析模型和最小二乘法解答

应用最小二乘法求解了模拟天然气管道裂纹稳定扩展问题的等效染模型，。在Ｗｉｌｌｉａｍｓ等效梁模型的基础上，通过管道裂纹后部截面开裂后变形协调条件。补充了弹性地基对于等效梁的作用，推导出梁的分析模型，并应用最小二乘法求得了裂纹驱动力和位移的解答。其结果与Ｌｅｅｖｅｒｓ所给出的不考虑弹性地基对于梁的作用结果比较，改进了梁模型的计算精度，使裂纹驱动力成为合理的有限值，与已有的试验结果吻合较好。     

血液入口段两相流动的层流NS方程数值模拟

本文采用球形刚性颗粒悬浮于牛顿流体的血液模型，基于有相间滑移的两相流动层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，采用相间滑移算法（Ｉｎｔｅｒ－ＰｈａｓｅＳｌｉｐＡｌｇｏｒｉｔｈｍ）对圆管内定常轴对称血液入口两相流动问题进行了计算，获得了与实验结果吻合的计算结果，较好地模拟血液流动中红细胞的径向迁移现象，结果表明，采用二相流动模型研究血液流动是一种有前途的方法。     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

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

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

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

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