确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

Effect of crack orientation statistics on effective stiffness of mircocracked solid

This paper addresses the problem of calculating effective elastic properties of a solid containing multiple cracks with prescribed orientation statistics. To do so, the representative unit cell approach has been used. The microgeometry of a cracked solid is modeled by a periodic structure with a unit cell containing multiple cracks; a sufficient number is taken to account for the microstructure statistics. The developed method combines the superposition principle, the technique of complex potentials and certain new results in the theory of special functions. A proper choice of potentials provides reducing the boundary-value problem to an ordinary, well-posed set of linear algebraic equations. The exact finite form expression of the effective stiffness tensor has been obtained by analytical averaging the strain and stress fields. The convergence study has been performed; the statistically meaningful results obtained show dependence of the effective elastic stiffness on angular scattering of cracks. Comparison has been made with the selected simple micromechanical models, namely, non-interaction approximation, differential scheme and modified differential scheme. It is found that, among these models, the differential scheme provides the best fit of the numerical data.     

混凝土三维裂纹应力强度因子的简化算法

根据混凝土特性,本文提出了一种计算脆性材料三维裂缝应力强度因子的有限元方法,用以分析混凝土结构物中裂缝的稳定性及危害程度.其特点是不需要在裂缝尖端附近区域划分过细的网络,也不用引入特殊尖端单元,便于在工程中使用.本方法曾用于某拱坝坝肩裂缝稳定判别及应力重分布计算,效果良好.     

A conservative local interface sharpening scheme for the constrained interpolation profile method

A conservative local interface sharpening scheme has been developed for the constrained interpolation profile method with the conservative semi‐Lagrangian scheme, because the conservative semi‐Lagrangian scheme does not feature a mechanism to control the interface thickness, thus causing an increase of numerical error with the advance of the time step. The proposed sharpening scheme is based on the conservative level set method proposed by Olsson and Kreiss. However, because their method can cause excessive deformation of the free‐surface in certain circumstances, we propose an improvement of the method by developing a local sharpening technique. Several advection tests are presented to assess the correctness of the advection and the improved interface sharpening scheme. This is followed by the validations of dam‐breaking flow and the rising bubble flows. The mass of the fluid is exactly conserved and the computed terminal velocity of the rising bubble agrees well with the experiments compared with other numerical methods such as the volume of fluid method (VOF), the front tracking method, and the level set method. Copyright © 2011 John Wiley & Sons, Ltd.     

A hybrid FVM–LBM method for single and multi‐fluid compressible flow problems

The lattice Boltzmann method (LBM) has established itself as an alternative approach to solve the fluid flow equations. In this work we combine LBM with the conventional finite volume method (FVM), and propose a non‐iterative hybrid method for the simulation of compressible flows. LBM is used to calculate the inter‐cell face fluxes and FVM is used to calculate the node parameters. The hybrid method is benchmarked for several one‐dimensional and two‐dimensional test cases. The results obtained by the hybrid method show a steeper and more accurate shock profile as compared with the results obtained by the widely used Godunov scheme or by a representative flux vector splitting scheme. Additional features of the proposed scheme are that it can be implemented on a non‐uniform grid, study of multi‐fluid problems is possible, and it is easily extendable to multi‐dimensions. These features have been demonstrated in this work. The proposed method is therefore robust and can possibly be applied to a variety of compressible flow situations. Copyright © 2009 John Wiley & Sons, Ltd.     

Use of characteristics method with cubic interpolation for unsteady-flow computation

The specified-time-interval (STI) scheme has been used commonly in applying the method of characteristics (MOC) to unsteady open-channel flow problems. However, with the use of STI scheme, the numerical error for the simulation results can always be induced due to the interpolation used to approximate the characteristics trajectory. Hence, in order to remedy the numerical errors caused by the interpolation, one needs to seek some kind of interpolation technique with higher-order accuracy. Instead of the linear interpolation technique, which has been used very commonly and can induce serious numerical diffusion, the Holly--Preissmann two-point, method, which is a cubic interpolation technique with fourth-order of accuracy, is proposed here to integrate with the method of characteristics for the computation of one-dimensional unsteady flow in open channel. The concept of reachback and reachout in space and time directions for the characteristics is also introduced to assure the model stability. The computed results from this new model are compared with those computed by using the Preissmann four-point scheme and the multimode method of characteristics with linear interpolation.     

High‐accuracy upwind method using improved characteristics speeds for incompressible flows

