文章快速检索    
  同济大学学报(自然科学版)  2019, Vol. 47 Issue (7): 906-913.  DOI: 10.11908/j.issn.0253-374x.2019.07.002
0

引用本文  

吴振丽, 陆洲导, 李凌志, 张佳旭. 梁侧锚固钢板法加固混凝土梁的非线性有限元分析[J]. 同济大学学报(自然科学版), 2019, 47(7): 906-913. DOI: 10.11908/j.issn.0253-374x.2019.07.002.
WU Zhenli, LU Zhoudao, LI Lingzhi, ZHANG Jiaxu. Non-linear Finite Element Analysis of Bolted Side-plated Beams[J]. Journal of Tongji University (Natural Science), 2019, 47(7): 906-913. DOI: 10.11908/j.issn.0253-374x.2019.07.002

基金项目

国家自然科学基金(51378397)

第一作者

吴振丽(1988—),女,博士生,主要研究方向为结构抗火及加固.E-mail:wzlsdd@tongji.edu.cn

通信作者

李凌志(1980—),男,副教授,工学博士,主要研究方向为结构抗火及加固.E-mail:lilingzhi@tongji.edu.cn

文章历史

收稿日期:2018-07-11
梁侧锚固钢板法加固混凝土梁的非线性有限元分析
吴振丽 , 陆洲导 , 李凌志 , 张佳旭     
同济大学 土木工程学院, 上海 200092
摘要:以4根梁侧锚固钢板加固混凝土梁(BSP梁)受剪性能试验研究为基础,基于弹塑性有限元软件OpenSees,建立了BSP梁的非线性有限元模型,模拟了构件加载全过程和受剪破坏时的受力性能.模型中考虑了材料的非线性、钢板和混凝土梁交界面上滑移等.计算结果与试验结果吻合良好.并对影响BSP梁性能的主要参数进行分析,包括锚栓行数、纵向锚固间距、钢板高度、钢板厚度、钢板屈服强度、锚栓直径及混凝土强度等,得出了各参数对BSP梁受剪承载力和纵横向滑移的影响规律,为BSP法加固设计提供了依据.
关键词梁侧锚固钢板    滑移    受剪承载力    有限元分析    
Non-linear Finite Element Analysis of Bolted Side-plated Beams
WU Zhenli , LU Zhoudao , LI Lingzhi , ZHANG Jiaxu     
College of Civil Engineering, Tongji University, Shanghai 200092, China
Abstract: Based on the experimental study of 4 bolted-side plates (BSP) beams, a numerical finite element analysis using OpenSees was conducted on the entire loading process and the shear failure of BSP beams.The nonlinear behaviors of materials and slip model on the interface between steel plate and concrete beam were taken into account. The results of the finite element analysis were found to be in good agreement with the experimental results. The main parameters affecting the behavior of BSP beams were analyzed, including the number of rows of anchor bolts, anchor bolt spacing along the beam, steel plate depth, steel plate thickness, yield strength of steel plate, anchor bolt diameter, concrete strength and so on. The results show the influence of the parameters on the shear bearing capacity, longitudinal and transverse slip, which provides basis for the strengthening design of BSP beams.
Key words: bolted side-plating (BSP)    slip    shear bearing capacity    finite element analysis    

对现有建筑功能的多样化需求,导致结构布置改变或者出现超载,从而需对其进行加固处理.对于钢筋混凝土梁而言,常用加固方法有采用混凝土或其他水泥基材料增大截面[1-2]、黏钢[3]或黏碳纤维加固[4-5]等.大量研究表明黏钢和黏碳纤维[3, 6-7]能显著提高梁的抗弯承载力和延性,但是在剪应力作用下由于锚固不足容易发生剥离破坏.梁侧锚固钢板法[8] (简称BSP法)是通过锚栓将钢板固定在梁侧面从而使其达到加固的技术.Barnes等人[9]对“强弯弱剪”梁分别在其梁侧采用黏钢法和BSP法进行加固,对比发现两者均能显著提高梁承载力,但是黏钢法加固效果低于BSP法,主要原因为黏贴钢板端部锚固不足,在剪应力作用下发生了剥离.

