冷弯不等肢卷边角钢轴压试验及承载力设计方法

姚行友; 胡成立; 刘亚菲; 郭彦利

doi:10.3969/j.issn.0258-2724.20230010

冷弯不等肢卷边角钢轴压试验及承载力设计方法

doi: 10.3969/j.issn.0258-2724.20230010

1.
南昌工程学院土木与建筑工程学院，江西南昌 330099
2.
南昌工程学院江西省水利土木工程基础设施安全重点实验室，江西南昌 330099

基金项目: 国家自然科学基金项目（51868049）；江西省自然科学基金重点项目（20242BAB26074）

详细信息

作者简介:
姚行友（1978—），男，教授，研究方向为冷弯型钢结构，E-mail：yaoxingyoujd@163.com

中图分类号: TU392.1
计量
- 文章访问数: 292
- HTML全文浏览量: 77
- PDF下载量: 30
- 被引次数: 1
出版历程
- 收稿日期: 2022-10-10
- 修回日期: 2023-05-09
- 网络出版日期: 2024-08-05
- 刊出日期: 2023-07-07

Axial Compression Test and Bearing Capacity Design Method of Cold-Formed Steel with Unequal-Leg Lipped Angles

1.
School of Civil and Architectural Engineering, Nanchang Institute of Technology, Nanchang 330099, China
2.
Jiangxi Provincial Key Laboratory of Hydraulic & Civil Engineering Infrastructure Security, Nanchang Institute of Technology, Nanchang 330099, China

摘要

摘要:
为研究轴心受压冷弯薄壁不等肢卷边角钢屈曲性能和承载力设计方法，采用试验和有限元程序分析其屈曲性能和极限承载力，并基于直接强度法提出承载力设计建议方法. 首先，开展32根不同截面、长细比和宽厚比的冷弯薄壁不等肢卷边角钢轴压试验；随后，采用ABAQUS有限元软件对冷弯薄壁不等肢卷边角钢在不同宽厚比、肢宽比、长细比等条件下的屈曲性能和承载力进行参数化分析；最后，基于试验和有限元分析结果提出轴心受压冷弯薄壁不等肢卷边角钢构件承载力计算的修正直接强度法公式. 结果表明：宽厚比较小，试件易发生弯扭屈曲，宽厚比较大，试件易发生局部屈曲（或弯扭）和局部相关屈曲；构件极限承载力随长细比增大而降低，随着宽厚比的增大，承载力增长趋势逐渐减缓；提出的修正直接强度法计算轴心受压冷弯薄壁不等肢卷边角钢构件承载力的可靠度指标均值均大于3.2，建议方法准确且安全可靠.
- 冷弯薄壁型钢 /
- 不等肢卷边角钢 /
- 宽厚比 /
- 长细比 /
- 局部屈曲
Abstract:
To study the buckling performance and bearing capacity design method of the cold-formed thin-walled steel (CFTWS) with unequal-leg lipped angles under axial compression, the tests were conducted, and a finite element program was used to analyze the buckling performance and ultimate bearing capacity. Based on the direct strength method, the bearing capacity design method was presented. Firstly, the 32 CFTWS with unequal-leg lipped angles with different sections, slenderness ratios, and width-to-thickness ratios were tested under axial compression. Then, ABAQUS finite element software was used to perform parametric analysis on the buckling performance and bearing capacity of CFTWS with unequal-leg lipped angles with different width-to-thickness ratios, width ratios of legs, and slenderness ratios. Finally, the modified direct strength method for calculating the bearing capacity of CFTWS with unequal-leg lipped angles under axial compression was proposed based on the results of the test and finite element analysis. The results show that the specimens with a small width-to-thickness ratio are prone to flexural-torsional buckling, while the specimens with a large width-to-thickness ratio are prone to the local buckling or buckling related to flexural-torsional buckling and local buckling. The ultimate bearing capacity of the specimens decreases with the increase in the slenderness ratio, and the increase amplitude of the bearing capacity slows down with the increase in the width-to-thickness ratio. The reliability analysis shows that the reliability indexes of the bearing capacity of the CFTWS with unequal-leg lipped angles under axial compression calculated using the modified direct strength method are greater than 3.2, indicating that the proposed method is accurate and reliable.
- cold-formed thin-walled steel /
- steel with unequal-leg lipped angles /
- width-to-thickness ratio /
- slenderness ratio /
- local buckling

HTML全文

隧址区岩体初始应力场是隧道工程设计的重要指标，也是引起隧道开挖应力、位移释放的荷载来源. 近年来，国家建设重心逐步西移，西部地区隧道朝着长、大、深方向发展，该区域地质构造运动频繁，断层、节理较为发育，其分布的岩体原位地应力场更加复杂^[1].

国内外学者在原位地应力反演方面做了大量研究，大致思路可归结为利用地应力实测资料^[2-4]同时采用数值模拟方法进行求解^[5]：汪波等^[6]利用苍岭隧道地应力实测值，采用多元线性回归法获得了整个工程区地应力分布情况并进行了岩爆预测，但其地应力反演需经大量试算才能确定边界条件；许传华等^[7]根据紫金山矿区地应力资料，提出利用支持向量机和模拟退火算法进行反演分析，提高了反演效率，但对于侧压力系数样本量的选定仍不够充分，且计算量较大；王庆武等^[8]基于RBF (radical basis function)神经网络和地层剥蚀原理进行了拉林铁路桑日至加查段地应力反演，并对河谷岸坡地应力场分布特征进行了分析，得出了岸坡浅表部的最大主应力与岸坡走向相近并随埋深加大向竖向偏转的规律；徐卫中等^[9]依托重庆蟠龙抽水蓄能电站工程，采用多元回归方法、神经网络和遗传算法分别反演了地下厂房区域原位地应力，得出多元回归方法误差更小、计算时间更短的结论.

许多学者还利用反演手段对特殊复杂地质区域地应力场分布规律进行了研究：李科等^[10]探讨了断裂扰动区地应力场分布情况，提出受断裂扰动区影响，该区域附近地应力呈现不均匀特性，其量值、方向均发生较大改变；赵辰等^[11]结合地层剥蚀原理及侧压力系数法对托巴水电站地下厂房区域地应力进行了反演分析，同时考虑了复杂地质条件下地形地貌、断层等对地应力的影响，得出了河谷段由深至浅，最大主应力方向朝平行于山坡地形面方向偏转的规律；张敏等^[12]基于川藏铁路复杂的地质条件，采用应力试算法对桑珠岭隧道区域初始地应力进行了反演分析，其结论表明断层破碎带区域内竖向主应力有骤减的规律.