In this study, the Nervier–Stokes equations for incompressible flows, modified by the artificial compressibility method, are investigated numerically. To calculate the convective fluxes, a new high‐accuracy characteristics‐based (HACB) scheme is presented in this paper. Comparing the HACB scheme with the original characteristic‐based method, it is found that the new proposed scheme is more accurate and has faster convergence rate than the older one. The second order averaging scheme is used for estimating the viscose fluxes, and spatially discretized equations are integrated in time by an explicit fourth‐order Runge–Kutta scheme. The lid driven cavity flow and flow in channel with a backward facing step have been used as benchmark problems. It is shown that the obtained results using HACB scheme are in good agreement with the standard solutions. Copyright © 2015 John Wiley & Sons, Ltd.     

THE NMM-BASED EFG METHOD AND SIMULATION OF CRACK PROPAGATION

数值流形方法（numerucal manifold method,NMM）通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题. 通过用移动最小二乘插值（moving least squares interpolation,MLS）中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法（element free Galerkin,EFG）. 该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题. 建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹（kinked crack）的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性. 大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.     

Cures for expansion shock and shock instability of Roe scheme based on momentum interpolation mechanism

The common defects of the Roe scheme are the non-physical expansion shock and shock instability. By removing the momentum interpolation mechanism (MIM), an improved method with several advantages has been presented to suppress the shock instability. However, it cannot prevent the expansion shock and is incompatible with the traditional curing method for expansion shock. To solve the problem, the traditional curing mechanism is analyzed. Effectiveness of the traditional curing method is discussed, and several defects are identified, one of which leads to incompatibility between curing shock instability and expansion shock. Consequently, an improved Roe scheme is proposed, which is with low computational costs, concise, easy to implement, and robust. More importantly, the proposed scheme can simultaneously solve the problem of shock instability and expansion shock without additional costs.     

竖向裂缝对功能梯度板的热响应影响

采用条单元法分析了竖向裂缝对功能梯度板的热响应影响。功能梯度板的材料属性沿板厚方向连续线性变化,温度变化热源作用于带竖向裂缝的功能梯度板上。通过研究得出了裂缝距离热源不同位置、不同长度裂缝及不同频率温度变化热源等状态时,功能梯度板上表面的热响应情况。采用条单元法研究带竖向裂缝的功能梯度板的热响应问题,为带裂缝的功能梯度板的热研究提供了一定的指导。     

Numerical schemes with high order of accuracy for the computation of shock waves

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｍｐｕｔａｔｉｏｎｏｆｆｌｏｗｆｉｅｌｄｗｉｔｈｓｈｏｃｋｗａｖｅｓｈａｓｂｅｅｎｔｈｅｓｕｂｊｅｃｔｏｆｒｅｓｅａｒｃｈｆｏｒｍａｎｙｙｅａｒｓ.Ｔｈｅｒｅａｒｅｂａｓｉｃａｌｌｙｔｗｏｃａｔｅｇｏｒｉｅｓｏｆｍｅｔｈｏｄｓ,ｎａｍｅｌｙ,ｓｈｏｃｋｆｉｔｔｉｎｇｍｅｔｈｏｄａｎｄｓｈｏｃｋｃａｐｔｕｒｉｎｇｍｅｔｈｏｄ.Ｔｈｅｆｏｒｍｅｒｄｉｖｉｄｅｓｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｉｎｔｏｓｕｂ＿ｄｏｍａｉｎｓｂｙｔｈｅｓｈｏｃｋｗａｖｅｓ.Ｉｎｅ…     

A finite difference method for elastic wave scattering by a planar crack with contacting faces

This paper presents a finite difference time-domain technique for 2D problems of elastic wave scattering by cracks with interacting faces. The proposed technique introduces cracks into the finite difference model using a set of split computational nodes. The split-node pair is bound together when the crack is closed while the nodes move freely when open, thereby a unilateral contact condition is considered. The development of the open/close status is determined by solving the equation of motion so as to yield a non-negative crack opening displacement. To check validity of the proposed scheme, 1D and 2D scattering problems for which exact solutions are known are solved numerically. The 1D problem demonstrates accuracy and stability of the scheme in the presence of the crack-face interaction. The 2D problem, in which the crack-face interaction is not considered, shows that the proposed scheme can properly reproduce the stress singularity at the tip of the crack. Finally, scattered fields from cracks with interacting faces are investigated assuming a stick and a frictionless contact conditions. In particular, the directivity and higher-harmonics are investigated in conjunction with the pre-stress since those are the basic information required for a successful ultrasonic testing of closed cracks.     

An efficient and accurate algebraic interface capturing method for unstructured grids in 2 and 3 dimensions: The THINC method with quadratic surface representation