国内外学者对BSP梁的受力性能进行了一系列的研究[10-21],发现由于螺栓的剪切变形,钢板与混凝土间存在纵向和横向相对滑移,其量值由钢板与混凝土梁的刚度比及螺栓的滑移响应决定[14-16].由于相对滑移的存在,BSP梁抗弯刚度降低、延性变大,而受弯承载力基本上保持不变[17-18].而由于BSP梁受剪时,主斜裂缝接近梁端部,钢板端部螺栓数量不足将使螺栓发生较大滑移甚至钢板端锚失效,导致受剪承载力严重降低,甚至受剪破坏模式的转变[19-21], 因此在数值分析中必须考虑锚栓滑移效应.

目前对BSP梁受剪性能的研究基本为试验研究,但由于进行试验所需成本较高,且受到试验设备、材料等外界因素的限制,尚未对其受剪性能进行全面参数分析,且实验过程中很难测得钢板和混凝土梁之间的相对滑移.而数值分析作为重要的研究手段,结合试验研究,建立准确的数值模型,能为其提供详细的研究依据.

有限元分析软件OpenSees是由加州大学伯克利分校主导研发的一款开源分析软件.国内外众多学者已基于此软件对混凝土构件进行了力学响应分析.陆新征[22-23]等人开发了分层壳单元,该单元基于复合材料力学原理,将一个壳单元划分为很多层,各层根据需要设置不同的厚度和材料性质,各层材料之间满足平截面假定,且考虑钢筋混凝土面内弯剪共同作用效应和面外弯曲效应.

为了准确全面了解BSP梁的受力性能和变形特性,本文采用OpenSees基于壳单元对BSP梁进行了非线性有限元分析,目的在于:模拟加载条件,采用合理的单元模型和本构关系对BSP梁进行数值模拟,并与试验结果进行对比分析,验证数值模型的准确性;对影响加固梁受剪性能的主要因素进行参数分析,研究其工作机理,得出影响其受剪性能的相关因素及其作用规律,提出BSP法加固的设计建议.

1 试验概况

试验模型[20]为5根钢筋混凝土梁,其中1根为未加固对比试件(Ctrl)另4根采用BSP法进行加固,加固方案见表 1, 其中Ctrl为参照组.试件截面尺寸、配筋情况、钢板及螺栓的布置方案、加载点布置见图 1.试件均采用4点加载,加载时利用千斤顶和分配梁对各点施加集中荷载(F).所有试件均按“强弯弱剪”设计,其混凝土强度等级为C60,纵向钢筋为Ⅲ级热轧钢筋(C),梁侧钢板用高强螺栓与混凝土连接,螺栓直径(D)为12mm,锚固长度取100mm,螺栓水平对称布置.

下载CSV 表 1 梁的加固方案 Tab.1 Strengthening details of bolted-side- plated specimens
图 1 BSP梁构造(单位:mm) Fig.1 Details of BSP beam (Unit: mm)
2 有限元分析 2.1 混凝土

本文采用弹塑性有限元分析软件OpenSees对BSP梁进行有限元模拟.混凝土采用分层壳单元模拟,该单元采用的多轴受力混凝土模型基于损伤力学和弥散裂缝模型来模拟混凝土的二维受力行为.弥散裂缝是将开裂的混凝土单元仍然当做一个连续的单元考虑[22],不考虑开裂处的单元离散化,认为裂缝弥散于混凝土单元内.

2.2 钢筋

在分层壳单元中,可以通过“弥散”钢筋模型来考虑钢筋布置.但是由于剪力作用下BSP梁中钢筋受力复杂,为了提供模型准确性,采用“离散”钢筋模型,使用Truss单元来模拟.该单元只考虑轴向力和轴向变形,忽略截面方向几何变形的影响,并赋予由Mengotto与Pinto提出的steel02本构模型[24].

2.3 钢板

由于BSP梁中钢板主要在平面内受力,因此采用平板壳单元进行模拟.该单元是一种高性能四边形平板壳单元,由平面膜单元和板弯曲元耦合而成,能够模拟壳体的复杂受力状态[24].钢材作为各向同性材料且破坏准则符合Von Mises弹塑性模型,因此可以使用OpenSees预设的基于偏二张量的J2 Plasticity Material本构关系来模拟.

2.4 锚栓连接

