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.     

Elastodynamic Green’s function for reinforced concrete beams

A semi-analytical solution procedure for three dimensional wave propagation in reinforced concrete (RC) beams has been presented in this paper. Elastodynamic Green’s function has been derived by employing the compatibility conditions and utilizing the symmetry conditions at the loaded cross section. Numerical procedure developed for the Green’s function has been validated using results available in the literature for an infinite laminated composite plate. Three-dimensional wave propagation analysis has been performed for reinforced concrete beam sections of T and L shapes which are common forms of structural elements. Steel reinforcement has been modeled in the finite element mesh. Effect of corrosion has also been included in the finite element model. Green’s function for reinforced concrete sections affected by corrosion of steel unit normalized frequency has been evaluated for illustration. Accuracy of the solution technique has been evaluated in terms of the percentage error in energy balance between the input energy of the applied unit load and the output energy carried by the propagating wave modes. The percentage error has been found to be negligible in all the cases considered here. A simple and accurate numerical method has been presented here as a tool to evaluate Green’s function for RC beams and can be used to detect corrosion.     

Isogeometric analysis of nonlinear Euler–Bernoulli beam vibrations

In this paper we analyze the vibrations of nonlinear structures by means of the novel approach of isogeometric finite elements. The fundamental idea of isogeometric finite elements is to apply the same functions, namely B-Splines and NURBS (Non-Uniform Rational B-Splines), for describing the geometry and for representing the numerical solution. In case of linear vibrational analysis, this approach has already been shown to possess substantial advantages over classical finite elements, and we extend it here to a nonlinear framework based on the harmonic balance principle. As application, the straight nonlinear Euler–Bernoulli beam is used, and overall, it is demonstrated that isogeometric finite elements with B-Splines in combination with the harmonic balance method are a powerful means for the analysis of nonlinear structural vibrations. In particular, the smoother k-method provides higher accuracy than the p-method for isogeometric nonlinear vibration analysis.     

Dispersion characteristics of wave propagation in layered piezoelectric structures: Exact and simplified models

This article presents a study of the dispersion characteristics of wave propagation in layered piezoelectric structures under plane strain and open-loop conditions. The exact dispersion relation is first determined based on an electro-elastodynamic analysis. The dispersion equation is complicated and can be solved only by numerical methods. Since the piezoelectric layer is very thin and can be modeled as an electro-elastic film, a simplified model of the piezoelectric layer reduces this complex problem to a non-trivial solution of a series of quadratic equations of wave numbers. The model is simple, yet captures the main phenomena of wave propagation. This model determines the dispersion curves of PZT4-Aluminum layered structures and identifies the two lowest modes of waves: the generalized longitudinal mode and the generalized Rayleigh mode. The model is validated by comparing with exact solutions, indicating that the results are accurate when the thickness of the layer is smaller or comparable to the typical wavelength. The effect of the piezoelectricity is examined, showing a significant influence on the generalized longitudinal wave but a very limited effect on the generalized Rayleigh wave. Typical examples are provided to illustrate the wave modes and the effects of layer thickness in the simplified model and the effects of the material combinations.     

The numerical solution of Stokes flow in a domain with re-entrant boundaries by the boundary element method

Numerical solutions are presented for two-dimensional low Reynolds number flow in a rotating tank with stationary barriers. The boundary element method is employed, assuming straight panels and quadratic source distribution. The feasibility of repositioning the nodes as a way to minimize the error is explored. A stretching parameter places smaller elements near the re-entrant regions. Elementary error analysis shows uniform improvement in the solution with stretching. The changing eddy pattern for different numbers and sizes of the barriers is compared with experimental results.     

Selection of finite-element mesh parameters in modeling the growth of hydraulic fracturing cracks