由此可见，针对复杂地质条件下地应力反演研究主要集中于探讨如断层、褶皱、河谷岸坡等特殊地质、地貌，鲜少有对于侵入岩地层初始地应力场分布规律的研究. 本文以宁南至攀枝花段高速公路工程火山隧道为依托，以水压致裂法实测地应力及隧址区实际地形地貌为基础，建立了三维数值长、大模型，采用最小二乘法寻优准则，结合多元线性回归法拟合精度、计算效率较高的特点，针对隧址区花岗质侵入岩地层附近初始地应力场进行了反演分析，探讨了侵入岩地层初始地应力分布规律，为火山隧道及类似工程设计、施工提供理论支持.

1. 工程概况

1.1 地形地貌及地层岩性

G4216线宁南至攀枝花段高速公路火山隧道全长7 145 m，最大埋深466.9 m. 纵断面如图1所示，SX-1~SX-5为竖向应力提取线；ZX-1~ZX-5为水平方向应力提取线. ZK390+755～ZK394+500段穿越地层岩性主要为泥岩、砂岩；ZK394+500～ZK397+900段围岩为花岗闪长岩. 隧道轴线穿越花岗质侵入岩附近，岩性变化剧烈，为硬、软岩交界带，且在高地应力作用下，坚硬的花岗闪长岩具备发生岩爆的条件，因此，十分必要进行隧址区域初始地应力场反演分析.

图 1 火山隧道纵断面

Figure 1. Longitudinal section of Huoshan tunnel

下载: 全尺寸图片幻灯片

1.2 研究区域地应力特征

线路区域位于川滇地区，以极为发育的南北向构造占主导，与之配套的北西、北东向构造和东西向构造均有发育，按照Anderson理论和摩尔库伦理论分析，研究区域应力场控制方向应为NNW～NWW向. 崔效锋等^[13]利用多种应力数据分析得到的结论是该区域整体的应力场方向为NNW向，约为N10°W，震源机制解类型主要为走滑型.

1.3 隧址区初始地应力测量

根据火山隧道地应力测试资料，在里程ZK395+320附近采用水压致裂法进行了钻孔火隧SK-02原位地应力测试，钻孔位置如图1所示，测试结果如表1所示，表中： ${\sigma _{\text{v}}}$ 为竖向主应力； ${\sigma _{\text{H}}}$ 为最大水平主应力； ${\sigma _{\text{h}}}$ 为最小水平主应力. 钻孔位置距本文重点研究的侵入体前沿区域较近，在该区域内反演误差相对较小，因此认为可以利用该钻孔数据进行花岗闪长岩侵入面附近地应力的反演分析.

表 1 SK-02钻孔水压致裂原位地应力测量结果

Table 1. Initial geostress measurement results of hydraulic fracturing method in SK-02 borehole

测点编号	测段深度/m	${\sigma _{\text{H}}}/{\rm MPa}$	${\sigma _{\text{h}}}/{\rm MPa}$	${\sigma _{\text{v}}}/{\rm MPa}$	${\sigma _{\text{H}}}$ 方位
1	165.6	6.14	4.07	4.22
2	220.7	7.43	5.01	5.62
3	261.3	8.15	5.65	6.66	N15°W
4	285.6	8.54	5.99	7.28	N9°W

下载: 导出CSV

| 显示表格

由表1可知：钻孔附近区域初始地应力场以水平主应力为主，构造应力占主导地位；4个孔深位置的地应力分布均呈现 ${\sigma _{\text{H}}}$ > ${\sigma _{\text{v}}}$ > ${\sigma _{\text{h}}}$ 的特征，且随着埋深的增加，各应力分量也相应增大；该孔两段的印模结果表明，最大水平主应力方向为N9°W～N15°W，测孔附近的地应力场以NNW向为主，测试结果与前述区域地应力分布特征吻合，验证了测量结果的可靠性.

2. 火山隧道隧址区地应力反演分析

2.1 三维数值模型的建立

以隧道线路平面约8000 m × 1000 m的长方形区域作为计算域（考虑隧道轴线附近几何信息的录入），模型底面取至隧道轴线以下300 m，基于地勘资料提供的等高线平面图建立顶面地表，导入Surfer软件并将克里金插值法得到的高精度地形坐标值输入有限元软件ANSYS中求得，具体数值模型如图2所示. 岩体材料参数如表2所示，在模型以砂岩、泥岩为主的区域，沿隧道轴线砂岩、泥岩以50～150 m不等的厚度大量交替分布，且无断层、褶皱等不连续带，同时考虑到花岗闪长岩弹性模量几乎是另外两种岩体的5倍，在实际材料参数的设定中按照砂岩、泥岩两种岩体所占权重，将该区域岩体材料参数均设定为表2中砂岩、泥岩.

图 2 有限元计算模型

Figure 2. Finite element calculation model

下载: 全尺寸图片幻灯片

表 2 岩体物理力学性质参数

Table 2. Physico-mechanical parameters of rock masses

岩体类型	弹性模量/GPa	泊松比	密度/（kg•m⁻³）
花岗闪长岩	13	0.28	2610
砂岩	3.5	0.32	2670
泥岩	1.5	0.40	2640
砂岩、泥岩	3	0.34	2660

下载: 导出CSV

| 显示表格

2.2 模型边界条件

由于水压致裂法所测得原位地应力未包含铅垂面内剪切应力，因此仅将对工程区域岩体初始地应力场影响较大的4种基本因素作为模型边界条件（所受荷载如图3所示，其中：G为自重荷载；F为水平构造荷载；U、V分别为X、Y方向平面剪切构造荷载）：1）自重应力场，对模型底面及侧面施加法向约束，并对模型施加Z向重力荷载；2） X向水平构造应力场，对模型两侧（左、右面）施加X方向水平均布荷载模拟该方向上的构造应力，同时对底面及未加载的两个侧面施加法向约束；3） Y向水平构造应力场，对模型两侧（前、后面）施加Y方向水平均布荷载，底面及未加载的两个侧面施加法向约束；4） XOY平面剪切构造应力场：通过对模型四个侧面施加边界位移用以等效剪应力^[14]，具体在前、后面施加0.01 m位移，左、右面施加0.08 m位移.

2.3 原位地应力坐标系的转化

现场实测原位地应力所采用的坐标系以主应力实际方向建立，而数值模型建立的坐标系显然与其不同. 为方便验证反演效果，需将地应力实测值进行坐标转化，变为与模型相同的坐标系. 按照弹性力学应力分量坐标变换理论，新坐标系下实测地应力分量坐标转化公式如式（1），其中， $v,w = x,y,{\textit{z}}$ .

