横观各向同性轴对称问题的加权残数法及其工程应用

利用文[1,5]的理论,本文提出用两个调和函数的方法求解横观各向同性的轴对称问题,将它引入加权残数法,推得了一系列简洁的边界型最小二乘配法公式,就结构工程和岩土工程的两个例题作了计算,得到了有用的结论。本文的方法比列赫尼茨基的重调和函数法更为简便。     

非线性马休方程的亚谐分叉解及欧拉动弯曲问题

本文利用动力系统分叉理论,给了非线性马休方程的一种新的解法,并得到了整个系统参数平面上的不同参数域中分叉图各种可能的拓扑结构,利用本文所给的方法研究了欧拉动弯曲问题,得到了一些新的结果。     

超越摄动:同伦分析方法基本思想及其应用

介绍一种新的、求解强非线性问题解析近似的一般方法------同伦分析方法.该方法从根本上克服了摄动理论对小参数的过分依赖, 其有效性与所研究的非线性问题是否含有小参数无关, 因此,适用范围广.此外, 不同于所有其他解析近似方法,同伦分析方法提供了一个简单的途径, 确保所得到的级数解收敛, 从而获得足够精确的解析近似.而且, 不同于所有其他解析近似方法, 同伦分析方法(HAM)提供了选取基函数之自由, 从而可以选择较好的基函数, 更有效地逼近问题的解.同伦分析方法为非线性问题的解析近似求解提供了一个全新的思路, 为非线性问题(特别是不含小参数的强非线性问题)的求解开辟了一个全新的途径.简要描述同伦分析方法的基本思想, 其在非线性力学、物理、化学、生物、金融、工程和计算数学等领域的应用举例, 以及与摄动方法、Lyapunov 人工小参数法、$\delta$展开法、Adomian 分解法、同伦摄动方法之区别和联系.      

7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力．计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同．为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.     

涂层结构中温度场的边界元解

涂层结构由于其优良的物理化学性能而备受人们关注,但受其厚度尺寸的影响,涂层材料中物理量的数值计算一直是工程中的难点。边界元法分析涂层结构时,难点在于涂层子域的数值分析,边界量计算既涉及奇异积分又涉及几乎奇异积分。本文基于间接规则化边界积分方程,准确高效地计算奇异边界积分。针对计算边界量及内点物理量时涉及的几乎奇异积分,采用一类非线性变量替换法,有效地改善了被积函数的震荡特性,从而消除了积分核的几乎奇异性。通过采用二次单元逼近几何边界,使得高效准确地计算超薄的涂层结构成为可能。     

弹塑性接触问题数值计算的二次规划算法

针对弹塑性接触问题所推得的数值求解式子,运用二次规划法具体设计了算法,该算法采用有限的基底交换运算就可得到收敛的数值解,具有较好的收敛性及较小的计算工作量.工程计算算例结果表明文章所提出的接触问题的求解方法是有效的,针对结构中接触问题所建立的数值计算模型能真实反映实际工作状况是可靠的.文章中还详细给出了接触单元相关矩阵和向量的具体推导形式.     

MODE III 2-D FRACTURE ANALYSIS BY THE SCALED BOUNDARY FINITE ELEMENT METHOD

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.     

The design of characterizing-integral schemes and the application to the shallow water wave problems

In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn ’t only limited to conservation laws unlike other methods and maybe easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems.The paper shows us now to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem.     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

Model reduction using L1‐norm minimization as an application to nonlinear hyperbolic problems

We are interested in the model reduction techniques for hyperbolic problems, particularly in fluids. This paper, which is a continuation of an earlier paper of Abgrall et al, proposes a dictionary approach coupled with an L¹ minimization approach. We develop the method and analyze it in simplified 1‐dimensional cases. We show in this case that error bounds with the full model can be obtained provided that a suitable minimization approach is chosen. The capability of the algorithm is then shown on nonlinear scalar problems, 1‐dimensional unsteady fluid problems, and 2‐dimensional steady compressible problems. A short discussion on the cost of the method is also included in this paper.     