The effect of the mesh geometry on the accuracy of solutions obtained by the finite-element method for problems of linear fracture mechanics is investigated. The guidelines have been formulated for constructing an optimum mesh for several routine problems involving elements with linear and quadratic approximation of displacements. The accuracy of finite-element solutions is estimated based on the degree of the difference between the calculated stress-intensity factor (SIF) and its value obtained analytically. In problems of hydrofracturing of oil-bearing formation, the pump-in pressure of injected water produces a distributed load on crack flanks as opposed to standard fracture mechanics problems that have analytical solutions, where a load is applied to the external boundaries of the computational region and the cracks themselves are kept free from stresses. Some model pressure profiles, as well as pressure profiles taken from real hydrodynamic computations, have been considered. Computer models of cracks with allowance for the pre-stressed state, fracture toughness, and elastic properties of materials are developed in the MSC.Marc 2012 finite-element analysis software. The Irwin force criterion is used as a criterion of brittle fracture and the SIFs are computed using the Cherepanov–Rice invariant J-integral. The process of crack propagation in a linearly elastic isotropic body is described in terms of the elastic energy release rate G and modeled using the VCCT (Virtual Crack Closure Technique) approach. It has been found that the solution accuracy is sensitive to the mesh configuration. Several parameters that are decisive in constructing effective finite-element meshes, namely, the minimum element size, the distance between mesh nodes in the vicinity of a crack tip, and the ratio of the height of an element to its length, have been established. It has been shown that a mesh that consists of only small elements does not improve the accuracy of the solution.     

大直径SHPB弥散效应的二维数值分析

采用轴对称动态有限元HONDO程序对大直径SHPB装置中压杆横向泊松效应引起的应力波弥散进行二维数值分析,并从以下三个方面讨论波形弥散的影响：（1）SHPB装置中压杆直径和杆长对弥散结果（主要是升时）的影响;（2）压杆中的波形弥散对试件应力-应变曲线的影响;（3）弥散对试件应变率的影响。分析表明,在直径SHPB弥散效应对实验结果的影响很大,必须考虑。     

Acoustic simulation using a novel approach for reducing dispersion error

It is well‐known that the traditional finite element method (FEM) fails to provide accurate results to the Helmholtz equation with the increase of wave number because of the ‘pollution error’ caused by numerical dispersion. In order to overcome this deficiency, a gradient‐weighted finite element method (GW‐FEM) that combines Shepard interpolation and linear shape functions is proposed in this work. Three‐node triangular and four‐node tetrahedral elements that can be generated automatically are first used to discretize the problem domain in 2D and 3D spaces, respectively. For each independent element, a compacted support domain is then formed based on the element itself and its adjacent elements sharing common edges (or faces). With the aid of Shepard interpolation, a weighted acoustic gradient field is then formulated, which will be further used to construct the discretized system equations through the generalized Galerkin weak form. Numerical examples demonstrate that the present algorithm can significantly reduces the dispersion error in computational acoustics. Copyright © 2016 John Wiley & Sons, Ltd.     

Synthesis of highly convergent 2D and 3D finite elements for short acoustic wave simulation

New higher-order finite elements of enhanced convergence properties for acoustic wave simulation are presented in the paper. The element matrices are obtained by combining modal synthesis and optimization techniques in order to achieve minimum errors of higher modes of the computational domain. As a result, simulation models of propagating wave pulses require a smaller number of finite element divisions per wavelength compared to the conventional element model thus significantly reducing computational costs. Though finite element matrices are obtained in optimization, the resulting patterns of the matrices are versatile and further can be used in any wave propagation model. The mass matrices of the elements are diagonal, so explicit time integration schemes are applicable. The usage of new elements is especially efficient in situations where wavelengths of the simulated signal are much shorter than the dimensions of the computational domain. This is referred to as short wave propagation analysis. The results of wave propagation simulation for ultrasonic measurements are presented as application examples. The B-scans and computed dispersion curves are provided for visual interpretation of the results.     

透射边界条件在波动谱元模拟中的实现:二维波动