在实际工作中,BSP梁横向剪切变形和纵向弯曲变形会使锚栓受到横向和纵向两个方向的合力,从而产生滑移,而该锚栓滑移会导致BSP梁受剪承载力严重降低,甚至受剪破坏模式的转变,因此在数值分析中必须考虑锚栓滑移效应.本数值模拟中,引用了零长度单元作为非线性弹簧,赋予其Hysteretic单轴本构模型来近似模拟,并利用矢量加法的原理将赋值的单轴本构分解到横向和纵向两个方向上去,从而可考虑纵横向滑移的综合效应,如图 2所示.

图 2 BSP梁单元(一半) Fig.2 Finite element model of BSP beam(symmetry)
2.5 网格划分、边界、加载及分析设置

混凝土及钢板模型网格沿梁高度、长度均按50 mm划分,如图 2所示.

试验中加固梁设置为简支梁,故在建模过程中在支座中心设置约束,一端设置平面内纵向和横向的约束,其余4个方向无约束;另一端仅设置横向约束,其余5个方向无约束.

为了使模型中加载方式尽量与试验加载相同,在荷载定义时需要注意避免出现应力集中,采用多节点等效加载的方式.收敛准则采用了位移控制方式,采用了Krylovnewton算法进行迭代计算.

3 分析结果 3.1 跨中位移-荷载曲线

将上述有限元模型计算所得结果与试验测得的荷载-跨中挠度曲线对比如图 3所示,其中FEM表示模拟,TEST表示实验.试件极限荷载(Pu)及跨中位移模拟与试验的对比列于表 2.

图 3 试验结果与模拟结果对比 Fig.3 Comparison between test results and simulation results
下载CSV 表 2 试验和模拟结果对比 Tab.2 Comparison between test results and simulation results

图 3表 2可看出:试验数据和模拟结果具有良好的一致性,宏观上说明本文建立的有限元模型、选用的材料本构及滑移本构等可以较好地模拟BSP梁的受剪承载力和变形.

3.2 钢筋应变对比

图 4为数值模拟所得的试件P2B1钢筋应变曲线与试验结果的对比.由于混凝土离散性和数值模型使用的弥散裂缝模型无法准确模拟出实际裂缝情况等客观原因,导致模拟试件结果稍大,但相对误差不大,且没有改变钢筋总体受力趋势,表明该模型能较合理地模拟BSP梁的箍筋和纵筋应变jn.

图 4 P2B1梁试验和模拟钢筋应变对比 Fig.4 Strain curves of steel reinforcement in P2B1
3.3 剪跨区钢板应变对比

图 5为数值模拟所得的钢板应变(P2B2)与试验结果比较.其钢板上下边缘应变测点Sn和钢板中部主拉压应变测点hn布置如图 5a所示.由图 5b可知:模拟曲线在加载初期和试验曲线相差极小,当加载到0.7Pu左右时,和试验曲线存在一定误差,主要由于平板壳单元无法模拟钢板平面外屈曲.由图 5c可知,钢板中部主压应变在加载初期呈线性增长,模拟曲线与试验曲线基本吻合;当加载到0.4Pu左右,模拟结果与试验结果开始出现离散,模拟结果增长变缓,但是直到加载至Pu,其相对误差仍然很小.可以看出,数值模拟可很好地模拟BSP梁钢板应变变化.

图 5 P2B2梁钢板应变曲线对比 Fig.5 Strain curves of steel plate in P2B2
4 参数分析

为了弥补试验中试件数量有限的不足,在验证有限元模型的准确性之后,采用有限元分析对影响BSP梁受剪性能的主要参数进行分析.

4.1 锚栓行数

图 6变化的参数是锚栓行数, 其中P3B1-nn表示行数.从图 6a中可以看出,在加载初期,曲线斜率没有明显变化,即增加锚栓行数对BSP梁的初始刚度无影响;加载后期,两行锚栓试件比单行锚栓试件有明显的上升,主要由于两行锚栓可以在钢板上传递力偶,但是多行锚栓曲线与两行曲线基本吻合,即加固设计中仅通过增加锚栓行数无法增加试件受剪承载力.

图 6 不同锚栓行数变形对比 Fig.6 Different row number of anchor bolt

图 6b为垂直于梁轴线的横向滑移,在集中力作用处出现了正向和负向滑移的最大值.图 6c为沿梁轴线的纵向滑移,从图中可看出锚栓行数越多,相对纵横向滑移越小.可知增加锚栓行数可有效减小相对滑移,提高BSP梁承载力和延性.

