《计算长度系数数值计算》由会员分享，可在线阅读，更多相关《计算长度系数数值计算（11页珍藏版）》请在金锄头文库上搜索。
1、C:.DESKTOP计算长度系数数值计算方法研究.DOCPage 1Arup International Consultants (Shanghai) Co Ltd Shenzhen Branch Draft 1 13 August 20091概述单层网壳在壳体平面内、外的屈曲模态不同，因此其杆件在壳体曲面内、外的计算长 度不同。关于该部分内容，空间网格结构技术规程已给出杆件的计算长度的取值 依据并在条文说明中做出解释。但是目前空间钢结构体型越来越复杂，完全按照规范 取值可能导致不安全。2数值方法的提出对于刚性节点的网壳，当处于最不利状态的压杆失稳时，其他杆件也随同弯曲，从而 对失稳杆件起到约
2、束作用。与失稳杆件直接相连的杆件约束作用大，相距较远的杆件 约束作用小，相连杆件约束作用的大小取决于它的线刚度和受力状态。同时，荷载的 分布模式对杆件屈曲也有非常大的影响。所以，合理确定荷载分布模式和周围杆件对 目标杆件的约束程度是较为精确的数值计算的前提。 基于弹性理论的计算杆件的计算长度稀疏的数值方法主要有以下两种：1， 基于整体模型中对目标杆件施加单位力的求解方法，以下简称为整体模型单位力 法（）；2， 基于目标杆件单体模型的弹簧约束刚度系数法，以下简称为目标杆件弹簧系数法 （）；2.1整体模型单位力法（）的求解步骤1， 在整体模型中，沿构件两端施加一对轴向单位压力；2， 在此工况下，对
3、整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临界荷 载系数和静力分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹陷失稳的 临界荷载值；crP3， 根据欧拉临界荷载的计算公式便可以反算出杆件的计算长度22LEIPcr系数。2.2目标杆件弹簧系数法（）的求解步骤利用该方法求解计算长度系数的步骤如下：1， 在整体模型中删除目标杆件，在删除的目标杆件两端节点上分别施加 Fx1、Fy1、Fz1、Mx1、My1、Mz1、Fx2、Fy2、Fz2、Mx2、My2、Mz2 共 12 组力， 每组力作为一种单独工况；2， 计算分析得到每组力下相对应方向下的位移，通过矩阵运算，从而可. KF 得到目标杆件的
4、约束弹簧钢度系数；3， 建立单个目标杆件模型，保证杆件的局部坐标系和整体坐标系的对应关系，将 上述所得的弹簧系数作为边界条件施加在目标杆件上；4， 沿目标杆件两端施加一对轴向单位压力；5， 在此工况下，对整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临界 荷载系数和静力分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹陷失 稳的临界荷载值；crP6， 根据欧拉临界荷载的计算公式便可以反算出杆件的计算长度22LEIPcr系数。3算例分析及比较本章将基于深圳湾体育中心项目对此进行研究。该项目为大跨度空间结构其中单层网 壳部分杆件空间扭曲特征明显，体育场周边钢结构网壳悬挑，具有代表性。C:.DE
5、SKTOP计算长度系数数值计算方法研究.DOCPage 2Arup International Consultants (Shanghai) Co Ltd Shenzhen Branch Draft 1 13 August 20093.1整体模型单位力法算例 在整体模型中，沿构件两端施加一对轴向单位压力，如图 4-1 所示。图 4-1 施加的单位荷载示意图 在此工况下，对整体模型进行弹性屈曲分析，可以得到杆件局部屈曲时的临 界荷载系数和静力分析时对应的杆件轴向压力，两者的乘积即为杆件局部凹 陷失稳的临界荷载值。屈曲分析对应的第一模态云图如图 4-2 所示：crP图 4-2 杆件的第一屈曲模态
6、根据欧拉临界荷载的计算公式便可以反算出杆件的计算长22LEIPcr度系数。根据以上方法，分别选取了体育场悬挑屋面的边缘杆件和中间杆件，单层网壳的 其他部分杆件和树形柱支撑进行数值计算，选取的杆件部位如图 4-3 中红色标示 杆件所示。计算结果如表 4-4 所示。需要说明的是，以上方法是一种近似算法，其荷载分布模式的简化对结果有一定 的影响。由于仅对计算杆件施加轴向压力，其轴向相邻杆件处于受拉状态，可以 认为考虑了杆件的刚度约束条件，但是荷载条件和由荷载条件决定的几何刚度条 件考虑不充分。但是由于该方法操作简单，易于分析，对于我们初步判断杆件的 计算长度系数还是比较有效的。C:.DESKTOP计
7、算长度系数数值计算方法研究.DOCPage 3Arup International Consultants (Shanghai) Co Ltd Shenzhen Branch Draft 1 13 August 2009图 4-3 计算长度系数数值计算选取杆件部位图示表 4-4 杆件的计算长度系数单位：mKN,杆件 编号所处 部位杆件 长度施加 荷载EI1N临界荷 载系数crP计算 长度14.321.01135904.60.96117133.22112946.889.96 2.31 25.681.0972237.60.9486907.0281304.1310.86 1.91 35.491.09
8、72237.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 1.35 6体 育 场 悬 挑 边 缘5.351.0626178.20.9687973.0684315.148.56 1.60 73.141.0409610.40.9298226.1090417.136.69 2.13 84.661.0409610.40.9289232.9482167.477.01 1.50 94.611.0409610.
9、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.071.0409610.40.93126479.83117636.365.86 1.44 13体 育 场 悬 挑 中 间4.081.0409610.40.92184068.73168437.614.90 1.20 1423.861.02062678.00.74105594.2378501.9216.10 0.67 15树形 柱支撑
10、17.601.0316045.20.9521819.0420656.3012.29 0.70 3.2目标杆件弹簧系数法按单自由度方法计算弹簧刚度系数按照前面对该方法操作步骤地论述，对一个杆件，在两端节点上共需加载 12 组力， 即 12 种荷载工况，每种工况在每个节点上又将产生 6 个自由度方向上的力，因此位 移矩阵为 6x6 的对称阵。但是实际上进行矩阵运算较为繁琐，工程应用不是很方便， 所以尝试采用单自由的计算方法，即每种力作用的结果仅考虑该方向上的位移，用数 学表达式表示为：（ 分别代表 6 个自由度方向） iiiKFi按上述描述方法，在整体模型中删除目标杆件，在删除的目标杆件两端节点上
11、分别施 加 Fx1、Fy1、Fz1、Mx1、My1、Mz1、Fx2、Fy2、Fz2、Mx2、My2、Mz2 共 12 组力，每C:.DESKTOP计算长度系数数值计算方法研究.DOCPage 4Arup International Consultants (Shanghai) Co Ltd Shenzhen Branch Draft 1 13 August 2009组力作为一种单独工况，大小取为 100000KN 或 100000KN.m，计算目标杆件的弹簧系 数如表 4-6 所示。表 4-6 目标杆件的弹簧刚度系数节点荷载 工况施加力DX (m)DY(m)DZ (m)RX (rad)RY (
12、rad)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 8760
13、Mz2100000-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
14、.0221.413-0.018K5570791 8761Mz11000000.3090.1560.0660.039-0.0180.333K66300586 目标杆件弹性屈曲分析模型如图 4-38 所示。目标杆件的截面特性为，,，4040450700BmL32. 426/102mKNE ，施加的杆端压力为，分析后所得杆件内力为，第40.0055141mIKN1KN9851. 0一阶屈曲模态对应的临界荷载系数为 87144.02，那么，那么：KNPcr57.85845，,计算值要大于 4.3.3 节中描述的方法mPcrEILe43.11265. 2LLe计算值 2.31。图 4-38 目标杆件弹性
15、屈曲分析模型图目标杆件第一阶屈曲模态如图 4-39 所示。C:.DESKTOP计算长度系数数值计算方法研究.DOCPage 5Arup International Consultants (Shanghai) Co Ltd Shenzhen Branch Draft 1 13 August 2009图 4-39 目标杆件第一阶屈曲模态图3.3目标杆件弹簧系数法考虑完整的 6x6 刚度矩阵通过上面分析可知，按照单自由度体系的方法计算刚度矩阵是有缺陷的，其忽略了各 个自由度方向上刚度的相互影响。为验证其具体影响，下面补充采用 6x6 刚度矩阵的 计算方法的分析。计算方法仍然同 4.6.3 节中相关
16、内容，只是该处要求得完整的 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 -11194 -23433 8760Fz210000021070 15255 9130 -17881 9229 11596 8760Mx2100000-19695 -46684 -17881 253561 -1049 -51645 8760My