$\left\{ { \begin{array}{l} {\sigma _{1'w'}} = {\alpha _{1'v}}{\alpha _{w'w}}{\sigma _{vw}},\\ {\sigma _{2'w'}} = {\alpha _{2'v}}{\alpha _{w'w}}{\sigma _{vw}},\\ {\sigma _{3'w'}} = {\alpha _{3'v}}{\alpha _{w'w}}{\sigma _{vw}}, \end{array} } \right.$

(1)

式中：σ_vw为原坐标系应力； ${\alpha _{1'v}}、{\alpha _{2'v}}、{\alpha _{3'v}}、{\alpha _{w'w}}$ 为坐标转化系数.

运用式（1）将表1实测原位地应力进行坐标转化，得到与数值模型相同坐标系下的应力分量，如表3所示.

图 3 模型荷载示意

Figure 3. Model load indication

下载: 全尺寸图片幻灯片

2.4 多元线性回归法反演原理

将第k观测点的原位地应力回归计算值 ${\hat \sigma _k}$ 作为因变量，将地应力实测点处有限元数值模拟求得的自重应力场和构造应力场作用下的计算值 $\sigma _{ki}$ （第k观测点，第i种工况）作为自变量，则多元线性回归方程为

${\hat {\boldsymbol{\sigma}} _k} = {b_0} + \sum\limits_{i = 1}^n {{b_i}{\boldsymbol{\sigma}} _{ki}} ,$

(2)

式中： ${\hat {\boldsymbol{\sigma}} _k}$ 、 ${\boldsymbol{\sigma}} _{ki}$ 分别为原位地应力回归计算值单列矩阵、有限元计算的单列矩阵； ${b_0}$ 为自由项， ${b_i}$ 为多元回归系数， $n$ 为工况数（共计4种工况， $n = 4$ ）.

统计学中，利用最小二乘估计的方法求得回归系数 ${b_i}$ 是求解多元线性回归方程的关键. 假定有 $m$ 个观测点，每一个测点有6个观测分量，则最小二乘法的残差平方和为

${Q_{{\rm{can}}}} = \sum\limits_{k = 1}^m {{{\sum\limits_{j = 1}^6 {\left( {{\sigma _{jk}} - {b_0} - \sum\limits_{i = 1}^n {{b_i}} \sigma _{jki}} \right)} }^2}} ,$

(3)

式中： $\sigma _{jk}^{}$ 为第 $k$ 观测点 $j$ 应力分量的观测值，应力分量 $j = 1,2, \cdots ,6$ ，对应6个初始应力分量； $\sigma _{jki}$ 为第 $i$ 种工况下第 $k$ 观测点 $j$ 应力分量的有限元数值模拟计算值.

依据最小二乘法原理,观测值与回归值之间的残差平方和将达到最小使得残差平方和为最小值的方程式如式（4）所示. 求解该方程式，得到4+1个待定回归系数 ${\boldsymbol{b}} = {({b_0},{b_1},{b_2},{b_3},{b_4})^{\text{T}}}$ ，则隧址区域内所有点的回归初始地应力可通过式（2）得到.

表 3 实测原位地应力值与回归值对比

Table 3. Comparison of measured initial geostress values and regression values MPa

测点编号	${\sigma _X}$				${\sigma _Y}$				${\sigma _{\textit{Z}}}$				${\tau _{XY}}$
测点编号	实测值	回归值	绝对误差	相对误差/%	实测值	回归值	绝对误差	相对误差/%	实测值	回归值	绝对误差	相对误差/%	实测值	回归值	绝对误差	相对误差/%
1	−4.15	−4.59	0.44	10.6	−6.06	−7.19	1.13	18.6	−4.22	−3.96	0.26	6.2	−0.41	−0.53	0.12	29.3
2	−5.11	−5.17	0.06	1.2	−7.33	−7.29	0.04	0.5	−5.62	−5.66	0.04	0.7	−0.48	−0.48	0	0
3	−5.75	−5.60	0.15	2.6	−8.05	−7.41	0.64	8.0	−6.66	−6.88	0.22	3.3	−0.50	−0.44	0.06	12.0
4	−6.10	−5.84	0.26	4.3	−8.44	−7.50	0.94	11.1	−7.28	−7.60	0.32	4.4	−0.51	−0.41	0.10	19.6
注： ${\sigma _X}、{\sigma _Y}、{\sigma _{\textit{Z}}}$ 分别为 $X、Y、Z$ 方向正应力值； ${\tau _{XY}}$ 为 $XOY$ 平面剪应力.

下载: 导出CSV

| 显示表格

$\left[ {\begin{array}{*{20}{c}} {6m}&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk1}} } }&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk2}} } }& \cdots &{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jkn}} } }\\ {}&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {{{(\sigma _{jk1})}^2}} } }&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk1}\sigma _{jk2}} } }& \cdots &{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk1}\sigma _{jkn}} } }\\ {}&&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {{{(\sigma _{jk2})}^2}} } }& \cdots &{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk2}\sigma _{jkn}} } }\\ {}&{}&& \ddots & \vdots \\ {}&{}&{}&{}&{\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {{{(\sigma _{jkn})}^2}} } } \end{array}} \right] \left[ {\begin{array}{*{20}{c}} {{b_0}}\\ {{b_1}}\\ {{b_2}}\\ \vdots \\ {{b_n}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {{\sigma _{jk}}} } }\\ {\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk1}{\sigma _{jk}}} } }\\ {\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jk2}{\sigma _{jk}}} } }\\ \vdots \\ {\displaystyle\sum\limits_{k = 1}^m {\displaystyle\sum\limits_{j = 1}^6 {\sigma _{jkn}{\sigma _{jk}}} } } \end{array}} \right] .$

(4)

3. 隧址区地应力反演结果

3.1 多元回归系数的确定

基于多元线性回归法原理，利用式（4）解得多元回归系数：自由项 ${b_0}{\text{ = 1}}{\text{.57}}$ ，自重应力场回归系数 ${b_{\text{1}}}{\text{ = 1}}{\text{.08}}$ ，X向构造应力场 ${b_{\text{2}}}{\text{ = 5}}{\text{.0}}$ ，Y向构造应力场 ${b_{\text{3}}}{\text{ = 4}}{\text{.23}}$ ，水平剪切构造应力场 ${b_{\text{4}}}{\text{ = 19}}{\text{.33}}$ . 依据统计学原理，对多元回归效果进行显著性检验，计算得到负相关系数 $R{\text{ = 0}}{\text{.987}}$ ，接近于1，说明自变量与因变量之间相关程度很高，且 ${{F}}$ 显著性统计量 ${{F = }} {\text{1.2}} \times {\text{1}}{{\text{0}}^{{{ - 8}}}}$ ，远小于显著性水平0.01，同样说明回归效果显著. 则火山隧道隧址区岩体原位地应力场回归反演值为