4.2 纵向锚栓间距

图 7变化参数为纵向锚栓间距,其中P3Bnn表示锚栓间距.如图 7a所示,四组荷载-位移曲线均表现出明显的弹性阶段和塑性阶段,锚栓间距越小,弹性阶段刚度越大;且BSP梁最大位移和极限荷载也随锚栓间距的减小而增大.对比图 7b7c发现,锚栓间距越小,相对滑移越小,BSP梁承载力越高,延性越好.

图 7 不同纵向锚栓间距 Fig.7 Different anchor bolt spacing along the beam
4.3 钢板高度

图 8变化的参数是钢板高度,其中PnB2中n表示钢板高度为n00.从图 8a中看出,当钢板高度小于300 mm时,荷载-位移曲线仅有线性增长阶段,此时BSP梁表现出明显的剪切脆性破坏;当钢板高度达到300 mm后,荷载-位移曲线表现出明显的弹性阶段和塑性阶段,钢板高度越大,弹性阶段刚度越大;且BSP梁最大位移也随钢板高度增大而增大.从图 8b8c可看出,钢板高度越大,相对滑移越大,BSP梁承载力越高,延性越好.

图 8 不同钢板高度 Fig.8 Different steel plate depth
4.4 钢板厚度

图 9变化的参数为钢板厚度,其中P3B2-tntn表示厚度为n.从图 9a中可看出四组BSP梁的荷载-位移曲线均表现出明显的弹性阶段和塑性阶段,钢板厚度越大,弹性阶段刚度越大;且BSP梁最大位移和极限荷载也随钢板厚度增大而增大.从图 9b9c可知,钢板厚度越大,相对滑移越大,BSP梁受剪承载力越高,延性越好.这与前面讨论的结果不同,这是因为钢板厚度越大,钢板发生相同形变时,钢板承受的荷载越大.

图 9 不同钢板厚度 Fig.9 Different steel plate thickness
4.5 钢板屈服强度

图 10变化的是钢板屈服强度,其中P3B1-SnSn表示屈服强度为n00.从图 10a中可看出,四组荷载-位移曲线均表现出明显的弹性阶段和塑性阶段,钢板屈服强度对BSP梁弹性阶段刚度的影响可以忽略不计.当钢板屈服强度从150 MPa增加到400 MPa,BSP梁的受剪承载力有明显提升.但是当其超过400 MPa后,继续增大钢板屈服强度对BSP法加固梁的极限荷载影响极小,可以忽略不计.这主要由于当钢板屈服强度大于一定数值后,梁在加载过程中钢板应力未达到钢板屈服强度就已屈曲,并未发挥钢板全部的受力性能.

图 10 不同钢板屈服强度 Fig.10 Different yield strengths of steel plates

图 10b10c可知,当钢板屈服强度小于400 MPa时,相对纵横向滑移随钢板强度增加而增加;但当钢板屈服强度大于400 MPa后,其对相对滑移、加固梁受剪承载力及延性的影响不大.

4.6 锚栓直径

图 11变化参数为锚栓直径,其中P3B2-dndn表示直径为n.从图a可以看出,BSP梁荷载-位移曲线均表现出明显的弹性阶段和塑性阶段,当锚栓直径小于16 mm时,随着直径的增大,BSP梁的承载能力和最大位移稍有提高;当锚栓直径增大至16 mm后继续增加,曲线与16 mm时基本吻合.从图 11b11c可以看出,锚栓直径越大,钢板和混凝土的相对滑移越小,钢板与混凝土梁连接越紧密.

图 11 不同锚栓直径 Fig.11 Different anchor bolt diameters
4.7 剪跨比

图 12变化参数为BSP梁的剪跨比,其中P3B2-nn表示剪跨比为n.从图a可知,五组荷载-位移曲线均表现出明显的弹性阶段和塑性阶段.随着剪跨比增加,BSP梁的受剪承载力和刚度均减少.从图 12b12c可看出,横向滑移在弯剪区的零点位置、最大值位置与剪跨比有关,即与荷载作用位置有关,具体位置和随荷载作用点位置呈正相关;纵向荷载在弯剪区的变化率也与剪跨比有关,剪跨比越小,变化率越大.

图 12 不同剪跨比 Fig.12 Different shear span ratios
4.8 混凝土强度等级

