基于边界元法的非均质薄板弯曲问题的解

本文用边界元法研究非均质无限域弹性薄板弯曲问题.在数值实施过程中,对于夹杂和基体分别形成边界积分方程.通过离散边界积分方程,得到相应的方程组,然后结合界面条件,最终获得问题的求解方程组.在界面的相关量求得之后,可以根据需要来求解基体和夹杂中的有关位置的弯矩.数值结果与已有的解做了对比.     

瞬态热传导问题的精细积分-双重互易边界元法

采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。     

流体晃动时域双协边界元法分析

利用双协边界元法在时域内对流体晃动问题进行分析,推导出边界积分方程及相应的边界条件。分析过程中考虑流体的粘性,自由面上的动力学条件为法向正应力和切向剪应力为零。固壁面上采用流体质点与固壁质点速度相等的条件。时域的离散采用差分法,并利用时间步迭代,逐步追踪流体自由面,在流体的不断变动的边界上考虑其边界条件。数值结果表明本文的双协边界元法是可行的。     

边界元法计算切口多重应力奇性指数

提出采用边界元法直接计算V形切口的多重应力奇性指数。首先在切口尖端挖出一微小扇形域,在该域边界列常规边界积分方程,后将扇形域内的位移场和应力场表示成关于切口尖端距离ρ的渐近级数展开式,回代入切口边界积分方程,离散后得到关于切口奇性指数的代数特征方程,从而求解获得V形切口的应力奇性指数。该法避免了常规边界元法和有限元法在切口尖端附近布置细密单元的缺陷,并可同时求得多阶应力奇性指数。     

扩展边界元法分析切口和裂纹结构应力场的准确性

在线弹性理论中,切口/裂纹结构尖端区域存在奇异应力场,数值方法不易求解。本文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离 的渐近级数展开式表达,其级数项的幅值系数作为基本未知量,而外部区域采用常规边界元法离散方程。两者方程联立求解可获得切口和裂纹结构完整的位移和应力场。扩展边界元法具有半解析法特征,适用于一般的切口和裂纹结构应力场分析,其解可精细描述从尖端区域到整体结构区域的应力场。作者研制了扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明扩展边界元法求解切口和裂纹结构应力场的准确性和有效性。     

直接边界元法三维自由面兴波时域模拟

应用直接边界元法在时域中求解稳定航速运动的三维自由面兴波问题.基于格林定理,在所有边界面上划分网格,对边界积分方程进行数值离散,采用线性自由面边界条件,随时间步进更新自由面势.由于物体空间位置移动辐射条件不需要单独表述,迭代过程中自由面计算域保持不变.以割划水面NACA0024为例,计算模拟了自由面兴波稳定波形;提出了求解矩阵方程组奇异性的处理方法和解决割划问题的动网格技术.本文计算结果和有限体积法及有关试验结果对比表明,该方法是可靠的.     

边界元法在流—固相互作用分析中的应用

本文采用边界元法研究水与粘弹性介质在地震波入射条件下的相互作用问题。首先,采用半空间格林函数建立液体介质区域的边界元方程,然后考虑液—固界面的连续条件,获得水对固体介质的影响矩阵,最后用此影响矩阵修正固体介质的边界元方程。通过对简单问题的计算表明,采用本文方法求解流—固相互作用问题,其计算量与单独求解固体介质的波动问题时相当。     

物性值随温度变化热弹性问题的摄动—边界元分析

本文将摄动法和边界元法相结合求解物性值随温度变化的热弹性问题,简述了基本方程和积分方程的建立,导出了有关计算公式。算例表明本文方法简便、有效。     

一种新型的边界元法——边界轮廓法

利用传统边界元积分方程的被积函数的散度等于零的特性，提出一种新型的边界元法——边界轮廓法，使求解问题的维数降低两维。对线弹性平面问题，选择二次位移形函数，求得相应的位移和应力势函数，使二维问题的求解转化为边界点的数值计算，给出了边界点的位移和面力及域内点的应力和位移的计算公式。实例计算表明，该方法具有较高的精度。     

平面热弹性问题的边界元分析

本文利用位移法由平面热弹性问题的基本方程出发,简要地叙述了边界积分方程的建立及离散化手法,导出了由边界上的位移和表面力直接计算边界应力的公式。作为数值计算例,计算了圆形区域,同心圆区域和具有偏心圆孔的圆形区域的热应力。计算结果与解析解或实验结果进行了比较,两者相当吻合。计算表明,边界元法对求解平面热弹性问题十分有效.本文也适用于有体积力的平面弹性问题.     

