混凝土密肋式锥面网壳的动力特性分析

杨水艳 张华刚 马克俭

摘 要：混凝土密肋式锥面网壳是由相等的密肋平板在脊线处交汇形成的一种新型网壳结构。为了解網壳结构的动力特性，采用子空间迭代法求解结构的自振频率及振型，并通过讨论影响结构基频的主要因素，对基频进行公式拟合。结果表明：结构的振型以竖向振动为主;因结构刚度分布均匀，频率多次“重频”且有明显的跳跃点;网格数不变时，增大结构矢跨比对提高结构整体刚度是不利的，建议屋盖矢跨比的选取不宜过大;边梁刚度的增大会加强对密肋板的约束，进而提高结构整体刚度，建议边梁截面高度的选取宜按跨度的1/40～1/30确定;增大密肋梁刚度能有效提高结构整体刚度，但过大的密肋梁高度会加大结构自重，建议密肋梁截面高度的取值宜按跨度的1/100～1/75确定;结构基频的公式拟合值误差满足工程精度设计要求。

关键词：混凝土密肋式锥面网壳;动力特性;基频;结构刚度

中图分类号：TU399

文献标志码：A

我国薄壁结构最早使用于20世纪40年代后期，如1948年在常州建成的一个仓库便是用了圆柱面壳体结构[1]。因薄壁结构曲面施工困难和理论计算复杂等问题，至60年代后期，薄壁结构的使用逐渐萧条。近年来，数值计算技术的发展使薄壁结构的计算分析变得较为容易，但施工困难是尚未有效解决的问题，因此，我国许多学者致力于壳体结构的形式创新研究。马克俭等[2]提出了将空腹夹层板弯曲成曲面形式的空腹网壳结构，并于2008年在中国电建集团贵阳勘测设计院会议中心工程中成功应用。常玉珍等[3]提出了外包U型钢的组合肋壳结构。金杰等[4]在钢管中注入混凝土而提出了钢管混凝土网壳结构。上述结构形式推动了壳体结构的发展。

为有效解决混凝土曲面支模难和降低工程造价等问题，张华刚等[5-8]提出由密肋平板在脊线处交汇形成新型混凝土折板式密肋网壳。其中，混凝土人字形折板式密肋网壳已成功运用到实际工程中。在关岭美食城工程中，用钢量仅为43.0 kg/m2，可见这类新型结构的用钢量较低，具有良好的经济技术指标。

为了解混凝土密肋式锥面网壳结构的动力特性，基于数值模拟分析，采用子空间迭代法[9]求解结构自振频率及振型，以期为这种结构的抗震分析提供参考。

1 结构形式及算例基本情况

1.1 结构形式

结构形式如图1所示。通过将圆锥沿曲面等分切割成三角形密肋平板，再由密肋平板在脊线处交汇形成锥面网壳。密肋梁网格采用正交正方的形式布置，在边梁与脊线交汇处设置支座，并约束其全部自由度。

1.2 算例情况

如图2所示，结构跨度为30 m，矢高为7.5 m，经承载力设计后，屋面板厚为60 mm，边梁截面尺寸为400 mm×800 mm，脊线截面尺寸为300 mm×700 mm，密肋梁截面尺寸为150 mm×400 mm。结构自振分析时，屋面板采用板壳单元，其余构件均采用空间梁单元，混凝土材料的弹性模量Ec=3.25×104 N/mm2，泊松比v=0.2，钢筋混凝土密度为2.42×103 kg/m3;结构不计自重的恒载取5.0 kN/m2，上述均为本文所有算例的共性参数。

2 自振特性分析

2.1 自振频率

结构前50阶自振频率如图3所示。频谱分布较为密集，且随振型阶数的增大有明显的跳跃性，可见结构刚度分布较为均匀，结构成对出现的振型较多，这是结构有多条对称轴的缘故[10-12]。

2.2 振型

结构前9阶振型如图4所示。低阶振型主要以竖向振动为主，模态坐标较大区域位于边梁附近的密肋平板上。第1阶振型沿环向关于脊线各有6个半波反对称振动，第2～5阶振型呈现2个半波振动，第6阶振型沿环向关于脊线各有6个半波正对称振动，第7～9阶振型在每块密肋平板上均出现多个振动区域。低阶振型的振动节线大体为脊线，可见脊线对结构具有拱向支承的作用。

综上分析，结构第1阶竖向振动频率可近似反映其整体刚度。在下文参数化分析时，主要考察不同因素对结构基频的影响[13]。

3 自振频率的参数化分析

3.1 参数化分析算例

