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.

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.

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.

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.

Concurrent fatigue crack growth simulation using extended finite element method

In this paper,a concurrent simulation framework for fatigue crack growth analysis is proposed using a novel small time scale model for fatigue mechanism analysis and the extended finite element method(X-FEM)for fatigue crack growth simulation.The proposed small time scale fatigue model does not require the cycle counting as those using the classical fatigue analysis methods and can be performed concurrently with structural/mechanical analysis.The X-FEM greatly facilitates crack growth simulation without remeshing requirements ahead of the crack tip as in the classical finite element method.The basic concept and theory of X-FEM was briefly introduced and numerical predictions of stress intensity factors are verified with reference solutions under both uniaxial and multiaxial loadings.The small time scale fatigue model is integrated into the numerical simulation algorithm for concurrent fatigue crack growth analysis.Model predictions are compared with available experimental observations for model validation.

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

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

径向井辅助前置液酸压裂缝扩展数值模拟

基于扩展有限元及双尺度连续模型理论建立径向井辅助前置液酸压裂缝扩展模型,引入偏离系数作为量化评价指标,分析各因素对裂缝形态的影响,并利用灰色关联分析方法明确主控因素。结果表明:径向井对前置液酸压裂缝扩展具有显著引导作用,裂缝优先沿径向井方向起裂并扩展一定距离,而后逐渐偏向最大水平主应力方向;径向井方位角、水平地应力差和地层弹性模量越小,径向井长度越长,偏离系数越小,径向井引导效果越好,其中水平地应力差为主控因素;当径向井方位角为15°时,前置液酸压裂缝沿着径向井方向扩展14.32 m后偏向最大水平主应力方向,径向井的引导效果最好,偏离系数为0.005。

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

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

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

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

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

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

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

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

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

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

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

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

基于扩展有限元的现场尺度水力裂缝扩展机制模拟研究

干热岩储层水力压裂形成复杂缝网对于建设增强型地热系统(EGS)至关重要。为了深入认识干热岩储层水力裂缝扩展机制,揭示工程参数和地质参数对水力裂缝扩展的影响规律,本文以青海共和盆地某干热岩储层为研究对象,基于扩展有限元(XFEM)数值模拟软件ABAQUS展开研究,建立现场尺度的干热岩储层模型,实现热流固耦合分析水力裂缝的扩展机制。结果显示:正断型应力状态下,裂缝发生了很明显的朝向Z轴的转向,走滑型应力状态下,裂缝扩展主要沿水平方向,有轻微倾向于Z轴转向的趋势,逆冲型应力状态下,裂缝沿水平方向扩展;随着弹性模量的增大,裂缝的起裂压力减小,裂缝延伸的注水压力也减小,弹性模量越大裂缝起裂的时间越早,泊松比越大水力裂缝的破裂压力越大,破裂压力整体呈现出正断型>走滑型>逆冲型;随着温度的升高,水力裂缝破裂压力减小,大约降低2~3 MPa左右,3个应力状态下都呈现出相同的趋势。本文丰富了干热岩水力裂缝扩展的数值模拟手段,相关研究成果可为干热岩储层中水力压裂扩展预测与不同参数对其的影响分析提供技术支撑。

基于XFEM的大体积结构波动传播规律及裂纹反演方法

大体积混凝土结构被广泛应用于土木、水利等领域的重大工程中,而混凝土抗拉强度低的力学特性决定了其易产生裂纹,因此,发展高效的检测方法,识别大体积混凝土结构中的裂纹信息十分必要.论文提出了一种新的方法,通过提取响应信号频谱中特定频率的幅值特征,基于BP人工神经网络建立幅值特征与裂纹信息间的映射关系,从而有效识别出裂纹信息.首先采用扩展有限元法(eXtended Finite Element Methods, XFEM)和人工吸收边界模型,分别模拟了单裂纹和双裂纹情形下,大量不同裂纹信息下特定位置传感器的响应,分析其频谱曲线并提取特征,建立频谱特征—裂尖位置数据集,以训练人工神经网络,测试集的反演效果显示,该方法具有较好的准确度,可有效识别出裂纹信息.

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

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