正交异性体某些断裂动力学问题

首先对具有任意自相似指数的正交异性弹性动力学问题推导了解的一般表示,给出了一般解法,然后用这一方法对若干具体问题进行求解,利用本文方法可以迅速将所论问题化为半平面上的Ｒｉｍａｎｎ－Ｈｉｌｂｅｒｔ问题,并可以相当简单地得到问题的闭合解。     

论粘性激波层理论的理论基础及其进一步发展

本文将Davis的量级分析方法改用匹配渐近展开方法,作为一级近似推导出了高超音速化学反应粘性激波层方程。证明了粘性激波层方程是NS方程在匹配渐近意义上的一级近似方程。进一步讨论了这一方程的基本假设条件。本文首次推导出的二级近似方程,是对Davis粘性激波层方程的修正。这种修正可以提高数值解的精度,有助于对问题获得更全面的了解,对进一步发展与完善高超音速钝头体绕流问题的数值求解方法起一定的作用。     

Stress singularity analysis of anisotropic multi-material wedges under antiplane shear deformation using the symplectic approach

Symplectic approach has emerged a popular tool in dealing with elasticity problems especially for those with stress singularities. However, anisotropic material problem under polar coordinate system is still a bottleneck. This paper presents a subfield method coupled with the symplectic approach to study the anisotropic material under antiplane shear deformation. Anisotropic material around wedge tip is considered to be consisted of many subfields with constant material properties which can be handled by the symplectic approach individually. In this way, approximate solutions of the stress and displacement can be obtained. Numerical examples show that the present method is very accurate and efficient for such wedge problems. Besides, this paper has extended the application of the symplectic approach and provides a new idea for wedge problems of anisotropic material.     

基于多模型的参数空间预测推断

当建立多个模型对工程结构进行数值仿真时,为了得到更可靠的预测结果,需要综合考虑模型选择不确定性和模型形式不确定性对预测结果的影响. 联合贝叶斯方法与实验数据计算不同模型的可信度,采用调节因子方法传播模型选择不确定性得到系统响应置信区间,并叠加模型形式不确定性的影响获得综合模型计算结果的置信区间,再通过插值得到关心量在预测点的置信区间. 最后通过某飞行器气动力系数的预测推断检验了该方法的可行性.     

THERMOELASTIC PROBLEMS IN THE HALF SPACE—AN APPLICATION OF THE GENERAL SOLUTION IN ELASTICITY

In this paper,some thermoelastic problems in the half space are studied by using thegeneral solutions of the elastic equations.The method presented here is extremely effectivefor the axisymmetric problems of the half space as well as the half plane problems.     

无限元与相同子结构串的一致性及其算法

讨论了无限元方法￣［１］与相同子结构串分析￣［２］的一致性，综合文［１］和［２］得到一种通用算法，依据文［３］尚可应用于ＬＱ控制问题，且提出一种直接求出转移矩阵的算法，具有简便和较高精度的优点．     

The perturbation finite element method for solving problems with nonlinear materials

The perturbation method is one of the effective methods for so-lving problems in nonlinear continuum mechanics.It has been de-veloped on the basis of the linear analytical solutions for the o-riginal problems.If a simple analytical solution cannot be ob-tained.we would encounter difficulties in applying this method tosolving certain complicated nonlinear problems.The finite ele-ment method appears to be in its turn a very useful means for sol-ving nonlinear problems,but generally it takes too much time incomputation.In the present paper a mixed approach,namely,theperturbation finite element method,is introduced,which incorpo-rates the advantages of the two above-mentioned methods and enablesus to solve more complicated nonlinear problems with great savingin computing time.Problems in the elastoplastic region have been discussed anda numerical solution for a plate with a central hole under tensionis given in this paper.     

卷积型加权残值法求解梁的动力学问题

 利用卷积型加权残值法推导出Gurtin变分原理, 并应用卷积型加 权残值法计算了简支梁的振动、冲击问题, 算例表明此方法是成功的.