图 13变化参数为混凝土强度等级,其中P3B2-nn表示混凝土强度等级为n.从图 13a中可以看出,随着混凝土强度等级增加,BSP梁的初始刚度、最大位移和极限荷载有明显的增加;当混凝土强度低于40 MPa时,荷载位移曲线基本呈直线,呈现明显的脆性破坏;当强度增加至40 MPa后,荷载-位移曲线开始出现塑性阶段,BSP梁的延性和变形能力有显著提高.从图 13b13c中可知,混凝土强度越大,钢板与混凝土之间的相对滑移越大.

图 13 不同混凝土强度等级 Fig.13 Different concrete strength grades
5 设计建议

根据上述分析结果可知,影响BSP梁受剪性能的主要因素有:锚栓间距、钢板高度、钢板厚度,而锚栓排数、锚栓直径、钢板屈服强度对加固效果有一定影响.因此,在实际工程应用中,可选钢板屈服强度为400 MPa,锚栓直径为16 mm,且钢板高度不大于300 mm时,可只设置2排螺栓,并在允许范围内减少锚栓间距,增大钢板高度、厚度,以获得较大的受剪承载力、刚度和较好的延性.

6 结论

采用OpenSees软件对BSP梁进行了非线性有限元分析.通过对模拟结果的分析,可得出下列结论:

(1) 有限元模拟结果与试验的宏观、微观数据均保持较高的一致性,说明OpenSees用于模拟BSP梁受剪性能试验具有较高的准确性和可靠性.

(2) 当加强连接件,即增大锚栓排数和直径,减少锚栓间距时,钢板与混凝土间相对滑移减少,增加了BSP梁的整体性,从而提高了受剪承载力和延性;当加强混凝土梁或钢板,即增加钢板高度、厚度和屈服强度或混凝土强度时,BSP梁的相对滑移变形显著增加,钢板与混凝土梁的整体协作能力减弱,但是钢板或混凝土的增强增加了BSP梁的变形能力,因而其受剪承载力有显著提高.

