Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

Dynamic Crack Propagation Analysis Using Scaled Boundary Finite Element Method

The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.

A New Formulation of the Scaled Boundary Finite Element Method for Heterogeneous Media:Application to Heat Transfer Problems

The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.

ALE FINITE ELEMENT ANALYSIS OF NONLINEAR SLOSHING PROBLEMS BASED ON F LUID VELOCITY POTENTIAL

Base d on fluid velocity potential, an ALE finite element formulation for the analysi s of nonlinear sloshing problems has been developed. The ALE kinemat ical description is introduced to move the computational mesh independently of f luid motion, and the container fixed noninertial coordinate system is employed to establish the governing equations so that the mesh is needed to be updated in this coordinate system only. This leads to a very simple mesh moving algorithm which makes it easy to trace the motion of the moving boundaries and the free su rface without producing undesirable distortion of the computational mesh. The fi nite element method and finite difference method are used spacewise and timewise , respectively. A numerical example involving either forced horizontal oscillati on or forced pitching oscillation of the fluid filled container is presented to illustrate the effectiveness and the robustness of the method. In additi on, this work can be extended for the fluid structure interaction problems.

基于SBFEM的剪切型裂纹动态开裂模拟

剪切型断裂是岩土工程中常见的破坏模式,了解剪切破坏机理并准确预测剪切型裂纹的萌生、扩展过程对保障工程结构的安全性与稳定性具有重要意义.文章建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)和非局部宏-微观损伤模型的剪切型裂纹动态开裂模拟方法,定义了基于偏应变概念的物质点对的正伸长量,可作为预测剪切型裂纹扩展行为的动态开裂准则,一点的损伤定义为该点影响域范围内连接的物质键损伤的加权平均值,而物质键的损伤则与基于偏应变概念的物质点对的正伸长量相关联,并引入能量退化函数建立结构域几何拓扑损伤与能量损失之间的关系,将拓扑损伤与应力应变联系起来,通过能量退化函数修正了SBFEM的刚度系数矩阵,得到了子域在损伤状态下的刚度矩阵,推导了考虑结构损伤的SBFEM动力控制方程,采用Newmark隐式算法对控制方程进行时间离散.最后,通过3个典型算例验证了建议的模型可较好地模拟剪切型断裂问题,能够很好地捕捉剪切型裂纹的扩展路径,并得到较为准确的载荷-位移曲线.

考虑材料参数空间变异性的SBFEM开裂模拟

论文建立了基于比例边界有限元法(scaled boundary finite element methods,SBFEM)框架的非局部宏微观损伤模型,考虑材料细观物理参数的空间变异性,探讨了材料参数的空间变异性对结构开裂过程的影响。结果表明:考虑材料参数空间变异性后,裂纹扩展路径具有不确定性,建议的模型能够很好地反应材料内在的随机性;随着结构受力情况的复杂化和结构本体缺陷的增多,裂纹开裂模式的变异性也会增大。自相关长度和参数变异系数对结构开裂分析结果有重要影响。

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Modeling Technology in Finite Element Analysis of Electrostatic Proximity Fuze Problem

In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.

基于ABAQUS的SBFEM及在混凝土坝非线性地震响应分析中的应用

通过用户子程序接口,将比例边界有限元法嵌入到商用程序ABAQUS中.采用蒙皮单元技术完成了接触界面属性定义,将比例边界有限元法与ABAQUS中的接触非线性模型相结合.通过耦合比例边界有限元法-经典有限元法,采用黏弹性边界建立了考虑横缝接触非线性的坝体-库水-地基动力相互作用模型.将接触非线性混凝土坝模型应用到NG5拱坝中,计算结果均与ABAQUS自带单元的结果吻合良好,并具有较高的计算效率,从而为工程结构动力分析提供了一种新手段.

NONLINEAR GALERKIN METHOD FOR THE EXTERIORNONSTATIONARY NAVIER-STOKES EQUATIONS