$\begin{split} &{\sigma }_{{\rm{hui}}}={1.57+1}{.08}{\sigma }_{{\rm{zi}}}+{5}{.0}{\sigma }_{{{\rm{gou1}}}}+\\ &\quad { 4}{.23}{\sigma }_{{{\rm{gou2}}}}+{19}{.33}{\sigma }_{{{\rm{gou3}}}},\end{split}$

(5)

式中： ${\sigma }_{{\rm{zi}}}$ 为自重应力场； ${\sigma }_{{{\rm{gou1}}}}$ 、 ${\sigma }_{{{\rm{gou2}}}}$ 分别为X、Y方向均匀挤压构造运动引起的应力场； ${\sigma }_{\rm{gou3}}$ 为水平面内纯剪构造运动引起的应力场.

利用有限元数值模拟得到相应于地应力实测点处的计算值，代入式（5）后即可求出该点处的回归反演值，表3中列出了火隧SK-02钻孔不同深度实测值与有限元回归值. 其中正应力反演结果较为接近于实测值，最小绝对误差为0.04 MPa，最大绝对误差为1.13 MPa，最大相对误差为18.6%，最小相对误差仅0.5%. 剪切应力反演结果的绝对误差很小，但由于剪切应力量值本就较小，导致相对误差偏大. 方位角的反演值在测点3处为N12.9°W，测点4为N13.1°W，对比表1实测值可知误差最大不超过4°. 总体而言反演效果较好，且结果可靠.

3.2 隧道轴线初始地应力场分布探究

首先，对花岗质侵入岩内各重要部位进行名称的约定，方便后文讨论. 将花岗闪长岩侵入体内右侧临近砂、泥岩区域称为“侵入体前沿区域”（如图1所示），将花岗闪长岩与砂、泥岩交接面称为“侵入面”. 受限于火山隧道钻孔数量较少，虽距钻孔位置较远处存在一定的误差影响，在重点讨论侵入体前沿区域附近地应力特征之前，仍有必要对隧道整体轴线的初始地应力场分布进行探讨.

基于区域地应力场反演结果，进一步对隧道轴线（图1中ZX-2）初始地应力场分布规律进行探讨，考虑到模型边界效应带来的误差影响，在隧道轴线初始地应力提取时，剔除了左、右洞口段约800 m范围内的应力值，其分布见图4所示.

图 4 隧道轴线初始地应力分布

Figure 4. Initial geostress distribution of the tunnel axis

下载: 全尺寸图片幻灯片

从图中可以看出：竖向应力的分布与隧道埋深规律相似，即随着埋深的增加，竖向主应力增大，这是由于竖向主应力主要源自岩体自重应力；在两个山峰最高点附近区域ZK396+200～ZK395+400、ZK394+500～393+500内初始应力场以竖向主应力占主导；轴线上其他位置竖向应力的分布基本符合山体轮廓线的起伏.

水平主应力的分布规律较竖向应力更加复杂，在ZK397+000～ZK394+600花岗闪长岩区域，水平主应力分布规律与埋深关联性不大，先小幅度减小后急剧增大. 结合图1，花岗质侵入岩横截面中间部位较宽，两侧相对较窄，而较窄的横截面区域内侵入体受挤压构造作用更大，应力集中现象明显，这是造成花岗闪长岩区域水平主应力值中间小、两侧大的原因. 同时花岗质侵入岩产状为岩床型，其经过地壳断裂张力持续作用下挤入围岩层间空隙中，在隧道轴线接近侵入体前沿区域，在地质构造作用下所产生的构造应力值也就越大，这是造成邻近侵入面处应力急剧增加的原因. 图1中右侧砂、泥岩范围内，由于岩性更加单一，水平主应力分布较为稳定，这与该区域埋深变化不大相符合，其中最大水平主应力随埋深近乎线性变化，最小水平主应力则在小范围内浮动.

总体来看，隧道轴线初始应力场由于花岗质侵入岩的存在呈较为复杂的分布：1）在隧道桩号ZK396+200～ZK395+400内 ${\sigma _{\text{v}}} > {\sigma _{\text{H}}} > {\sigma _{\text{h}}}$ ，在隧道桩号ZK393+200～ZK392+100范围内 ${\sigma _{\text{H}}} > {\sigma _{\text{v}}} > {\sigma _{\text{h}}}$ ，其余区段由于受埋深及岩性变化的影响，分布规律较为离散. 2）花岗质侵入岩一侧（图1中左侧区域）主应力量值均整体大于泥岩、砂岩一侧（图1中右侧区域），这是由于受岩性变化的影响，花岗闪长岩弹性模量较大，造成区域内主应力量值较高. 3）主应力值在花岗质侵入岩边界ZK394+600附近均急剧降低，是由于岩性突变引起，造成花岗闪长岩一侧应力集中现象明显，其中最大水平主应力受其影响较大，减幅5.30 MPa；而竖向应力及最小水平主应力受影响较小，减幅分别为2.70、2.00 MPa.

4. 花岗质侵入岩区域初始地应力分布

4.1 花岗质侵入岩区域初始地应力量值分布规律

综合已有研究，大多分析了埋深对初始地应力场分布的影响，鉴于本文工程背景较为特殊，考虑到地层中花岗质侵入岩对初始地应力场的扰动，其分布规律的影响因素显然不单只有埋深，侵入体分布形态往往决定了围岩的受挤压程度，因此也应从侵入体分布形态的角度进行讨论.

基于隧址区初始地应力场反演结果，在桩号ZK395+400～ZK393+500区域内，以竖向间隔每60 m提取与隧道轴线ZX-2平行的ZX-1、ZX-3、ZX-4、ZX-5沿线（如图1）的地应力结果，如图5所示.

