三维弹塑性问题的比例边界有限元法

为了研究三维复杂结构的非线性硬化问题,在考虑关联流动准则的Von Mises塑性屈服准则和各向同性/随动硬化准则的本构模型中,引入一致切线模量,并基于比例边界有限元方法和平衡八叉树算法对本构模型进行求解。采用Newton-Raphson迭代法求解位移和应力的非线性方程,为了减少计算量,仅在比例中心点执行并储存迭代过程的变量,同时采用Fortran语言编制了求解程序。两个算例计算结果表明,该模型的计算结果具有较高的准确性。

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

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

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

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

第五届比例边界有限元法最新进展国际研讨会在我校成功召开

2023年11月24日至26日,由我校、澳大利亚新南威尔士大学和河海大学联合主办,我校土木工程防灾减灾湖北省引智创新示范基地和土木与建筑学院承办,大连理工大学、宁波大学协办的“第五届比例边界有限元法最新进展国际研讨会(SBFEM-2023)”在我校成功召开。本次研讨会由大连理工大学林皋院士担任名誉主席,澳大利亚新南威尔士大学宋崇民教授、三峡大学陈灯红教授、河海大学杜成斌教授担任主席。

稳态非线性热传导问题的比例边界有限元法

提出一种基于比例边界有限元法和Kirchhoff变换求解稳态非线性热传导问题的数值方法。为了消除导热系数随温度变化引起的非线性,引入Kirchhoff变换将非线性的偏微分控制方程简化为Laplace方程,因而可消除迭代计算。作为一种兼顾了有限元法和边界元法优势的半解析数值方法,比例边界有限元法只需在计算域的边界上划分网格且不需基本解。在变换空间采用比例边界有限元法求得高精度的数值解之后,借助Kirchhoff反变换可求得温度场。数值结果表明,所提的稳态非线性热传导问题分析方法是行之有效的。

基于比例边界有限元法的成层半空间中斜坡场地地震响应计算模型

由于左右侧不对称且不规则的边界条件延伸到远场,斜坡地形场地的地震波动场求解存在一定困难。为得到成层半空间中斜坡地形场地的地震动力响应,基于波场理论提出一种以相似拼接线为缩放中心的比例边界有限元法计算模型:首先将待求波动场分解为具有规则边界条件的已知波场和由已知波场在真实不规则边界处引起的散射波场,然后通过在斜坡场地的不规则边界处施加等效地震荷载将散射波场的求解由外源输入问题转化为内源辐射问题,最后推导了适用于水平及倾斜成层、左右岸不对称等复杂地基的比例边界有限元法,可沿径向解析地求解斜坡场地的内源辐射问题。通过与文献中均质和成层弹性半空间中凹陷地形场地在SH波入射条件下地表动力响应进行对比,表明了计算模型的准确性,并通过分析成层半空间中斜坡地形场地的波动响应验证了模型的有效性。提出的计算模型,为复杂地基-结构相互作用分析的输入波动场计算提供了新的技术手段。

静电场分析的比例边界有限元法

比例边界有限元法（SBFEM）是一种半解析数值分析的新方法，集合了有限元法和边界元法的优点，又具有独特的优点．在其辐射坐标上保持了解析性，因此其模拟精度较高，另外可以自动满足无限远的边界条件．从拉普拉斯方程出发，利用加权余量法并通过比例坐标和笛卡儿坐标变换，推导出静电场分析的比例边界有限元方程、电位求解公式以及电场求解公式．算例计算结果与解析解和其他数值方法比较结果表明，此方法具有精度高、计算工作量小的优点．

应用比例边界有限元法求解狭缝对双箱水动力的影响

比例边界有限元法(SBFEM)是一种半解析数值分析的新方法,既融合了有限元法和边界元法的优点,又有其特有的优点。用该方法可求解有限水深下狭缝对双箱水动力作用的影响,为波浪与多浮体超大型结构的相互作用探索一些规律。整个计算域划分成2个无限子域和4个有限子域,并利用加权余量法在各个子域上推导了SBFEM的积分方程;计算了4个数值算例并与边界元等其它数值方法进行了比较,验证了该方法是一种用很少单元便能得到精确结果的高效方法。应用SBFEM对不同箱体宽度、不同狭缝宽度、不同吃水深度条件的双箱作了计算,得出了狭缝对双箱水动力干涉影响的一些规律,对超大型浮体水动力分析和结构设计具有一定的参考价值。

应用比例边界有限元法研究波浪与带狭缝三箱作用的共振现象

