微分求积法分析平面接头应力奇异性北大核心CSCD

对于双材料平面接头问题提出了一个分析应力奇性指数的新方法:微分求积法(DQM).首先,将平面接头连接点处位移场的径向渐近展开格式代入平面弹性力学控制方程,获得了关于应力奇性指数的常微分方程组(ODEs)特征值问题.然后,基于DQM理论,将ODEs的特征值问题转化为标准型广义代数方程组特征值问题,求解之可一次性地计算出双材料平面接头连接点处应力奇性指数,同时,一并求出了接头连接点处相应的位移和应力特征函数.数值计算结果说明该文DQM计算平面接头连接点处应力奇性指数的结果是正确的.     

基于MPS-FEM耦合方法对比研究刚性与弹性挡板对液舱晃荡的抑制作用

由流体冲击载荷引起的流固耦合问题广泛存在于船舶与海洋工程领域.例如:在特定激励频率下载液货舱内流体的非线性运动引起对舱壁的砰击作用,进而可能影响液舱围护系统的安全性.由于此类流固耦合问题通常涉及多学科知识,且流体自由面的变化具有强非线性特征,对研究人员带来较大挑战.考虑到Lagrange类方法在处理结构和流体自由面大变形问题上的优势,基于MPS-FEM耦合方法开发了流固耦合求解器.其中,采用MPS方法来数值模拟流体场瞬态变化,FEM方法来分析结构场的变形问题.此外,该求解器采用了弱耦合的方式来实现流体场和结构场之间的数据传递.为了验证该方法在处理流固耦合问题上的可靠性,首先数值研究了溃坝泄洪流与弹性挡板之间的流固耦合标准算例,数值结果与实验标准结果能够较好地吻合.此后,采用该求解器数值研究了带刚性挡板和弹性挡板的液舱晃荡问题,对比分析了多种激励频率下两种挡板对液舱内流体运动及舱壁上冲击压力的抑制效果.     

伴有磁场和纳米固体颗粒时的Jeffery-Hamel流动解析研究——Adomian分解法

用一种强有力的解析方法,称为Adomian分解法(ADM),来研究磁场和纳米颗粒对Jeffery-Hamel流动的影响.将该问题模型的控制方程,即将传统的流体力学Navier-Stokes方程和Maxwell电磁方程,简化为非线性的常微分方程.该方法得到的结果与Runge-Kutta方法得到的数值结果相一致,结果用表格列出.不同α,Ha和Re数下的图形表明,本方法可以得到高精度的结果.首先对不同的Hartmann数和管壁倾角,研究喇叭形管道中的流场;最后在没有磁场作用时,研究纳米固体颗粒体积率的影响.     

二维非稳态对流扩散边界控制问题的简化算法

针对二维非稳态对流扩散边界控制问题计算量大的问题,提出了基于降阶模型的最优实时控制方法.利用POD(the Proper Orthogonal Decomposition)和奇异值分解以及Galerkin投影方法得到了具有高精度离散形式的状态空间降阶模型.在所得的降阶状态空间模型中,利用离散时间线性二次调节器方法设计出了最优控制器.对流-扩散过程的控制模拟结果说明了所提方法的有效性和准确性.     

灰色时序组合模型及其在地下水埋深预测中的应用

地下水埋深的变化过程是一个复杂的非线性过程,这种具有复杂的非线性组合特征的序列,使用某一种模型进行预测,结果往往不理想.在分析了灰色GM(1,1)模型、灰色GM(1,1)周期性修正模型和时序AR(n)模型的优点和缺点基础上,提出了一种新的灰色时序组合预报模型.该方法利用了GM预测所需原始数据少、方法简单的优点,用周期修正方法反映其地下水位埋深周期性波动的特征,用AR(n)模型预报其地下水位埋深的随机变化.实例研究表明,这种方法方便简洁实用且预测结果接近于实际观测值,为其它地区的地下水位埋深和相关时间序列的分析研究提供参考与借鉴作用.     

利用场发射电子显微镜研究具有原子分辨的单壁碳纳米管的场发射

利用场发射电子显微镜(FEM)研究了组装在金属钨(W)针尖上的单壁碳纳米管(SWCNTs)的电子发射特性. 透射电子显微镜(TEM)结果表明在W针尖上存在一单壁碳纳米管束. 在适当的实验条件下得到了一组具有原子分辨率的FEM像, 该像可能是从单壁碳纳米管束中突出的一根(16,0)锯齿型单壁碳纳米管开口端的发射形成的. 理论计算了FEM图像放大倍数和分辨率. 当压缩因子β的值设定为β= 1.76时, FEM图像放大倍数的理论计算值与实际测量值相一致. 利用考虑镜像力的Gomer分辨率公式计算了FEM的分辨本领. 在典型电场强度ε = 5.0×10⁷ V/cm)情况下, FEM的分辨本领可以达到0.277 nm, 这接近(16,0)型单壁碳纳米管开口端碳原子沿锯齿边的间距(0.246 nm). 这些进一步支持了实验得到的结果.     

