Mechanical behavior of concrete-filled steel tubular columns subjected to coupled fire and impact loading
-
摘要: 为研究火灾高温与撞击联合作用下钢管混凝土柱的力学性能,基于ABAQUS建立了高温作用下考虑轴力影响的钢管混凝土柱侧向撞击有限元模型。首先,对高温与撞击联合作用下考虑轴力影响的钢管混凝土柱的破坏模式与受力全过程进行了分析,探讨了高温下钢管混凝土柱的抗撞性能与工作机理;其次,重点研究了受火时间、材料强度、含钢率以及撞击能量对抗撞性能的影响,并给出了相关设计建议。研究结果表明:高温与撞击联合作用下,钢管混凝土柱主要发生受弯破坏;受火15 min后,构件抗撞性能明显降低。轴压力对构件抗撞性能产生不利影响,轴压比从0增加到0.2,受火60 min构件抗撞性能下降了7.8%;混凝土强度对高温下钢管混凝土柱抗撞性能有显著影响,受火90 min后,混凝土强度由30 MPa增加到50 MPa,构件抗撞性能提高约85%;外钢管强度与含钢率对高温下抗撞性能影响不大。Abstract: To investigate the mechanical behavior of concrete-filled steel tubular (CFST) columns under coupled fire and impact loads, a finite element (FE) model was established with ABAQUS software to describe the impact resistance of axial-loaded CFST columns under elevated temperatures. The commonly used ISO 834 standard fire and rigid-body impact were employed to model the fire and impact loads, respectively. In the model, the static implicit and dynamic explicit analysis were coupled by using “Restart” and “Import” commands and the strain-rate effect was taken into account. The numerical model was validated by comparing the impact force time history curves obtained from relevant tests at different temperatures. Based on the validated FE models, the impact responses of axial-loaded CFST columns under coupled fire and impact loads were analyzed. Then, the influences of the fire duration, concrete and steel strength, steel ratio and impact energy on the mechanical behavior were investigated, and some design suggestions were proposed. The platform impact force and the maximum mid-span deflection were employed to quantitatively analyze the impact resistance of the CFST columns. The results show that the CFST column presents a flexural failure mode when it is exposed to coupled fire and impact loads. With the increase of fire duration, the proportion of energy consumption of the steel tube reduces. The impact resistance of the column decreases obviously when it is subjected to fire for a duration of 15 min. Axial compression load has an adverse influence on the impact performance. When the axial-load level increases from 0 to 0.2, the platform value of the impact force reduces by 7.8% at a fire duration of 60 min. The concrete strength has a significant effect on the impact resistance. When the cubic compressive strength of the concrete increases from 30 MPa to 50 MPa, the impact resistance is improved by approximately 85% when the fire lasts for 90 min. The steel ratio and steel strength have marginal influences on the impact resistance of the CFST columns at elevated temperatures.
-
Key words:
- impact-resistance /
- fire /
- concrete-filled steel tube /
- mechanisms analysis /
- axial-load level
-
薄壁金属构件由于其低成本和高吸能性而在汽车结构中作为一种性能良好的缓冲吸能元件得到广泛的应用。为了确保碰撞过程中乘员和汽车主要部件的安全, 薄壁构件作为受撞时主要的承载和吸能元件[1], 在耐撞性设计中, 薄壁构件应能够以稳定可控的变形方式吸收尽可能多的能量; 同时最大限度的降低碰撞过程中薄壁构件所承受的压溃载荷, 较高的压溃载荷会使汽车在碰撞过程中产生较大的加速度, 使车内乘员受到严重的伤害, 故在碰撞过程最大压溃载荷应尽可能的小[2]。
诱导结构是薄壁构件应力集中的地方, 可以控制薄壁构件变形形式和降低碰撞时的压溃载荷, 从而使乘员免受高压溃力, 同时也降低了其它部件所承受的压溃载荷。Y.B.Cho等[3]采用数值方法对具有孔洞型诱导结构的前梁装置进行了抗撞性模拟研究, 计算表明设定合理的诱导结构能提高能量吸收同时降低压溃载荷; 张涛等[4]研究了薄壁组合结构及其设计缺陷结构的吸能特性, 提出设置一定的诱导缺陷可以降低碰撞过程中的峰值载荷。
本文中提出在薄壁圆管上增加圆弧形凹槽的诱导变形原则, 研究其对抗撞性的影响。在碰撞分析中大多采用最大峰值压溃力为评价指标, 为了提高结构的抗撞性, 同时满足轻量化的设计要求, 综合考虑比吸能、最大峰值压溃力及压溃力的平均值, 建立多目标薄壁构件优化任务, 并结合径向基函数法构造优化目标的近似函数, 采用理想点法求解多目标优化问题, 详细分析凹槽数量与其半径对薄壁构件抗撞性的影响, 以期为吸能构件抗撞性设计提供依据。
1. 问题描述
(1) 最大峰值压溃力Fmax:在压溃过程中沿轴向产生的压溃载荷的最大值[5]。最大峰值压溃力会引起较高的加速度, 是碰撞过程中需要考虑的重要参数。为了减小对乘员的伤害, 故希望在整个压溃过程中最大峰值压溃力尽可能的小。
(2) 平均压溃力Fm:薄壁构件在碰撞过程中的平均压溃力为:
Fm=∫δmax0F dδδmax (1) 式中:F为压溃力, δ为压溃位移, δmax最大压溃位移。在碰撞过程中压溃力一般围绕平均压溃力波动。
(3) 压溃力效率η, 即平均压溃力与最大压溃力之比, 其表达式为:
η=FmFmax (2) 压溃力效率η综合考虑了平均压溃力与最大压溃力对薄壁构件抗撞性的影响, η越大说明最大压溃力越小越接近平均压溃力, 压溃力曲线平稳, 理想的吸能构件η应为1。
(4) 比吸能χ, 单位质量的薄壁结构所吸收的能量[6]为:
χ=EW (3) 式中:E为薄壁结构所吸收的总能量, W为薄壁结构的总质量。为了提高结构的抗撞性, 同时满足轻量化的设计要求, 应使薄壁构件的比吸能最大化。
以比吸能和压溃力效率为优化设计指标, 在优化设计过程中, 有2个设计参数即凹槽数量m和凹槽半径r, 它们在一定范围内变化。该多目标优化问题可表示为
{maxη(m,r)maxχ(m,r) s. t. mmin⩽m⩽mmaxrmin⩽r⩽rmax (4) 事实上, 目标函数很难用数学方法精确的表达, 本文中采用径向基函数模型构造其近似表达式, 并采用理想点法寻找多目标问题最优解。
2. 结构抗撞性设计
2.1 有限元分析
建立薄壁金属圆管模型, 如图 1所示, 圆管下端固定, 上端自由并承受质量25kg, 冲击速度v=15.49m/s的刚性体冲击作用, 圆管厚度t=0.8mm, 长度L=90mm, 直径d=31mm。
薄壁构件的材料选为高强度钢, 密度ρ=7.82×103 kg/m3, 弹性模量E=207.2GPa, 泊松比ν=0.3, 屈服应力σ0=446MPa。高强度钢材料的动态变形受材料应变率的影响较大, 需要在材料模型中考虑应变率的影响。一般采用Cowper-Symonds塑性材料模型[7]:
σy=σ0(1+ε/c)1/p (5) 式中:σy为考虑应变率之后的动态屈服应力, σ0为静态屈服应力, ε为应变率, c=40s-1和p=5为该模型中与应变率相关的参数, 静态应力应变变化曲线如图 2所示。
为了验证有限元分析的精确性, 与文献[7]中的实验数据进行比较, 表 1给出了薄壁构件原模型由LS-DYNA计算得到的最大峰值压溃力、总吸能、最大压溃位移及其与文献[7]中对应的实验值。通过对比可以看出有限元结果和实验值之间的相对误差很小, 数值模拟的结果与实验值能够较好地吻合。
2.2 模型改进
诱导结构是薄壁构件应力集中的地方, 使结构的屈曲变形稳定。为了提高原模型在碰撞过程中的抗撞性, 本文对原模型结构做了改进, 如图 3所示, 沿光滑表面的薄壁结构(原模型结构尺寸不变)添加均布的圆弧形凹槽。现分别施加1~4个半径为0.8mm诱导槽, 在前文的工况下对4个模型进行实验及有限元数值模拟计算得到其碰撞响应, 如表 2所示。由表 2可知, 改进后的模型与原模型相比η值增大, 最大峰值压溃力显著降低了, 平均压溃力却没有明显减少, 压溃力曲线趋于平稳; 而吸能几乎没有受到影响, 由于增加凹槽导致结构质量的增加, 因而比吸能χ有所减小, 可见添加圆弧形凹槽, 显著提高了薄壁结构的抗撞性。
表 2 改进模型的有限元分析结果及实验值对比Table 2. Comparison of the improved model between experiment and simulationm Fmax/kN Fm/kN η E/kJ χ/(kJ·kg-1) 实验 有限元 实验 有限元 0 120.30 120.09 50.36 0.420 0 2.95 2.94 55.60 1 80.26 77.82 46.14 0.592 9 2.93 2.91 54.29 2 91.63 88.39 42.48 0.480 6 2.87 2.90 53.90 3 94.32 89.26 42.32 0.474 1 2.90 2.93 54.26 4 90.53 87.12 40.92 0.469 7 2.91 2.94 54.24 为了进一步研究凹槽结构对抗撞性的影响, 取凹槽数m=1~4, 半径r=0.8~3.2mm。采用全因子实验设计, 在m、r的设计域内, 选用了28个样本点研究其抗撞性, 并对凹槽结构进行优化设计。表 2所示为原模型和4个改进模型的有限元分析结果与实验值对比, 表 2中m=0表示未加凹槽, 为原模型, m=1~4表示加入0.8mm凹槽的数量, 也即是4个改进模型, 从数据对比看有限元分析结果与实验值较接近, 从而验证了改进模型有限元分析的可靠性。为了节约资源, 对余下的24个设计样本点只进行有限元分析计算从而得到其碰撞响应。
该多目标优化问题数学表达式可进一步表示为:
{min−η(m,r)min−χ(m,r) s. t. 1⩽m⩽40.8 mm⩽r⩽3.2 mm (6) 3. 径向基函数法优化设计
3.1 径向基函数模型
径向函数[8-11]是以预测点与样本点间的欧氏距离为自变量的一族函数。以径向函数为基函数, 通过线性组合构造出来的近似模型即为径向基函数模型。
在设计空间中, 关于设计变量响应函数的近似表达式可以定义为:
˜y(x)=n∑i=1λiφ(‖x−xi‖) (7) 式中:λT=(λ1, λ2, …, λn)为权系数, ‖x-xi‖为预测点和样本点间的欧氏距离, φ(‖x-xi‖)为径向函数。将n个样本点及每个样本点对应的响应函数值代入到式(7), 可以得到方程组矩阵形式为:
Y=Φλ (8) 式中: Y=[˜y(x1),˜y(x2),⋯,˜y(xn)]T,Φ=[Φij]=[φ(‖xi−xj‖)](i,j=1,2,⋯,n)
在得到响应函数的径向基函数模型后, 需要对近似模型精度进行验证, 通过额外测试点(不包括样本点)决定系数R2和误差平方和方均根 √¯R2来验证径向基函数模型对设计变量的拟合情况。决定系数R2和误差平方和方均根 √¯R2分别定义为:
R2=1−n∑i=1(yi−˜yi)2n∑i=1(yi−ˉyi)2 (9) √¯R2=√n∑i=1(yi−˜yi)2n (10) 式中:n是测试点的个数, 其中yi是i个测试点的有限元分析结果yi的平均值, ˜yi是第i个测试点的径向基函数近似值。
3.2 理想点法求解多目标优化问题
多目标优化问题在数学上一般可表示为
{minf1(x1,⋯,xn)⋮minfp(x1,⋯,xn) s. t. gi(x)⩾0,i=1,2,⋯,mhi(x)=0,i=1,2,⋯,l (11) 先求解p个单目标问题
minx∈Dfj(x),j=1,2,⋯,p (12) 设其最优值为fj*, 称f*=(f1*, …, fp*)T为值域中的一个理想点, 因为一般很难达到, 故在期望的某种度量下, 寻求距离f*最近的f作为近似值。最直接的方法是构造评价函数并极小化, 既求解:
minx∈Dφ[f(x)]=√p∑i=1[fi(x)−f⋆i]2 (13) 并将它的最优解x*作为式(8)在这种意义下的最优解。
于是多目标优化问题可进一步表示为
{minf1(x)=−χ(x)minf2(x)=−η(x) s. t. 1⩽m⩽40.8 mm⩽r⩽3.2 mm (14) 3.3 优化结果
对于带有圆弧形凹槽的薄壁构件, 凹槽数量m和半径r是影响结构抗撞性能的2个重要参数, 通过有限元计算得到了28个样本点的碰撞响应, 从而得到η和χ关于设计变量m、r的响应函数和响应面, 如图 4~5所示。从图中可见凹槽数量半径对结构的比吸能和压溃力效率有很大影响, 呈现出一定的非线性关系, 不具备单调性。
表 3给出了χ和η这2个单目标函数的优化值, 可见2个单目标函数最优解不同, 2个目标函数很难同时达到最优。采用理想点法对式(14)的多目标优化问题进行求解, 得到带有圆弧形凹槽薄壁构件的最优值为m=3和r=1.224mm, χ=53.723kJ/kg, η=0.673。图 6分别给出了采用有限元方法计算得到的最大的比吸能和压溃力效率的结构, 以及采用理想点法计算得到多目标问题的优化结构。
表 3 单目标函数优化结果Table 3. Optimums of single objective functions单目标函数 m r/mm χ/(kJ·kg-1) η max χ 1 0.800 54.29 0.592 9 max η 4 1.198 51.94 0.682 5 优化后构件的η=0.673, 比表 2中m=0的原模型的η=0.420 0提高了52.90%, 而χ=53.723kJ/kg, 与55.06kJ/kg相比仅降低了4.05%。图 7~8分别给出了该多目标问题达到最优时, 在碰撞过程中压溃力和吸收能量随压溃位移和压溃时间变化情况, 并将其与相应条件的原模型进行比较。从压溃载荷图 7可以看出, 与原模型相比优化后构件最大峰值压溃力降低了52.1%, 最大峰值压溃力更接近平均压溃力, 压溃力曲线平稳, 可见以η为优化设计指标是合理的, η综合考虑了最大峰值压溃力和压溃力曲线的平稳性。图 8所示为总吸能与压溃时间的关系, 可以看出, 优化后构件所吸收的能量并没有很大的减少, 在碰撞时间为10ms时, 优化后构件所吸收的能量与同时间原模型吸收的能量比较相近了, 可见优化后构件的抗撞性得到了显著提高。
优化后构件在碰撞过程中发生叠缩变形如图 9所示, 这种变形模式有利于结构在变形过程中吸收更多的碰撞动能。带有圆弧形诱导凹槽结构的薄壁构件变形首先在凹槽处开始, 最大的塑性弯矩和塑性铰发生在凹槽处, 薄壁结构在该处容易压溃。凹槽能控制塑性变形, 塑性应变在该处达到最大值, 这也就意味着沿着变形的凹槽处能消散大量的压溃力, 大部分的冲击能在该处以塑性变形的方式被吸收。可见增加圆弧形凹槽诱导结构后, 可大大提高薄壁构件的抗撞性能。
4. 结论
对薄壁圆管结构进行改进引入圆弧形诱导凹槽, 以比吸能和压溃力效率为评价指标, 建立了兼顾最大峰值压溃力、平均压溃力及比吸能的多目标优化问题的数学模型, 研究凹槽数量和半径对抗撞性能的影响。数值分析的结果表明:引入圆弧形诱导凹槽结构在保证不减少薄壁构件吸能能力的情况下达到了减小最大峰值压溃力、使压溃力曲线平稳的目的。优化结果表明最优构件的最大峰值压溃力比原模型减小了52.1%, 压溃力曲线与原模型相比更趋于平稳, 从而有效的保证了碰撞过程中乘员的安全性; 从吸能角度来看, 最优构件所吸能的能量并没有很大减少, 在碰撞时间为10ms时, 最优构件所吸收的能量达到了与同一时间原模型吸收的能量; 而最优构件的比吸能只比原模型减少了4.05%, 可见引入圆弧形诱导凹槽提高了薄壁结构的抗撞性, 为进一步研究吸能构件的抗撞性设计提供了依据。
-
表 1 构件详细参数
Table 1. Detailed parameters of specimens
组别 构件编号 D0×ts/mm n t0/min m0/kg v0/(m·s−1) fcu/MPa fy/MPa 不加轴力 W0 400×8 0 0 2000 7 40 345 W3 400×8 0 30 2000 7 40 345 W6 400×8 0 60 2000 7 40 345 W9 400×8 0 90 2000 7 40 345 加轴力 Z0 400×8 0.2 0 2000 7 40 345 Z3 400×8 0.2 30 2000 7 40 345 Z6 400×8 0.2 60 2000 7 40 345 -
[1] HUO J S, ZHENG Q, CHEN B, et al. Tests on impact behaviour of micro-concrete-filled steel tubes at elevated temperatures up to 400℃ [J]. Materials and Structures, 2009, 42(10): 1325–1334. DOI: 10.1617/s11527-008-9452-0. [2] 何远明. 高温下钢管混凝土SHPB抗冲击性能试验研究和理论分析[D]. 长沙: 湖南大学, 2011.HE Y M. Experimental and numerical study on SHPB impact behaviors of comcrete-filled steel tube at elevate temperatures [D]. Changsha: Hunan University, 2011. [3] 霍静思, 任晓虎, 肖岩. 标准火灾作用下钢管混凝土短柱落锤动态冲击试验研究 [J]. 土木工程学报, 2012, 45(4): 9–20. DOI: 10.15951/j.tmgcxb.2012.04.009.HUO J S, REN X H, XIAO Y. Impact behavior of concrete-filled steel tubular stub columns under ISO-834 standard fire [J]. China Civil Engineering Journal, 2012, 45(4): 9–20. DOI: 10.15951/j.tmgcxb.2012.04.009. [4] JIN L, ZHANG R B, LI L, et al. Impact behavior of SFRC beams at elevated temperatures: experimental and analytical studies [J]. Engineering Structures, 2019, 197: 109401. DOI: 10.1016/j.engstruct.2019.109401. [5] 丛珊. 高温下钢框架抗冲击性能及梁柱连接性能的实验和数值模拟分析[D]. 济南: 山东建筑大学, 2018.CONG S. Experimental and numerical simulation analysis on impact resistance and beam-column connection performance of steel frames at elevated temperatures [D]. Jinan: Shandong Jianzhu University, 2018. [6] 史艳莉, 纪孙航, 王文达, 等. 高温作用下钢管混凝土构件侧向撞击性能 [J]. 爆炸与冲击, 2020, 40(4): 043303. DOI: 10.11883/bzycj-2019-0293.SHI Y L, JI S H, WANG W D, et al. The lateral impact performance of concrete-filled steel tubular (CFST) members at high temperatures [J]. Explosion and Shock Waves, 2020, 40(4): 043303. DOI: 10.11883/bzycj-2019-0293. [7] International Organization for Standardization. Fire-resistance tests: elements of building construction: Part 1: general requirements: ISO 834-1—1999 [S]. 1999. [8] 王蕊, 李珠, 任够平, 等. 钢管混凝土梁在侧向冲击荷载作用下动力响应的试验研究和数值模拟 [J]. 土木工程学报, 2007(10): 34–40. DOI: 10.15951/j.tmgcxb.2007.10.010.WANG R, LI Z, REN G P, et al. Experimental study and numerical simulation of the dynamic response of concrete filled steel tubes under lateral impact load [J]. China Civil Engineering Journal, 2007(10): 34–40. DOI: 10.15951/j.tmgcxb.2007.10.010. [9] 姜珊, 王蕊. 中空夹层不锈钢钢管混凝土构件的侧向撞击试验及有限元分析 [J]. 工业建筑, 2016, 46(11): 161–167. DOI: 10.13204/j.gyjz201611031.JIANG S, WANG R. Experiment study and finite element analysis of concrete filled stainless and steel double skin tubes member under lateral impact [J]. Industrial Construction, 2016, 46(11): 161–167. DOI: 10.13204/j.gyjz201611031. [10] 任够平. 预加轴力的钢管混凝土构件受刚性块侧向冲击的动力响应分析[D]. 太原: 太原理工大学, 2009.REN G P. On dynamic responses of axially preloaded concrete filled steel tubular members under lateral impact by a rigid body [D]. Taiyuan: Taiyuan University of Technology, 2009. [11] WANG R, HANG L H, ZHAO X L, et al. Experimental behavior of concrete filled double steel tubular (CFDST) members under low velocity drop weight impact [J]. Thin-Walled Structures, 2015, 97: 279–295. DOI: 10.1016/j.tws.2015.09.009. [12] ZHAO H, WANG R, HOU C C, et al. Performance of circular CFDST members with external stainless-steel tube under transverse impact loading [J]. Thin-Walled Structures, 2019, 145: 10638. DOI: 10.1016/j.tws.2019.106380. [13] LIE T T, KODUR V K R. Fire resistance of steel columns filled with bar-reinforced concrete [J]. Journal of Structural Engineering, 1996, 122(1): 30–36. DOI: 10. 1061/(ASCE)0733-9445(1996)122: 1(30). [14] HONG S, VARMA A H. Analytical modeling of the standard fire behavior of loaded CFT columns [J]. Journal of Applied Physics, 2009, 65(1): 54–69. DOI: 10.1016/j.jcsr.2008.04.008. [15] British Standard Institution. Design of steel structures: part 1-2: general rules: structural fire design: EN 1993-1-2—2005[S]. 2005. [16] 刘发起. 三面受火的矩形钢管混凝土柱抗火性能分析[D]. 哈尔滨: 哈尔滨工业大学, 2010.LIU F Q. Fire resistance of concrete filled RHS columns under three-surface fire loading [D]. Harbin: Harbin Institute of Technology, 2010. [17] 韩林海. 钢管混凝土结构理论与实践[M]. 北京: 科学出版社, 2016.HAN L H. Concrete filled steel tubular structures theory and practice [M]. Beijing: Science Press, 2016. [18] DING J, WANG Y C. Realistic modelling of thermal and structural behaviour of unprotected concrete filled tubular columns in fire [J]. Journal of Constructional Steel Research, 2008, 64(10): 1086–1092. DOI: 10.1016/j.jcsr.2007.09.014. [19] 侯川川. 低速横向冲击荷载下圆钢管混凝土构件的力学性能研究[D]. 北京: 清华大学, 2012.HOU C C. Study on performance of circular concrete-filled steel tubular (CFST) members under low velocity transverse impact [D]. Beijing: Tsinghua University, 2012. [20] CHEN H, LIEW J Y. Explosion and fire analysis of steel frames using mixed element approach [J]. Journal of Engineering Mechanics, 2005, 131(6): 606–616. DOI: 10.1061/(ASCE)0733-9399(2005)131:6(606). [21] CHEN L, FANG Q, JIANG X Q, et al. Combined effects of high temperature and high strain rate on normal weight concrete [J]. International Journal of Impact Engineering, 2015, 86: 4056. DOI: 10.1016/j.ijimpeng.2015.07.002. [22] ZHAO H, WANG R, LI Q M, et al. Experimental and numerical investigation on impact and post-impact behaviours of H-shaped steel members [J]. Engineering Structures, 2020, 216: 110750. DOI: 10.1016/j.engstruct.2020.110750. [23] WANG Y, QIAN X D, LIEW J Y, et al. Experimental behavior of cement filled pipe-in-pipe composite structures under transverse impact [J]. International Journal of Impact Engineering, 2014, 72: 1–16. DOI: 10.1016/j.ijimpeng.2014.05.004. 期刊类型引用(7)
1. 雷成,李福胜,缪得祥,朱涛. 基于回归分析法的矩形槽圆管吸能特性研究. 机械强度. 2024(02): 292-300 . 百度学术
2. 王春华,杨雨泽,安达,唐治. 液压支架凹角圆管式薄壁构件吸能特性分析与优化. 机械设计与制造. 2024(10): 266-271+276 . 百度学术
3. 邓敏杰,刘志芳. 仿马尾草薄壁结构的设计与耐撞性研究. 高压物理学报. 2022(03): 111-120 . 百度学术
4. 张欣玥,惠旭龙,葛宇静,舒挽,白春玉,刘小川. 中低速压缩加载下不同截面构型复合材料薄壁结构吸能特性及失效分析. 爆炸与冲击. 2022(06): 36-49 . 本站查看
5. 霍鹏,许述财,范晓文,李建平,杨欣,黄晗. 鹿角骨单位仿生薄壁管斜向冲击耐撞性研究. 爆炸与冲击. 2020(11): 127-138 . 本站查看
6. 谭丽辉,谭洪武,鲁帅. 诱导缺陷结构在薄壁构件耐撞性中的应用. 河南科技. 2016(21): 74-75 . 百度学术
7. 李玉如,柳忠彬,肖守讷,王欢. 预制螺旋槽薄壁管的吸能研究. 四川理工学院学报(自然科学版). 2015(02): 21-25 . 百度学术
其他类型引用(16)
-