比例边界有限元法(SBFEM)是一种半解析的新数值分析方法,比完全数值方法具有更高的分析精度。本文应用该方法求解波浪与有限水深中三个相同,但之间存在狭缝的二维固定方箱的作用问题,为波浪与多模块超大型浮体的相互作用提供一些规律。文中给出了各箱的波浪力及反射、透射系数的数值结果,与BEM结果进行了比较,吻合良好。发现波浪与三箱的作用存在双共振现象而且中间箱的共振形式与两边箱不同。进一步分析了箱体吃水深度、狭缝宽度和箱体宽度对三箱共振频率的影响。这对超大型浮体水动力分析和结构设计具有一定的参考价值。数值结果证实该方法用较少单元便能得到精确结果,这为超大型浮体的水动力分析提供了一个可行的分析方法。

饱和土体固结3D比例边界有限元法分析

基于饱和土体Biot固结理论,考虑土体体力与表面作用力,结合饱和土体中土骨架动力平衡方程与孔隙流体的连续方程,在计算域内采用Galerkin加权余量法,推导了3D无限、有限域全耦合饱和土体固结的比例边界有限元方程.理论推导表明:与单相弹性土体不同的是,由于孔隙水的存在,饱和土体固结比例边界有限元法方程中不仅有位移、应力矩阵,而且还存在孔隙水压力的影响耦合矩阵;该方程不仅能够满足无限域无穷远边界辐射条件,而且在径向(与饱和土体-结构接触面垂直方向,指向无限远)上具有严格的精确性,环向上(与饱和土体-结构接触面平行方向)具有有限元单元无限收敛的准确性.

比例边界有限元法在电机电磁场问题中的应用