将邢浩洁和李鸿晶提出的多次透射公式(multi-transmitting formula,MTF)的谱元格式应用于均匀介质中线弹性SH波动问题的谱元模拟.假定紧邻人工边界的一层谱单元为具有直线边界的四边形单元,以保证每个人工边界节点都唯一对应一条指向内域的离散网格线.人工边界节点在某时刻的位移由该离散网格线上的节点在前若干时刻的位移确定,按照MTF谱元格式进行计算.通过平面波以一定角度传播的外源问题算例和点源脉冲自由扩散的内源问题算例,验证了方法的可行性以及对实际复杂波动问题的适用性.通过不同类型初值问题算例,在时域内分析了插值多项式阶次、人工波速和透射阶次三个参数对反射误差的影响.结果表明:插值多项式阶次较高的格式会表现出更好的精度,但总体上对反射误差的影响较小;人工波速对反射误差具有显著影响,当人工波速小于介质物理波速时反射误差较大,而当人工波速等于或稍大于介质物理波速时反射误差处于较低水平;透射阶次对反射误差具有决定性影响,表现在不失稳的情形下提高透射阶次能够迅速降低反射误差,但内源问题从三阶MTF开始出现飘移失稳,外源问题从二阶MTF开始出现轻微的飘移失稳.     

用二次形函数薄层法分析弹性层状地基中的动力问题

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法。本文在土层垂直方向离散中利用Galerkin加权残值法推导出二次形函数薄层元的计算公式。采用薄层单元模拟半空间上的层状场地,模型底面用阻尼器边界或傍轴边界代替半空间。利用点源简谐荷载作用下的土层反应与其它数值分析方法的对比讨论薄层模型的设置指标。并用一次和二次两种形函数离散方法计算了层状地基中的面波频散曲线、圆形均布简谐荷载作用下弹性半空间的位移反应和半无限地基中单桩的竖向阻抗函数,分别讨论其计算精度。     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

梯形应力脉冲在弹性杆中的传播过程和几何弥散

在应力波传播过程中,几何弥散效应往往难以避免.对应力波在弹性杆中传播的几何弥散效应进行解析分析,对于基础波动问题研究以及材料动态力学行为表征等课题,显得至关重要.本文简单说明了弹性杆中考虑横向惯性修正的一维 Rayleigh-Love应力波理论,概述了其波动控制方程的变分法推导过程;针对 Hopkinson杆实验中常用的梯形应力加载脉冲,建立了相应的偏微分方程初边值问题的求解模型,并运用 Laplace变换方法研究了脉冲在杆中传播的几何弥散现象;根据留数定理进行 Laplace反变换,给出了杆中不同位置和时刻的应力波的级数形式解析解,分析了计算项数对结果收敛性的影响;将解析计算结果与采用三维有限元数值模拟的计算结果进行对比,两者吻合程度良好,从而证明 Rayleigh-Love横向惯性修正理论可以有效地表征典型 Hopkinson杆实验中的几何弥散效应.在此基础上围绕梯形加载脉冲的弥散效应进行参数研究,定量描述了传播距离、泊松比、脉冲斜率等参数的影响.本文给出的 Rayleigh-Love杆在梯形加载条件下的解析解,揭示了几何弥散效应的本质规律,可以用于实际实验的弥散修正过程.      

On Taylor weak statement finite element methods for computational fluid dynamics

A Taylor series augmentation of a weak statement (a ‘Taylor weak statement’ or ‘Taylor-Galerkin’ method) is used to systematically reduce the dispersion error in a finite element approximation of the one-dimensional transient advection equation. A frequency analysis is applied to determine the phase velocity of semi-implicit linear, quadratic and cubic basis one-dimensional finite element methods and of several comparative finite difference/finite volume algorithms. The finite element methods analysed include both Galerkin and Taylor weak statements. The frequency analysis is used to obtain an improved linear basis Taylor weak statement finite element algorithm. Solutions are reported for verification problems in one and two dimensions and are compared with finite volume solutions. The improved finite element algorithms have sufficient phase accuracy to achieve highly accurate linear transient solutions with little or no artificial diffusion.     

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q ₁ ^rot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H ¹-norm and the pressure in the L ²-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.     

Transient responses of repeated impact of a beam against a stop

