Bearing Capacity of Concrete-Filled Double Skin Stub Columns with Square outer Stainless Steel tube Under Axial Compression
-
摘要:
为促进方中空不锈钢管混凝土构件在土木工程中的应用,以不锈钢外管厚度和混凝土强度为变量的6组试件为研究对象,首先,进行轴压试验,得到了不同试件在轴压荷载作用下的破坏模式、荷载-位移曲线、荷载-应变曲线,并进一步分析了不锈钢方管宽厚比、核心混凝土强度以及不锈钢方管约束效应系数对方中空不锈钢管混凝土短柱极限承载力的影响;然后,初步讨论了倒角对强度和延性的影响,提出了避免内管先于外管屈曲的最小厚度计算方法;最后,基于试验结果以及已有文献数据,采用拟合方法推导了方中空不锈钢管混凝土短柱的抗压承载力计算式,并与已有文献的简化模型及国外主要规范的计算结果进行对比. 研究结果表明:试件宽厚比由34.9降至20.9,极限承载力的提升率平均为98.5%,核心混凝土强度由C40提升至C60时,试件极限承载力的提升率平均为7.3%;短柱的轴压极限承载力随约束效应系数近似呈线性增加,约束效应系数
$ \xi $ 越大,短柱的承载力越高;本文得到的计算式可以较好地预测方中空不锈钢管混凝土短柱的轴压承载力.-
关键词:
- 方中空不锈钢管混凝土短柱 /
- 宽厚比 /
- 轴压承载力 /
- 破坏模式 /
- 约束效应
Abstract:In order to promote the application of square concrete-filled double skin tube (SCFDST) columns with outer stainless steel tube in civil engineering, six groups of specimens with different outer tube thicknesses and different core concrete strength were tested under axial compression, and the failure modes, load–displacement curves, and load–strain curves were obtained. The influences of the width-to-thickness ratio of the stainless steel square tube, the strength of core concrete, and the restraint effect coefficient of stainless steel square tube on the ultimate bearing capacity of SCFDST short columns were further analyzed. Meanwhile, the influence of chamfer on the strength and ductility of specimens was discussed preliminarily, and the minimum thickness calculation method to avoid the inner tube buckling before the outer tube buckling was proposed. Finally, based on the test results and the data in the existing literature, a fitted formula for calculating the compressive capacity of SCFDST short columns was derived; the calculated results by the proposed method were compared with those by the simplified model in the literature and the main foreign specifications, to verify its effectiveness. The results show that with the width-to-thickness ratio decreasing from 34.9 to 20.9, the ultimate bearing capacity was increased by 98.5% on average. When the core concrete strength was increased from C40 to C60, the ultimate bearing capacity of the specimen was increased by 7.3% on average. Besides, the axial ultimate bearing capacity of SCFDST short columns increases linearly with the constraint effect coefficient. Compared with the simplified model in the literature, the formula obtained in this study can predict the bearing capacity of SCFDST short columns.
-
近年来,振动半主动控制系统在各工程领域受到广泛关注,其能提供接近主动控制装置的减振效果,且能耗低,同时还能提供如同被动装置的可靠性. 在众多半主动控制系统中,磁流变(MR)阻尼器在应用的简单性、可靠性和鲁棒性方面具有显著优势,应用磁流变阻尼器的半主动控制悬架系统受到车辆振动舒适性研究人员的重视.
对磁流变半主动悬架系统的研究内容主要包括磁流变阻尼器动力学建模[1-2]、悬架系统建模及振动传递特性研究[3]、半主动悬架控制系统仿真分析及试验验证[4-5]等. 磁流变阻尼器作为半主动悬架的执行元件,其力学特性分析及正逆模型建立一直是半主动悬架的研究热点问题. MR阻尼器的阻尼力表现出明显的非线性及滞回特性,且具有显著的场强依赖性及激励特性(振幅及频率)依赖性[6]. 阻尼力具有多值特性,不仅取决于阻尼器活塞速度,同时受活塞运动历史状态及控制电流的影响.
用于描述磁流变阻尼力的模型众多,如:包含微分环节的Bouc-Wen模型、Duhem模型、LuGre模型等;包含代数模型的sigmoid模型、双曲正切模型、魔术公式模型等;由非参数模型表示的多项式模型、神经网络模型等[7-9]. 上述模型各自具有准确性、复杂性等方面的优缺点,目前还没有一种通用的应用于半主动悬架的磁流变阻尼器力学模型.
建立准确的控制电流逆模型是使阻尼器产生期望阻尼力的重要基础,传统的逆模型建立方法包括不基于阻尼力正模型的闭环控制逆模型(如阶梯函数逆模型、Sigmoid函数逆模型、状态连续反馈逆模型等)、基于阻尼力正向模型直接逆推的逆动力学模型(如基于Bingham模型、Bouc-Wen模型、魔术公式模型模型等)[10-12]. 基于正向模型直接逆推逆动力学模型的主要难点在于容易产生难以直接计算解析解的超越方程,不可避免地需要对正向模型进行简化处理,影响建模精度. 上述建模方法无法准确描述不同工况下控制电流与阻尼器运动状态及期望阻尼力的关系,因此,基于神经网络及模糊推理系统的逆模型建模方法受到研究者们的关注.
本文在试验获取磁流变阻尼器力学性能的基础上,对磁流变阻尼器正逆动力学模型进行建模及参数识别,并将建立的正逆模型引入半主动座椅悬架控制系统模型中进行分析,对所提出正逆模型的效果进行验证.
1. 磁流变阻尼器建模及参数识别
1.1 磁流变阻尼器滞回分离模型
建立结构简单、拟合精度高且通用的磁流变阻尼器模型是研究磁流变阻尼器力学特性及半主动悬架控制系统的重点和基础. 基于魔术公式的磁流变阻尼器力学模型[13]与Bouc-Wen模型等包含微分方程的模型相比,计算复杂度明显降低,进而能够有效提高模型的参数辨识效率,如式(1)所示.
{Fd=Fk+Fc+Fm+Fb+f,Fk=k(x−x0),Fc=c˙x,Fm=m¨x,Ff=f0sign˙x,Fb=Dsin[Carctan(B˙x−E(B˙x−arctanB˙x))], (1) 式中:
$ {F}_{{\rm{d}}} $ 、$ {F}_{{\rm{k}}} $ 、$ {F}_{{\rm{c}}} $ 、$ {F}_{{\rm{m}}} $ 、$ {F}_{{\rm{f}}} $ 及$ {F}_{{\rm{b}}} $ 分别为模型输出阻尼力、弹性项、黏性项、惯性项、摩擦项及剪切项,该模型滞回特性主要由$ {F}_{{\rm{k}}} $ 及$ {F}_{{\rm{m}}} $ 表征;B、C、D、E为调整$ {F}_{{\rm{b}}} $ 曲线形状的参数,$ B={b}_{1}I + {b}_{2} $ ,$ D={d}_{1}I + {d}_{2} $ ,I为电流,b1、b2、d1、d2为与电流相关的参数;$ x $ 和$ \dot{x} $ 分别为阻尼器活塞位移及速度;$ k $ 、$ c $ 、$ {f}_{0} $ 分别为弹性项、黏性项、摩擦项系数;$ {x}_{0} $ 为弹性项初始位移;f为频率.在建立模型参数与电流的关系中,魔术公式内部参数B包含与电流相关参数,无法直接逆向推导出控制电流与期望力的关系式,对控制电流的计算精度产生影响.
本文在魔术公式模型的基础上,参考于建强等[14-15]建立的基于骨架曲线与滞回分离的建模方法,建立基于正弦、余弦型魔术公式的磁流变阻尼器滞回分离模型,如式(2). 滞回分离模型建立的原理如图1所示.
{Fd=FmaxFnorm +f0,Fnorm =Fpf+Fhyssignx+c˙x,Fpf=Dpsin[Cparctan(Bp˙x−Ep(Bp˙x−arctanBp˙x))],Fhys=Dhcos[Charctan(Bh˙x−Eh(Bh˙x−arctanBh˙x))], (2) 式中:
$ {F}_{{\rm{max}}} $ 和$ {F}_{{\rm{norm}}} $ 分别为最大阻尼力及归一化后的阻尼力;$ {F}_{{\rm{pf}}} $ 和$ {F}_{{\rm{hys}}} $ 分别对应骨架曲线与滞回分离曲线;Bp、Cp、Dp、Ep及Bh、Ch、Dh、Eh分别为影响Fpf及Fhys曲线形状的参数.1.2 滞回模型参数辨识
本文采用德国SCHENCK公司液压作动器及Lord RD-8040-1磁流变阻尼器(图2),试验工况包括振幅5~20 mm、频率0.5~9.0 Hz的正弦位移激励,阻尼器输入电流为0~1 A.
磁流变减振器力学模型参数识别问题本质上是寻找模型计算结果和减振器实际力学特性误差最小值的数学优化问题,如式(3)所示.
Fobj=min{RMS(Fcalcu−Ftest)}, (3) 式中:
$ {F}_{{\rm{obj}}} $ 、$ {F}_{{\rm{test}}} $ 和$ {F}_{{\rm{calcu}}} $ 分别为目标函数、试验值及模型计算值.为考虑磁流变阻尼力动态力不对称性,通过处理试验数据,得到不同工况下滞回环速度中心,其随电流及激励的变化规律如图3所示. 由图3可知:滞回环速度中心随施加电流增加而逐渐增加;随激励幅值及频率增加均呈现出增加趋势,当频率大于2 Hz时,该趋势更加明显.
因此,可建立滞回环速度中心
$ {\dot{x}}_{{\rm{center}}} $ 与电流及激励特性的关系,如式(4)所示.˙xcenter=(v1I+v2)(v3lnvm+v4), (4) 式中:
$ {v}_{{\rm{m}}}\mathrm{为} $ 最大速度,对于正弦激励信号,${v}_{{\rm{m}}}= 2{\text{π}} Af=\sqrt{{\dot{x}}^{2}-x\ddot{x}}$ ,$\ddot{x} $ 为阻尼器加速度,A 为激励幅值;$ {v}_{1} $ 、$ {v}_{2} $ 、$ {v}_{3} $ 及$ {v}_{4} $ 为拟合公式系数.将阻尼力非对称特性考虑进模型中,式(2)中所有速度项
$ \dot{x} $ 变为$ \dot{x}-{\dot{x}}_{{\rm{center}}} $ .图4为不同激励及电流条件下最大阻尼力变化曲线,本文试验所得阻尼力随施加电流增加而逐渐增加,并未达到完全磁饱和状态. 在相同电流条件下,最大阻尼力均随阻尼器最大活塞速度增加而增加,利用式(5)形式函数可以较准确地描述最大阻尼力和电流及激励特性的关系.
Fmax(I,vm)=(k1I2+k2I+k3)(k4vk5m), (5) 式中:
$ {k}_{1} $ ~$ {k}_{5} $ 为拟合公式系数.为对考虑阻尼力非对称性及最大阻尼力与激励特性及电流变化关系后的模型参数进行辨识,本文采用一种新型的混合优化算法——基于Sobol 序列的差分-禁忌混合算法(SS-DTHA),该算法基于差分进化算法较强的全局搜索能力,引入禁忌搜索增强算法局部搜索能力,并通过调节算法对全局搜索和局部搜索的侧重比例,进一步改进算法的探索和挖掘能力,对解决非线性系统参数辨识问题具有积极作用,如图5所示.
图6为本文提出的改进魔术公式模型与tanh模型及双Sigmod模型计算阻尼力的辨识结果对比(激励工况为15 mm-4.0 Hz-1.0 A). 由图可知:本文提出的模型在滞回环过渡区及屈服后区域均具有更好的拟合精度;3种模型与试验数据的误差均方根值分别为41.07、96.83、90.26 N.
图7为不同输入电流及激励频率条件下,模型计算阻尼力的辨识结果,并与试验测试值进行对比,图中实线和虚线分别代表试验值及模型计算值. 由图可知,在不同激励频率及电流条件下,本文提出的阻尼力模型均能准确地表征磁流变阻尼器动态阻尼力的非线性及滞回特性.
表1为不同激励幅值、频率及输入电流条件下平均百分比误差
$ {E}_{{\rm{average}}} $ 的计算结果,如式(6). 可以看出,本文提出的模型在不同工况下均能得到较低的误差,模型精度较高. 各工况平均百分比误差在3.4%附近,误差变化范围为2.7%~4.8%. 模型对试验阻尼力拟合效果较好,满足半主动控制系统建模要求.表 1 不同工况模型平均百分比误差Table 1. Comparison of average percentage error of models under different working conditions电流/A 幅值 5 mm 幅值 15 mm 3 Hz 4 Hz 3 Hz 4 Hz 0 0.0447 0.0355 0.0457 0.0452 0.2 0.0474 0.0355 0.0289 0.0396 0.4 0.0316 0.0321 0.0302 0.0319 0.6 0.0310 0.0322 0.0269 0.0284 0.8 0.0344 0.0349 0.0305 0.0315 1.0 0.0308 0.0359 0.0282 0.0273 Eaverage=RMS(Fcalcu−Ftest)/RMS(Ftest). (6) 为建立适用于阻尼器工作范围的通用模型,通过分析模型各参数的识别结果可以发现:参数Bp、Bh、c随电流及激励特性变化具有明显的规律,其他参数与电流、激励无明显变化关系,设置为常数. 参数Bp、Bh、c随电流及激励特性变化关系如图8所示. 由图可知:参数Bp、Bh随激励频率增加逐渐降低,在相同激励工况下,电流变化对参数Bp、Bh影响不大;参数c随电流变化无明显增降规律,随频率增加而逐渐降低.
上述参数与激励和电流的关系拟合为
{Bp=p1vp2m,Bh=h1vh2m,c=c1vc2m, (7) 式中:
$ {p}_{1} $ 、$ {p}_{2} $ 、$ {h}_{1} $ 、$ {h}_{2} $ 、$ {c}_{1} $ 及$ {c}_{2} $ 分别为各参数的拟合系数.基于以上参数的变化规律,对幅值为5 mm正弦激励下1.0~9.0 Hz及0~1.0 A电流共47个工况模型参数进行重新辨识,可得到不同工况通用模型的辨识结果,如表2所示.
表 2 滞回模型参数辨识结果Table 2. Parameter identification of hysteretic model变量 辨识结果 ${F}_{{\rm{max}}}$ $\left(-9.36{I}^{2} + 123.97I + 23.05\right)\left(4.29{ {v}^{0.19}_{{\rm{m}}} }\right)$ ${\dot{x} }_{{\rm{center}}}$ $\left(-3.47I + 1.91\right)(-2.35{\rm{ln} }\;{v}_{ {\rm{m} } } + 9.21)$ ${B}_{{\rm{p}}}$ $1.97{v}^{-1.09}_{{\rm{m}}}$ ${B}_{{\rm{h}}}$ $1.01{v}_{{\rm{m}}}^{-0.86}$ $ c $ $0.32{ {v}_{ {\rm{m} } } ^{-1.01}}$ Cp 1.32 Dp 0.70 Ep −7.64 Ch 1.10 Dh 0.42 Eh −11.50 f0 −80.34 图9为所建立的通用模型对10 mm正弦激励下1.0~5.0 Hz、0.2~1.0 A电流不同工况的拟合结果与试验数据对比曲线,图中实线代表试验数据,虚线代表模型计算结果. 本文所建通用模型对不同工况下的阻尼器动态力学特性均具有较好的拟合能力.
2. 基于ANFIS建立阻尼器逆模型
半主动悬架控制系统需要根据座椅运动状态及控制策略计算期望阻尼力,需建立控制电流逆模型计算获得期望阻尼力所需的电流. 本文选用自适应神经模糊系统(ANFIS)建立磁流变逆动力学模型.
为构建准确的ANFIS逆模型,建立准确全面的训练数据集是建模的重点和基础. 本文将阻尼器动态力学特性试验中所有工况案例中均匀抽取数据点整合在起,包含激励幅值5、10、15、20 mm,激励频率0.5~9.0 Hz,输入电流0~1.0 A工况,形成16400个数据点的总集合. 为避免阻尼器运动状态突变对模型的影响,选取每组数据开始及结束速度均接近0,并在每组数据后加10个全零数组,以避免先前数据集合对下一集合的输出预测的影响.
为增加ANFIS模型对变电流工况的推理能力,本文在试验获取的测试数据的基础上,根据磁流变减振器工作范围,建立激励幅值范围为2~14 mm、频率范围为0.5~10.0 Hz的扫频激励. 通过对位移信号进行微分来获得速度信号,时间步长为0.01 s,通过30 s模拟生成总共3000个数据集. 输出电流训练数据由幅值 0~1.0 A 的随机信号产生. 在激励数据不变的基础上,生成5组随机电流数组共15 000个数据点.
本文逆模型训练数据在试验恒定电流数据集及通过正向模型计算得到的随机电流数据集中各抽取7500组数据组成15000个数据点的训练集合,其余数据作为验证数据. 选取
$x\left(t\right)、\dot{x}\left(t\right)、\ddot{x}\left(t\right)、F\left(t\right)、 F\left(t-1\right)$ 作为模型输入,建立包含32 个“if-then”模糊规则的五输入单输出ANFIS模型.图10(a)为训练后逆模型的预测电流与验证目标恒定电流对比,由图可知,预测电流在目标电流值附近波动,两者之间误差均方根值为0.0869 A,电流在0附近误差较大. 图10(b)为随机电流序列模型验证结果,预测电流能够跟随目标电流的变化趋势,两者之间误差均方根值为0.1171 A.
将阻尼器逆向动力学模型与阻尼力模型结合,使该系统根据期望阻尼力产生预测电流进而产生预测阻尼力,模型验证流程如图11所示. 预测阻尼力与目标阻尼力对比结果见图12(a). 结果表明,由 ANFIS 逆模型得到的预测电流以及产生的预测阻尼力能够较好地跟随目标阻尼力. 图12(b)为一组在阻尼力活塞运动速度范围内的验证数据,可以发现,在阻尼器运动范围内,预测阻尼力与目标阻尼力变化范围基本一致. 阻尼器全工况范围内预测阻尼力误差均方根值为135.79 N,为阻尼器最大阻尼力的5.6%.
3. 半主动座椅悬架控制系统仿真
建立包含上述建立的阻尼器正逆模型的单自由度半主动座椅悬架系统,簧上质量、刚度及被动阻尼系数分为75 kg、47.37 kN/m及1200 Ns/m. 控制策略采用on-off控制、连续天棚控制及LMS控制算法进行对比分析,on-off控制根据座椅运动状态在电流为最小值(0)和最大值(1.0 A)之间切换,天棚阻尼系数设置为2000 Ns/m,LMS算法步长因子设置为0.005,除控制算法外的半主动座椅悬架模型剩余部分均设置相同. 激励信号选取幅值5 mm、频率1.0~9.0 Hz正弦扫频激励. 图13为半主动座椅悬架控制系统示意.
图14为被动和半主动座椅悬架位移传递率. 可以发现:被动座椅悬架在共振频段将座椅下地板振动激励放大;与被动座椅悬架相比,3种控制算法在低频共振频段均能够有效地衰减振动. 通过计算,在4.0 Hz处3种算法的座椅悬架位移传递率较原被动座椅分别降低了47.62%、48.72%、42.09%. 天棚控制算法及最小均方算法(least mean square, LMS)在共振频率以上仍能保持较好的振动传递特性,on-off控制对高频激励的振动衰减性能较差.
半主动控制系统的实时性对控制效果具有显著影响. 基于磁流变阻尼器的半主动座椅悬架系统各环节时间延迟如图15所示,包括传感器信号采集时延
$ \Delta {t}_{3} $ 、 控制器时延$ \Delta {t}_{2} $ (数据处理、计算期望阻尼力、输出驱动电流时延)及磁流变阻尼器响应时延$ \Delta {t}_{1} $ .本文用经典控制算法(on-off控制、天棚控制)对磁流变阻尼器响应时滞的鲁棒性进行分析,计算结果如图16所示. 可以发现:对于两种算法,增加阻尼器响应时滞均会降低振动控制的改善效果;对于采用连续控制模型的天棚控制算法,对于阻尼器响应时滞具有一定鲁棒性,时延18 ms时在仿真频域范围内仍可以实现优于被动座椅悬架的减振效果;on-off控制算法对阻尼器响应时滞较敏感,特别是超过
$ \sqrt{2} $ 倍固有频率后的高频段,当阻尼器时延18 ms时,on-off控制座椅悬架振动传递特性变差. 因此,尽可能减小半主动座椅悬架系统各个环节的时滞,对于改善半主动控制系统减振效果具有积极作用.通过上述分析可知,本文建立的阻尼力模型及控制电流逆模型能够产生满足控制策略需求的期望阻尼力,有效地衰减座椅下底板振动,对改善车辆乘坐舒适性具有积极作用.
4. 结 论
针对半主动座椅悬架控制系统,分析磁流变阻尼器非线性滞回特性,提出了一种改进魔术公式模型,并用改进差分进化算法进行模型辨识;利用自适应神经模糊系统建立磁流变阻尼器控制电流逆模型,并对包含阻尼器正逆模型的半主动座椅悬架进行验证计算. 结果表明,本文建立的正逆模性能够较好地表征磁流变阻尼器非线性滞回特性,满足座椅悬架的半主动控制需求.
1) Magic Formula 模型作为一种经典的轮胎模型,具有广泛的应用,对于非线性系统具有良好的拟合能力. 本文探索正弦型及余弦型魔术公式在磁流变减振器力学模型中的应用.
2) 在不同激励特性及电流条件下,本文建立的阻尼力模型均能较好地描述磁流变阻尼器的非线性滞回特性,对应于5 mm、3 Hz激励条件下,电流0~1.0 A变化范围内的模型平均百分比误差分别为0.0447、0.0474、0.0316、0.0310、 0.0344及0.0308.
3) 选取
$ x\left(t\right),\dot{x}\left(t\right),\ddot{x}\left(t\right),F\left(t\right),F\left(t-1\right) $ 作为模型输入,建立包含32 个“if-then”模糊规则的五输入单输出ANFIS磁流变逆动力学模型. 本文建立的逆模型对恒定电流工况及随机电流工况均可以得到较好的拟合效果,误差均方根值分别为0.0869 A及0.1171 A.4) 阻尼器逆向动力学模型与阻尼力模型串联模型对期望阻尼力的预测能力较好,阻尼器全工况范围内预测阻尼力误差均方根值与阻尼器最大阻尼力的百分比可达到5.6%.
5) 与被动座椅悬架相比,磁流变半主动座椅悬架在共振频率范围内均能够有效地衰减振动. on-off控制、连续天棚控制及LMS控制算法在4 Hz处的座椅悬架位移传递率较原被动座椅分别降低了47.62%、48.72%、42.09%. 阻尼器响应时滞增加会降低半主动控制的减振效果. 因此,应尽可能减小半主动座椅悬架系统各个环节的时滞.
-
表 1 试件参数
Table 1. Measured parameters of all specimens
组 试件编号 D/mm d/mm to/mm ti/mm L/mm r/mm ti.min/mm fsyi,min/MPa A S120-3-40 119.6 57.2 3.43 3.41 359.1 5 0.66 185 S120-4-40 120.7 57.6 4.18 3.51 358.4 2 0.90 336 S120-5-40 121.2 57.9 5.81 3.62 359.0 5 1.12 497 S120-3-60 119.8 57.8 3.42 3.47 359.2 5 0.66 183 B S120-4-60 120.6 57.6 4.14 3.51 359.0 2 0.90 332 S120-5-60 120.6 57.8 5.80 3.59 358.8 5 1.12 499 表 2 金属材料性能
Table 2. Measured material properties
类型 名称 σ0.2/MPa E/GPa f0.2,c/MPa μ n 不锈钢 S120-3 363 203 579 0.3 8 S120-4 550 200 1010 0.3 4 S120-5 586 204 999 0.3 4 碳钢 ϕ57-3 430 220 0.3 -
[1] HAN L H, TAO Z, HUANG H, et al. Concrete-filled double skin (SHS outer and CHS inner) steel tubular beam-columns[J]. Thin-Walled Structures, 2004, 42(9): 1329-1355. doi: 10.1016/j.tws.2004.03.017 [2] PAGOULATOU M, SHEEHAN T, DAI X H, et al. Finite element analysis on the capacity of circular concrete-filled double-skin steel tubular (CFDST) stub columns[J]. Engineering Structures, 2014, 72: 102-112. doi: 10.1016/j.engstruct.2014.04.039 [3] UENAKA K, KITOH H, SONODA K. Concrete filled double skin circular stub columns under compression[J]. Thin-Walled Structures, 2010, 48(1): 19-24. doi: 10.1016/j.tws.2009.08.001 [4] HASSANEIN M F, KHAROOB O F, LIANG Q Q. Circular concrete-filled double skin tubular short columns with external stainless steel tubes under axial compression[J]. Thin-Walled Structures, 2013, 73: 252-263. doi: 10.1016/j.tws.2013.08.017 [5] LI W, HAN L H, ZHAO X L. Axial strength of concrete-filled double skin steel tubular (CFDST) columns with preload on steel tubes[J]. Thin-Walled Structures, 2012, 56: 9-20. doi: 10.1016/j.tws.2012.03.004 [6] GUO Z, CHEN Y, WANG Y, et al. Experimental study on square concrete-filled double skin steel tubular short columns[J]. Thin-Walled Structures, 2020, 156: 107017.1-107017.20. [7] AYOUGH P, SULONG N H R, IBRAHIM Z, et al. Nonlinear analysis of square concrete-filled double-skin steel tubular columns under axial compression[J]. Engineering Structures, 2020, 216: 110678.1-110678.26. [8] HASSANEIN M F, KHAROOB O F, GARGNER L. Behaviour and design of square concrete-filled double skin tubular columns with inner circular tubes[J]. Engineering Structures, 2015, 100: 410-424. doi: 10.1016/j.engstruct.2015.06.022 [9] HASSANEIN M F, ELCHALAKANI M, KARRECH A, et al. Behaviour of concrete-filled double-skin short columns under compression through finite element modelling: SHS outer and SHS inner tubes[J]. Structures, 2018, 14: 358-375. doi: 10.1016/j.istruc.2018.04.006 [10] 丛术平,彭敏,王继升,等. 方中空夹层钢管混凝土短柱轴压性能试验[J]. 中国科技论文,2019,14(10): 1085-1089.CONG Shuping, PENG Min, WANG Jisheng, et al. Experimental research on axial compression performance of square double skin steel tube short columns filled with concrete[J]. China Sciencepaper, 2019, 14(10): 1085-1089. [11] HUANG H, HAN L H, TAO Z, et al. Analytical behaviour of concrete-filled double skin steel tubular (CFDST) stub columns[J]. Journal of Constructional Steel Research, 2010, 66(4): 542-555. doi: 10.1016/j.jcsr.2009.09.014 [12] HAN L H, LI Y J, LIAO F Y, et al. Concrete-filled double skin steel tubular (CFDST) columns subjected to long-term sustained loading[J]. Thin-Walled Structures, 2011, 49(12): 1534-1543. doi: 10.1016/j.tws.2011.08.001 [13] HAN L H, HUANG H, TAO Z, et al. Concrete-filled double skin steel tubular (CFDST) beam-columns subjected to cyclic bending[J]. Engineering Structures, 2006, 28(12): 1698-1714. doi: 10.1016/j.engstruct.2006.03.004 [14] WANG F Y, YOUNG B, GARDNER L, et al. CFDST sections with square stainless steel outer tubes under axial compression: experimental investigation, numerical modelling and design[J]. Engineering Structures, 2020, 207: 110189.1-110189.13. doi: 10.1016/j.engstruct.2020.110189 [15] WANG F C, HAN L H, LI W, et al. Analytical behavior of CFDST stub columns with external stainless steel tubes under axial compression[J]. Thin-Walled Structures, 2018, 127: 756-768. doi: 10.1016/j.tws.2018.02.021 [16] HAN L H, REN Q X, LI W, et al. Tests on stub stainless steel-concrete-carbon steel double-skin tubular (DST) columns[J]. Journal of Constructional Steel Research, 2011, 67(3): 437-452. doi: 10.1016/j.jcsr.2010.09.010 [17] TANG H Y, CHEN J L, FAN L Y, et al. Experimental investigation of FRP-confined concrete-filled stainless steel tube stub columns under axial compression[J]. Thin-Walled Structures, 2020, 146: 106483.1-106483.14. [18] 唐红元,李政周,范璐瑶,等. 矩形不锈钢管混凝土短柱轴压性能试验研究[J]. 西南交通大学学报,2022,57(4): 855-864.TANG Hongyuan, LI Zhengzhou, FAN Luyao, et al. Experimental investigation on behavior of rectangular concrete-filled stainless steel tubular stub columns under axial loading[J]. Journal of Southwest Jiaotong University, 2022, 57(4): 855-864. [19] 中国建筑科学研究院, 中华人民共和国建设部, 国家质量监督检验检疫总局. 普通混凝土力学性能试验方法标准: GB/T 50081—2002[S]. 北京: 中国建筑工业出版社, 2002. [20] 朱浩川,姚谏. 不锈钢材料的应力-应变模型[J]. 空间结构,2011,17(1): 62-69.ZHU Haochuan, YAO Jian. Stress-strain model for stainless steel[J]. Spatial Structures, 2011, 17(1): 62-69. [21] CRUISE R B, GARDNER L. Strength enhancements induced during cold forming of stainless steel sections[J]. Journal of Constructional Steel Research, 2008, 64(11): 1310-1316. doi: 10.1016/j.jcsr.2008.04.014 [22] WANG Y L, CAI G, LARBI A S, et al. Monotonic axial compressive behaviour and confinement mechanism of square CFRP-steel tube confined concrete[J]. Engineering Structures, 2020, 217: 110802.1-110802.16. [23] PHAM T M, HADI M N S. Stress prediction model for FRP confined rectangular concrete columns with rounded corners[J]. Journal of Composites for Construction, 2014, 18(1): 04013019.1-04013019.10. [24] HAN T H, STALLINGS J M, KANG Y J, et al. Nonlinear concrete model for double-skinned composite tubular columns[J]. Construction and Building Materials, 2010, 24(12): 2542-2553. doi: 10.1016/j.conbuildmat.2010.06.001 [25] KERR A D, SOIFER M T. The linearization of the prebuckling state and its effect on the determined instability loads[J]. Journal of Applied Mechanics, 1969, 36(4): 775-783. doi: 10.1115/1.3564770 [26] European Committee for Standardization. Eurocode 4: design of composite steel and concrete structures part1-1: general rules and rules for buildings[S]. London: British Standards Institution, 1994. [27] American Concrete Institute (ACI). Building code requirements for structural concrete and commentary: ACI 318-99[S]. Detroit: American Concrete Institute, 1999. [28] American Institute of Steel Construction. Specification for structural steel buildings: AISC 360 [S]. Chicago: [s.n.], 2016 [29] Standards Australia. Bridge design, part 6: steel and composite construction: AS5100.6—2004[S]. Sydney: Standards Australia International Ltd., 2004. [30] HAN L H, ZHAO X L, TAO Z. Tests and mechanics model for concrete-filled SHS stub columns, columns and beam-columns[J]. Steel and Composite Structures, 2001, 1(1): 51-74. doi: 10.12989/scs.2001.1.1.051 [31] 韩林海. 钢管混凝土结构——理论与实践[M]. 3版. 北京: 科学出版社, 2016. -