圆柱分层媒质中并矢格林函数的快速计算

提出了一种适于在 PC机上计算圆柱分层媒质中格林函数的数值方法. 运用迭代方法和渐近方法分区域计算了复宗量圆柱组合函数, 得到了谱域格林函数的精确值. 用高精度的数值积分法将谱域格林函数经逆傅里叶变换, 求得了格林函数的空域结果. 计算结果与现有文献比较, 取得了良好的一致性, 证明了数值模型的正确性和精确性.     

圆柱分层煤质中并知格林函数的快速计算

提出了一种适于在PC机上计算圆柱分层煤质中和林函数的数值方法，运用迭代方法和渐近方法分区域计算了复宗量圆柱组合函数，得到了谱域和式函数的精确度，用高精度的数值积分法将谱域格林函数经逆傅里叶变换，求得出了格林函数的空域结构，计算结果与现有文献比较，取得了良好的一致性，证明了数值模型的正确性和精确性。     

黏弹性松弛函数的积分表示

讨论了松弛函数的积分表示.用Laplace变换方法得到了Maxwell模型的应力应变关系方程,方程中的松弛函数可以用Mittag-Leffler函数表示.由于大的负变元的存在,计算起来非常困难,运用连续松驰谱法将Mittag-Leffler函数用积分的形式来表示,从而解决了这个问题.通过数值算例说明了该结果的有效性.     

黏弹性松弛函数的积分表示（英文）

讨论了松弛函数的积分表示.用Laplace变换方法得到了Maxwell模型的应力应变关系方程,方程中的松弛函数可以用Mittag-Leffler函数表示.由于大的负变元的存在,计算起来非常困难,运用连续松驰谱法将Mittag-Leffler函数用积分的形式来表示,从而解决了这个问题.通过数值算例说明了该结果的有效性.     

三层媒质二维谱域Green函数的路径追踪算法

为了搜索三层媒质二维谱域Green函数中的全部模式极点,提出了一种路径追踪算法.该算法首先利用数学变换,将三层媒质二维谱域Green函数的特征方程替换为复平面上的4个超越方程;然后,应用路径追踪法和Newton-Raphson迭代法求解各超越方程的根,得到谱域Green函数的全部模式极点,即表面波极点、漏波极点和反常极点.该算法克服了围线积分法中积分路径穿过分支割线进入不同Riemann面时容易造成失根的缺点,可应用于离散复镜像方法和最陡下降路径积分法中,以实现三层媒质二维空域Green函数近场和远场的精确、快速计算.最后,利用数值算例验证了该算法的正确性和有效性.     

平面波谱-表面积分法的快速算法

针对平面波谱-表面积分法计算效率严重依赖于天线口径的问题,提出谱域积分预设方案,提取出谱域积分插值点,建立天线口径平面波角谱密度函数数据库,避免了大量重复计算,并将平面波谱理论的近场计算由四重积分弱化为二重积分.然后,根据平面波谱-表面积分法关于带罩天线远场电性能计算的思想,提出了基于等效源区域分解的并行平面波谱-表面积分算法.开发并在集群系统中实测了基于MPI环境的并行程序,作为应用计算了弹载天线-罩系统远场方向图,结果表明该程序高效准确.     

离散复镜像法求取层状介质的格林函数

对分层地层进行井间电磁成像时,需要求取分层介质中的格林函数.利用频率域中电磁场的边界条件,求出频率域分层介质的格林函数,再利用傅里叶逆变换,将格林函数变换到空间域,得到以Sommerfeld积分形式表达的解.为了避免积分中的奇异性,利用离散复镜像法将积分核用复镜像的指数求和式表示,引入广义函数束方法,可以在不提取积分核中表面波项的条件下,采用数值方法提取准静态项,给出复镜像点的数目、位置和强度,使得该方法在多层介质情况下对格林函数的计算更为有效.对一个7层介质中的垂直磁偶极子的矢量势和标量势进行了计算,离散复镜像法与Sommerfeld精确积分的结果吻合较好,说明离散复镜像法是比较准确的.     

谱域法分析微带线色散特性

研究了微带线的色散特性.采用谱域法进行数值计算得到微带线表切纵向和切向电流分布特性.利用电流分布特性,基函数选用切比雪夫多项式,在不同的频段选取适当数量的基函数,对相位常数及有效相对介电常数进行数值计算.计算结果表明:基函数的选择非常关键,采用切比雪夫多项式作为基函数可以得到非常精确的结果,提高计算效率.计算结果与HFSS软件的仿真结果非常吻合.     

M函数的数值计算

M函数是由δ冲击模型寿命分布的计算中引出一类多重积分,它的计算对精确计算冲击模型寿命分布有很大帮助.文中利用蒙特卡罗法给出了M函数的近似数值计算算法,根据算法用Matlab编程计算,给出不同参数的数值表.     

A window-based accelerating algorithm for SDA integration

给出了一个求解片上多层互连线的复数传播常数的谱域分析法(SDA).由于SDA是一种全波分析法,其计算量通常很大,为了提高积分方程的求解速度,引入了一个新颖的钟形窗口函数.该函数在谱域内具有低通滤波特性,并且拥有简单闭合的计算公式,能够有效地降低计算复杂度.基于该窗口函数的特性研究,推导出了一个新的用于提高SDA积分求解速度的加速算法.实验证明,该算法具有较高的精度和速度.     

An element by element spectral element method for elastic wave modeling