Radial integral boundary element method for simulating phase change problem with mushy zone

A radial integral boundary element method(BEM) is used to simulate the phase change problem with a mushy zone in this paper. Three phases, including the solid phase, the liquid phase, and the mushy zone, are considered in the phase change problem. First, according to the continuity conditions of temperature and its gradient on the liquid-mushy interface, the mushy zone and the liquid phase in the simulation can be considered as a whole part, namely, the non-solid phase, and the change of latent heat is approximated by heat source which is dependent on temperature. Then, the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone, respectively. Moreover, an iterative predictor-corrector precise integration method(PIM) is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary. According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone, the front-tracking method is used to track the move of the interface. The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally, four numerical examples are provided to assess the performance of the proposed numerical method.     

Transient heat conduction analysis using the NURBS-enhanced scaled boundary finite element method and modified precise integration method

The Non-uniform rational B-spline(NURBS)enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper.The scaled boundary finite element method is a semi-analytical technique,which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction.In this method,only the boundary is discretized in the finite element sense leading to a reduction of the spatial dimension by one with no fundamental solution required.Nevertheless,in case of the complex geometry,a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often unavoidable in the conventional finite element approach,which leads to huge computational efforts and loss of accuracy.NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape.In the proposed methodology,the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions,while the straight part of the boundary is discretized by the conventional Lagrange shape functions.Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analysis and the solution is obtained using the modified precise integration method.The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion.Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method.The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

The intersection marker method for 3D interface tracking of deformable surfaces in finite volumes

Currently, the majority of computational fluid dynamics (CFD) codes use the finite volume method to spatially discretise the computational domain, sometimes as an array of cubic control volumes. The Finite volume method works well with single‐phase flow simulations, but two‐phase flow simulations are more challenging because of the need to track the surface interface traversing and deforming within the 3D grid. Surface area and volume fraction details of each interface cell must be accurately accounted for, in order to calculate for the momentum exchange and rates of heat and mass transfer across the interface. To attain a higher accuracy in two‐phase flow CFD calculations, the intersection marker (ISM) method is developed. The ISM method is a hybrid Lagrangian–Eulerian front‐tracking algorithm that can model an arbitrary 3D surface within an array of cubic control volumes. The ISM method has a cell‐by‐cell remeshing capability that is volume conservative and is suitable for the tracking of complex interface deformation in transient two‐phase CFD simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

Boundary Integral Analysis of the Symmetric Dynamic Problem for an Infinite Bimaterial Solid with an Embedded Crack

The symmetric frequency domain problem for two ideally bonded elastic half-spaces with a perpendicular plane crack is considered. It is reduced to the boundary integral equation (BIE) with integration over the limited crack region. The contact conditions on the bimaterial interface are satisfied identically in the initial stage of obtaining the equation. After boundary element solution of the equation, the stress concentration in the vicinity of a penny-shaped crack under time-harmonic loading of constant amplitude is studied. The mode I stress intensity factors as functions of angular coordinate of a crack front point and wave number for various relations between the material parameters are computed. The crack depth relative to the bimaterial interface is determined, when the effect of the material dissimilarity on the crack can be neglected.     

Numerical simulation of solidification processes in enclosures

The present work deals with the development and application of numerical models for the simulation of solidification problems liquid/solid taking diffusion and convection into account. For the calculation of the thermal coupled flow process the finite element method is applied. In order to improve the numerical stability of the free convection problems, the streamline-upwind/Petrov–Galerkin method is used. Solidification processes are moving boundary problems. Three different models are set up which consider latent heat at the solidification front respectively in the mixed zone during the phase transition. Moreover, numerical methods are investigated in order to describe the behaviour of the flow at the boundary of the moving phase. Three examples serve illustrations; the technical example – casting of a transport and storage container – was provided by the company Siempelkamp Gießerei GmbH.     

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

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

Analyses of dynamic characteristics of a fluid-filled thin rectangular porous plate with various boundary conditions

Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are derived in the frequency domain. By considering the coupling effect between the solid phase and the fluid phase and without any hypothesis for the fluid displacement, the model presented here is rigorous and close to the real materials. Owing to the use of extended homogeneous capacity precision integration method and precise element method, the model can be applied in higher frequency range than pure numerical methods. This model also easily adapts to various boundary conditions. Numerical results are given for two different porous plates under different excitations and boundary conditions.     

Timoshenko梁谐响应求解新方法探索

在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.