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.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

A scaled boundary node method applied to two-dimensional crack problems

A boundary-type meshless method called the scaled boundary node method （SBNM） is developed to directly evaluate mixed mode stress intensity factors （SIFs） without extra post-processing. The SBNM combines the scaled boundary equations with the moving Kriging （MK） interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter. As a result, the SBNM requires only a set of scattered nodes on the boundary, and the displacement field is approximated by using the MK interpolation technique, which possesses the 5 function property. This makes the developed method efficient and straightforward in imposing the essential boundary conditions, and no special treatment techniques are required. Besides, the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction. Therefore, the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip. Numerical examples using the SBNM for computing the SIFs are presented. Good agreements with available results in the literature are obtained.

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.

The Calculation of Stress-Strain State of Anisotropic Composite Finite-Element Area with Different Boundary Conditions on the Surface

The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.

STRAIGHTFORWARD MULTI-SCALE BOUNDARY ELEMENT METHOD FOR GLOBAL/LOCAL MECHANICAL ANALYSIS OF ELASTIC HETEROGENEOUS MATERIAL

A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.

COMPUTER COMPUTATION OF THE METHOD OF MULTIPLE SCALES-DIRICHLET PROBLEM FOR A CLASS OF SYSTEM OF NONLINEAR DIFFERENTIAL EQUATIONS

The method of boundary layer with multiple scales and computer algebra were applied to study the asymptotic behavior of solution of boundary value problems for a class of system of nonlinear differential equations . The asymptotic expansions of solution were constructed. The remainders were estimated. And an example was analysed. It provides a new foreground for the application of the method of boundary layer with multiple scales .

APPLICATION OF THE MODIFIED METHOD OF MULTIPLE SCALES TO THE BENDING PROBLEMS FOR CIRCULAR THIN PLATE AT VERY LARGE DEFLECTION ANDTHE ASYMPTOTICS OF SOLUTIONS (Ⅰ)

In this paper, the modified method of multiple scales is applied to study the bending problems for circular thin plate with large deflection under the hinged and simply supported edge conditions. Theseries solutions are constructed, the boundary layer effects are analysed and their asymptotics are proved.

Statistical second-order two-scale analysis and computation for heat conduction problem with radiation boundary condition in porous materials

This paper discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configuration for the structure with random distribution is briefly characterized.Secondly,the SSOTS formulae for computing the heat transfer problem are derived successively by means of the construction way for each cell.Then,the statistical prediction algorithm based on the proposed two-scale model is described in detail.Finally,some numerical experiments are proposed,which show that the SSOTS method developed in this paper is effective for predicting the heat transfer performance of porous materials and demonstrating its significant applications in actual engineering computation.

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

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

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

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

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

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

基于改进四叉树和比例边界有限元法的自适应设计域拓扑优化方法

针对大型结构拓扑优化计算成本高和固体各向同性材料惩罚模型(SIMP)在优化后结构边界处求解精度低的问题,提出了一种基于SIMP法的自适应设计域(ADD)拓扑优化方法。将改进四叉树法应用在拓扑优化的过程中,通过自动划分不同等级的网格单元来减少网格数量、减轻计算负担并提高边界处求解精度;采用比例边界有限元法(SBFEM)实时计算划分后结构的有限元信息,解决了不同等级网格间悬挂节点的问题。所提方法可在初始网格相对较少的情况下得到更加精确的结果,大幅度地降低了计算成本。数值算例结果表明,所提方法在最终结构边界处精度相同的情况下,计算时间最快可缩短为原来的1/100,可以为后续降低大型结构拓扑优化的计算成本提供参考。

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

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

基于图像八叉树的三维比例边界有限元多面体网格生成算法

基于图像八叉树方法,提出了平衡八叉树和多面体网格修剪相结合的三维比例边界有限元多面体网格算法,该算法根据结构尺寸建立恰好完全包含整个结构的立方体网格,再建立结构图像像素信息,根据像素信息,按照2∶1的平衡分割原则递归地进行等分,完成平衡八叉树网格生成。结构内单元网格完全保留,结构外单元网格删除,对于结构边界单元网格,提出采用面-立方体相交判断方法进行边界单元与结构边界相交面的筛选,搜寻结构单元与结构边界表面相交点,通过有序连接相交点形成边界单元切割面,再结合边界单元其他几个面,构成裁剪后的多面体单元。数值算例结果表明,基于本文算法生成的比例边界有限元网格计算结果具有较好的精度和边界适应性。

基于快速多极杂交边界点法的三维弹性力学问题求解

杂交边界点法是近年来快速发展的一种无网格方法,国内外学者已成功将快速多极算法与边界元方法相结合,使得边界元的应用范围进一步扩大。为此,开展了将快速多极算法与杂交边界点法相结合计算三维弹性力学问题的研究。三维弹性力学的基本解利用球谐函数展开为级数,并通过自适应八叉树结构将求解域分解为分级不相邻的区域。数值算例表明,快速多极杂交边界点法具有很好的计算精度和计算效率,该算法在普通个人电脑上,可以在令人接受的时间内完成40万以上自由度的计算。因此,该算法在大规模问题的计算中具有显著优势,对于复合材料、加筋混凝土等更为复杂材料的力学问题具有实用性。

基于比例边界有限元的复合梁自由振动频率计算

将比例边界有限元方法(SBFEM)拓展用于计算复合梁的自由振动频率.该方法将梁简化为一维模型,并且仅选用x和z方向的弹性线位移作为基本未知量.从弹性力学基本方程出发,通过比例边界坐标、虚功原理和对偶变量技术推导得到了复合梁的一阶常微分比例边界有限元动力控制方程,其通解为解析的矩阵指数函数.利用Padé级数求解矩阵指数函数可得各个梁层的动力刚度矩阵,根据自由度匹配原则组装得到复合梁的整体刚度和质量矩阵.求解特征值方程,最终可得复合梁的自由振动频率.该方法对复合梁的层数和边界条件均无限制,具有广泛的适用性.将该文的解与三层、四层和十层复合梁振动频率的数值参考解以及阶梯型悬臂梁固有频率的实验实测值进行对比,验证了比例边界有限元算法的准确性、高效性和快速收敛性.