将比例边界有限元法应用到电机磁场分析计算中,采用加权余量法推导无旋磁场的比例边界有限元方程及其离散格式,给出相邻子域公共边上内部节点磁位求解方法.以具有解析解的电机磁场问题为例进行计算和对比分析,验证离散格式的正确性和比例边界有限元法求解电机磁场的可行性,为该方法在电机磁场问题中的应用打下基础.

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite element methods,XFEM)和比例边界有限元法(scaled boundary finite element methods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improved extended scaled boundary finite element methods,i XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的i XSBFEM的有效性,相比于常规XFEM,i XSBFEM的基于位移范数的相对误差收敛性较好;采用i XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合较好.

移动荷载作用下半无限弹性空间中地铁隧道动力响应的频域—波数域比例边界有限元法分析

基于比例边界有限元法与傅里叶积分变换,分析移动荷载作用下半无限弹性空间中地铁隧道动力响应。考虑到移动荷载运动方向与地铁隧道轴线方向的一致性,将半无限弹性空间动力控制方程进行频域—波数域的傅里叶积分变换;在变换域内,隧道径向采用比例坐标系,环向采用有限元意义离散,在计算域内利用Galerkin法,建立频域—波数域内的比例边界有限元方程,进而推导出半无限空间动力刚度的一阶微分矩阵方程;再利用高频渐近展开和傅里叶积分逆变换,得到时间—空间域系统的动力响应。该方法极大地减小了计算分析量,同时避免了无限边界的计算误差。利用该方法进行实例计算的结果表明:随着荷载移动速度的增大,地铁隧道振动波传播到土体表面,引起的土体振动有放大增强的现象,将会对地铁隧道上部结构的安全性造成一定影响。

半无限弹性空间中移动荷载动力响应的频域-波数域比例边界有限元法分析

基于Fourier积分变换和虚功原理,形成了频域-波数域比例边界有限元法,分析了移动荷载作用下半空间域弹性空间动力响应。首先对半无限域弹性体的动力控制方程进行时间到频域,荷载移动方向的空间域到波数的Fourier积分变换,然后选择比例中心,利用虚功原理,在地铁隧道孔洞横截面环向上采用有限元法意义离散,建立了频域-波数域比例边界有限元方程,进而形成了一阶微分矩阵方程形式的半无限空间动力刚度。文中理论推导表明:利用文中方法分析半无限域中沿地铁隧道结构纵轴向的移动荷载动力响应问题,不仅可避免无穷边界计算处理误差,而且可极大减小计算分析量。计算结果表明:半无限弹性地基的振动响应随移动荷载速度增大而增大,尤其是当荷载速度增大到土体剪切波速后,振动波传播到土体表面引起土体振动显著增大,土体振动性增大,将会对土体及表面结构的安全性形成一定影响,另一方面土体的振动在沿地铁隧道纵轴向的衰减比竖向慢。

基于图像四叉树的改进型比例边界有限元法研究

为提高数值计算的精度,断裂力学问题的数值模拟需要在裂纹扩展的局部区域采用较密的网格,而远离裂纹扩展的区域可采用较疏的网格,且对于裂纹扩展问题的数值模拟,大多数数值方法又存在局部网格重剖分的问题.论文提出了一种基于图像四叉树的改进型比例边界有限元法用于模拟裂纹扩展问题,该方法可根据结构域几何外边界的图像全自动进行四叉树网格剖分,无需任何人工干预,网格剖分效率极高,由于比例边界有限元法本身的优势,四叉树网格的悬挂节点可以直接地视为新的节点,无需任何特殊处理.通过引入虚节点的思想,将裂纹与四叉树单元边界交叉点作为虚节点,虚节点的自由度作为附加自由度处理,并采用水平集函数表征材料内部的裂纹面,含不连续裂纹面的子域可通过节点水平集函数识别,使得裂纹扩展时无需进行网格重剖分,界面的几何特征通过比例边界有限元子域的附加自由度表征.最后,通过若干算例验证了该方法的性能,建议的改进型比例边界有限元法在求解复合型应力强度因子和模拟材料内部裂纹扩展路径时均具有较高的精度.

基于比例边界有限元法的结构-地基动力相互作用时域算法的改进

本文基于比例边界有限元法(Scaled Boundary Finite Element Method)的缩减基函数解法,对结构-无限地基动力相互作用的时域算法进行了改进。通过选择适当的基函数数目,缩减结构-地基接触面上的自由度,以减小卷积积分所带来的计算工作量,推导了缩减自由度后运动方程的表达式。通过重力坝和拱坝加速度频响曲线的算例,对比了不同程度缩减的基函数和全部基函数对计算精度和效率的影响。结果表明,基函数的缩减可使计算效率明显提高,但精度损失不大。当采用60%的基函数时,计算效率提高5倍,而峰值频响曲线的精度损失却不超过4%。因此,该算法比较适合于大型结构-地基动力相互作用问题的时域分析。

应用比例边界有限元法求解多种二维浮体水动力特性

比例边界有限元法(SBFEM)是一种半解析的数值方法,比完全数值方法具有更高的精度,该方法结合了有限元和边界元的优点,采用相对少的剖分单元就可以得到较高精度的模拟结果。通过改变有限子域内部比例中心的位置,使这种方法可以应用到多种形式浮体在波浪作用下的水动力特性的计算中。同时还给出了各种形式浮体的波浪力及反射、透射系数的数值结果,并与边界元方法(BEM)计算结果和特征函数展开方法得到解析解进行了比较,均吻合良好。研究表明比例边界有限元不仅可以计算矩形的浮体结构,而且对于多种结构形式的浮体都可以计算,这为多种结构形式浮体的水动力分析提供了一个可行的方法。

基于比例边界有限元法和灰狼优化算法的裂纹尖端位置识别

基于比例边界有限元法(SBFEM)和灰狼优化(GWO)算法,提出了一种裂纹尖端识别方法.首先,借助SBFEM解决断裂力学问题特有的优势,快速准确地计算出反演所需的测点位移,并验证了正问题求解的正确性.其次,建立与裂纹尖端位置有关的目标函数,将求解裂纹尖端位置转换为求解目标函数最小值的优化问题.最后,采用GWO算法对目标函数进行了优化,进而搜索裂纹尖端的最佳位置.数值算例结果表明:利用SBFEM的高精度、半解析的优点,在反演过程中采用其求解正问题是非常有效的;GWO算法具有良好的全局收敛性,且相比经典的粒子群算法,能够更快速准确地搜索出裂纹尖端的位置;GWO算法具有较好的抗噪性.

瞬态分析中的比例边界有限元法递推公式

为改进现有比例边界有限元法卷积积分公式的计算效率,推导了一种新的递推公式。本文求解了比例边界有限元法的动态刚度矩阵,并将其分解成高频和低频两部分;对低频部分利用无约束高斯-牛顿迭代算法建立其对应的有理式公式,确定递推系数;将高频和低频两部分通过傅里叶逆变换建立时域公式,其中高频部分对应阻尼矩阵,低频部分对应刚度矩阵与递推项。结果表明:当递推阶数m=15时,递推公式的解与解析解曲线形状相似,在最大值处的相对误差为0.8%;当分析步数为3183时,卷积积分公式的计算时间为递推公式的14倍。此外,在瞬态分析中该递推公式可用常规有限元法求解;通过具有解析解的坝库耦合算例,验证了该递推公式的正确性和有效性。

基于比例边界有限元法的高阶透射边界应用

基于比例边界有限元法和连分式展开推导了无限域弹性动力分析的求解方程,实现了一种局部的高阶透射边界.采用改进的连分式法求解无限域的动力刚度矩阵,克服了原连分式算法可能会造成矩阵运算病态的问题.该局部高阶透射边界在时域里表示为一阶常微分方程组,其稳定性取决于其系数矩阵的广义特征值问题.如果出现虚假模态,采用移谱法来校正系数矩阵以消除虚假模态.通过两个算例验证了该高阶透射边界的精确性、鲁棒性.