从图5（a）可知：在花岗闪长岩侵入面处，ZX方向5条轴线最大水平主应力有明显的衰减，由于不同埋深的5条轴线与花岗闪长岩侵入面交点不同，因此其发生衰减的里程桩号也不一样，按照ZX-5～ZX-1的顺序先后骤减，其正好符合花岗闪长岩边界线的分布；每条曲线衰减幅度并不相同，每条曲线应力量值的减小量如表4所示，其中ZX-1曲线在最大水平主应力衰减处埋深约为372 m，应力减小幅度为4.5 MPa，ZX-5曲线应力衰减处埋深约为128 m，应力减小幅度为5.7 MPa；除ZX-4外，埋深较浅的ZX-3、ZX-5曲线应力减小量大于埋深较深的ZX-1及ZX-2曲线应力减小量，由此可见在花岗质侵入岩边界处最大水平主应力减小幅度在埋深较浅处大于埋深较深处；埋深越大，最大水平主应力分量受地层岩性改变的影响越小，量值越趋于稳定；在最大水平主应力发生骤减之前，即花岗闪长岩区域内，ZX-1～ZX-5应力值相差不大，反映出最大水平主应力在花岗闪长岩内部分布较为稳定，而最大水平主应力发生骤减之后，最大水平主应力分布随埋深的增大而增大.

由图5（b）、（c）可知：竖向应力和最小水平主应力主要受埋深影响明显，无论是花岗闪长岩内部，还是外部，其主应力量值均随埋深的增加而增加；在花岗闪长岩边界，主应力量值减小幅度较最大水平主应力更小，且变化趋势更加平缓.

图 5 ZX-1～ZX-5初始地应力分布

Figure 5. Initial geostress distribution (ZX-1−ZX-5 )

下载: 全尺寸图片幻灯片

表 4 最大水平主应力减小量值

Table 4. Maximum horizontal principal stress reduction magnitude

轴线编号	应力减小处埋深/m	减小幅度/MPa
ZX-1	372	4.5
ZX-2	324	5.3
ZX-3	296	6.3
ZX-4	214	4.5
ZX-5	128	5.7

下载: 导出CSV

| 显示表格

花岗闪长岩区域5条竖向应力提取线SX-1～SX-5的初始地应力分布如图6所示.

图 6 SX-1～SX-5初始地应力分布

Figure 6. Initial geostress distribution (SX-1−SX-5)

下载: 全尺寸图片幻灯片

由图6可知，越接近花岗质侵入岩前沿，侵入面斜率越大，主应力量值越大. 1）其中最大水平主应力更为明显，图6（a）中在花岗闪长岩内部（高程约1.90 km以下），SX-5的最大水平主应力平均比其余曲线大约3.0 MPa，而其余曲线最大水平主应力量值相差不大. 高程在花岗闪长岩以上时，其量值分布符合埋深的变化情况，埋深较浅的SX-3、SX-4沿线最大水平主应力值更小. 2）竖向应力（图6（b）），在边界衰减幅度变弱，量值随高程几乎线性分布，符合竖向应力分布特征. 3）最小水平主应力的分布将其背后机理呈现得更为明显，首先从距离花岗质侵入岩前沿区域远近而言，图6（c）中越接近花岗质侵入岩前沿，受构造挤压作用越大，最小水平主应力量值也越大，如在花岗闪长岩内部，最小水平主应力位置从大到小依次为SX-5、SX-4、SX-3、SX-1、SX-2，仅在SX-1及SX-2处出现例外. 结合图1可知，形成这一例外是由于SX-1的埋深远大于SX-2对应位置的埋深. 这说明在侵入体内部，埋深和侵入体分布形态均是影响初始地应力分布的重要因素.

4.2 侵入面区域初始地应力方位角分布规律

在花岗闪长岩侵入面，由于岩性的突变及侵入构造作用的影响，主应力方向也会发生较大改变.

提取图1中ZX-1～ZX-5及SX-1～SX-5沿线与花岗闪长岩边界线交点两侧的回归计算结果，为便于讨论，以下对方位角的叙述均采用主应力方向与X轴（隧道轴线）的夹角代替，各交点两侧不同岩性的最大水平主应力与X轴夹角如表5所示. 图7给出了表5中各点夹角的示意图，图中将各点最大水平主应力及最小水平主应力方向投影到坐标系的水平面上，并引入了X轴以方便讨论.

表 5 最大水平主应力与X轴夹角

Table 5. Angles between maximum horizontal principal stress and X-axis

测点编号	侵入体内侧 ${\sigma _{\text{H}}}$ 与 X 轴夹角/（°）	侵入体外侧 ${\sigma _{\text{H}}}$ 与 X 轴夹角/（°）
1	88.1	−86.7
2	89.6	−74.8
3	89.0	−74.6
4	83.6	−54.9
5	53.5	−60.8
6	79.7	−86.2
7	77.9	−54.9
8	88.5	−90.0
9	89.9	−87.4
10	89.4	−89.4
11	86.7	−56.1
12	15.2	−35.5
注：最大主应力与X轴夹角以顺时针为正，以逆时针为负.

下载: 导出CSV

| 显示表格

从表5可以发现，最大水平主应力方向在花岗闪长岩及砂、泥岩中发生了大角度变化. 图7中测点1在花岗闪长岩内最大水平主应力方向为88.1°，而到了花岗闪长岩外侧（砂、泥岩内侧），方向剧烈变化为 −86.7°. 另一方面，除个别点外，大致存在着埋深越大，最大水平主应力与X轴夹角越大，如埋深较浅的测点12和5，花岗闪长岩内最大水平主应力与X轴夹角分别为15.2° 及53.5°，而在埋深较深的测点10和1，夹角分别达到了89.4° 和88.1°. 从花岗闪长岩分布形态对最大水平主应力方向的影响可知，越接近侵入体前沿，测点4、3、2、1最大水平主应力夹角逐渐变大，同时花岗闪长岩内外两侧最大水平主应力之间的夹角也逐渐变大，这与4.1节中讨论量值变化时的规律一致，即越接近侵入体前沿，主应力受构造作用影响更为明显. 而在远离花岗闪长岩前沿区域该规律并不明显，大部分区域夹角的分布主要还是受埋深的影响.

图 7 最大水平主应力与X轴夹角示意

注：图中长线表示最大水平主应力方向，短线表示最小水平主应力方向.

Figure 7. Diagram of the angles between maximum horizontal principal stress and X-axis

下载: 全尺寸图片幻灯片

5. 结　论

1）隧址区初始地应力场受花岗闪长岩侵入体地层的影响，其分布规律较一般工况更加复杂. 在侵入体内部，埋深和侵入体分布形态均是影响初始地应力场的重要因素.

2）主应力量值在侵入体内侧到外侧出现急剧降低的特性，其中最大水平主应力减小幅度最大. 同时侵入体一侧的主应力量值整体大于砂岩、泥岩一侧. 隧道轴线方向：竖向应力主要受埋深的影响，其随着埋深的增大而增大；最大水平主应力和最小水平主应力在侵入体内部的分布与埋深关联性较弱，两者与对应位置处侵入体横截面大小及距侵入体侵入面位置有关，侵入体横截面越小、距侵入面越近，水平主应力量值越大. 侵入体前沿区域挤压构造作用更加明显，其应力量值明显大于侵入体其余区域；除侵入体前沿区域外，最大水平主应力量值变化不大；同时经过侵入面前后应力量值的减小量与埋深成反比.

