金属拱型波纹屋盖结构的简化设计方法

金属拱型波纹屋盖结构的简化设计方法张勇1刘锡良2王元清1石永久1清华大学土木工程系北京1000842.天津大学土木工程系天津300072)摘要：金属拱型波纹屋盖结构在力学性能上具有很强的几何非线性，结构设计必须考虑几何非线性的影响。本文通过对这种结构进行线性和几何非线性对比分析，找到了两种计算结果之间的对应规律，从而提出了一套基于线性分析的简化设计方法。本文运用这套方法对若干实际工程进行了试设计，试设计结果证明了简化设计方法的合理性和可操作性。关键词：金属拱型波纹屋盖；简化设计方法；几何非线性分析；线性分析中文分类号:TU318.1文献标识码：AAsimplifieddesignmethodforarchedcorrugatedmetalroofZhangYong1LiuXiliang2WangYuanqing1ShiYongjiu1DepartmentofCivilEngineering,TsinghuaUniversity,Beijing100084,China2.DepartmentofCivilEngineering,TianjinUniversity,Tianjin300072,China)Abstract:Archedcorrugatedmetalroofexhibitsstronggeometricalnonlinearitywhenloaded,whichcannotbeignoredinstructuraldesign.Throughthecomparisonanalysisoflinearanalysisandnonlinearanalysis,therulesoftherelationshipbetweentheresultsofthetwomethodswasfound,andasimplifieddesignmethodbasedonlinearanalysisforthiskindofstructurewasproposedbytheauthorsofthispaper.Thereasonabilityandfeasibilityofthismethodwasverifiedthroughthetrialdesignofseveralpracticalprojects.Keywords:archedcorrugatedmetalroof;simplifieddesignmethod;geometricalnonlinearAnalysis;LinearAnalysis1引言金属拱型波纹屋盖结构的刚度较小，在荷载作用下具有显著几何软化特征，因而其稳定问题实质上是一个弹塑性分析问题，其临界荷载是结构几何非线性和物理非线性同时作用的结果[2]。而结构几何、物理双非线性分析方法以结构非线性分析理论为基础，具有一定的理论难度，难以作为通用的结构设计方法而为广大工程技术人员所接受。目前，在世界范围内金属拱型波纹屋盖结构尚未形成普遍接受的设计计算方法。在我国这种结构的设计通常由施工单位来完成，施工单位确定结构的设计参数后向专业设计单位提供屋盖对下部结构的支承反力，设计单位只进行下部支承结构的设计。事实上绝大多数施工单位对这种结构的力学性能也缺乏清楚的认识，结构参数往往根据主观经验来确定，其结果自然难以保证结构设计“安全、合理”，故此这种结构的恶性工程事故在我国时有发生[1],[6]。探索一种简便、适用的简化设计方法，对于这种结构的健康发展具有非常紧迫的现实意义。本文作者在对金属拱型波纹屋盖结构进行系统的理论分析与试验研究的基础上，提出了一套设计方法，只需对结构进行线弹性分析，不仅易于理解，而且简便易行。经过大量工程试设计后，该方法已被编制中的国家标准《金属拱型波纹屋盖结构技术规程》（目前处于报批阶段）采用。本文将详细介绍这套简化设计方法的推导过程，明确采用该方法时的设计步骤。本文给出的部分实际工程试设计结果也在一定程度上验证了该方法的合理性和可操作性。2金属拱型波纹屋盖结构受力机理分析金属拱型波纹屋盖结构所用材料与一般的压型钢板相同，结构成型工艺也与压型钢板相似，但两者在力学性能上却差别很大：作为一种板架合一的屋盖结构，金属拱型波纹屋盖目前在我国最大跨度已做到38m；而压型钢板在建筑中只能配合梁檩使用，作为围护结构，跨度一般在3m以内。这种现象的根本原因是两者的受力模式不同：金属拱型波纹屋盖结构呈拱型（或筒壳）结构受力机理，内力分布较均匀，是一种靠合理的几何外形来提供较大承载能力或跨越能力的结构。压型钢板的受力机理类似于梁，而实腹梁是一种最不合理的结构形式。另外金属拱型波纹屋盖结构由于成型的需要，结构板件上均匀地压有许多横向小波纹，这些小波纹虽然在一定程度上降低了结构在跨度方向的刚度，但却大大增强了板件的局部稳定性，突出了结构的壳体受力特征，更为施工组装过程中的结构构件提供必要的抗扭刚度。实际工程中金属拱型波纹屋盖结构所受荷载类型比较单一，设计时要考虑的荷载主要有：自重、保温荷载、吊顶荷载、雪荷载（活荷载）、积灰荷载、风荷载等。文献[2]根据这几种荷载的特点将它们分成自重类荷载、类雪荷载、半跨三角形分布荷载以及类风荷载等四类，并运用非线性有限元法对在这四类荷载作用下的不同支承形式、不同矢跨比屋盖结构进行了几何及物理双非线性分析，从中明确了这种结构的受力机理和破坏模态。对结构受力机理的把握是探讨结构简化设计方法的基础。图1轴力为压力时的弯矩-荷载关系曲线Fig.1Moment-loadcurvewhenaxialforceiscompressiveforce图2轴力为拉力时的弯矩-荷载关系曲线Fig.2Moment-loadcurvewhenaxialforceistensileforce由于金属拱型波纹屋盖结构的板件较薄，因此截面刚度小，在荷载尤其是非对称荷载作用下结构变形较大。拱型结构的轴力较大，而轴力对于变形的作用也就是常说的P-图1轴力为压力时的弯矩-荷载关系曲线Fig.1Moment-loadcurvewhenaxialforceiscompressiveforce图2轴力为拉力时的弯矩-荷载关系曲线Fig.2Moment-loadcurvewhenaxialforceistensileforce3金属拱型波纹屋盖结构简化设计公式从非线性分析的角度来考察，结构的稳定性问题和强度问题是联系在一起的。结构失稳前的大变形导致弯矩内力迅速增加，因此失稳前结构一般都要进入弹塑性变形阶段。文献[2]的分析证实，作为一种薄壁钢结构，一旦金属拱型波纹屋盖结构上出现塑性区，结构刚度将受到很大削弱，因而变形增加更快，因变形而引起的附加弯矩也相应增大，塑性区将迅速发展，最终导致整个结构的破坏。因此对这种结构要提出一种基于线性分析的简化设计方法，必须解决的是如何通过线性弯矩内力求出非线性弯矩内力，从而将稳定问题转化为强度问题并与强度问题一起计算。对于薄壁结构的强度问题有一点已是共识：控制结构于弹性工作阶段，即采用边缘屈服准则作为薄壁结构设计时的强度评判标准。在对这种结构进行非线性分析过程中发现结构二阶弯矩内力相对于一阶弯矩内力的变化幅度主要与荷载类型及实际荷载和临界荷载之间的比率有关，这使得很自然想到参照压弯杆件的相关公式来表达这种弯矩放大作用。因此本文采用了下面一组公式作为这种结构的简化设计公式：（1）式中Aeq、Weq分别为考虑了波纹效应的等效截面面积及等效截面抵抗矩；、分别为计算工况下截面的一阶轴力和二阶弯矩。可分别通过下列公式求得： （2） （3）其中：计算工况中第种设计荷载作用下的单位宽度截面的一阶轴力；为计算工况中第种设计荷载作用下的单位宽度截面的一阶弯矩；为与计算工况中第类设计荷载相对应的弯矩放大系数。前面已经提到对于封闭的拱型结构，在风荷载作用下弯矩内力与外荷载的关系曲线是上凹的，即一阶弯矩比二阶弯矩（即实际弯矩）大，因此为了简便，不再缩小风荷载作用下的一阶弯矩，此时弯矩放大系数，其他设计荷载作用下的弯矩放大系数可通过下式求得： （4）其中的为弯矩调整系数，与支承条件、荷载类型及结构的矢跨比有关。为计算工况中第类荷载设计值；为拱型屋盖在计算工况中第类荷载作用下的弹性临界荷载：（5）式中的为考虑了波纹效应的结构等效截面惯性矩；为临界荷载系数，本文作者在文献[4]中推导出（5）式，并给出了不同条件下的值。如前所述，这种结构的板件上分布着很多轧制小波纹，小波纹对板件的局部刚度和结构的整体刚度都有很大的影响，文献[3]分析了这种影响，并按等刚度原则给出了等效截面刚度、、的计算方法。到目前为止以上公式中只有一个参数尚未确定，这就是弯矩调整系数。根据（2）式及（3）式，可通过下式求出：（6）这里的即为二阶分析得到的弯矩。为了求解上式，必须对结构进行二阶与一阶计算结果的对比分析。首先对结构进行一阶分析，求得各种荷载水平下的一阶最大弯矩，然后进行二阶弹性分析，求得作用下结构的二阶最大弯矩。由于值是通过结构的最大弯矩求出的，因此严格意义上说公式（4）中的应为最大弯矩放大系数，分析表明采用最大弯矩放大系数来放大结构各截面上的一阶弯矩得到的结果与采用各截面相应弯矩放大系数得到的结果差别不大，且偏于安全。本文大量的参数分析表明，在支承条件一定的情况下的取值仅与荷载类型、的比值及结构矢跨比有关。接下来本文将求出各种荷载类型下的弯矩调整系数。为了节约篇幅，本文只给出部分结果的分析过程。3.1铰支结构在自重类荷载作用下的弯矩调整系数图3为矢跨比为0.25的铰支结构稳定承载力系数随变化的关系曲线。图中包含了18条曲线，每条曲线代表一特定跨度及板厚结构的计算结果。基本重合的几条曲线代表的结构跨度相同而板厚不同。从此图中可以看出在自重类荷载作用下随着跨度的变化有所改变，但变化的幅度不大（变化幅度在0.0025之间，约占最小取值的0.65%）。当结构跨度一定、板厚改变时几乎不发生变化。本文认为随结构跨度也不应改变，这是这种结构的内部规律。计算时之所以出现值随跨度改变而波动的现象，主要和计算模型引起的系统偏差及公式（6）中各参数存在误差有关。图4为当矢跨比为0.4时的分析结果。此图中同样包含18条曲线，每条曲线代表一跨度及板厚确定的结构，可以看出这18条曲线几乎完全重合在一起。这证明了与结构跨度及所用板厚无关。当结构矢跨比分别取0.15、0.2、0.30、0.35及0.45时也可得出类似的结果。总的说来结构的矢跨比越小，随结构跨度的改变波动的幅度越大。但各种矢跨比下，波动的绝对幅度都不大。图3矢跨比0.25Fig.3Therisetospanratiois0.25图4矢跨比0.40Fig.4Therisetospanratiois0.40图5自重类荷载下的弯矩调整系数Fig.5图3矢跨比0.25Fig.3Therisetospanratiois0.25图4矢跨比0.40Fig.4Therisetospanratiois0.40图5自重类荷载下的弯矩调整系数Fig.5Momentadjustingfactorunderself-weighttypeload图6类雪荷载下的弯矩调整系数Fig.6Momentadjustingfactorundersnow-loadtypeload3.2铰支结构在类雪荷载作用下的弯矩调整系数金属拱型波纹屋盖结构在类雪荷载作用下的变化情况类似于自重类荷载作用下的情况，因此本文仅将最后的结果绘于图6中。从图中可以看出值随矢跨比改变而产生的波动也不大（约9.5%），另外各条曲线也较平直。为简化设计，本文统一取类雪荷载作用下结构的弯矩调整系数＝0.41。该取值含盖了图中所有曲线，因此偏于安全的。.3.1.3铰支结构在半跨三角形荷载作用下的弯矩调整系数图7半跨荷载下的弯矩调整系数Fig.7Momentadjustingfactorunderhalfspanload半跨三角形荷载作用下各种矢跨比结构的曲线绘于图7中。从图中可以看出随着矢跨比的增减值变化非常显著，此时若只选取一个值作为该工况下的弯矩调整系数，其设计结果偏差较大。故此将图中曲线分成两类，矢跨比不大于0.3时，曲线比较平缓，可以归为一类，取＝0.6；矢跨比大于0.3时，与之间明显呈非线性关系，且各条曲线的距离较近，因而归为另一类，取＝0.58503+0.07952×(q/qcr)-0.00312×(q/qcr)2+0.02688×(q/qcr)3。分析表明在半跨荷载作用下，这种结构的弹性稳定承载力远大于弹塑性稳定承载力，因此图中值较大（>0.5）处的曲线已不具有实际意义。图7半跨荷载下的弯矩调整系数Fig.7Momentadjustingfactorunderhalfspanload同样的推导过程可求出固支结构在上述荷载类型下的弯矩调整系数。显然固支拱型结构的稳定承载力要显著大于铰支结构。由于实际工程中，这种结构在拱脚处很难做到理想固支，实际采用的支座形式基本都为铰支座，因此本文不再列出固支时的弯矩调整系数。4、金属拱型波纹屋盖结构设计步骤根据本文推导简化设计公式的过程及思路，可总结按照上述方法设计金属拱型波纹屋盖结构的步骤如下：对结构所受各类荷载进行组合，确定各种工况下各类荷载的设计值；对于多雪地区，须考虑半跨非均布雪荷载的作用。因此一般情况下设计时要考虑的荷载组合类型有：自重类荷载＋全跨均布雪荷载/活荷载自重类荷载＋全跨均布雪荷载/活荷载＋风荷载自重类荷载＋风荷载自重类荷载＋半跨三角形分布雪荷载自重类荷载＋半跨三角形分布雪荷载＋风荷载（2）按照文献[3]提供的方法计算结构单位宽度截面的各等效弹性常数；（3）按照线弹性理论计算出结构在各种设计荷载下的各截面的弯矩和轴力。（4）利用公式（5）求出各类荷载对应的结构弹性临界荷载：（5）利用公式（4）求出各类设计荷载对应的弯矩放大系数，并进而利用公式（2）和公式（3）求出特定荷载工况下结构各截面的一阶轴力及二阶弯矩。（6）按照公式（1）验算各工况下结构各截面的承载力。5、工程试设计本文运用前面推导的简化设计公式，并根据上一节的简化设计步骤编制了金属拱型波纹屋盖设计软件（ACMRSAP）。该软件采用VB语言编程，具有很好的用户界面，能够自动进行荷载组合并对各种荷载工况下的结构进行内力分析及承载力验算，可以图形显示并从图上查询各种计算结果，可以直接打印结构内力图及设计计算书。本文运用此软件进行了大量的工程试设计，均得到较满意结果。下面将列举出其中三个工程的设计结果。5．1工程一结构参数：跨度30m，拱高6m，板厚1.4mm，板型MIC-240，支座类型铰支座，钢板屈服强度280MPa，设计强度235MPa荷载标准值：自重0.20kN/m2，基本雪压0.75kN/m2，活荷载0.30kN/m2，基本风压0.60kN/m2计算结果：截面最大应力420.4MPa>235MPa结构不安全。该工程为一实际工程，工程地点在乌鲁木齐市。工程完工于1999年6月，于2000年1月在一场大雪后倒塌。对事故现场的勘测表明此工程施工质量很好。5．2工程二结构参数：跨度33m，拱高6.6m，板厚1.25mm，板型MIC-240，支座类型铰支座，板屈服强度280MPa，设计强度235MPa荷载标准值：自重0.18kN/m2，基本雪压0.40kN/m2，活荷载0.30kN/m2，基本风压0.55kN/m2计算结果：截面最大应力291.9MPa>235MPa结构不安全。此工程也是一事故工程，建于1995年并于1997年1月在一场大雪后塌落，工程地点在辽宁省鞍山市。在事故现场的足尺模型试验表明，结构在半跨雪荷载作用下的承载能力不足是导致事故的根本原因[5]。5．2工程三结构参数：跨度24m，拱高