一种求解三维热辐射输运方程的整体预处理迭代方法及并行计算

为三维灰体热辐射输运方程的隐式离散纵标方法发展一个整体预处理迭代方法并研制并行程序。该方法采用组装线性代数方程组策略,同时求出所有离散方向上的辐射强度。借助预处理的Krylov子空间迭代法,避免复杂网格上扫描算法可能遇到的死锁问题,能够提高健壮性和计算效率。空间离散上采用一阶迎风有限体积格式。数值实验测试变形六面体网格上的收敛率、评估预处理迭代方法的性能并计算辐射和物质的耦合问题,给出三维弯管和黑腔问题的模拟结果,验证程序的正确性和方法的适应性。     

三维有限体中平片裂纹的超奇异积分方程方法—Ⅱ.数值方法

超奇异积分方程方法的理论分析已在本文的第Ｉ部分中给出，这一部分是经的数值方法，及用此方法求解的若干典型的平片裂纹问题。     

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

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

沿板平面材料梯度变化功能材料板件分析

本文发展一种功能梯度构件分析的细观元法.细观元法在构件的常规有限元内部设置密集观细单元以反映材料特性变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算.这种细观元法既能充分反映材料功能梯度沿各方向任意变化特性;而其计算单元又和常规有限元一样,是一种针对功能梯度结构分析的有效数值方法.现有功能梯度板件分析中无论对不同形状还是不同边界的功能梯度构件,其材料特性均沿板厚度方向梯度变化,本文用细观元法进行计算与分析,给出了目前尚未得到的沿板平面方向功能梯度变化构件的力学量三维分布形态.     

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

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

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

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

考虑法向压力梯度的三维旋转边界层计算

本文以量级分析为基础,建立了一般曲线坐标系上的三维旋转边界层方程。对旋转在边界层中的影响进行分析之后,提出了一个能够处理壁面法向压力梯度不为零问题的压力梯度迭代方法。在传统的Box法的基础上发展了一套完整的求解三维旋转边界层的方法和程序,并对螺旋面、压气机转子叶面以及圆柱面上的旋转边界层进行了计算,与他人的计算和实验的对比分析表明,该方法和程序是正确的,可用于求解任意几何物面上的三维旋转边界层。     

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) 弹粘塑性有限元与弹性边界元的耦合及数值计算稳定性条件。