“吕”字形巷道布置的卸压机理分析

利用有限元法分析了综放采场瓦斯巷和回风巷或运输巷采用“吕”字形布置后的卸压机理，得到了上部瓦斯巷宽度对下部回风巷或运输巷卸压效果的定量关系，为改革综放采场巷道布置并提高瓦斯巷抽放效果提供了新的设计思路和可靠依据。     

某边坡变形破坏机理的数值分析

某花岗岩矿岩土体主要有填土、粉质粘土、淤泥及淤泥质土、含泥中粗砂、粘土、残积砂质粘土,其下为花岗岩基岩。边坡最大高度75m,其中残积土层以上的土质边坡高度约50m。本文运用有限元强度折减法对该边坡的稳定性进行了研究,得出了潜在的破坏模式及各个滑动部分的安全系数。分析结果表明,本边坡存在3个潜在滑动面,其中两个为浅层滑动面,另一个为深层滑动面。两个浅层滑动面的稳定系数在1.13～1.17,处于临界状态。深层滑动面稳定性系数在1.29,目前尚处在相对稳定状态。根据分析结果,提出了如下加固措施:（1）对滑坡体进行削坡处理;（2）在其滑动面下方设置抗滑台阶;（3）采用预应力锚索框架梁结合肋柱进行加固处理。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

压电传感器探测层合板分层的数值分析

用修正的准三维有限元法分析了受单轴拉伸表面覆盖有压电薄膜的复合材料层合板。文中首先讨论了压电薄膜层对层合板层间应力的影响，然后在存在分层裂纹的情况下，计算了压电层的电压分布。结果表明将压电薄膜作为传感器来进行分层探测是可行的。文中最后考虑了多种因素对传感灵敏度的影响。     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

三维非线性有限元与弹性边界元耦合数值方法

本文系统地讨论了以下三个问题:(1) 有限元与边界元耦合中的几个数值问题,其中包括:边界积分方程的凝聚、等效刚度矩阵的对称化及面力不连续的处理;(2) 弹塑性有限元与弹性边界元的耦合;(3) 弹粘塑性有限元与弹性边界元的耦合及数值计算稳定性条件。     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

用有限元自适网格控制结构应力分析的精度

本文涉及有限元结构应务分析精度控制的三个基本方面：１）自适应的有限元网格自动生成技术；２）有限元应力分析的误差估计：３）用单元尺寸场将误差分布用于形成新网格。     

物理参数变化对短脉冲激光激励温度场的影响

 为研究多物理参数（耦合系数、电子热导率、电子热容、晶格热容）同时随温度变化对短脉冲激光辐照金属材料产生温度场分布的影响,基于双温耦合理论,建立了短脉冲激光辐照金属材料金的加热过程的有限元求解模型。在同时考虑脉冲激光的空间、时间分布和多参数同时随温度变化的情况下,得到短脉冲激光辐照金属材料金激励产生的温度场二维瞬态分布,并进一步比较了多物理参数同时随温度变化和采用室温物理参数两种情况下温度场分布的区别。数值结果表明：多物理参数同时随温度变化使电子温度和晶格温度的上升变快,最大值变大,而且使得材料中激光穿透直接辐照到的区域温度变高。     

Computational performance of Free Mesh Method applied to continuum mechanics problems

The free mesh method (FMM) is a kind of the meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, or a node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm. The aim of the present paper is to review some unique numerical solutions of fluid and solid mechanics by employing FMM as well as the Enriched Free Mesh Method (EFMM), which is a new version of FMM, including compressible flow and sounding mechanism in air-reed instruments as applications to fluid mechanics, and automatic remeshing for slow crack growth, dynamic behavior of solid as well as large-scale Eigen-frequency of engine block as applications to solid mechanics.