Analysis of numerical integration error for Bessel integral identity in fast multipole method for 2D Helmholtz equation

In 2D fast multipole method for scattering problems, square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel, and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule, and the relationship between trapezoidal and square quadrature rule, sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.     

超对策模型中结局偏好模糊认知信息的融合

提出了一种新的超对策结局偏好模糊认知信息的融合方法.用一种模糊集结算子来表示专家群体对其他局中人结局偏好的模糊认知,通过求解线性规划模型得到结局偏好群体模糊认知的水平截集;以与模糊数的截集相关的解模糊函数和Newton-Cotes积分公式,得到专家群体对其他局中人结局偏好清晰认知的关系矩阵;运用大多数一致性优胜者模糊集来确定其他局中人的结局偏好向量.给出的数值例子说明了文中提出方法的有效性.     

含有裂纹水泥混凝土路面的理论分析

以断裂力学理论为基础,将水泥混凝土路面视为Winkler地基上的弹性板,利用傅立叶积分变换,并引入位错密度函数,建立了奇异积分方程,推导出水泥混凝土路面表面存在裂纹时裂纹尖端应力强度因子的解析表达式,然后应用Lobatto-Chebyshev数值方法求解奇异积分方程,得到了应力强度因子的数值解.以实际路面为例,计算了裂...     

基于SD模型的省域城市群物流联盟体系优化研究

根据省域城市群经济的发展趋势,将系统动力学基本原理与省域物流联盟体系有机结合,探讨了省域物流系统动力学的建模方法;依据省域物流系统动力学模型的特点,推导出了数学解析式.通过对京津冀省域物流联盟系统的分析,论证了该模型的实用性与实践指导意义,为省域物流系统的整合与调配提供了一个新的研究方向.     

含裂纹水泥混凝土路面的解析算法

水泥混凝土路面在使用过程中出现裂纹是一种常见的损坏形式。以断裂力学理论为基础,对含裂纹水泥混凝土路面进行了理论分析。将水泥混凝土路面视为Winkler地基上的弹性板,利用傅立叶积分变换,并通过引入位错密度函数建立奇异积分方程,推导出水泥混凝土路面含有垂直裂纹时裂纹尖端的应力强度因子的解析表达式。然后,再应用Lobatto—Chebyshev法求解奇异积分方程,得到应力强度因子的数值解。为了说明问题,以实际路面为例,给出水泥混凝土路面内部存在裂纹时裂纹尖端应力强度因子的计算结果,并讨论了影响应力强度因子大小的因素。     

路面坡度对水膜厚度的影响分析

利用满宁公式,推导出水膜厚度的计算公式。通过水膜厚度公式,分别求出水膜厚度关于横坡和纵坡的偏导函数。根据所得偏导函数进行分析,得出结论:关于横坡坡度的偏导函数小于0恒成立,即水膜厚度随着横坡坡度的增大而减小;关于纵坡坡度的偏导函数大于0恒成立,即水膜厚度随着纵坡坡度的增大而增大;横坡与纵坡坡度的偏导函数绝对值之间的比值大于2槡2,即横坡坡度对水膜厚度的影响要大于纵坡坡度。     

含分数阶导数微分方程的数值求解

基于精细积分方法,结合差分格式,提出了含有分数阶导数微分方程的数值求解方法,并对含有分数导数的一阶和二阶微分方程进行求解.所论方法首先引入差分格式,将含有分数阶导数的一阶和二阶微分方程变为一阶的常微分方程,然后再用精细积分方法逐步积分进行求解.文中不仅给出了详细的理论推导,而且还给出了相关的数值算例.数值算例对不同的分数阶微分方程进行了讨论,探讨了不同的分数阶和时间步长对计算结果的影响,并将计算结果与文献中相关算法的计算结果进行了比较.数值结果表明了所提的求解策略在求解分数阶导数模型时的可行性,以及较高的计算精度和较好的稳定性.     

Calculation Method of Wave Forces on Large Round-Ended Caisson Foundation