在基本算例的基础上，分别考虑矢跨比、边梁刚度、脊线刚度、密肋梁刚度和屋面板厚等因素对结构基频的影响。算例取值情况如下：

1）仅改变屋盖矢跨比计算6个算例，矢跨比分别取1/8、1/7、1/6、1/5、1/4和1/3。

2）考虑边梁刚度的影响时，主要通过改变其截面高度实现，脊线截面高度取650 mm，边梁截面高度分别取700、750、800、850、900、950 mm。

3）考虑脊线刚度的影响时，边梁截面高度取900 mm，脊线截面高度分别取600、650、700、750、800、850 mm。

4）考虑密肋梁刚度的影响时，密肋梁截面高度分别取300、350、400、450、500、550 mm。

5）考虑板厚的影响时，屋面板厚分别取50、60、70、80、90、100 m。

3.2 矢跨比的影响

改变屋盖矢跨比对基频的影响如图5所示。基频随着矢跨比的增大而持续减小。由图5可见：当矢跨比为1/8时，基频为8.19 Hz;当矢跨比为1/3时，基频为6.57 Hz，基频降幅约为19.8%;矢跨比在1/8～1/5范围内基频变化率较小，这是因为随着矢跨比的增加，结构的展开面积增大，在网格数不变时，将降低结构整体刚度，并增大结构质量。在考虑结构网格数不变的情况下，相当于间接稀疏了斜板的部分网格，从而导致结构刚度的减小。因此，矢跨比的改变对结构刚度影响较为明显，建议结构矢跨比的选取不宜大于1/5。

3.3 边梁刚度的影响

边梁刚度的改变对基频的影响如图6所示。随着边梁截面高度的增大，结构的基频不断增大。由图6可见：当边梁截面高度为700 mm时，基频为6.72 Hz;当边梁截面高度为950 mm时，基频为7.88 Hz，基频增幅约为17.3%。边梁作为结构的主要传力构件，增大边梁的刚度会加强对密肋梁的约束，进而延缓结构变形，提高结构整体刚度。因此，改变边梁刚度对提高结构整体刚度有一定贡献，建议边梁截面高度的取值宜按跨度的1/40～1/30确定。

3.4 脊线刚度的影响

脊线刚度对结构基频的影响如图7所示。提高脊线刚度对结构的基频影响较小。由图7可见：当脊线截面高度为600 mm时，基频为7.67 Hz;当脊线截面高度增到850 mm时，基频为7.75 Hz，基频增幅仅为1.0%。因为脊线成拱后自身刚度较大，而密肋平板才是结构变形的薄弱环节，因此，增大脊线截面高度对结构整体刚度的提高效果不明显，脊线截面高度的选取满足强度设计要求即可。

3.5 密肋梁刚度的影响

密肋梁刚度的改变对基频的影响如图8所示。基频随着密肋梁刚度的增大而增大。由图8可见：当密肋梁截面高度取300 mm时，基频为6.51 Hz;当密肋梁截面高度取 550 mm时，基频为8.05 Hz，基频增幅约为23.7%。因结构整体刚度受密肋平板控制，加大密肋梁刚度可直接增大每块密肋板的刚度，进而明显提高屋盖整体刚度，但过大的密肋梁刚度会加大结构自重。密肋梁截面高度的选取宜按跨度的1/100～1/75确定。

3.6 板厚的影响

屋面板厚的改变对基频的影响如图9所示。随着屋面板厚的增大，基频呈持续下降的趋势。由图9可见：当屋面板厚为50 mm时，基频为7.46 Hz;当屋面板厚为100 mm时，基频为6.67 Hz，基频降幅约为10.6%。

综上分析可得，屋面板厚度增加对结构刚度的贡献不足以抵消其质量增加对基频的影响，而质量增加又会加大地震作用，因此，建议屋盖板厚选取满足构造要求即可。

4 基频公式拟合

4.1 等效刚度的计算

由上述数值模拟分析可知，结构的刚度由边梁、脊线、密肋梁和屋面板厚提供，可将其连续化进行等效刚度的计算[14]。根据图10分别计算密肋梁的等效薄膜刚度B1和等效抗弯刚度D1分别为：

5 结语

1）结构的振型以竖向振动为主，且在密肋折板中下区域的竖向振动出现较早。因结构刚度分布均匀，频率多次“重频”且有明显的跳跃点。

2）网格数不变时，增大结构矢跨比对结构整体刚度的控制是不利的，建议选取屋盖矢跨比不宜过大。

3）边梁刚度的增大会加强对密肋板的约束，进而提高结构整体刚度。建议边梁截面高度的选取宜按跨度的1/40～1/30确定。