A new algorithm combining nonlinear Galerkin method and coupling method of finite element and boundary element is introduced to solve the exterior nonstationary Navier-Stokes equations. The regularity of the coupling variational formulation and the convergence of the approximate solution corresponding to the algorithm are proved. If the fine mesh h is chosen as coarse mesh H-square, the nonlinear Galerkin method, nonlinearity is only treated on the coarse grid and linearity is treated on the fine grid. Hence, the new algorithm can save a large amount of computational time.

数据驱动的半无限介质裂纹识别模型研究

缺陷识别是结构健康监测的重要研究内容,对评估工程结构的安全性具有重要的指导意义,然而,准确确定结构缺陷的尺寸十分困难.论文提出了一种创新的数据驱动算法,将比例边界有限元法(scaled boundary finite element methods,SBFEM)与自编码器(autoencoder,AE)、因果膨胀卷积神经网络(causal dilated convolutional neural network,CDCNN)相结合用于半无限介质中的裂纹识别.在该模型中,SBFEM用于模拟波在含不同裂纹状缺陷半无限介质中的传播过程,对于不同的裂纹状缺陷,仅需改变裂纹尖端的比例中心和裂纹开口处节点的位置,避免了复杂的重网格过程,可高效地生成足够的训练数据.模拟波在半无限介质中传播时,建立了基于瑞利阻尼的吸收边界模型,避免了对结构全域模型进行计算.搭建了CDCNN,确保了时序数据的有序性,并获得更大的感受野而不增加神经网络的复杂性,可捕捉更多的历史信息,AE具有较强的非线性特征提取能力,可将高维的原始输入特征向量空间映射到低维潜在特征向量空间,以获得低维潜在特征用于网络模型训练,有效提升了网络模型的学习效率.数值算例表明:提出的模型能够高效且准确地识别半无限介质中裂纹的量化信息,且AE-CDCNN模型的识别效率较单CDCNN模型提高了约2.7倍.

基于SBFEM研究T应力对岩石起裂特性的影响

传统断裂判据在分析岩石裂纹尖端起裂特性时忽略了非奇异项参数的影响,从而导致理论结果与试验结果产生较大的偏差。比例边界有限元法(Scaled boundary finite element method-SBFEM)是一种半解析的数值计算方法,可以较精确的计算实数、复数以及幂-对数奇异性问题同时不需要对裂纹尖端进行特殊处理,应力强度因子等裂尖奇异参数都可以根据定义直接进行提取。基于此方法对复合型岩石裂纹尖端奇异参数I、II型应力强度因子和T应力进行了提取,并建立了考虑T应力的断裂准则,对不同裂纹形式下岩石的起裂特性进行研究,并与现有试验结果进行对比。结果表明,T应力为正值时会加剧裂纹扩展,T应力为负值时会抑制裂纹扩展,可见裂纹尖端的非奇异参数即T应力的影响不容忽视。

An adaptive scaled boundary finite element method by subdividing subdomains for elastodynamic problems

The scaled boundary finite element method(SBFEM) is a semi-analytical numerical method,which models an analysis domain by a small number of large-sized subdomains and discretises subdomain boundaries only.In a subdomain,all fields of state variables including displacement,stress,velocity and acceleration are semi-analytical,and the kinetic energy,strain energy and energy error are all integrated semi-analytically.These advantages are taken in this study to develop a posteriori h-hierarchical adaptive SBFEM for transient elastodynamic problems using a mesh refinement procedure which subdivides subdomains.Because only a small number of subdomains are subdivided,mesh refinement is very simple and efficient,and mesh mapping to transfer state variables from an old mesh to a new one is also very simple but accurate.Two 2D examples with stress wave propagation were modelled.The results show that the developed method is capable of capturing propagation of steep stress regions and calculating accurate dynamic responses,using only a fraction of degrees of freedom required by adaptive finite element method.

