Restrained Torsion Analysis of Box Girders with Corrugated Steel Webs Based on Reissner’s Principle
-
摘要:
为更加合理地分析波形钢腹板箱梁约束扭转效应,考虑波形钢腹板的褶皱效应推演了翘曲正应力和剪应力计算式,应用Reissner原理建立了波形钢腹板箱梁约束扭转控制微分方程,给出了不同于乌曼斯基第二理论的翘曲系数公式. 通过简支梁数值算例验证了所推导公式的正确性,并分析了腹板厚度和悬臂板宽度变化对箱梁横截面应力的影响. 研究结果表明:相对于乌曼斯基第二理论,基于Reissner原理计算的应力与有限元解吻合更好;按乌曼斯基第二理论与按Reissner原理计算的翘曲系数的比值可达到4.70;波形钢腹板主要承担剪应力,几乎不承担翘曲正应力,而顶底板既承担翘曲正应力也承担剪应力,应对顶底板予以重视,防止斜裂缝的产生;腹板厚度增大能减小翘曲正应力;随着悬臂板宽度的增大,当悬臂板宽度比大于0.10时,翘曲正应力减小,而当悬臂板宽度比大于0.30时,总剪应力几乎无变化.
Abstract:To analyze the restrained torsion effect of the box girder with corrugated steel webs (CSWs) more reasonably, the calculation formulas of warping normal stress and shear stress were deduced considering the accordion effect of CSWs. Moreover, the governing differential equation for analyzing the restrained torsion was established using Reissner’s principle, and the warping coefficient formula differing from Umanskii’s second theory was presented. The calculation formulas were then verified through a numerical simulation of a simply supported box girder, and the influence of web thickness and cantilever slab width on cross-section stresses of the box girder was analyzed. The results show that the stresses calculated based on Reissner’s principle are in better agreement with the finite element solution than those based on Umanskii’s second theory. The ratio of warping coefficient calculated by Umanskii’s second theory to that calculated by Reissner’s principle can reach 4.70. CSWs mainly bear shear stress and almost no warping normal stress, while top and bottom slabs bear both warping normal stress and shear stress. Therefore, more attention should be paid to top and bottom slabs to prevent oblique cracks. In addition, growing web thickness can reduce warping normal stress. With the increase of the cantilever slab width, warping normal stress decreases when the cantilever slab width ratio is above 0.10, while the total shear stress hardly changes when the cantilever slab width ratio is above 0.30.
-
波形钢腹板箱梁因为能减轻结构自重、减小预应力损失、避免腹板斜裂缝而成为具有广阔应用前景和吸引力的新型桥梁结构[1-2]. 不同于混凝土箱梁,波形钢腹板箱梁用波形钢腹板替换混凝土腹板后,箱梁的抗扭刚度被显著削弱[2-3],因此其扭转和畸变效应也更加突出.
国内外学者对波形钢腹板箱梁的约束扭转效应开展了不少研究工作. 杨丙文等[4]分析了单箱单室波形钢腹板箱梁的约束扭转效应,但在悬臂板自由端出现了剪应力,与实际不符. 马磊等[5]推导了单箱双室波形钢腹板箱梁扭转控制微分方程,但仅分析了翘曲正应力. 邓文琴等[6]推导了单箱三室波形钢腹板悬臂梁扭转控制微分方程,并给出了翘曲正应力和剪应力计算式. Nie等[7-8]等在波形钢腹板箱梁的试验研究结果中靠近顶底板局部区域内出现了可观的正应变. 张元海等[9]指出分析波形钢腹板箱梁的约束扭转时应考虑顶底板对波形钢腹板的约束作用. 以上都是基于乌曼斯基第二理论(下文简称乌-Ⅱ理论)建立约束扭转控制微分方程并开展的理论分析,而乌-Ⅱ理论应用了剪应力的两种算式,出现了两种不一致的纵向翘曲位移[10]. 张元海等[11]指出约束扭转剪应力应按微元体平衡条件或根据自由扭矩和二次扭矩计算,不能直接按胡克定律计算. 此外,Kato等[12]采用矩阵法提出了波形钢腹板箱梁扭转效应计算理论. 李运生等[13]采用能量变分法推导了波形钢腹板曲线梁弯扭控制微分方程,并采用伽辽金法求得其解析解. 徐勋等[14]在广义坐标法的基础上基于混合变分原理建立了薄壁箱梁约束扭转分析新理论,与Reissner原理的结果是统一的,均考虑了全部次生剪应力对中面剪切变形的影响,计算精度更高. 鲍永方等[15]对适用于刚性扭转的乌曼斯基理论和Reissner原理进行比较,指出Reissner原理的计算精度高于乌曼斯基理论.
为了更加合理地分析波形钢腹板箱梁约束扭转效应,本文将位移和应力同时作为变分参量,应用Reissner原理推导出约束扭转控制微分方程,给出了不同于乌-Ⅱ理论的翘曲系数. 同时考虑波形钢腹板的褶皱效应,合理计算截面几何特性,推演出翘曲正应力和剪应力的计算公式. 结合波形钢腹板简支箱梁数值算例,研究了应力分布特征及截面几何参数变化对应力的影响规律.
1. 波形钢腹板的等效
1.1 等效剪切模量及等效厚度
图1为波形钢腹板的形状简图. 图中:a、c、d分别为直板段长、斜板段投影长、斜板段长;hw、tw分别为波高和板厚.
因组成波形钢腹板箱梁的各板件材料不同,约束扭转分析时需通过等效原则将材料统一. 本文根据波形钢腹板的力学性能先将其等效为正交异性平钢板,再等效为混凝土腹板. 将波形钢腹板等效为正交异性平钢板后的等效剪切模量[16]为
Gse=(a+c)/(a+d)Gs=ηGs, (1) 式中:Gs为钢板剪切模量;
η=(a+c)/(a+d) .波形钢腹板在箱梁的受力中主要承受剪力,由式(1)采用剪切模量比可计算波形钢腹板等效为混凝土腹板的厚度,波形钢腹板的等效厚度为
twe=Gse/Gctw=ηnGtw, (2) 式中:Gc为混凝土剪切模量;
nG 为钢板与混凝土的剪切模量比,nG=Gs/Gc .1.2 纵向等效弹性模量
在一个波长的波形钢腹板两端施加一对等值、反向的轴向力,采用卡氏第二定理计算出弯矩和轴力下的轴向位移,根据位移等效原则可求得波形钢腹板等效为平钢板的纵向弹性模量为
Ese=κEs, (3) 式中:Es为钢板弹性模量;
κ 为系数,κ=a+c(a+c2/d)(tw/hw)2+(3a+d)(twhw)2. 根据
Ecetwe=Esetw 及式(2)、(3)求出波形钢腹板等效为混凝土腹板的纵向等效弹性模量为Ece=κnE/(ηnG)Ec = λEc, (4) 式中:Ec为混凝土弹性模量;
nE 为钢板与混凝土的弹性模量比,nE=Es/Ec ;系数λ=κnE/(ηnG) .由于tw/hw一般很小,κ、λ则更小,波形钢腹板表现出明显的褶皱效应.
1.3 基本假定
根据薄壁箱梁的约束扭转分析理论,结合波形钢腹板的力学性能,对波形钢腹板箱梁约束扭转分析时作以下基本假定:
1) 刚性周边假设,横截面周边不变形;
2) 横截面上的翘曲正应力和剪应力沿壁厚均匀分布;
3) 混凝土顶底板与波形钢腹板间不产生界面滑移,波形钢腹板不发生剪切屈曲;
4) 波形钢腹板箱梁处于线性弹性工作状态.
2. 约束扭转应力
2.1 约束扭转翘曲正应力
分析波形钢腹板箱梁的约束扭转效应时需考虑二次剪切变形的影响,引入广义翘曲位移
β′(x) . 壁厚中面上任一点的纵向翘曲位移u(x,s)和切向位移v(x,s)分别为u(x,s)=−β′(x)ˉω(s), v(x,s)=φ(x)ρ(s), 式中:
β′(x) 为广义翘曲位移,x为纵向坐标;ˉω(s) 为广义主扇性坐标,s为沿箱壁的曲线坐标,本文s原点取在竖向对称轴与顶板中面的交点处;ρ(s) 为扭心至各板壁厚中面的垂直距离;φ (x) 为横截面扭转角.由式(5)可得壁厚中面任一点的翘曲正应变和剪应变分别为
ε(x,s)=∂u/∂x=−β″ˉω, γ(x,s)=∂u/∂s+∂v/∂x=φ′ρ−β′ˉω′. 根据刚性周边假设得各板件横向正应变
εs=0 ,将式(6a)代入平面问题的物理方程求得翘曲正应力为σ={−Eoβ″ˉω, 顶底板,−λEoβ″ˉω,等效腹板, (7) 式中:
Eo=Ec/(1−μ2c) ,μc为混凝土泊松比.引起翘曲正应力的内力定义为扭转双力矩
Bˉω ,借助式(7),可得Bˉω=∫s1σˉωtds=−EoIˉωβ″, (8) 式中:
Iˉω 为波形钢腹板箱梁的广义主扇性惯性矩,Iˉω=∫bt+2bfˉω2ttds+∫bbˉω2tbds+2λ∫bwˉω2tweds, ∫s1ds、 ∫bt+2bfds、 ∫bbds、 ∫bwds 分别为在全截面(s1)、顶板宽度(bt)(含2倍悬臂板宽bf)、底板宽度(bb)、腹板宽度(bw)的积分;tt、tb分别为顶、底板的厚度.联立式(7)和式(8)可得翘曲正应力与扭转双力矩的关系式为
σ={Bˉωˉω/Iˉω, 顶底板,λBˉωˉω/Iˉω,等效腹板. (9) 2.2 约束扭转剪应力
约束扭转剪应力应按微元体平衡条件或根据自由扭矩和二次扭矩计算. 约束扭转总扭矩Mz分解为自由扭转扭矩M1和二次扭矩
M2 ,对应总剪应力τz包含自由扭转剪应力τ1和二次剪应力τ2 . 在闭口箱壁和悬臂板内同时存在τ1和τ2 ,但是在悬臂板内的τ1沿壁厚呈线性分布且在壁厚中面上为0,故悬臂板承担的自由扭转扭矩和提供的抗扭刚度较小.闭口箱壁中面上任一点的自由扭转剪应力为
τ1=Msc/(Ωt)=M1ψ/(Idt), (10) 式中:Msc为闭口箱壁承担的自由扭矩;
ψ=Ω/∮1/tds ,Ω为箱梁横截面闭口箱壁中面所围面积的2倍;Id=Ω2/∮1/tds+2bft3t/3 ,为全截面抗扭惯性矩;t为壁厚.根据波形钢腹板箱梁的抗扭特性需对抗扭惯性矩进行修正,考虑悬臂板的贡献,抗扭惯性矩[17]为
Id=Ω2/Ω2[bttt(1−ζ)+bbtb(1−ζ)+2bwtwe(1+ζ)][bttt(1−ζ)+bbtb(1−ζ)+2bwtwe(1+ζ)]+2bft3t3, 式中:ζ为修正系数,对矩形波形钢腹板箱梁,
ζ=0.400bw/bt−0.060⩾ .将式(9)代入各板件微元体平衡方程
{{\partial}} {{q_{\bar \omega }}}/{{\partial}} s + ({\partial} {\sigma}/{\partial} x)t = 0 ,根据截面翘曲位移连续性条件\displaystyle \oint {({{q_{\bar \omega }}}/t)} {\text{d}}s = 0 ,求得各板件的二次剪力流{q_{\bar \omega }} ,进而可得壁厚中面上任一点的二次剪应力{\tau _{2 }} 为{\tau _{2}} = {q_{\bar \omega }}/t = - ({M_{2}}{\bar S_{\bar \omega }})/({I_{\bar \omega }}t), (11) 式中:
{M_{2}} = {B' _{\bar \omega }} ;{\bar S_{\bar \omega }} 为广义扇性静面矩.各板壁厚中面上任一点的
{\bar S_{\bar \omega }} 为{\bar{S}}_{\bar{\omega }}=\left\{\begin{array}{l}{S}_{\bar{\omega }}, \quad 悬臂板,\\ {S}_{\bar{\omega }}-{\displaystyle \oint ({S}_{\bar{\omega }}/t\text{)d}s}\Big/{\displaystyle \oint \;1/t\text{d}s},\quad 闭口箱壁,\end{array}\right. 式中:
{S_{\bar \omega }} 为壁厚中面上任一点的广义主扇性静面矩,{S_{\bar \omega }} = {\bar \lambda} \displaystyle \int_0^s {\bar \omega t} {\text{d}}s ,对于顶底板,\bar \lambda = 1 ,对于等效腹板,\bar \lambda = \lambda ;由式(10)、(11)可得各板壁厚中面上任一点的总剪应力为
{\tau }_{{\text{z}}}=\left\{\begin{array}{l}-{M}_{2}{\bar{S}}_{\bar{\omega }}/({I}_{\bar{\omega }}t) , \quad \ 悬臂板 , \\ {M}_{1}\psi/({I}_{{\rm{d}}}t)-{M}_{2}{\bar{S}}_{\bar{\omega }}/({I}_{\bar{\omega }}t) , \quad 闭口箱壁. \end{array}\right. (12) 3. 基于Reissner原理建立控制微分方程
基于Reissner原理建立波形钢腹板箱梁的约束扭转控制微分方程时,将位移和应力同时作为了变分参量,建立的能量泛函[18]为
\begin{split} &{\varPi_{\rm{R}}}({\boldsymbol{\sigma}} ,{\boldsymbol{u}}) = \int_V {\left({{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} ({\boldsymbol{u}}) - {\boldsymbol{V}}({\boldsymbol{\sigma}} ) - {\boldsymbol{Fu}} \right)} {\text{d}}V -\\ &\quad \int_l {{\boldsymbol{Pu}}} {\text{d}}x, \end{split} (13) 式中:
{\boldsymbol{\sigma}} 为应力向量;{\boldsymbol{u}} 为位移向量;{\boldsymbol{\varepsilon}} ({\boldsymbol{u}}) 为应变向量;{\boldsymbol{V}}({\boldsymbol{\sigma}} ) 为余能密度;{\boldsymbol{F}} 为体力向量,本文中为0;{\boldsymbol{P}} 为面力向量;V为体积;l为计算跨径.能量泛函中的应变能为
\int_V {{{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} ({\boldsymbol{u}})} {\text{d}}V = \int_V {({\sigma }\varepsilon + {\tau _{\text{z}}}\gamma )} {\text{d}}V. (14) 将式(6)、(9)和式(12)代入式(14),得
\int_V {{{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} } ({\boldsymbol{u}}){\text{d}}V = \int_0^l {( - {B_{\bar \omega }}\beta '' + {M_{\text{z}}}\varphi ' - {M_{2}}\beta '{\text{)d}}x} . (15) 余能函数为
\int_V {{\boldsymbol{V}}({\boldsymbol{\sigma}} )} {\text{d}}V = \int_V {\left(\dfrac{{\sigma ^2}}{{2E_{\rm{o}}}} + \dfrac{{\tau _{\text{z}}^2}}{{2{G_{\rm{c}}}}}\right)} {\text{d}}V. (16) 将式(9)和式(12)代入式(16),得
\int_V {{\boldsymbol{V}}({\boldsymbol{\sigma}})} {\text{d}}V = \int_0^l \left({\dfrac{{{B_{\bar \omega }^2}}}{{2E_{\rm{o}}{I_{\bar \omega }}}} + \dfrac{{M_{\rm{1}}^2}}{{2{G_{\rm{c}}}{I_{\rm{d}}}}} + } \dfrac{{\chi M_{2 }^2}}{{2{G_{\rm{c}}}{I_{\rm{d}}}}}\right){\text{d}}x , (17) 式中:
\chi = \dfrac{{{I_{\rm{d}}}}}{{I_{\bar \omega }^2}}\displaystyle \int_{s_{1}} {\left(\bar S_{\bar \omega }^2/t\right)} {\text{d}}s .当箱梁上仅作用分布扭矩mt时,外力势能为
- \int_0^l {{\boldsymbol{Pu{}}\text{d}}x} = - \int_0^l {{m_{\rm{t}}}\varphi {\text{d}}x} . (18) 将式(15)、(17)和式(18)代入式(13),并对能量泛函数
{\varPi _R} 进行一阶变分运算,可得\begin{split} &\delta {\varPi_{\rm{R}}} = \int_0^l \left[{\left( - \beta '' - \dfrac{{{B_{\bar \omega }}}}{{E_{\rm{o}}{I_{\bar \omega }}}}\right)\delta {B_{\bar \omega }} + \left(\varphi ' - \dfrac{{{M_{\rm{s}}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}\right)\delta {M_{1}}} + \right.\\ &\quad \left.\left(\varphi ' - \beta ' - \dfrac{{\chi {M_{2 }}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}\right)\delta {M_{2}}{\text{ + }}\left({B'_{\bar \omega }} - {M_{2}}\right)\delta \beta' - \right. \\ &\quad \Bigg.({M'_1}+ {M'_{2}} + {m_{\rm{t}}})\delta \varphi \Bigg]{\text{d}}x . \\[-12pt] \end{split} (19) 由
\delta {\varPi _R} = 0 可得能量泛函的约束条件:1) 力与位移的关系条件
\left\{\begin{array} {l} - \beta '' = \dfrac{{{B_{\bar \omega }}}}{{E_{\rm{o}}{I_{\bar \omega }}}}, \\ \varphi ' = \dfrac{{{M_{\rm{1}}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}, \\ \varphi ' - \beta ' = \dfrac{{\chi {M_{2}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}. \end{array} \right. (20) 2) 力的平衡条件
\left\{ \begin{array}{l} {B'_{\bar \omega }} = {M_{2}}, \\ {M'_1} + {M'_{2}} = - {m_{\rm{t}}} . \end{array} \right. (21) 借助式(20)和式(21),并忽略mt的高阶微分,可得波形钢腹板箱梁约束扭转控制微分方程为
\varphi^ {(4)} - k_{\text{R}}^2\varphi '' = {\xi _{\text{R}}}{m_{\rm{t}}}/(E_{\rm{o}}{I_{\bar \omega }}), (22) 式中:
{k_{\text{R}}} = \sqrt {({\xi _{\text{R}}}{G_{\rm{c}}}{I_{\rm{d}}})/(E_{\rm{o}}{I_{\bar \omega }})} ;{\xi _{\text{R}}} 为按Reissner原理计算的翘曲系数,{\xi _{\text{R}}} = 1/(1 + \chi ) .不同于Reissner原理,采用乌Ⅱ-理论推导波形钢腹板箱梁约束扭转控制微分方程时,首先要确定广义翘曲位移
\beta ' 与扭率\varphi ' 之间的关系. 由剪切胡克定律及式(6b)可得总剪力流qz为q_{\text{z}} = {\tau _{\text{z}}}t = {G_{\text{c}}}t(\varphi '\rho - \beta '\bar \omega '). (23) 将式(23)代入静力平衡方程
{M_{\text{z}}} =\displaystyle \int_{s_1} {{q_{\textit{z}}}\rho } {\text{d}}s ,可得\varphi ' - {\xi _{\text{w}}}\beta = {M_{\text{z}}}/({G_{\rm{c}}}{I_\rho }). (24) 式中:
{\xi _{\text{w}}} 为按乌-Ⅱ理论计算的翘曲系数,{\xi _{\text{w}}} = 1 - {I_{\text{d}}}/{I_\rho } ;{I_\rho } 为极惯性矩,{I_\rho } =\displaystyle \int_{s_1} {{\rho ^2}t{\text{d}}s} .由微元体平衡方程
\partial{{q_{\bar \omega }}} / {\partial}s + ({\partial} {\sigma}/{\partial} x)t = 0 及截面翘曲位移连续性条件\displaystyle \oint {({{q_{\bar \omega }}}/t)} {\text{d}}s = 0 ,可求得二次剪力流{{q_{\bar \omega }}} . 将其代入静力平衡方程{M_{2}} = \displaystyle \int _{s_1}{{q_{\bar \omega }}}\rho {\text{d}}s ,得{M_{2}} = - E_{\rm{o}}{I_{\bar \omega }}\beta ''' .由功能原理及自由扭转剪力流qs满足的静力平衡方程
{M_1} = \displaystyle \oint {{q_{\rm{s}}}\rho {\text{d}}s} 得{M_1} = {G_{\rm{c}}}{I_{\text{d}}}\varphi ' .M1和M2叠加得Mz,即
{M_{\text{z}}} = {G_{\rm{c}}}{I_{\text{d}}}\varphi ' - E_{\rm{o}}{I_{\bar \omega }}\beta ''' . (25) 将式(24)及
{\text{d}}{M_{\text{z}}} = - {m_{\rm{t}}}{\text{d}}\textit{z} 代入式(25),忽略mt的高阶微分,可得约束扭转控制微分方程为\varphi^ {(4)} - k_{\text{w}}^2\varphi '' = {\xi _{\text{w}}}{m_{\rm{t}}}/(E_{\rm{o}}{I_{\bar \omega }}) , (26) 式中:
{k_{\text{w}}} = \sqrt {({\xi _{\text{w}}}{G_{\rm{c}}}{I_{\rm{d}}})/(E_{\rm{o}}{I_{\bar \omega }})} .两种方法推导的约束扭转控制微分方程表达式相同,且均可采用初参数法求解[19],但翘曲系数公式不同. 采用乌-Ⅱ理论推导约束扭转控制微分方程时应用了由胡克定律和微元体平衡条件得到的剪力流的两种算式及两种不一致的纵向翘曲位移[10],且未考虑全部次生剪应力对中面剪切变形的影响. 应用Reissner原理建立扭转控制微分方程时将位移和应力同时作为变分参量,并考虑了全部次生剪应力对中面剪切变形的影响,通过建立能量泛函并对其变分,得出力与位移的关系条件及力的平衡条件,与刚性周边假设下的力与位移的关系条件及力的平衡条件是相符的,因此该理论更准确.
4. 算例分析
4.1 算例概况
算例采用计算跨径为40 m的简支波形钢腹板箱梁,梁端设置1.5 m厚混凝土横隔板,跨中作用集中扭矩荷载
\bar M = 1\;000 \;{\text{kN}}{\text{•}} {\text{m}} ,箱梁横截面尺寸见图2. 图2中还示出了计算点Ⅰ ~ Ⅵ 的位置. 波形钢腹板采用1600型,其几何尺寸:a=d=430 mm,c=370 mm,hw=220 mm,tw=16 mm. 混凝土的材料特性:Ec=34.5 GPa,Gc=13.8 GPa,μc=0.2;钢板的材料特性:Es=206 GPa,Gs=79 GPa,μs=0.31. 求得扭心至顶板中面的距离为0.595 0 m,极惯性矩Iρ=4.8236 m4,抗扭惯性矩{I_{\rm{d}}} = 4.059\;8\;{{\text{m}}^4} ,广义主扇性惯性矩{I_{\bar \omega }} =0.7555 m6,按乌-Ⅱ理论和本文方法计算的翘曲系数分别为ξw=0.1583,ξR=0.069 3.4.2 理论验证
图3为按本文方法和乌-Ⅱ理论计算的扭转内力分布曲线. 由图3可以看出:按不同方法计算的双力矩和二次扭矩在跨中截面有明显差异,按乌-Ⅱ理论计算的双力矩和二次扭矩分别比本文方法计算的约增大了51%和128%. 由于两种方法计算的翘曲系数差异较大导致扭转内力的差异较大,也说明不同方法计算的翘曲系数对扭转内力有显著影响.
为验证本文方法合理性,按本文方法、乌-Ⅱ理论计算横截面翘曲正应力和总剪应力,并与有限元解对比. 有限元解是采用ANSYS软件建立空间有限元模型求得,有限元模型如图4所示. 混凝土板和波形钢腹板分别采用SOLID45实体单元和SHELL63壳单元模拟,由于分析时假设波形钢腹板与混凝土顶底板之间无界面滑移[20],故在模型中壳单元和实体单元间用主从约束法作刚性连接处理;每侧梁端底板上设置2个支座,左侧梁端约束为3个坐标轴方向的平动位移Ux 、Uy和Uz,右侧梁端约束为横向和竖向平动位移Uy和Uz;集中扭矩荷载等效成刚性扭转荷载的力流形式,均匀施加于跨中截面闭合箱壁各单元的节点上,顶底板施加反向水平荷载
\bar M/(2h) ,其中,h为梁高,波形钢腹板施加反向竖直荷载\bar M/(2{b_{\text{t}}}) . 刚性扭转荷载下的箱梁有限元模型在纵向仅产生扭转翘曲正应力.为避免有限元模型中集中扭矩施加位置的局部应力影响,根据扭转内力分布及有限元模型的应力分布,选取距跨中1.6 m的左截面正面(外法线方向与x轴一致的横截面为正面)为分析截面. 以该截面上的Ⅰ、Ⅱ、Ⅴ 3个计算点进行翘曲正应力对比,Ⅲ、Ⅳ、Ⅵ 3个计算点进行总剪应力对比,各计算点的应力结果见表1. 表1中的相对误差为本文方法解或乌-Ⅱ理论解与有限元解的差除以有限元解求得. 由表1可以看出,本文方法计算的应力与有限元解吻合更好. 但是由于刚性周边等假设,理论解与有限元解之间还是存在一定的误差.
表 1 距跨中1.6 m左截面应力比较Table 1. Comparison of stresses at left section 1.6 m from mid-span名称 计算点 有限元
解/kPa本文方
法/kPa乌-Ⅱ
理论/kPa相对误差/% 本文
方法乌-Ⅱ
理论{\sigma} Ⅰ 54.06 57.35 63.61 6.09 17.67 Ⅱ −27.83 −25.42 −28.18 −8.66 1.28 Ⅴ −75.83 −81.24 −90.10 7.13 18.82 {\tau _{\text{z}}} Ⅲ 82.37 84.59 77.37 2.70 −6.07 Ⅳ 1670.20 1480.32 1354.03 −11.37 −18.93 Ⅵ 154.55 154.50 176.98 −0.03 15.51 4.3 应力分布
为了解波形钢腹板箱梁的应力在横截面上的分布规律,图5为跨中左截面翘曲正应力、二次剪应力、总剪应力的分布图,其中箭头方向为剪应力的方向. 由图5可以看出:点Ⅰ和点Ⅴ的翘曲正应力绝对值较大,而波形钢腹板几乎不产生翘曲正应力;点Ⅴ和点Ⅵ的二次剪应力绝对值较大,腹板的二次剪应力与自由扭转剪应力叠加后使得总剪应力减小且沿梁高呈均匀分布,而底板的二次剪应力与自由扭转剪应力叠加后使得总剪应力增大,悬臂板壁厚中面上不产生自由扭转剪应力,但存在二次剪应力.
4.4 参数影响分析
为进一步探索截面几何参数变化对波形钢腹板箱梁约束扭转效应的影响,选取不同高宽比的梁为分析对象. 梁高h分别取值为2.15、4.30 m,即高宽比h/bt分别为0.5、1.0. 根据图2分别改变tw、bf,腹板厚度tw以2 mm为步长从8 mm增大至22 mm,悬臂板单侧宽度bf从0开始以0.215 m为步长增大至4.300 m,即悬臂板宽度与闭合箱壁宽度的比值(简称悬臂板宽度比,bf /bt)从0开始以0.05为步长增大至1.00.
为比较乌-Ⅱ理论与本文方法计算的翘曲系数的差异,图6绘制了乌-Ⅱ理论与本文方法计算的翘曲系数比值
{\xi _{\text{w}}}/{\xi _{\rm{R}}} 随tw、bf /bt的变化曲线. 由图6可以看出:翘曲系数比值受腹板厚度和悬臂板宽度比的影响显著,且随腹板厚度和悬臂板宽度比的增大呈增长趋势,不同高宽比梁的翘曲系数比值增长趋势有差异;当tw>10 mm和bf /bt>0.20时,翘曲系数比值大于1;翘曲系数比值偏离1越远,则两种方法的计算结果差异越大;当腹板厚度为22 mm时,高宽比较小的梁的翘曲系数比值达到了4.70;当悬臂板宽度比为1.00时,高宽比较大的梁的翘曲系数比值达到了3.10.由于点Ⅰ和点Ⅴ的翘曲正应力绝对值较大,图7绘制了跨中左截面点Ⅰ和点Ⅴ翘曲正应力随tw、bf /bt的变化曲线. 由图7可以看出:点Ⅰ和点Ⅴ的翘曲正应力绝对值均随腹板厚度的增大而减少,随悬臂板宽度比的增大先增大后减少,bf /bt=0.10时开始减小;当bf /bt>0.10时,悬臂板宽度增大可以减少翘曲正应力.
图8绘制了跨中左截面点 Ⅳ 和点 Ⅵ 的总剪应力随tw、bf/bt的变化曲线. 由图8可以看出:高宽比较小的梁的总剪应力明显大于高宽比较大的梁;对高宽比较大的梁,当tw=8 mm和bf/bt=0.30时,点Ⅵ总剪应力分别达到点 Ⅳ总剪应力的34.4%、30.2%,而混凝土与钢板的剪切模量比为0.17,底板总剪应力显然不容忽视,经分析,悬臂板根部的总剪应力同样不容忽视. 因此,在波形钢腹板箱梁设计时,应对顶底板进行精细化设计,防止斜裂缝的产生;随着悬臂板宽度的增大,当bf /bt>0.30时,总剪应力几乎不变;对高宽比较小的梁的总剪应力随腹板厚度增大而减小,而对高宽比较大的梁,腹板厚度对总剪应力的影响不显著.
5. 结 论
1) Reissner原理将位移和应力同时作为变分参量,考虑约束扭转全部次生剪应力对中面剪切变形的影响来推导波形钢腹板箱梁约束扭转控制微分方程,理论更准确. 本文推导的约束扭转控制微分方程表达式与乌-Ⅱ理论的相同,但翘曲系数公式不同. 数值算例的结果表明,相对于乌-Ⅱ理论,本文方法的计算结果与有限元解吻合更好,验证了本文方法的合理性.
2) 按乌-Ⅱ理论与本文方法计算的翘曲系数的比值表征了两种理论计算结果的差异,翘曲系数比值偏离1越远,则计算结果的差异越大. 两种方法计算的翘曲系数受腹板厚度和悬臂板宽度的影响显著,翘曲系数比值最大可达到4.70.
3) 波形钢腹板箱梁发生约束扭转时,波形钢腹板几乎不承担翘曲正应力,主要承担剪应力,混凝土顶底板既承担翘曲正应力又承担剪应力. 底板总剪应力可达到腹板总剪应力的34.4%,而混凝土相对钢板的抗剪性能较差,应对顶底板予以重视,并对其精细化设计,防止主拉应力超限引起斜裂缝.
4) 波形钢腹板厚度增大可以减小约束扭转翘曲正应力;随着悬臂板宽度的增大,悬臂板宽度比大于0.10时,翘曲正应力减小,而当悬臂板宽度比大于0.30时,总剪应力几乎无变化.
-
表 1 距跨中1.6 m左截面应力比较
Table 1. Comparison of stresses at left section 1.6 m from mid-span
名称 计算点 有限元
解/kPa本文方
法/kPa乌-Ⅱ
理论/kPa相对误差/% 本文
方法乌-Ⅱ
理论{\sigma} Ⅰ 54.06 57.35 63.61 6.09 17.67 Ⅱ −27.83 −25.42 −28.18 −8.66 1.28 Ⅴ −75.83 −81.24 −90.10 7.13 18.82 {\tau _{\text{z}}} Ⅲ 82.37 84.59 77.37 2.70 −6.07 Ⅳ 1670.20 1480.32 1354.03 −11.37 −18.93 Ⅵ 154.55 154.50 176.98 −0.03 15.51 -
[1] 李立峰,周聪,王连华,等. 基于纽玛克法的变截面波形钢腹板组合箱梁畸变效应分析[J]. 中国公路学报,2018,31(6): 217-226.LI Lifeng, ZHOU Cong, WANG Lianhua, et al. Analysis on distortion effect of non-prismatic composite box girders with corrugated steel webs based on Newmark method[J]. China Journal of Highway and Transport, 2018, 31(6): 217-226. [2] ZHU Y B, WAN S, SHEN K J, et al. Theoretical study on the nonlinear performance of single-box multi-cell composite box-girder with corrugated steel webs under pure torsion[J]. Journal of Constructional Steel Research, 2021, 178: 106487.1-106487.20. [3] JIANG R J, AUF T K, XIAO Y F. Prestressed concrete girder bridges with corrugated steel webs: review[J]. Journal of Structural Engineering, 2015, 141(2): 04014108.1-04014108.9. [4] 杨丙文,黎雅乐,万水,等. 波形钢腹板箱梁的扭转应力分析[J]. 华南理工大学学报(自然科学版),2012,40(2): 19-22,34.YANG Bingwen, LI Yale, WAN Shui, et al. Stress analysis of box girder with corrugate steel webs under torsion[J]. Journal of South China University of Technology (Natural Science Edition), 2012, 40(2): 19-22,34. [5] 马磊,万水,蒋正文. 单箱双室波形钢腹板箱梁扭转与畸变性能研究[J]. 中国公路学报,2016,29(10): 77-85.MA Lei, WAN Shui, JIANG Zhengwen. Research on torsion and distortion performance of single box double-cell girder with corrugated steel webs[J]. China Journal of Highway and Transport, 2016, 29(10): 77-85. [6] 邓文琴,毛泽亮,刘朵,等. 单箱三室波形钢腹板悬臂梁扭转与畸变分析及试验研究[J]. 建筑结构学报,2020,41(2): 173-181.DENG Wenqin, MAO Zeliang, LIU Duo, et al. Analysis and experimental study on torsion and distortion of single box three-cell cantilever girder with corrugated steel webs[J]. Journal of Building Structures, 2020, 41(2): 173-181. [7] NIE J G, ZHU Y J, TAO M X, et al. Optimized prestressed continuous composite girder bridges with corrugated steel webs[J]. Journal of Bridge Engineering, 2017, 22(2): 04016121.1-04016121.17. [8] HE J, LIU Y Q, WANG S H, et al. Experimental study on structural performance of prefabricated composite box girder with corrugated webs and steel tube slab[J]. Journal of Bridge Engineering, 2019, 24(6): 04019047.1-04019047.17. [9] 张元海,孙成成. 波形钢腹板箱梁考虑腹板局部纵向刚度影响的扭转效应分析[J]. 东南大学学报(自然科学版),2021,51(2): 195-201.ZHANG Yuanhai, SUN Chengcheng. Analysis on the torsion effect of box girders with corrugated steel webs considering local longitudinal stiffness of webs[J]. Journal of Southeast University (Natural Science Edition), 2021, 51(2): 195-201. [10] 强士中,李乔. 关于闭口薄壁杆件约束扭转的周边不变形理论[J]. 桥梁建设,1985(1): 63-75. [11] 张元海,刘泽翔,王晨光. 带悬臂板薄壁箱梁极惯性矩的合理计算方法[J]. 铁道学报,2019,41(8): 94-99.ZHANG Yuanhai, LIU Zexiang, WANG Chenguang. Rational method for computing polar inertia moment of thin-walled box girders with cantilever slab[J]. Journal of the China Railway Society, 2019, 41(8): 94-99. [12] KATO H, NISHIMURA N. Practical analysis methods for continuous girder and cable stayed bridges composed of beams with corrugated steel webs[J]. Structural Engineering/Earthquake Engineering, 2004, 21(2): 207-222. doi: 10.2208/jsceseee.21.207s [13] 李运生,陈留剑,刘蓓,等. 波形钢腹板曲线结合梁弯扭效应的解析解推导及参数分析[J]. 铁道学报,2019,41(1): 101-108.LI Yunsheng, CHEN Liujian, LIU Bei, et al. Analytical solution derivation and parametrical analysis of bending-torsional effects of curved composite beam with corrugated steel webs[J]. Journal of the China Railway Society, 2019, 41(1): 101-108. [14] 徐勋,叶华文,强士中. 带悬臂板薄壁箱梁的扭转和畸变分析[J]. 铁道学报,2015,37(10): 83-91.XU Xun, YE Huawen, QIANG Shizhong. Torsion and distortion analysis of thin-walled box girder with cantilever flanges[J]. Journal of the China Railway Society, 2015, 37(10): 83-91. [15] 鲍永方,黄文彬. 矩形箱梁约束扭转理论的分析与比较[J]. 工程力学,1997,14(3): 132-137.BAO Yongfang, HUANG Wenbin. Analysis and comparison of restrained torsion theories for rectangular box-shaped beam[J]. Engineering Mechanics, 1997, 14(3): 132-137. [16] JOHNSON R P, CAFOLLA J. Corrugated webs in plate girders for bridges[J]. Proceedings of the Institution of Civil Engineers, Structures and Buildings, 1997, 122(2): 157-164. doi: 10.1680/istbu.1997.29305 [17] 李宏江. 波形钢腹板预应力混凝土组合箱梁扭转与畸变研究进展[J]. 建筑结构学报,2017,38(7): 59-67.LI Hongjiang. Review on torsion and distortion in prestressed concrete box girders with corrugated steel webs[J]. Journal of Building Structures, 2017, 38(7): 59-67. [18] 杜国华, 毛昌时, 司徒妙龄. 桥梁结构分析[M]. 上海: 同济大学出版社, 1994: 62-79. [19] 郭金琼, 房贞政, 郑振. 箱形梁设计理论[M]. 2版. 北京: 人民交通出版社, 2008: 59-71. [20] ZHOU C, LI L F, WANG L H. Improved softened membrane model for prestressed composite box girders with corrugated steel webs under pure torsion[J]. Journal of Constructional Steel Research, 2019, 153: 372-384. doi: 10.1016/j.jcsr.2018.10.023 -