This paper presents an analytical formula for estimating the longitudinal wave forces on a large roundended caisson foundation. The establishment of the formula is based on the superposition of the theoretical formula of wave forces on a large circular cylinder and the empirical formula of wave forces on a large rectangular cylinder. With the formula transformed into an inertial force form, a specific inertia coefficient with an exact expression is extracted from the formula. The numerical calculations of the wave forces on round-ended cylinders are carried out by the boundary-element method. The undetermined coefficients in the expression of the inertial coefficient are determined by the numerical results. It is obvious that the numerical values can be well expressed by the computation values from the established formula. By a model experiment carried out in laboratory wave flume, the correctness of the analytical formula is further verified by the measured wave forces on a test model of a round-ended bridge caisson foundation. The comparison shows that the experimental forces can be approximately estimated by this simple calculation method.     

一个新数值积分公式的加速算法

本文讨论了具有二阶代数精度的插值型数值积分公式J。的余项表达式,给出了J。的 Bulir鸽h一Stoer外推算法和对偶算法,并且分别将它们在奇异积分和振荡函数中进行应用。 结果表明:J。的这两种加速算法具有使用方便,计算简单,稳定性好,收敛快的特点。      

Initial pre-stress finding and structural behaviors analysis of cable net based on linear adjustment theory

The tensile cable-strut structure is a self-equilibrate pre-stressed system. The initial pre-stress calculation is the fundamental structural analysis. A new numerical procedure was developed. The force density method is the cornerstone of analytical formula, and then introduced into linear adjustment theory;the least square least norm solution, the optimized initial pre-stress, is yielded. The initial pre-stress and structural performances of a particular single-layer saddle-shaped cable-net structure were analyzed with the developed method, which is proved to be efficient and correct. The modal analyses were performed with respect to various pre-stress levels. Finally, the structural performances were investigated comprehensively.     

潜水器水动力系数计算方法研究

为了保证潜水器具有良好的操纵性能,构建操纵运动方程的水动力系数必须具有相当的准确性.文中以某长航程潜水器为研究对象求解其水动力系数,分别采用了近似公式估算、基于势流理论的面元法和提出的应用FLUENT模拟循环水槽试验3种方法.通过计算结果与模型循环水槽操纵性能试验结果的对比分析,总结了上述三种计算方法的适用范围和效果评价.     

双曲面本构中单轴棘轮描述的塑性模量研究

开展了两种材料在不同应力水平和不同应力幅值下的单轴棘轮效应实验,进行了实验结果的塑性流动特性分析;揭示了单轴应力循环下的棘轮效应特征和塑性模量在正、反向加载时所遵循的不同演化规律,分析并指出了现有双曲面本构模型无法正常描述棘轮效应的实质。在此基础上,引入了新的塑性模量演化公式,并开展了数值模拟,通过与实验结果比较发现,两者吻合较好,从而验证了塑性模量演化公式的合理性。     

列车动力学的非线性数值分析方法

列车动力学方程是个复杂非线性方程,必须采用数值积分方法求解。本 文从工程实用的观点出发,提供了二类简单有效的数值积分方法—Ne物ark显式积分法及其显一隐式预刚一校正积分法。这类积分方法,积 分过程简捷,.程序编制容易;由于积分过程中无须求解高阶线性代数方 程组,这就显著地提高了积分速度,并有效地节约了计算内存,从而在 普通微机上实现了对大型列车动力学问题的模拟分析。 文中建议了一种判定工程实际问题之非线性数值积分稳定性的实用办 法,并成协地应用子重载列车动力学分析石在此基拙上进一步实现了变 步长积分计算,大大提高了计算速度。通过算例校验及与试验结果的对 照,证实了方法的非线性数值分析有效性。      

压缩波惯性作用对其波形演变的影响