关于方程φ(abc)=2(φ(a)+φ(b)+φ(c))

设n为任意正整数,φ(n)是Euler函数.主要研究了方程φ(abc)=2(φ(a)+φ(b)+φ(c))的可解性问题,利用数论中的理论和方法,获得了该方程的所有正整数解.     

基于群智能算法的预防性维修周期优化

分析将蚁群优化算法应用于预防性维修周期工程寻优问题时遇到的算法参数选择困难等问题,提出将粒子群优化算法和空间划分方法引入该过程以改进原蚁群算法的寻优规则和历程.建立混合粒子群和蚁群算法的群智能优化策略:PS_ACO(Particle Swarm and Ant Colony Optimization),并将其应用于混联系统预防性维修周期优化过程中,以解决由于蚁群算法中参数选择不当和随机产生维修周期解值带来的求解精度差、寻优效率低等问题.算法的寻优结果对比分析表明:该PS_ACO算法应用于预防性维修周期优化问题,在寻优效率及寻优精度上有部分改进,且可相对削弱算法参数选择对优化结果的影响.     

一类偏微分方程(组)非古典对称存在性的判定方法

基于微分特征集理论和算法,提出在一定条件下判定偏微分方程(组)非古典对称存在性的机械化方法.该方法对Clarkson P A提出的关于偏微分方程(组)的非古典对称的公开问题给出了部分回答,为完全解决该问题提供了一个思路.通过若干个发展方程的非古典对称的确定说明了该方法的有效性.     

基于三次Bézier基函数插值的GM(1,1)模型背景值优化研究

本文研究了传统灰色GM(1,1)模型存在模型精度不高的问题.利用带形状参数的三次Bézier基函数,给出插值函数的表达式,并结合复化梯形公式,给定误差限的方法,获得了比传统灰色GM(1,1)模型更高精度的结果.推广了传统灰色GM(1,1)预测模型的结果.     

Band structure computation of two-dimensional and three-dimensional phononic crystals using a finite element-least square point interpolation method

In the present study, a finite element-least square point interpolation method (FE-LSPIM) is proposed for calculating the band structures of in-plane elastic waves in two-dimensional (2D) and three-dimensional (3D) phononic crystals (PCs). This method utilizes new shape functions by combining mesh-free shape functions and finite element shape functions to exploit the specific advantages of the mesh-free method and finite element method (FEM). As a result, FE-LSPIM inherits the completeness properties of the mesh-free method and the compatibility properties of FEM, and thus the solutions obtained tend to be more accurate. Indeed, according to our previous research, the present method obtains excellent accuracy, especially in the high-frequency domain. The proposed FE-LSPIM was combined with Bloch's theory and applied to compute the band gaps (BGs) for 2D PCs and 3D PCs in the present study, where several PCs were investigated to verify the high accuracy when computing the BGs. Numerical analysis showed that the proposed method can predict the BGs more precisely compared with the FEM and modified FEM.     

Free vibration and buckling analysis of plates using B-spline wavelet on the interval Mindlin element

A finite element method (FEM) of B-spline wavelet on the interval (BSWI) is used in this paper to solve the free vibration and buckling problems of plates based on Reissner–Mindlin theory. By aid of the high accuracy of B-spline functions approximation for structural analysis, the proposed method could obtain a fast convergence and a satisfying numerical accuracy with fewer degrees of freedoms (DOF). The numerical examples demonstrate that the present BSWI method achieves the high accuracy compared to the exact solution and others existing approaches in the literatures. The BSWI finite element has potential to be used as a numerical method in analysis and design.     

HIGH ACCURACY NUMERICAL METHOD OF THIN－FILM PROBLEMS IN MICROMAGENETICS

In this paper,a new high accuracy numerical method for the thin-film problems of micron and submicron size ferromagnetic elements is proposed,For the computaion of stray field,we use the finite element method(FEM) by introducing a semi-discrete artificial boundary condition [1,2],In our numerical experiments about the domain patterns and their movement,we can see that the results are accordant to that of experments and other numerical methods.Our method are very conveient to deal with arbitrary shape of thin films such as a polygon with high accuracy.     

An efficient mixed methodology for free vibration and buckling analysis of orthotropic rectangular plates

In this paper, a simple and efficient mixed Ritz-differential quadrature (DQ) method is presented for free vibration and buckling analysis of orthotropic rectangular plates. The mixed scheme combines the simplicity of the Ritz method and high accuracy and efficiency of the DQ method. The accuracy of the proposed method is demonstrated by comparing the calculated results with those available in the literature. It is shown that highly accurate results can be obtained using a small number of Ritz terms and DQ sample points. The proposed method is suitable for the problem considered due to its simplicity and potential for further development.     