3）侵入体内外侧水平主应力方向发生大角度变化，最大水平主应力与X轴夹角随埋深的增大而增大.

4）探明侵入岩地层隧址区地应力分布规律，对预测隧道未开挖区域应力状态提供了依据，可对类似工程提供参考. 施工中应加强现场地应力实测，在修正反演结果的基础上实现动态预测、动态设计.

图 1 不等肢卷边角钢截面

Figure 1. Section of steel with unequal-leg lipped angles

下载: 全尺寸图片幻灯片

图 2 应力-应变曲线

Figure 2. Stress-strain curves

下载: 全尺寸图片幻灯片

图 3 纵向初始缺陷量测位置

Figure 3. Measuring locations of longitudinal initial defect

下载: 全尺寸图片幻灯片

图 4 试件初始缺陷

Figure 4. Initial defect of specimens

下载: 全尺寸图片幻灯片

图 5 试验装置

Figure 5. Test setup

下载: 全尺寸图片幻灯片

图 6 位移计布置

Figure 6. Layout of displacement meters

下载: 全尺寸图片幻灯片

图 7 长度为400 mm的角钢柱屈曲模式

Figure 7. Buckling mode of steel angle with length of 400 mm

下载: 全尺寸图片幻灯片

图 8 长度为900、1500 mm的角钢柱屈曲模式

Figure 8. Buckling mode of steel angle with lengths of 900 mm and 1 500 mm

下载: 全尺寸图片幻灯片

图 9 长度为2100 mm的角钢柱屈曲模式

Figure 9. Buckling mode of steel angle with length of 2 100 mm

下载: 全尺寸图片幻灯片

图 10 荷载-位移曲线

Figure 10. Load-displacement curves

下载: 全尺寸图片幻灯片

图 11 有限元分析模型

Figure 11. Finite element analysis model

下载: 全尺寸图片幻灯片

图 12 有限元分析一阶屈曲模态

Figure 12. First-order buckling mode of finite element analysis

下载: 全尺寸图片幻灯片

图 13 试验与有限元分析屈曲模式对比

Figure 13. Comparison of buckling modes between test and finite element analysis

下载: 全尺寸图片幻灯片

图 14 试验与有限元分析荷载-位移曲线对比

Figure 14. Comparison of load-displacement curves between test and finite element analysis

下载: 全尺寸图片幻灯片

图 15 极限承载力和长细比关系曲线

Figure 15. Relationship between ultimate bearing capacity and slenderness ratio

下载: 全尺寸图片幻灯片

图 16 试验和有限元分析角钢屈曲模式和承载力对比

Figure 16. Comparisons of buckling modes and bearing capacities of steel angles between test and finite element analysis

下载: 全尺寸图片幻灯片

表 1 试件名义几何尺寸

Table 1. Nominal geometric dimension of specimens

试件编号	b_l/mm	b₂/mm	a_l/mm	a₂/mm	t/mm	b₁/t	b₂/t
LL4030	40	30	15	15	2	20	15
LL6040	60	40	20	20	2	30	20
LL9060	90	60	20	20	2	45	30
LL12080	120	80	24	24	2	60	40

下载: 导出CSV

表 2 试件实测几何尺寸

Table 2. Measured geometric dimensions of specimens

试件编号	b_l/mm	b₂/mm	a_l/mm	a₂/mm	t/mm	长度 L/mm	长细比 λ
LL4030-400-1	42.13	31.43	15.97	15.54	1.98	400.00	39.19
LL4030-400-2	41.51	31.45	15.56	16.34	1.98	399.27	39.88
LL4030-900-1	42.00	31.61	15.96	14.84	1.99	900.00	87.28
LL4030-900-2	42.25	31.05	15.14	15.79	2.00	900.00	89.45
LL4030-1500-1	41.98	31.48	15.99	15.28	1.97	1498.37	147.14
LL4030-1500-2	41.93	31.42	15.93	15.95	2.00	1499.33	146.65
LL4030-2100-1	42.07	31.32	15.60	15.96	1.98	2101.17	206.51
LL4030-2100-2	42.20	31.59	15.52	15.66	1.98	2101.27	206.93
LL6040-400-1	61.68	41.48	21.63	19.48	1.98	400.00	29.54
LL6040-400-2	61.35	41.52	21.17	21.00	1.97	399.00	29.40
LL6040-900-1	62.05	41.60	21.45	19.37	1.98	900.00	66.46
LL6040-900-2	62.16	41.63	21.54	19.61	2.09	900.00	68.25
LL6040-1500-1	62.10	41.64	21.06	20.48	1.99	1499.43	110.40
LL6040-1500-2	61.80	41.76	21.35	20.98	2.07	1499.27	111.54
LL6040-2100-1	62.24	41.37	21.32	19.86	1.91	2101.10	152.65
LL6040-2100-2	61.44	41.91	20.92	20.47	1.98	2100.00	154.12
LL9060-400-1	91.80	61.41	20.25	20.02	1.97	400.00	22.15
LL9060-400-2	91.92	61.31	20.12	20.50	1.97	399.50	22.16
LL9060-900-1	91.51	61.43	19.90	21.05	1.99	900.00	49.79
LL9060-900-2	91.90	61.89	20.67	20.19	2.00	900.00	49.39
LL9060-1500-1	91.50	62.59	20.51	20.88	2.00	1498.90	81.56
LL9060-1500-2	91.43	62.42	21.08	19.88	1.97	1499.10	81.69
LL9060-2100-1	92.07	61.52	20.99	19.69	1.98	2100.00	115.66
LL9060-2100-2	91.60	61.94	21.83	21.01	2.00	2101.03	113.94
LL12080-400-1	121.53	81.77	24.92	23.98	1.98	400.00	16.95
LL12080-400-2	121.30	81.82	24.19	24.22	1.98	399.27	17.00
LL12080-900-1	122.30	81.71	25.11	24.05	1.97	900.00	38.10
LL12080-900-2	121.87	81.95	23.68	24.90	1.97	900.00	38.23
LL12080-1500-1	121.97	82.20	24.71	24.17	2.00	1498.22	63.34
LL12080-1500-2	122.48	81.87	24.95	24.70	1.97	1499.13	63.33
LL12080-2100-1	122.79	81.89	23.76	24.03	1.98	2101.17	89.38
LL12080-2100-2	121.04	82.24	24.71	24.49	1.97	2100.14	88.61