LOCALLY STABILIZED FINITE ELEMENT METHOD FOR STOKES PROBLEM WITH NONLINEAR SLIP BOUNDARY CONDITIONS

Based on the low-order conforming finite element subspace （Vh, Mh） such as the P1-P0 triangle element or the Q1-P0 quadrilateral element, the locally stabilized finite element method for the Stokes problem with nonlinear slip boundary conditions is investigated in this paper. For this class of nonlinear slip boundary conditions including the subdifferential property, the weak variational formulation associated with the Stokes problem is an variational inequality. Since （Vh, Mh） does not satisfy the discrete inf-sup conditions, a macroelement condition is introduced for constructing the locally stabilized formulation such that the stability of （Vh, Mh） is established. Under these conditions, we obtain the H1 and L2 error estimates for the numerical solutions.

MODE III 2-D FRACTURE ANALYSIS BY THE SCALED BOUNDARY FINITE ELEMENT METHOD

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.

基于比例边界有限元法计算应力强度因子的不确定量化分析

应力强度因子是预测荷载作用下结构中裂纹产生和扩展的重要参数。半解析的比例边界有限元法结合了有限元和边界元法的优势,在裂纹尖端或存在奇异应力的区域不需要局部网格细化,可以直接提取应力强度因子。在比例边界有限元法计算应力强度因子的框架下,引入随机参数进行蒙特卡罗模拟(Monte Carlo simulation, MCS),并提出一种新颖的基于MCS的不确定量化分析。与直接的MCS不同,采用奇异值分解构造低阶的子空间,降低系统的自由度,并使用径向基函数对子空间进行近似,通过子空间的线性组合获得新的结构响应,实现基于MCS的快速不确定量化分析。考虑不同荷载状况下,结构形状参数和材料属性参数对应力强度因子的影响,使用改进的MCS计算应力强度因子的统计特征,量化不确定参数对结构的影响。最后通过若干算例验证了该算法的准确性和有效性。

重力坝地震断裂的多边形比例边界有限元模型研究

为研究重力坝的地震断裂破坏机理,基于线弹性断裂力学和多边形比例边界有限元法(Polygon SBFEM)提出一种全自动的重力坝动态断裂分析模型。该模型继承了SBFEM在断裂分析中高精度和高效率的优势,通过建立多边形比例边界有限单元广义动态应力强度因子的时域分析方法实现任意时刻断裂参数的求解,并结合裂缝扩展准则可实时判定裂缝稳定性;对达到临界状态的裂缝采用多边形网格局部重剖分技术,开发了考虑裂缝张开-闭合行为的动态接触模拟算法,进而实现了裂缝动态扩展高效自动化模拟。以Koyna重力坝为研究对象,考虑大坝库水动力相互作用,模拟了地震作用下坝体裂缝扩展过程,获得了断裂路径,验证了模型的正确性。其后,进一步探讨了数值模型中网格密度和裂缝扩展步长等因素的影响,结果显示仅用粗网格即可获得满意的计算结果,且三种裂缝扩展步长模拟的大坝裂缝扩展轨迹接近,但随着裂缝扩展步长增大,大坝更早失效破坏。研究成果可为混凝土坝的地震断裂分析提供有力的技术手段。

基于SBFEM的坝-库水相互作用分析

将W o lf和Song提出的比例边界有限元法(sca led boundary fin ite e lem en t m ethod,简称SBFEM)应用于坝-库水相互作用分析.在库水不可压缩假定的前提下推导了坝面动水压力与附加质量矩阵的基本方程.二维重力坝和三维拱坝坝面动水压力的算例表明,与有限元法比较,SBFEM计算精度明显提高,同时由于具有降维的特点,计算工作量也在很大程度上有所节约.还分析了库底和水面坡度变化以及水库边界形状适当变化对动水压力分布和数值大小的影响.结果表明,一般情况下水库边界的变化对动水压力的分布形状基本上不发生影响,对动水压力数值大小也影响不大.