HIGH ACCURACY NUMERICAL METHOD OF THIN-FILM PROBLEMS IN MICROMAGNETICS

In this paper, a new high accuracy numerical method for the thin-film problems of micron and submicron size ferromagnetic elements is proposed. For the computation of stray field, we use the finite element method(FEM) by introducing a semi-discrete artificial boundary condition [1, 2]. In our numerical experiments about the domain patterns and their movement, we can see that the results are accordant to that of experiments and other numerical methods. Our method are very convenient to deal with arbitrary shape of thin films such as a polygon with high accuracy.     

Static contact analysis of spiral bevel gear based on modified VFIFE (vector form intrinsic finite element) method

Since the intrinsic limitations of FEM (Finite element method) and lumped-mass method, we derive the formula of 8-node hexahedral element based on VFIFE (vector form intrinsic finite element method) method and applied it in contact analysis of gears. This paper proposed a new method to determine pure nodal deformation, which could simplify the computation compared to the traditional VFIFE method. Combining the VFIFE method and matching contact algorithm, we analyzed spiral bevel gear meshing problems. Spiral bevel models with two different mesh densities are calculated analyzed by the VFIFE method and FEM. Performance indicators of gears are extracted and compared, including contact forces, contact and bending stresses, contact stress patterns and loaded transmission errors. The results show that the VFIFE method has a stable performance and reliable accuracy under coarse or refined mesh conditions, while the FEM inaccurately calculates the contact stress of the coarse mesh model. The examples demonstrate that the proposed method could precisely analyze gear meshing problems with a coarse mesh model, which provides a new solution for gear mechanics.     

Accurate vibration analysis of skew plates by the new version of the differential quadrature method

An accurate free vibration analysis of skew plates is presented by using the new version of the differential quadrature method (DQM). Eight combinations of simply supported (S), clamped (C) and free (F) boundary conditions are considered. Detailed solution procedures are given and key points for success by using the DQM are emphasized. A way to simplifying the programming in using the DQM is proposed. Convergence study is made for the simply supported skew plate with a large skew angle. Good convergence of frequencies is observed. The DQ results agree very well with the existing first known accurate upper bound solutions, obtained by using Ritz method taking into considerations of the bending stress singularities occurred at corners having obtuse angles. Since slight discrepancy between the DQ data and the known accurate solutions is observed for plates with large skew angles, the DQ results are also compared with data obtained by using finite element method with very fine meshes to verify their accuracy.     

The construction of plane elastomechanics and Mindlin plate elements of B-spline wavelet on the interval

Based on two-dimensional tensor product B-spline wavelet on the interval (BSWI), a class of C0 type plate elements is constructed to solve plane elastomechanics and moderately thick plate problems. Instead of traditional polynomial interpolation, the scaling functions of two-dimensional tensor product BSWI are employed to form the shape functions and construct BSWI elements. Unlike the process of direct wavelets adding in the previous work, the elemental displacement field represented by the coefficients of wavelets expansions is transformed into edges and internal modes via the constructed transformation matrix in this paper. The method combines the versatility of the conventional finite element method (FEM) with the accuracy of B-spline functions approximation and various basis functions for structural analysis. Some numerical examples are studied to demonstrate the proposed method and the numerical results presented are in good agreement with the closed-form or traditional FEM solutions.     

FLIM - A Variant of FEM based on Line Integration

A variant of the finite element method (FEM) for modelling and solving partial differential equations based on triangular and tetrahedral meshes is proposed. While FEM is based on integration over finite elements, the new approach - briefly denoted as FLIM hereafter - uses integration along edges (finite lines). The stiffness matrix, which - for linear triangles and tetrahedra - is identical with the one obtained with FEM, as well as the load vector can solely be obtained by summing up the edge contributions. This new variant requires much lower storage than FEM, especially for three-dimensional problems, but yields the same approximation error and convergence rate as the finite element method. It is shown that its performance, when applied to linear problems, is in close agreement with the performance of the finite element method. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于时域精细积分算法的瞬态传热多宗量反演

基于有限元法和精细积分算法,提出了一种求解瞬态热传导多宗量反演问题的新方法．采用有限元法和精细积分算法分别对空间、时间变量进行离散,可以得到正演问题高精度的半解析数值模型,由此建立了多宗量反演的计算模式,并给出敏度分析的计算公式．对一维和二维的热物性参数、热源项、边界条件等进行了单宗量和多宗量的反演求解,初步考虑了初值和噪音等对反演结果的影响,数值算例验证了该方法的有效性．