下载: 导出CSV

表 3 试验和有限元分析角钢屈曲模式和承载力对比

Table 3. Comparisons of buckling modes and bearing capacities of steel angles between test and finite element analysis

试件编号	屈曲模式			Δ_Lmax/ mm	Δ_Gmax/ mm	P_t/kN	P_FE/kN	P_D/kN	P_MD/kN	P_Y/kN	P_t/P_Y	P_t/ P_FE	P_t/ P_D	P_t/P_MD
试件编号	试验	有限元	一阶模态	Δ_Lmax/ mm	Δ_Gmax/ mm	P_t/kN	P_FE/kN	P_D/kN	P_MD/kN	P_Y/kN	P_t/P_Y	P_t/ P_FE	P_t/ P_D	P_t/P_MD
LL4030-400-1	FT	FT	F	0.49	0.28	79.66	83.17	83.77	74.38	72.51	1.0986	0.9578	0.9509	1.0710
LL4030-400-2	FT	FT	F	0.41	0.34	79.38	82.66	83.67	73.51	73.41	1.0813	0.9603	0.9488	1.0799
LL4030-900-1	FT	FT	F	1.30	0.83	50.68	53.93	56.47	50.71	27.68	1.8312	0.9397	0.8974	0.9995
LL4030-900-2	FT	FT	F	1.54	0.77	52.72	54.70	56.90	50.97	27.97	1.8852	0.9638	0.9265	1.0342
LL4030-1500-1	FT	FT	F	1.30	0.38	40.70	40.93	40.36	42.15	23.01	1.7689	0.9944	1.0085	0.9655
LL4030-1500-2	FT	FT	F	1.14	0.17	40.98	41.52	41.89	43.75	23.88	1.7159	0.9869	0.9783	0.9366
LL4030-2100-1	FT	FT	F	1.84	0.08	30.52	31.19	24.48	25.57	13.96	2.1867	0.9786	1.2467	1.1936
LL4030-2100-2	FT	FT	F	1.10	0.09	29.26	31.35	24.50	25.59	13.97	2.0948	0.9335	1.2343	1.1434
LL6040-400-1	L	L	L	0.67	0.46	103.70	105.82	102.94	95.81	119.12	0.8706	0.9800	1.0074	1.0824
LL6040-400-2	L	L	L	0.75	0.31	102.46	105.26	105.33	96.41	121.73	0.8417	0.9734	0.9728	1.0628
LL6040-900-1	FT	FT	F	1.08	0.63	83.69	85.39	94.17	76.27	52.34	1.5989	0.9800	0.9886	1.0973
LL6040-900-2	FT	FT	F	1.17	0.74	83.62	84.70	98.45	80.28	54.25	1.5414	0.9872	0.9494	1.0416
LL6040-1500-1	FT	FT	F	1.15	0.60	59.96	61.98	74.21	61.37	42.31	1.4172	0.9673	0.9480	0.9769
LL6040-1500-2	FT	FT	F	1.42	0.43	60.82	63.31	78.23	64.30	44.60	1.3638	0.9607	0.9775	0.9459
LL6040-2100-1	FT	FT	F	1.17	0.24	48.66	50.69	41.59	43.44	23.71	2.0519	0.9600	1.1699	1.1201
LL6040-2100-2	FT	FT	F	1.42	1.08	49.94	50.96	42.99	44.90	24.51	2.0376	0.9800	1.1617	1.1122
LL9060-400-1	L	L	L	0.99	0.35	126.76	129.35	104.56	124.03	153.51	0.8257	0.9800	1.2123	1.0220
LL9060-400-2	L	L	L	0.57	0.57	125.71	129.60	106.03	124.39	155.32	0.8094	0.9700	1.1855	1.0106
LL9060-900-1	FT + L	FT + L	L	1.12	0.71	94.42	96.00	108.68	101.84	68.10	1.3865	0.9835	0.9688	0.9271
LL9060-900-2	FT + L	FT + L	L	1.02	0.92	95.60	97.20	107.99	102.07	67.07	1.4253	0.9835	0.9853	0.9366
LL9060-1500-1	FT + L	FT + L	L	0.79	0.71	77.34	80.36	90.86	94.90	51.80	1.4929	0.9625	0.9512	0.8149
LL9060-1500-2	FT + L	FT + L	L	0.88	0.83	77.04	78.87	85.54	89.34	48.77	1.5797	0.9768	0.9006	0.8623
LL9060-2100-1	D	D	F	1.22	0.90	59.95	61.37	48.74	50.90	27.79	2.1574	0.9768	1.2300	1.1776
LL9060-2100-2	D	D	F	1.08	0.77	61.51	63.42	54.77	57.20	31.22	1.9701	0.9700	1.2232	1.0754
LL12080-400-1	L	L	L	0.59	0.22	147.54	150.67	117.71	146.76	224.28	0.6578	0.9792	1.2535	1.0053
LL12080-400-2	L	L	L	0.58	0.14	146.78	150.41	117.29	146.49	223.06	0.6580	0.9759	1.2515	1.0020
LL12080-900-1	FT + L	FT + L	L	1.20	0.56	136.59	138.37	117.55	118.57	105.70	1.2923	0.9872	1.1620	1.1520
LL12080-900-2	FT + L	FT + L	L	1.52	1.11	131.36	138.86	118.35	118.60	106.10	1.2380	0.9460	1.1099	1.1075
LL12080-1500-1	FT + L	FT + L	L	1.20	1.56	97.14	98.99	116.57	104.35	81.31	1.1948	0.9813	0.8334	0.9309
LL12080-1500-2	FT + L	FT + L	L	1.52	1.11	98.56	98.19	116.67	102.58	83.41	1.1817	1.0038	0.8447	0.9608
LL12080-2100-1	FT + L	FT + L	L	1.44	1.42	85.43	85.91	73.43	76.89	41.97	2.0354	0.9948	1.2634	1.1110
LL12080-2100-2	FT + L	FT + L	L	1.45	1.54	85.82	86.48	76.31	80.96	44.19	1.9419	0.9926	1.2246	1.0600
均值											1.4760	0.9741	1.1246	1.0318
方差											0.2060	0.0002	0.0227	0.0077
变异系数											0.1395	0.0002	0.0201	0.0074

下载: 导出CSV

表 4 冷弯薄壁型不等肢卷边角钢轴压构件可靠指标

Table 4. Reliability indexes of CFTWS with unequal-leg lipped angles under axial compression

