计算长度系数数值计算方法研究.docx

1、计算长度系数数值计算方法研究1 概述单层网壳在壳体平面内、外的屈曲模态不同，因此其杆件在壳体曲面内、外的计算长度不同。关于该部分内容，空间网格结构技术规程已给出杆件的计算长度的取值依据并在条文说明中做出解释。但是目前空间钢结构体型越来越复杂，完全按照规范取值可能导致不安全。2 数值方法的提出对于刚性节点的网壳，当处于最不利状态的压杆失稳时，其他杆件也随同弯曲，从而对失稳杆件起到约束作用。与失稳杆件直接相连的杆件约束作用大，相距较远的杆件约束作用小，相连杆件约束作用的大小取决于它的线刚度和受力状态。同时，荷载的分布模式对杆件屈曲也有非常大的影响。所以，合理确定荷载分布模式和周围杆件对目标杆件的约

2、束程度是较为精确的数值计算的前提。 基于弹性理论的计算杆件的计算长度稀疏的数值方法主要有以下两种：1， 基于整体模型中对目标杆件施加单位力的求解方法，以下简称为整体模型单位力法（）；2， 基于目标杆件单体模型的弹簧约束刚度系数法，以下简称为目标杆件弹簧系数法（）；2.1 整体模型单位力法（）的求解步骤1， 在整体模型中，沿构件两端施加一对轴向单位压力；2， 在此工况下，对整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临界荷载系数和静力分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹陷失稳的临界荷载值；3， 根据欧拉临界荷载的计算公式便可以反算出杆件的计算长度系数。2.2 目标杆件弹簧系数

3、法（）的求解步骤利用该方法求解计算长度系数的步骤如下：1， 在整体模型中删除目标杆件，在删除的目标杆件两端节点上分别施加Fx1、Fy1、Fz1、Mx1、My1、Mz1、Fx2、Fy2、Fz2、Mx2、My2、Mz2共12组力，每组力作为一种单独工况；2， 计算分析得到每组力下相对应方向下的位移，通过矩阵运算，从而可得到目标杆件的约束弹簧钢度系数；3， 建立单个目标杆件模型，保证杆件的局部坐标系和整体坐标系的对应关系，将上述所得的弹簧系数作为边界条件施加在目标杆件上；4， 沿目标杆件两端施加一对轴向单位压力；5， 在此工况下，对整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临界荷载系数和静力

4、分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹陷失稳的临界荷载值；6， 根据欧拉临界荷载的计算公式便可以反算出杆件的计算长度系数。3 算例分析及比较本章将基于深圳湾体育中心项目对此进行研究。该项目为大跨度空间结构其中单层网壳部分杆件空间扭曲特征明显，体育场周边钢结构网壳悬挑，具有代表性。3.1 整体模型单位力法算例1 在整体模型中，沿构件两端施加一对轴向单位压力，如图4-1所示。图4-1 施加的单位荷载示意图2 在此工况下，对整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临界荷载系数和静力分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹陷失稳的临界荷载值。屈曲分析对应的第一模态云图如

5、图4-2所示：图4-2 杆件的第一屈曲模态3 根据欧拉临界荷载的计算公式便可以反算出杆件的计算长度系数。根据以上方法，分别选取了体育场悬挑屋面的边缘杆件和中间杆件，单层网壳的其他部分杆件和树形柱支撑进行数值计算，选取的杆件部位如图4-3中红色标示杆件所示。计算结果如表4-4所示。需要说明的是，以上方法是一种近似算法，其荷载分布模式的简化对结果有一定的影响。由于仅对计算杆件施加轴向压力，其轴向相邻杆件处于受拉状态，可以认为考虑了杆件的刚度约束条件，但是荷载条件和由荷载条件决定的几何刚度条件考虑不充分。但是由于该方法操作简单，易于分析，对于我们初步判断杆件的计算长度系数还是比较有效的。图4-3 计

6、算长度系数数值计算选取杆件部位图示表4-4杆件的计算长度系数单位： 杆件编号所处部位杆件长度施加荷载临界荷载系数计算长度1体育场悬挑边缘4.321.01135904.60.96117133.22112946.889.96 2.31 25.681.0972237.60.9486907.0281304.1310.86 1.91 35.491.0972237.60.9777580.4375638.5911.26 2.05 44.011.0626178.20.98128606.47125599.657.01 1.75 55.281.0626178.20.89136881.72122119.037.11

7、 1.35 65.351.0626178.20.9687973.0684315.148.56 1.60 7体育场悬挑中间3.141.0409610.40.9298226.1090417.136.69 2.13 84.661.0409610.40.9289232.9482167.477.01 1.50 94.611.0409610.40.92102061.9193751.016.57 1.42 104.561.0409610.40.92114636.10104983.746.21 1.36 114.531.0547218.40.93124855.84116374.396.81 1.50 124.