为了准确地评估隧道出口微压波,基于一维平面波方程和压缩波的特性,推导了压缩波在隧道内传播时压缩波压力梯度的理论公式.采用该理论公式,分析了传播距离,初始压缩波幅值以及初始压缩波波前长度等因素对隧道不同位置处压缩波压力波梯度的影响.结果表明:压缩波压力梯度随着传播距离的增加而増大,其变化幅度越来越大,当传播距离接近临界长度时,其压力梯度几乎发生突变;相同的传播距离下,压缩波压力梯度随初始压缩波幅值的增加而增加,随初始压缩波波前长度的增加而减小;当传播距离小于临界长度时,理论公式和数值计算结果的相对误差小于6%,数值计算结果与文献报道的实测结果吻合良好.      

盾构施工引起地表沉降的数值计算与实例分析

通过对盾构施工所引起的地表沉降进行了经验公式的计算,并利用有限元软件MIDAS(GTS)进行了施工阶段的数值模拟,将计算值、模拟值与现场实测值进行了比较,并对于经验公式与数值模拟的结果合理性进行了验证.提出盾构隧道施工的横向影响范围,对于周围建筑环境保护和施工安全有指导意义.     

蜂窝梁弯扭屈曲分析及整体稳定性的计算方法

蜂窝梁是在工字钢或H型钢腹板上按一定的线形进行切割后变换位置重新焊接组合形成的新型钢梁,具有节省材料、便于铺设管道、平面内刚度增大、承载能力高等优点.由于蜂窝梁腹板开孔,与相同截面的实腹梁相比抗侧刚度被削弱,整体稳定性降低.文中以实腹梁临界弯矩计算公式为基础,考虑蜂窝梁的抗侧刚度、翘曲刚度和扭转刚度,给出蜂窝梁弯扭屈曲临界弯矩计算公式,采用ANSYS对纯弯状态下的蜂窝梁进行了弯扭屈曲分析,以蜂窝梁的孔高比和距高比为变量,给出了不同情况下蜂窝梁弯扭屈曲临界弯矩值,并与当量实腹梁临界弯矩公式计算结果进行对比,得出蜂窝梁临界弯矩与当量实腹梁临界弯矩之差随孔高比和距高比之间的变化关系,对蜂窝梁整体稳定性计算公式进行修正,最后提出了蜂窝梁整体稳定性的实用计算方法.     

Reasonable compensation coefficient of maximum gradient in long railway tunnels

This paper deals with the theory and calculation methods for compensation of the gradient in railway tunnels through theoretical analysis,numerical calculation,and statistic regression methods.On the basis of the principle that the resultant force is zero,the formula of the maximum calculated gradient was derived for the freight and passenger line and high-speed passenger special line.The formula of aerodynamic drag in tunnel is provided using the domestic and foreign relevant experimental investigations,an...     

线分布激励下功能梯度材料矩形板声辐射研究

根据Rayleigh公式,建立了预测功能梯度矩形板在无限障板中的声辐射模型并进行了理论研究.求解时,根据功能梯度复合材料矩形板的基本动力方程,首先利用模态叠加法求出线分布激励载荷作用下功能梯度板横向振动位移,然后对振动位移进行傅里叶变换,根据Rayleigh公式求解,最终得到远场声压近似解析解.通过数值仿真计算了不同阻尼、板厚、板宽和梯度指数下功能梯度板的辐射声压.对近场声辐射亦进行了数值积分分析.     

基于理实一体化的《数控机床编程与操作》课程改革探索

《数控机床编程与操作》是数控技术、机电一体化专业的核心课程，实践性较强。本文针对当前教学中存在的理论教学与实践操作之间的脱节现象，提出实施“理实一体化”教学模式改革的新思路。     

Component-based model integration approach for computer numerical control system development

A component-based model integration framework for computer numerical control system design and development is presented.The model integrates modeling,simulation,verification and implementation in a uniform environment.The computer numerical control(CNC) modeling language with well defined syntax and unambiguous semantics is developed.Using the proposed CNC model integration method,a three axis milling system model is developed in the case study.The approach is an attempt to create an infrastructure to suppo...