The transient behavior of a cantilever beam, driven by periodic force and repeated impacting against a rod-like stop, is the subject of this investigation. As impact and separation phase take place alternately, the transient waves induced either by impacts or by separations will travel in more complicated ways. Thus the transient responses of both the beam and the rod during repeated impact become an important issue. In both impact phase and separation phase, the transient wave propagations are solved by the expansion of transient wave functions in a series of Eigenfunctions (wave modes). From the solutions, the answer of impact force is derived directly, so that the divergence problem, encountered in solving impact force numerically by a strongly non-linear equation coupled the unknown impact force with motions, has been avoided. Numerical results show the convergence of the time-step size and truncation number of wave modes in the calculations of impact force by the present method. As the transient wave effect is considered, the numerical results can show several transient phenomena involving the propagation of transient impact-induced waves, sub-impact phases, long-term impact motion, chatter, sticking motion, synchronous impact, non-synchronous impact (including asynchronous impact) and impact loss.     

Local time-stepping procedures for the space-time conservation element and solution element method

A local time-stepping procedure for the space-time conservation element and solution element (CESE) method has been developed. This new procedure allows for variation of time-step size in both space and time, and can also be extended to become multi-dimensional solvers with structured/unstructured spatial grids. Moreover, it differs substantially in concept and methodology from the existing approaches. By taking full advantage of key concepts of the CESE method, in a simple and efficient manner it can enforce flux conservation across an interface separating grid zones of different time-step sizes. In particular, no correction pass is needed. Numerical experiments show that, for a variety of flow problems involving moving shock and flame discontinuities, accurate and robust numerical simulations can be achieved even with a reduction in time-step size on the order of 10 or higher for grids across a single interface.     

Resonant waves in elastic structured media: Dynamic homogenisation versus Green’s functions

We address an important issue of dynamic homogenisation in vector elasticity for a doubly periodic mass-spring elastic lattice. The notion of logarithmically growing resonant waves is used in the analysis of star-shaped wave forms induced by an oscillating point force. We note that the dispersion surfaces for Floquet–Bloch waves in the elastic lattice may contain critical points of the saddle type. Based on the local quadratic approximations of a dispersion surface, where the radian frequency is considered as a function of wave vector components, we deduce properties of a transient asymptotic solution associated with the contribution of the point source to the wave form. The notion of local Green’s functions is used to describe localised wave forms corresponding to the resonant frequency. The special feature of the problem is that, at the same resonant frequency, the Taylor quadratic approximations for different groups of the critical points on the dispersion surfaces (and hence different Floquet–Bloch vectors) are different. Thus, it is shown that for the vector case of micro-structured elastic medium there is no uniformly defined dynamic homogenisation procedure for a given resonant frequency. Instead, the continuous approximation of the wave field can be obtained through the asymptotic analysis of the lattice Green’s functions, presented in this paper.     

几何精确的非协调等几何NURBS有限元

本文尝试将传统的非协调有限元技术推广到等几何有限元领域,建立了基于精确几何的非协调等几何分析方法,旨在拓展等几何分析应用范围,以便于等几何分析技术能真正实现CAD和FEA的融合,从而真正实现了无需划分网格的目的。我们定义了非协调的NURBS几何（类似非协调元）,给出了NURBS曲面之间几何弱连续的充分条件,进而定义了非协调的等几何分析,将之归纳为带约束驻值问题,并用拉格朗日方法进行求解。两个算例证明这种方法的有效性。未来的工作主要是证明这种方法在不同几何连续性条件下的收敛性以及将之应用到更广的领域。     

Mach Wave Effect on Laminar-Turbulent Transition in Supersonic Flow over a Flat Plate

The effect of a Mach wave (N wave) on laminar-turbulent transition induced by the first instability mode (Tollmien–Schlichting wave) in the flat-plate boundary layer is investigated on the basis of the numerical solution of Navier–Stokes equations at the freestream Mach number of 2.5. In accordance with the experiment, the N wave is generated by a two-dimensional roughness at the computation domain boundary corresponding to the side wall of the test section of a wind tunnel. It is shown that the disturbance induced by the backward front of the N wave in the boundary layer has no effect on the beginning of transition but displaces downstream the nonlinear stage of the first mode development. The disturbance induced by the forward front of the N wave displaces the beginning of transition upstream.