8、071.0409610.40.93126479.83117636.365.86 1.44 134.081.0409610.40.92184068.73168437.614.90 1.20 14树形柱支撑23.861.02062678.00.74105594.2378501.9216.10 0.67 1517.601.0316045.20.9521819.0420656.3012.29 0.70 3.2 目标杆件弹簧系数法按单自由度方法计算弹簧刚度系数按照前面对该方法操作步骤地论述，对一个杆件，在两端节点上共需加载12组力，即12种荷载工况，每种工况在每个节点上又将产生6个自由度方向上的力，因此

9、位移矩阵为6x6的对称阵。但是实际上进行矩阵运算较为繁琐，工程应用不是很方便，所以尝试采用单自由的计算方法，即每种力作用的结果仅考虑该方向上的位移，用数学表达式表示为：（分别代表6个自由度方向） 按上述描述方法，在整体模型中删除目标杆件，在删除的目标杆件两端节点上分别施加Fx1、Fy1、Fz1、Mx1、My1、Mz1、Fx2、Fy2、Fz2、Mx2、My2、Mz2共12组力，每组力作为一种单独工况，大小取为100000KN或100000KN.m，计算目标杆件的弹簧系数如表4-6所示。表4-6目标杆件的弹簧刚度系数节点荷载工况施加力DX(m)DY (m)DZ (m)RX (rad)RY(rad)

10、RZ (rad)K(KN/m)8760Fx21000005.6790.963-16.282-0.5660.695-0.192K1117608 8760Fy21000000.9630.865-4.134-0.0450.3050.055K22115662 8760Fz2100000-16.282-4.13463.7562.470-3.5350.071K331568 8760Mx2100000-0.566-0.0452.4700.527-0.1250.055K44189732 8760My21000000.6950.305-3.535-0.1251.1050.047K5590470 8760Mz21

11、00000-0.1920.0550.0710.0550.0470.268K66372734 8761Fx11000006.3990.383-17.5750.2301.1670.309K1115626 8761Fy11000000.3830.663-0.7810.1250.0090.156K22150762 8761Fz1100000-17.575-0.78165.418-1.209-5.0450.066K331529 8761Mx11000000.2300.15-1.2090.452-0.0220.039K44221138 8761My11000001.1670.009-5.045-0.022

12、1.413-0.018K5570791 8761Mz11000000.3090.1560.0660.039-0.0180.333K66300586 目标杆件弹性屈曲分析模型如图4-38所示。目标杆件的截面特性为，,，施加的杆端压力为，分析后所得杆件内力为，第一阶屈曲模态对应的临界荷载系数为87144.02，那么，那么：，,计算值要大于4.3.3节中描述的方法计算值2.31。图4-38目标杆件弹性屈曲分析模型图目标杆件第一阶屈曲模态如图4-39所示。图4-39目标杆件第一阶屈曲模态图3.3 目标杆件弹簧系数法考虑完整的6x6刚度矩阵通过上面分析可知，按照单自由度体系的方法计算刚度矩阵是有缺陷的，

13、其忽略了各个自由度方向上刚度的相互影响。为验证其具体影响，下面补充采用6x6刚度矩阵的计算方法的分析。计算方法仍然同4.6.3节中相关内容，只是该处要求得完整的6x6刚度矩阵。依据表4-6中相应的力和位移，求得完整的6x6刚度矩阵如表4-8所示。为方便说明，4.6.3节所述方法记为方法1，该节中所述方法为方法2，4.3.3节中描述的整体计算法记为方法。表4-8目标杆件的弹簧刚度系数节点荷载工况施加力6x6刚度系数8760Fx210000074461 8858 21070 -19695 13875 47554 8760Fy21000008858 181662 15255 -46684 -1119

14、4 -23433 8760Fz210000021070 15255 9130 -17881 9229 11596 8760Mx2100000-19695 -46684 -17881 253561 -1049 -51645 8760My210000013875 -11194 9229 -1049 114696 -10106 8760Mz210000047554 -23433 11596 -51645 -10106 421311 8761Fx110000079006 -6194 23088 30781 16710 -77688 8761Fy1100000-6194 179169 -4 -39961

15、 2402 -73377 8761Fz110000023088 -4 9054 15206 13184 -24284 8761Mx110000030781 -39961 15206 262500 32679 -41833 8761My110000016710 2402 13184 32679 104315 -17433 8761Mz1100000-77688 -73377 -24284 -41833 -17433 415534 方法1和方法2中求得对角线上刚度系数的比较如表4-9所示。表4-9两种计算方法下的刚度系数比较节点刚度序号方法方法比值8760K1117608744614.23K221

16、156621816621.57K33156891305.82K441897322535611.34K55904701146961.27K663727344213111.138761K1115626790065.06K221507621791691.19K33152990545.92K442211382625001.19K55707911043151.47K663005864155341.38目标杆件的截面特性为，,，施加的杆端压力为，分析后所得杆件内力为，第一阶屈曲模态对应的临界荷载系数为111503.00，那么，那么，。计算值要大于4.3.3节中描述的方法计算值2.31。三种方法下计算的杆件

