对流项占优问题的MLPG/SUPG方法数值模拟

在计算对流项占优问题时易产生假扩散,本文把流线型迎风格式应用于MLPG方法中可以减少对流项的影响,通过两个典型例子（旋转流场问题和Brezzi问题）验证该格式的精度与有效性,并与文献中的迎风格式的计算结果进行比较,计算结果表明,该方法能有效地克服假扩散现象,有较好的稳定性和较高的计算精度。     

时一空守恒元解元（CE／SE）方法综述

时一空守恒元解元方法是近年来兴起的一种全新的高分辨率守恒型方程计算方法．它具有物理概念清晰、计算精度高和格式构造简单等优点,是一种具有广阔发展前景的计算方法．本文详细地介绍了cE／ss方法的基本原理、发展历史、应用情况和最新进展、并指出了当前研究的不足和发展方向．     

Nonlinear study of the dynamic behavior of a string-beam coupled system under combined excitations

In this paper,the nonlinear dynamic behavior of a string-beam coupled system subjected to external,parametric and tuned excitations is presented.The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system which are described by a set of ordinary differential equations with two degrees of freedom.The case of 1:1 internal resonance between the modes of the beam and string,and the primary and combined resonance for the beam is considered.The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system and obtain approximate solutions up to and including the second-order approximations.All resonance cases are extracted and investigated.Stability of the system is studied using frequency response equations and the phase-plane method.Numerical solutions are carried out and the results are presented graphically and discussed.The effects of the different parameters on both response and stability of the system are investigated.The reported results are compared to the available published work.     

流形元法固定边界约束处理研究

在数值流形方法中,对于材料的固定边界,一般采用罚函数的方法进行处理,即在固定边界上设置刚性弹簧约束其位移来实现固定约束条件的近似满足。罚函数法在理论上不是严格的固定约束处理方法,罚弹簧的布置与弹簧刚度的大小对模拟的效果都会产生影响。基于流形单元上位移函数的组成提出了流形方法固定边界约束处理的新方法,在组成流形单元的物理覆盖上,通过取消相应的覆盖函数在流形单元位移函数中的组成来实现双向固定的约束条件,通过使用只包含单方向位移的覆盖函数使x向固定约束条件和y向固定约束条件得到实现,推导了相应固定约束条件下的流形单元刚度矩阵的数值计算格式。该方法严格满足固定约束的物理意义,简化了固定边界的处理,并经算例证明是有效和准确的,有利于数值流形方法的程序实现和工程应用。     

海冰动力学数值模拟中改进的PIC方法

为了准确地模拟海冰的动力过程,需要建立精确有效的数值方法。本文结合质点网格法(PIC)和光滑质点流体动力学方法(SPH)发展了一种改进的PIC方法。该方法在欧拉坐标下对海冰动量方程进行差分计算,在拉格朗日坐标下进行海冰质点位移、厚度和密集度计算,并采用Gauss函数进行欧拉网格点与拉格朗日质点间海冰参数的交互插值。采用改进的PIC方法对规则区域内的海冰堆积过程进行了数值试验,对渤海海冰的动力过程进行了72小时数值模拟。计算结果均表明改进的PIC方法具有计算量小,计算结果平稳精确的优点,可很好地适用于海冰动力作用过程的数值模拟。     

三阶非线性发展方程的线性化有理逼近方法

研究了三阶非线性发展方程的初边值问题的解。采用基于Sinc函数的微分求积法发展了线性化有理逼近方法。通常的配点法不适用于上述三阶问题的求解。本文把提出的方法用于求解KdV方程,取得了良好的效果。     

预应力钢桁架结构分析中的摩擦滑移索单元

为精确模拟预应力钢桁架中连续长索在支撑点的滑动,本文创建了一种考虑摩擦力影响的新单元。被称为摩擦滑移索单元的新单元有三个节点,中间节点为支撑点。本文首先利用弹性悬链线的解析解,建立了弹性悬链线单元,并推导了单元两端点的张拉刚度。摩擦滑移索单元由两个弹性悬链线单元组合而成,根据支撑点处索的滑动方向、索力差、滑动摩擦力和滑移刚度调整两索段的原长,使支撑点两侧的索力满足给定的摩擦关系。算例验证了新单元算法的正确性和高效性,对设计的预应力钢桁架的分析,显示了张拉过程中及使用荷载作用下新单元在结构分析中的应用。新单元可直接用于常规的有限元分析,研究在工作状态或施工中存在连续长索滑移的索结构。     

Hamilton-Jacobi方程的小波Galerkin方法

本文选择Daubechies小波尺度函数空间作为Galerkin方法的测试函数空间,并将其应用于Hamilton-Jacobi方程,得到了求解Hamilton-Jacobi方程的小波Galerkin方法的数值格式．由于小波在时间和频率上的局部性,本方法适用于处理具有奇异解的问题,可以有效地防止数值振荡．数值试验显示,本方法是有效的．     

DATA PREORDERING IN GENERALIZED PAV ALGORITHM FOR MONOTONIC REGRESSION

Monotonic regression （MR） is a least distance problem with monotonicity constraints induced by a partiaily ordered data set of observations. In our recent publication [In Ser. Nonconvex Optimization and Its Applications, Springer-Verlag, （2006） 83, pp. 25-33], the Pool-Adjazent-Violators algorithm （PAV） was generalized from completely to partially ordered data sets （posets）. The new algorithm, called CPAV, is characterized by the very low computational complexity, which is of second order in the number of observations. It treats the observations in a consecutive order, and it can follow any arbitrarily chosen topological order of the poset of observations. The CPAV algorithm produces a sufficiently accurate solution to the MR problem, but the accuracy depends on the chosen topological order. Here we prove that there exists a topological order for which the resulted CPAV solution is optimal. Furthermore, we present results of extensive numerical experiments, from which we draw conclusions about the most and the least preferable topological orders.     

磁性单层膜磁学性质的Monte Carlo模拟

采用类Ising模型,利用Monte Carlo方法研究了磁性单层膜中退磁性偶极作用和铁磁性交换作用对系统磁学性质的影响.结果显示,随着偶极相互作用的增加,系统在低温下的磁化出现平台现象,此时磁化曲线可分为2个阶段,在低外场下,温度升高,系统易磁化,在高外场下则反之.这种新奇的磁化行为导致系统的磁熵变在低温低外场下出现大于零的反常行为.在模拟过程中,对长程力作用采用了比较精确的处理方法.