参考文献
[1]
周建庭, 李志刚, 刘思孟. 二次受力RC梁增大截面加固配筋限值研究[J]. 重庆交通大学学报(自然科学版), 2015, 34(1): 30
ZHOU Jianting, LI Zhigang, LIU Simeng. Reinforcement boundary of steel in the reinforcing layer of RC beam under the secondary load[J]. Journal of Chongqing Jiaotong University (Natural Science), 2015, 34(1): 30
[2]
YU Kequan, LI Lingzhi, YU Jiangtao, et al. Feasibility of using ultra-high ductility cementitious composites for concrete structures withoutsteelrebar[J]. Engineering Structures, 2018, 170(9): 11
[3]
ADHIKARY B B, MUTSUYOSHI H, SANO M. Shear strengthening of reinforced concrete beams using steel plates bonded on beam web: Experiments and analysis[J]. Construction and Building Materials, 2000, 14(5): 237 DOI:10.1016/S0950-0618(00)00023-4
[4]
YU Jiangtao, LIU Keke, LI Lingzhi, et al. A simplified method to predict the fire resistance of RC Beams strengthened with near-surface mounted CFRP[J]. Composite Structures, 2018, 193(6): 1
[5]
YANG Xu, GAO Wanyang, DAI Jianguo, et al. Flexural strengthening of RC beams with CFRP grid-reinforced ECC matrix[J]. Composite Structures, 2018, 189: 9 DOI:10.1016/j.compstruct.2018.01.048
[6]
YU J T, WANG Y C, HU K X, et al. The performance of near-surface mounted CFRP strengthened RC beam in fire[J]. Fire Safety Journal, 2017, 90(C): 86
[7]
YU Jiangtao, XU Yanning, YU Kequan, et al. Preliminary study to enhance ductility of CFRP-strengthened RC beam[J]. Journal of Composites for Construction, 2016, 21(1): 4016070
[8]
LI Lingzhi. New partial interaction models for bolted-side-plated reinforced concrete beams[D].Hong Kong: The University of Hong Kong, 2013.
[9]
BARNES R A, BAGLIN P S, MAYS G C, et al. External steel plate systems for the shear strengthening of reinforced concrete beams[J]. Engineering Structures, 2001, 23(9): 1162 DOI:10.1016/S0141-0296(00)00124-3
[10]
LI Lingzhi, JIANG Changjiu, JIA Liangjiu, et al. Local buckling of bolted steel plates with different stiffener configuration[J]. Engineering Structures, 2016, 119(7): 186
[11]
LI L Z, JIANG C J, SU R K L, et al. A piecewise linear transverse shear transfer model for bolted side-plated beams[J]. Structural Engineering and Mechanics, 2017, 62(4): 443 DOI:10.12989/sem.2017.62.4.443
[12]
孙川东.直剪型锚栓钢板加固rc梁斜截面受剪承载能力研究[D].重庆: 重庆大学, 2014.
SUN Chuandong. Research on shear capacity of reinforced concrete beams strengthened with direct shear bolting plate [D]. Chongqing: Chongqing University, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10611-1014043440.htm
[13]
甘元初, 朱建新, 宋会莲. 锚贴钢板加固钢筋混凝土梁受剪性能的试验研究[J]. 南华大学学报(自然科学版), 2007, 21(2): 93
GAN Yuanchu, ZHU Jianxin, SONG Huilian. Shear-resistance property of RC beam strengthened with rivets anchoring steel pate[J]. Journal of University of South China (Science and Technology), 2007, 21(2): 93 DOI:10.3969/j.issn.1673-0062.2007.02.024
[14]
LI L Z, LO S H, SU R K L. Experimental study of moderately reinforced concrete beams strengthened with bolted-side steelplates[J]. Advances in Structural Engineering, 2013, 16(3): 499 DOI:10.1260/1369-4332.16.3.499
[15]
SU R K L, LI L Z, LO S H. Longitudinal partial interaction in bolted side-plated reinforced concrete beams[J]. Advances in Structural Engineering, 2014, 17(7): 921 DOI:10.1260/1369-4332.17.7.921
[16]
SU R K-L, LI L Z, LO S H. Shear transfer in bolted side-plated reinforced concrete beam[J]. Engineering Structures, 2013, 56(11): 1372
[17]
LO S H, LI L Z, SU R K L. Optimization of partial interaction in bolted side-plated reinforced concrete beams[J]. Computers & Structures, 2014, 131(7): 70
[18]
LI L Z, JIANG C J, SU R K L, et al. Design of bolted side-plated reinforced-concrete beams with partial interaction[J]. Proceedings of the Institution of Civil Engineers—Structures and Buildings, 2016, 169(2): 81 DOI:10.1680/jstbu.15.00037
[19]
JIANG Changjiu, LU Zhoudao, LI Lingzhi. Shear performance of fire-damaged reinforced concrete beams repaired by a bolted side-platingtechnique[J]. Journal of Structural Engineering —ASCE, 2017, 143(5): 4017007 DOI:10.1061/(ASCE)ST.1943-541X.0001726
[20]
LI Lingzhi, CAI Ziwei, LU Zhoudao, et al. Shear performance of bolted side-plated reinforced concrete beams[J]. Engineering Structures, 2017, 144(8): 73
[21]
LI Lingzhi, WU Zhenli, YU Jiangtao, et al. Numerical simulation of the shear capacity of bolted side-plated RC beams[J]. Engineering Structures, 2018, 171: 373 DOI:10.1016/j.engstruct.2018.06.003
[22]
林旭川, 陆新征, 缪志伟, 等. 基于分层壳单元的RC核心筒结构有限元分析和工程应用[J]. 土木工程学报, 2009, 43(3): 49
LIN Xuchuan, LU Xinzheng, MIAO Zhiwei, et al. Finite element analysis and engineering application of RC core-tube structures based on the multi-layer shell elements[J]. China Civil Engineering Journal, 2009, 43(3): 49 DOI:10.3321/j.issn:1000-131X.2009.03.008
[23]
王丽莎, 岑松, 解琳琳, 等. 基于新型大变形平板壳单元的剪力墙模型及其在OpenSees中的应用[J]. 工程力学, 2016, 33(3): 47
WANG Lisha, CEN Song, XIE Linlin, et al. Development of a shear wall model based on a new flat shell element for large deformation simulation and application in Opensees[J]. Engineering Mechanics, 2016, 33(3): 47
[24]
LONG Yuqiu, LONG Zhifei, CEN Song. Advanced finite element method in structural engineering[M]. Beijing: Tsinghua University Press, 2009