17、计算长度系数比较如表4-10所示。表4-10三种方法下杆件的计算长度系数比较单元编号方法方法方法比值比值2.652.342.311.151.01注：比值以方法为基本值。通过表4-9的比较可以发现，对于弹簧系数法，如果刚度系数取用正确，和整体计算方法的计算结果是非常接近的。4 所取杆件数的影响类似如深圳湾体育中心这类结构类型，由于复杂的单层网壳曲面外无明显的支撑，将荷载模式仅选为加在目标杆件两端明显和实际受力特性不符。所以以下尝试将荷载模式选为加在多个杆件上。4.1 基于整体模型单位力法下加载杆件数目的影响分析方法同上，选取目标杆件时沿杆件长度方向依次选取一个杆件、两个杆件、三个杆件、四个杆件和

18、五个杆件，对比结果如下：选取杆件在悬挑部位：CASE1:杆件编号杆件长度施加荷载EIN1临界荷载系数Pcr计算长度计算长度系数NO.13.1414.10E+051.00114159.53114159.535.951.90CASE2杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.13.1414.10E+051.0064765.8864765.887.902.52NO.24.7114.10E+051.0064765.8864765.887.901.68CASE3杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.13.1414.10E+051.005

19、7347.8257347.828.402.68NO.24.7114.10E+051.0057347.8257347.828.401.78NO.34.6714.10E+051.0057347.8257347.828.401.80CASE4杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.13.1414.10E+051.0057356.0057356.008.402.67NO.24.7114.10E+051.0057356.0057356.008.401.78NO.34.6714.10E+051.0057356.0057356.008.401.80NO.44.6114.10

20、E+051.0057356.0057356.008.401.82CASE5杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.13.1414.10E+051.0059138.5459138.548.272.63NO.24.7114.10E+051.0059138.5459138.548.271.75NO.34.6714.10E+051.0059138.5459138.548.271.77NO.44.6114.10E+051.0059138.5459138.548.271.79NO.54.5814.10E+051.0059138.5459138.548.271.81选取杆件

21、在非悬挑部位：CASE1:杆件编号杆件长度施加荷载EIN1临界荷载系数Pcr计算长度计算长度系数NO.14.32 11.14E+060.96 136667.31 131782.82 9.22 2.14 CASE2杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.14.32 11.14E+060.86 71332.79 61700.72 13.48 3.12 NO.26.22 11.14E+060.83 71332.79 59073.53 13.78 2.22 CASE3杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.14.32 11.14E+0

22、60.76 61508.08 46586.22 15.51 3.59 NO.26.22 11.14E+060.72 61508.08 44395.92 15.89 2.56 NO.36.01 11.14E+060.77 61508.08 47572.19 15.35 2.55 CASE4杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.14.32 11.14E+060.66 60741.72 40260.22 16.69 3.87 NO.26.22 11.14E+060.63 60741.72 38192.57 17.13 2.76 NO.36.01 11.14E+06

23、0.67 60741.72 40545.70 16.63 2.77 NO.45.84 11.14E+060.65 60741.72 39370.35 16.87 2.89 CASE5杆件编号杆件长度施加荷载EN1临界荷载系数Pcr计算长度计算长度系数NO.14.32 11.14E+060.62 62486.41 38599.73 17.04 3.95 NO.26.22 11.14E+060.59 62486.41 36615.16 17.50 2.81 NO.36.01 11.14E+060.63 62486.41 39510.16 16.84 2.80 NO.45.84 11.14E+060

24、.62 62486.41 38491.63 17.07 2.92 NO.52.10 11.14E+060.60 62486.41 37784.28 17.23 8.20 4.2 目标杆件弹簧系数法按单自由度方法计算弹簧刚度系数求解目标杆件的弹簧刚度系数方法同上，计算值如表4-7所示。表4-7目标杆件的弹簧刚度系数节点荷载工况施加力DX(m)DY (m)DZ (m)RX (rad)RY(rad)RZ (rad)K(KN/m)8758Fx21000005.626 1.416 -15.729 -0.393 0.732 -0.051 K1117774 8758Fy21000001.416 1.028

25、-4.747 -0.013 0.328 0.039 K2297230 8758Fz2100000-15.729 -4.747 56.544 1.490 -3.211 0.040 K331769 8758Mx2100000-0.393 -0.013 1.490 0.371 -0.123 0.037 K44269557 8758My21000000.732 0.328 -3.211 -0.123 0.891 0.033 K55112235 8758Mz2100000-0.051 0.039 0.040 0.037 0.033 0.201 K66496929 8760Fx11000009.061 -1.524 -16.347 -0.108 1.082 -2.105 K1111037 8760Fy1100000-1.524 2.794 -4.507 -0.330 0.26