4）脊线是振动节线，且其成拱后自身刚度较大，因此，增大脊线刚度对提高结构整体刚度效果不明显。

5）增大密肋梁刚度对提高结构整体刚度的影响较大，但大的密肋梁高度会加大结构自重，密肋梁截面高度的取值宜按跨度的1/100～1/75确定。

6）结构基频的拟合公式计算误差能满足工程设计精度。

参考文献：

[1]常玉珍， 吴敏哲. 组合肋壳结构静力特性分析[J]. 建筑结构， 2009， 39（4）： 54-56.

[2]马克俭， 张华刚， 肖建春， 等. 中国钢筋混凝土空间网格结构新体系的开拓与发展[J]. 中国工程科学， 2008， 10（8）： 22-34.

[3]常玉珍， 吴敏哲. 组合肋壳结构静力特性分析[J]. 建筑结构， 2009， 39（4）： 54-56.

[4]金杰， 蔡相娟. 钢管混凝土网壳的可行性分析[J]. 混凝土， 2010（3）： 51-53， 57.

[5]张华刚， 杨期柱， 马克俭， 等. 人字形密肋式折板拱壳的静力性能分析及其工程应用[J]. 建筑结构学报， 2009， 30（S2）： 29-35.

[6]柳勇斌， 张华刚， 吴琴， 等. 混凝土折板形柱面网壳结构静力性能的参数化分析[J]. 空间结构， 2019， 25（3）： 22-29，37.

[7]方强， 张华刚， 齐卫东， 等. 混凝土幕形密肋网壳的静力性能及参数化分析[J]. 广西大学学报（自然科学版）， 2017， 42（4）： 1307-1316.

[8]张华刚， 张钰， 吴琴， 等. 新型混凝土壳体结构的研究概况及工程实践[J]. 空间结构， 2019， 25（3）： 3-12， 21.

[9]包世华. 结构动力学[M]. 武汉： 武汉理工大学出版社， 2005.

[10]姜腾钊， 张华刚， 魏威， 等. 混凝土折板式拱网壳动力特性的参数化分析[J]. 贵州大学学报（自然科学版）， 2020， 37（2）： 73-78.

[11]朱锐， 张华刚， 陈寿延， 等. 斜放网格棱柱面密肋折板网壳的动力特性分析[J]. 贵州大学学报（自然科学版）， 2019， 36（6）： 83-91.

[12]唐攒辉， 张华刚， 魏威， 等. 混凝土棱柱面网壳的弹性稳定分析[J]. 貴州大学学报（自然科学版）， 2019， 36（6）： 92-97.

[13]姜岚， 张华刚. 大跨度空腹夹层板楼盖基于舒适度要求的动力特性分析[J]. 空间结构， 2014， 20（3）： 56-60， 81.

[14]中华人民共和国住房和城乡建设部. 空间网格结构技术规程： JGJ 7—2010 [S]. 北京： 中国建筑工业出版社， 2010.

[15]ZHANG L W， WU J， ZHANG D L. Shape optimization and stability analysis for kiewitt spherical reticulated shell of triangular pyramid system[J]. Mathematical Problems in Engineering， 2019（1）： 1-11.

（责任编辑：周晓南）

Abstract：

Concrete Multi-ribbed conical reticulated shell is a new type of reticulated shell structure formed by the intersection of equal ribbed flat plates at the ridgeline. In order to understand the dynamic characteristics of reticulated shell structures， the subspace iteration method is used to solve the natural frequency and mode of the structure. By discussing the main factors affecting the fundamental frequency of the structure， the fundamental frequency formula is fitted. The results show that： the vibration mode of the structure is mainly vertical vibration. Because the structure stiffness is uniformly distributed， the frequency is "repetitive frequencies" for many times and there are obvious jumping points. When the number of grids is constant， it is unfavorable to increase the rise-span ratio to span of the structure to improve the overall stiffness of the structure. It is suggested that the rise-span ratio to span of the roof should not be too large. The increase of the stiffness of the edge beam will strengthen the constraint of the dense floor and thus improve the overall stiffness of the structure. It is recommended that the cross-section height of the edge beam should be determined according to 1/40-1/30 of the span. Increasing the stiffness of the multi-ribbed beam has a greater impact on improving the overall rigidity of the structure. But a large multi-ribbed beam height will increase the self-weight of the structure. The selection of beam section height should be determined according to 1/100-1/75 of the span. The error of the formula fitting value of the structural fundamental frequency can meet the engineering design accuracy.

Key words：

concrete cone reticulated shell; dynamic characteristics; fundamental frequency; structural stiffness