We present in this paper an efficient and accurate volume of fluid (VOF) type scheme to compute moving interfaces on unstructured grids with arbitrary quadrilateral mesh elements in 2D and hexahedral elements in 3D. Being an extension of the multi‐dimensional tangent of hyperbola interface capturing (THINC) reconstruction proposed by the authors in Cartesian grid, an algebraic VOF scheme is devised for arbitrary quadrilateral and hexahedral elements. The interface is cell‐wisely approximated by a quadratic surface, which substantially improves the numerical accuracy. The same as the other THINC type schemes, the present method does not require the explicit geometric representation of the interface when computing numerical fluxes and thus is very computationally efficient and straightforward in implementation. The proposed scheme has been verified by benchmark tests, which reveal that this scheme is able to produce high‐quality numerical solutions of moving interfaces in unstructured grids and thus a practical method for interfacial multi‐phase flow simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

双周期平行四边形排列裂纹反平面问题的有限元分析

提出了双周期平行四边形排列裂纹反平面问题的有限元方法,通过对单位胞元引入周期边界条件,在裂纹尖端采用奇异单元,解决了有限元分析这类问题的效率和精度问题.利用Ansys软件计算,在各种有解析解对照的情形下,应力强度因子的相对误差都在0.2%以内.与现有通常限于对称阵列的双周期裂纹的研究相比,本文发展的方法适用于一般的非对称平行四边形裂纹阵列.算例揭示了行向裂纹间的相互干涉放大应力强度因子,而叠向裂纹间的作用相互屏蔽.对于平行四边形阵列的情形,这两种相反的干涉效应使应力强度因子与裂纹错动参数间呈现非单调依赖关系.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

Partitioned subiterative coupling schemes for aeroelasticity using combined interface boundary condition method

The combined interface boundary condition (CIBC) method has been recently proposed for fluid–structure interaction. The CIBC method employs a Gauss–Seidel-like procedure to transform traditional interface conditions into velocity and traction corrections whose effect is controlled by a dimensional parameter. However, the original CIBC method has to invoke the uncorrected traction when forming the traction correction. This process limits its application to fluid–rigid body interaction. To repair this drawback, a new formulation of the CIBC method has been developed by using a new coupling parameter. The reconstruction is simple and the structural traction is removed completely. Two partitioned subiterative coupling versions of the CIBC method are developed. The first scheme is an implicit strategy while the second one is a semi-implicit strategy. Iterative loops are actualised by the fixed-point algorithm with Aitken accelerator. The obtained results agree with the well-documented data, and some famous flow phenomena have been successfully detected.     

薄板结构的边界元优化设计

完成了一种用边界元法对弹性薄板结构进行优化设计的计算机分析及彩色绘图系统.本文论述了平板弯曲问题的边界元法及其优化设计.系统采用本文作者提出的三方程协调方案,并精确计算分布荷载的域内积分.本文完成了多个平板弯曲问题的边界元分析算例及板结构边界元优化设计的工程算例.结果表明,本系统方法先进、结果精确可靠,具有十分明显的工程实用价值.     

An improvement of the Prandtl-Reuss theory and the hybrid / mixed vari-stiffness method of finite element analysis

The common Prandtl-Reuss theory has been improved in this paper. A quasi-flow law of the isotropic hardening Mises materials has been proposed as well, on the basis of which, an efficient iterative algorithm of finite element analysis, hybrid / mixed vari-stiffness method, has been obtained. The numerical examples calculated by the plane stress / strain element model are given. Compared with the common initial stress method, the hybrid / mixed vari-stiffness method shows its advantages in the convergent speed, calculating accuracy and treatment scheme of the incompressibility of materials.     

Incompressible Roe‐like scheme for turbulent flows

In the current study, numerical investigation of incompressible turbulent flow is presented. By the artificial compressibility method, momentum and continuity equations are coupled. Considering Reynolds averaged Navier–Stokes equations, the Spalart–Allmaras turbulence model, which has accurate results in two‐dimensional problems, is used to calculate Reynolds stresses. For convective fluxes a Roe‐like scheme is proposed for the steady Reynolds averaged Navier–Stokes equations. Also, Jameson averaging method was implemented. In comparison, the proposed characteristics‐based upwind incompressible turbulent Roe‐like scheme, demonstrated very accurate results, high stability, and fast convergence. The fifth‐order Runge–Kutta scheme is used for time discretization. The local time stepping and implicit residual smoothing were applied as the convergence acceleration techniques. Suitable boundary conditions have been implemented considering flow behavior. The problem has been studied at high Reynolds numbers for cross flow around the horizontal circular cylinder and NACA0012 hydrofoil. Results were compared with those of others and a good agreement has been observed. Copyright © 2012 John Wiley & Sons, Ltd.