The spectral element method which combines the advantages of spectral method with those of finite element method, provides an efficient tool in simulating elastic wave equation in complex medium. Based on weak form of elastodynamic equations, mathematical formulations for Legendre spectral element method are presented. The wave field on an element is discretized using high?order Lagrange interpolation, and integration over the element is accomplished based upon the Gauss?Lobatto?Legendre integration rule. This results in a diagonal mass matrix which leads to a greatly simplified algorithm. In addition, the element by element technique is introduced in our method to reduce the memory sizes and improve the computation efficiency. Finally, some numerical examples are presented to demonstrate the spectral accuracy and the efficiency. Because of combinations of the finite element scheme and spectral algorithms, this method can be used for complex models, including free surface boundaries and strong heterogeneity.     

短管炮模块装药两相流内弹道模拟

大口径火炮模块发射装药通常由多个模块组成,流场结构复杂。该文对燃烧场进行分区模化,并在此基础上给出两相流内弹道理论模型和计算步骤,编制程序对短管试验炮膛内发射过程进行了数值模拟。结果表明,膛内的压力、压差及弹丸速度等参量的理论计算结果与实验吻合得较好,模拟结果有较高的置信度。该模型和计算程序可用于模块装药结构火炮内弹道模拟、装药设计和安全性能分析     

渗流分析的非正交网格有限体积法

将有限体积法的体积积分用Ｇｒｅｅｎ公式转化为单元边界积分，使之能应用于非正交网格的渗流分析计算；保留有限体积法计算简单实用、计算精度高、省机时省内存的特点，又有网格划分灵活方便，对复杂计算域适应性强的优点。     

一种基于网格服务的信息系统框架

提出了基于网格服务的信息系统框架,探讨了框架的形成和发展,指出基于网格服务信息 系统的特点和实现的关键技术。此框架使传统的信息系统开发过程分化为应用服务开发和集成过程,简 化并优化了系统的实施,使系统具有更强的易维护性、稳定性、可靠性和适应性。此框架使应用服务提 供商和集成商分化出来,分工更细,专业化和社会化程度更高,符合社会发展的总趋势。框架实现了多 个网格的集成,提出了多种相关技术在网格计算中的研究需求和方向。     

Stability and Accuracy Analysis for Taylor Series Numerical Method

The Taylor series numerical method (TSNM) is a time integration method for solving problems in structural dynamics. In this paper, a detailed analysis of the stability behavior and accuracy characteristics of this method is given. It is proven by a spectral decomposition method that TSNM is conditionally stable and belongs to the category of explicit time integration methods. By a similar analysis, the characteristic indicators of time integration methods, the percentage period elongation and the amplitude decay of TSNM, are derived in a closed form. The analysis plays an important role in implementing a procedure for automatic searching and finding convergence radii of TSNM. Finally, a linear single degree of freedom undamped system is analyzed to test the properties of the method.     

InP N^＋／P太阳电池光谱响应的数值模拟和分析

用数值模拟方法对ＩｎＰ Ｎ＾＋／Ｐ太阳电池在ＡＭ１．０下的光谱响应和光电流进行了计算和分析，探讨引起ＩｎＰ太阳电池光谱特性曲线变化的原因，从而获得Ｎ＾＋区和Ｐ区的最佳掺杂浓度，厚度，表面复合速度等参数，为设计制造ＩｎＰ太阳电池提供参考数据。     

新型割缝筛管有限元网格划分

针对现有割缝筛管存在的问题，研究开发出了一种新型割缝防砂筛管．针对其结构特点，利用非自由边界分割法，对其进行了网格划分，并依据等参单元拟合割缝边界的情况来控制网格密度，经网格质量评估，证明了此种网格划分方法的有效性及精确度，这种方法不仅适用于各种类型割缝筛管，也适用于开有矩形、方形等孔口的厚壁筒形构件以及旋转构件裂缝的有限元模拟，对于实际应用有很大的利用价值．     

磁轫致辐射辅助函数的特性和数值计算

磁轫致辐射函数在小宗量和大宗量区域中的展开特性。提出了一种渐近和多项式展开与数值积分相结合计算磁轫致辐射辅助函数的方法。以该算法编制的磁轫致辐射辅助函数子程序可以有效地应用于天体磁轫致辐射谱分布特性的数值计算。     

Integrating Automatic Patch Snatching System by the Concept of Generalized Mechanisms

An automatic patch snatching system is described which has three machine cells including dividing equipment, tackling equipment and patch removing equipment. The electromechanical integration design of in-piece tackling system based on the concept of generalized mechanisms is introduced. The design of the dividing equipment, the design of micro.feeding machine equipped with micro-level pressure sensor, and the realizing of computer path for patch snatching are investigated.     

三维非定常不可压N-S方程一个高精度求解格式

对三维、非定常、不可压Navier-Stokes方程提出了一种新的数值计算方法.在时间离散上应用3阶精度混合显-隐相结合的分裂格式,空间离散在x及y轴向采用非等间距网格的紧致有限差分格式与z向应用Fourier谱展开相杂交的数值方法逼近.经平板边界层流的验证表明,该算法具有计算精度高、稳定性好、收敛速度快等特点.同时也研究了三维、非定常流体运动下游边界问题,提出了无反射出流边界条件,以减少在有限计算区域内人工出流边界反射引起的数值误差,保证直接数值模拟的精度和准确性.该算法的提出对于求解边界层、射流及混合层等流动中的转捩与演化问题具有重要的理论意义.