类型	ρ₂	组合形式	ρ₁=0.5	ρ₁=1.0	ρ₁=2.0	ρ₁=3.0	可靠指标平均值
办公楼	0	1.3G + 1.5L₁	4.518	4.467	4.442	4.403	4.457
	0.5	1.3G + 0.9 × 1.5 （L₁ + W）	4.306	4.298	4.253	4.198	4.264
	1.0	1.3G + 0.9 × 1.5 （L₁ + W）	4.178	4.161	4.102	4.087	4.132
	2.0	1.3G + 0.9 × 1.5 （L₁ + W）	4.096	4.011	3.918	3.824	3.962
	3.0	1.3G + 0.9 × 1.5 （L₁ + W）	3.954	3.898	3.837	3.763	3.863
	4.0	1.3G + 0.9 × 1.5 （L₁ + W）	3.781	3.675	3.621	3.601	3.670
住宅	0	1.3G + 1.5L₂	4.924	4.869	4.841	4.799	4.858
	0.5	1.3G + 0.9 × 1.5 （L₂ + W）	4.693	4.684	4.635	4.575	4.647
	1.0	1.3G + 0.9 × 1.5 （L₂ + W）	4.554	4.535	4.471	4.454	4.503
	2.0	1.3G + 0.9 × 1.5 （L₂ + W）	4.464	4.372	4.270	4.168	4.318
	3.0	1.3G + 0.9 × 1.5 （L₂ + W）	4.309	4.248	4.182	4.101	4.210
	4.0	1.3G + 0.9 × 1.5 （L₂ + W）	4.121	4.005	3.947	3.925	3.999
注：ρ₁ 为活荷载标准值与风荷载标准值之和与恒荷载标准值的比值，ρ₂ 为风荷载标准值与活荷载标准值的比值.

下载: 导出CSV

参考文献(18)

[1]	陈明,鲁卫波,武志远,等. 双肢冷弯C型钢蒙古包刚架的平面内稳定性[J]. 西南交通大学学报,2022,57(1):215-222. CHEN Ming, LU Weibo, WU Zhiyuan,et al. Investigation on in-plane stability of double-limb cold-formed C-shaped steel rigid frame of yurt[J]. Journal of Southwest Jiaotong University, 2022, 57(1):215-222.
[2]	ADLURI S M R, MADUGULA M K S. Flexural buckling of steel angles: experimental investigation[J]. Journal of Structural Engineering, 1996, 122(3): 309-317. doi: 10.1061/(ASCE)0733-9445(1996)122:3(309)
[3]	DINIS P B, CAMOTIM D, SILVESTRE N. On the mechanics of thin-walled angle column instability[J]. Thin-Walled Structures, 2012, 52: 80-89. doi: 10.1016/j.tws.2011.12.007
[4]	American Iron and Steel Institute. North American specification for the design of cold-formed steel structural members: AISI S100-16[S]. Washington D. C.: American Ir-on and Steel Institute, 2016.
[5]	YOUNG B. Tests and design of fixed-ended cold-formed steel plain angle columns[J]. Journal of Structural Engineering, 2004, 130(12): 1931-1940. doi: 10.1061/(ASCE)0733-9445(2004)130:12(1931)
[6]	DE BARROS CHODRAUI G M, SHIFFERAW Y, MALITE M, et al, On the stability of cold-formed steel angles undercompression[J]. Rev Esc Minas, 2007, 60(2): 355-363.
[7]	SILVESTRE N, DINIS P B, CAMOTIM D. Developments on the design of cold-formed steel angles[J]. Journal of Structural Engineering, 2013, 139(5): 680-694. doi: 10.1061/(ASCE)ST.1943-541X.0000670
[8]	ELLOBODY E, YOUNG B. Behavior of cold-formed steel plain angle columns[J]. Journal of Structural Engineering, 2005, 131(3): 457-466. doi: 10.1061/(ASCE)0733-9445(2005)131:3(457)
[9]	YOUNG B, ELLOBODY E. Buckling analysis of cold-formed steel lipped angle columns[J]. Journal of Structural Engineering, 2005, 131(10): 1570-1579. doi: 10.1061/(ASCE)0733-9445(2005)131:10(1570)
[10]	YOUNG B, CHEN J. Column tests of cold-formed steel non-symmetric lipped angle sections[J]. Journal of Constructional Steel Research, 2008, 64(7): 808-815.
[11]	SHIFFERAW Y, SCHAFER B W. Cold-formed steel lipped and plain angle columns with fixed ends[J]. Thin-Walled Structures, 2014, 80: 142-152. doi: 10.1016/j.tws.2014.03.001
[12]	ZHANG L L, WANG F Y, LIANG Y T, et al. Press-braked S690 high strength steel equal-leg angle and plain channel section stub columns: testing, numerical simulation and design[J]. Engineering Structures, 2019, 201: 109764.1-109764.12. doi: 10.1016/j.engstruct.2019.109764
[13]	姚行友,胡成立,聂泽谕,等. 冷弯不等肢复杂卷边角钢轴压试验及承载力设计方法[J]. 建筑钢结构进展,2024,26(2): 34-44. YAO Xingyou1,HU Chengli,NIE Zeyu,et al. Experiment and design method of cold-formed thin-walled steel unequal-legged complex-lipped angles under axial Compressio[J]. Progress in Steel Building Structures, 2024, 26(9): 34-44.
[14]	周宜松. 复杂卷边不等边角形截面冷弯薄壁型钢柱轴压稳定承载力研究[D]. 郑州:郑州大学,2020.
[15]	徐步洲. 双卷边不等肢角钢冷弯薄壁型钢柱轴压稳定承载力研究[D]. 西安:西安建筑科技大学,2017.
[16]	中华人民共和国建设部,中华人民共和国国家质量监督检验检疫总局. 冷弯薄壁型钢结构技术规范:GB 50018—2002[S]. 北京:中国计划出版社,2002.
[17]	国家市场监督管理总局,国家标准化管理委员会. 金属材料拉伸试验第1部分:室温试验方法:GB/T 228.1—2021[S]. 北京:中国标准出版社,2021.
[18]	姚行友,李元齐. 考虑畸变屈曲的冷弯薄壁型钢卷边槽形截面构件设计可靠度分析[J]. 建筑钢结构进展,2017,19(6): 65-74. YAO Xingyou, LI Yuanqi. Design Reliability Analysis of Cold-Formed Thin-Walled Steel Members with Lipped Channel Sections Considering Distortional Buckling[J]. Progress in Steel Building Structures, 2024, 19(6): 65-74.