不同围压作用下非均匀岩石水压致裂过程的数值模拟

从岩石细观非均匀性的特点出发，提出一个描述非均匀材料渗流和破裂相互作用的数值模型。在这个数值模型中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流，损伤耦合迭代来实现。算例表明，该模型能较好地模拟出岩石类材料在水力压裂作用下，微结构非均匀分布和不同围压比对破裂模式、失稳压力的影响，非均匀性导致试件的开裂压力、失稳压力明显不同，裂纹扩展路径不规则发展，模拟结果和实验结果较为一致。     

汶川震区北川9.24暴雨泥石流特征研究

2008年9月24日汶川震区的北川县暴雨导致区域性泥石流发生,这次9.24暴雨泥石流灾害导致了42人死亡,对公路和其他基础设施造成严重损毁。本研究采用地面调查和遥感解译方法分析地震与暴雨共同作用下的泥石流特征,获取的气象数据用于分析泥石流起动的临界雨量条件。本文探讨了研究区泥石流起动和输移过程,并根据野外调查,分析了泥石流形成的降雨、岩石和断层作用,特别是强降雨过程与物源区对泥石流发生的作用。根据应急调查发现北川县境内暴雨诱发的泥石流72处,其分布受岩石类型、发震断层和河流等因素控制。根据对研究区震前和震后泥石流发生的临界雨量和雨强的初步分析,汶川地震后,该区域泥石流起动的前期累积雨量降低了14.8%～22.1%,小时雨强降低25.4 %～31.6%。震区泥石流起动方式主要有二种,一是由于暴雨过程形成的斜坡表层径流导致悬挂于斜坡上的滑坡体表面和前缘松散物质向下输移,进入沟道后转为泥石流过程;二是"消防水管效应"使沟道水流快速集中,并强烈冲刷沟床中松散固体物质,导致沟床物质起动并形成泥石流过程。调查和分析发现沟内堆积的滑坡坝对泥石流的阻塞明显,溃决后可导致瞬时洪峰流量特别大。研究结果表明了汶川震区已进入一个新的活跃期。因此,应该开展对汶川地震区的泥石流风险评估和监测、早期预警,采取有效的工程措施控制泥石流的发生和危害。     

锡基巴氏合金磨损表面的分形与磨损率

利用结构函数分析法研究了滑动摩擦学系统中金属磨损表面轮廓线的分形特性。结果表明：磨损表面轮廓线在小于Ｓｍ的尺度上具有分形结构。采用结构函数法可以方便地确定粗糙表面轮廓线的分形参数，即垂直于滑动方向上磨损表面轮廓线的分维Ｄ可作为磨损表面分维的特征值，它与金属磨损率的变化有着密切的关系，最佳分维Ｄｏｐｔ值对应于材料的最低磨损率。     

架设阶段悬索桥静力问题的随机有限元分析

架设阶段的悬索桥不仅具有显著的几何非线性行为，而且存在材料特性的变异。本文建立了以增量形式表达的非线性摄动法随机有限元列式，并编制了悬索桥静力分析程序ＳＮＡＰ。针对一座国内在建的悬索桥，计算了该桥在某一施工阶段考虑缆索弹性模量随机性时主梁跨中点位移的变化。结果表明，材料随机性的影响十分明显。     

径向轴承启动过程瞬态热效应的研究

联立求解广义雷诺方程、油膜瞬态三维能量方程、轴瓦瞬态三维固体热传导方程及轴颈的运动方程，并考虑粘度和密度随温度及压力的变化，在油膜与轴瓦界面使用热流连续性边界条件，确定了在承受稳态载荷时圆轴承在启动过程中的温度变化情况．介绍了一种有效的用于求解轴承瞬态性能的改进Newton—Raphson算法．结果表明：整个系统达到热平衡的时间与启动加速度的大小无关；启动加速度越高，温度升高越快；在启动过程中，油膜最高温度出现的位置沿膜厚方向而变化．同有关试验结果的对比证实所建模型和所用算法正确可行．     

Study on failure process of delaminated stiffened composite plates under compression

Failure behavior of the delaminated stiffened composite plates under compression is studied by the finite element method, based on a Global-Local variational model. A virtual crack closure technique and a self-adaptive grid moving scheme are proposed to predict the delamination growth process. The contact effect along the delamination front is considered. The numerical results show that the influences of the distribution and location of the stiffeners, the configuration and size of the delamination, the boundary condition and the contact upon the failure behavior of the plates are significant. The project supported by the National Natural Science Foundation of China (59975014)     

飞秒激光辐射纳米薄膜脱附的能量输运机理研究

从微观能量传输角度,对飞秒脉冲激光的辐射传热及其产生的金属薄膜吸附层分子脱附过程的能量输运机理进行了研究。以非平衡态传热和随机过程理论为基础,建立了描述吸附分子脱附过程的数学模型;并分别对单脉冲与双脉冲飞秒激光辐射下,金属Ru基体表面层的电子和声子的温度分布、吸附层分子CO的脱附概率分布进行了数值模拟。其结果与实验规律相符,为飞秒脉冲激光应用于微机电系统中杂质脱附提供了重要信息。     

On damping coefficient due to phase transformation

The damping coefficient of capillary waves due to the evaporation-condensation process at the interface of the two phases of a fluid is evaluated. To highlight the mechanism of the effect of heat and mass transfer across the interface between regions of liquid and vapor, potential flow of incompressible fluids are assumed. Thus other mechanisms of damping are neglected. To fascilitate the analysis, the method of multiple-scale is employed in the analysis, even though the problem is linear.     

桁架结构稳定分析的几何非线性欧拉稳定理论

目前,桁架结构稳定分析有两个理论:经典的特征值理论和几何非线性临界点理论。前者发展较早,在许多结构力学书中都有叙述。但是经过工程应用检验,它过高地估计了结构抗稳定的能力。后者是20年前提出的,提出者认为适用于所有的扁桁架。经过几年的研究,作者连续提出三个稳定理论:线性与非线性欧拉理论;非线性临界点-欧拉理论和两个计算方法,进行桁架结构的线性和非线性截面优化设计。经过理论研究和对前人文献中数值例题的计算和比较,发现特征值理论不符合实际,临界点理论只适用于大扁度的桁架,而不适用于一切扁桁架。研究结果提供了若干有用的结论,给出了各种稳定理论的适用范围,指出了国际文献中若干例题的正确和错误,补充了桁架结构的稳定分析理论,使之更为完整和实用了。     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法，它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点，该方法被有效应用于充填流动过程的数值模拟.