CN101712438B - 利用塔式起重机运行数据获取评估点应力时域值的方法 - Google Patents

Abstract

本发明涉及一种利用塔式起重机运行数据获取评估点应力时域值的方法，它属于一种利用塔机运行数据记录仪采集的大批量载荷数据获得评估点应力时域值的方法。本发明主要是解决现有的方法存在的采集数据样本容量小、采集结构部位有限等技术难点。本发明的技术方案是：利用塔式起重机运行数据获取评估点应力时域值的方法，包括下列步骤：1)使用数据记录仪采集塔式起重机的起重量等实时运行数据；2)对采集的数据进行处理；3)建立塔式起重机结构计算模型；4)计算某时刻结构杆件内力和应力；5)将上步所得的某时刻应力值加上或减去修正值，得到修正后的应力值；6)重复上述步骤，得到多个时刻所对应的结构内力和应力数据并整理成应力谱数据。

Description

利用塔式起重机运行数据获取评估点应力时域值的方法

技术领域

本发明涉及一种利用塔式起重机运行数据获取评估点应力时域值的方法，它属于一种利用塔机运行数据记录仪采集的大批量载荷数据、运动及动力学数据获得塔式起重机评估点应力时域值的方法。 

背景技术

塔式起重机是现代施工中必不可少的重要施工设备，广泛应用于建筑施工中，也可以用于码头、货场的货物装卸及搬运。 

随着科学技术的不断进步和塔机安全性能要求的提高，对其安全运行要求也越来越高，为此，人们从各个方面做了大量的工作，如加强塔机运行的使用管理、增加塔机安全运行的各种限位装置等。但是这些还不够，由于塔机的主要组成部分金属结构存在疲劳问题，因此即使其他方面做到位了，塔身本身存在问题，也将是安全隐患的极大事件。因此，塔机金属结构可靠性评估对于保障塔机安全运行十分重要。 

目前，塔机可靠性评估的数据来源主要是通过应变片测定指定点的应力，虽然其简便易行，但是存在着采集数据样本空间小、采集结构部位有限等缺点，若长期、多部位采集，成本较高。 

另外，现在由塔机载荷数据到结构应力数据的处理方法主要是有限元法。有限元方法分析的主要步骤是：建立塔机模型并划分单元；加载某一时刻载荷数据，输出某一时刻的内力及应力等计算结果。对于实际工作的塔机，其数据是大批量的，需要多次处理，才可得到随时间历程变化的内力、应力数据。这种处理方法主要存在的问题是：(1)批量加载时，一般需要人为对数据进行处理后加载，批量数据处理需时较长，容易出错。(2)计算结果数据输出格式较固定，用于可靠性评估时，需人为存放到所需的数据库中，耗费资源，且工作周期长。(3)购买商业化有限元分析软件，实际使用时仅用到极少一部分功能，造成资源浪费。 

发明内容

本发明的目的是解决现有的方法存在的采集数据样本容量小、采集结构部位有限等技术难点，并提供一种能实时采集大批量塔式起重机运行数据且能依据实时运行数据快速计算塔式起重机结构杆件内力、应力数据的利用塔式起重机运行数据获取评估点应力时域值的方法。 

本发明为解决上述技术难点而采用的技术方案是：利用塔式起重机运行数据获取评估点应力时域值的方法，其包括下列步骤： 

1)使用数据记录仪采集塔式起重机的起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据，并将采集到的上述实时运行数据存放到存储单元中，以获得塔式起重机较长时间的运行数据； 

2)对采集的较长时间的实时运行数据进行处理，并依据速度公式： 求出起升、变幅和回转速度，然后根据求出的两点的速度与加速度公式： 求出起升、变幅和回转加速度，以获得起升机构、变幅机构和回转机构的速度和加速度的参数； 

3)建立塔式起重机结构计算模型，并编制用于快速计算的计算机程序； 

4)读入结构模型计算所需参数和选择载荷组合形式的工况数据，计算某时刻塔式起重机结构杆件内力和应力，得到塔机杆件内力和应力数据，并将计算结果分类存放在数据库中； 

5)将上述第4)步骤所得的某时刻应力值加上或减去修正值，得到修正后的应力值，修正值的范围为计算值的0％～15％； 

6)重复第4)和第5)步骤，就可以得到塔式起重机多个时刻所对应的结构内力和应力数据，把数据整理成多个评估点的应力谱数据。 

所述选择载荷组合形式是指对起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据全部选用或选用其中的一部分。 

所述数据记录仪包括用于监控塔式起重机参数信息的传感器采集装置、微 处理单元、数据存储单元和通讯网络接口单元；塔式起重机参数信息的传感器采集装置的信号与微处理单元的信号输入接口连接；微处理单元的数据输出接口连接数据存储单元的输入接口；通讯网络接口单元的输入端连接微处理单元的通讯输出接口。 

所述塔式起重机参数信息的传感器采集装置由传感器和模数转换芯片构成，传感器的输出端与模数转换芯片的输入端连接。 

由于本发明采用了上述技术方案，实现了大批量塔机运行数据到结构应力数据的快速处理，为塔机寿命评估提供了准确可靠的大样本空间数据，解决了现有的方法存在的采集数据样本容量小、采集结构部位有限的技术难点。因此，与背景技术相比，本发明具有能实时采集大批量塔式起重机运行数据且能依据采集到的大批量实时运行数据快速计算得到塔式起重机结构杆件内力、应力数据时域值的优点。 

附图说明

图1是本发明的流程框图； 

图2是本发明塔身腹杆体系示意图； 

图3是双吊点臂架支反力计算模型图； 

图4是自重载荷作用下斜面桁架杆件计算模型图； 

图5是移动载荷重用下斜面桁架杆件计算模型图； 

图6是垂直载荷作用臂架水平桁架计算简图； 

图7是风载荷作用臂架水平桁架计算简图； 

图8是臂架回转惯性载荷作用臂架水平桁架计算简图； 

图9是塔式起重机数据记录仪结构框图； 

图10是本发明计算程序总体流程图； 

图11是参数输入模块流程图； 

图12是计算模块流程图； 

图13是结果显示及输出模块流程图。 

图中：L₁、L₂和L₃--吊臂根部铰点A到第一吊点B、第一吊点B到第二吊点C、第二吊点C到吊臂端部最大幅度处的距离；L_B和L_C--拉杆BF和CF 的长度；θ_B和θ_C--拉杆BF和CF与吊臂间的夹角；A_B、A_C--拉杆BF和CF的截面积；α--臂架斜面桁架与水平桁架夹角，β--两斜面桁架夹角，γ--臂架斜面桁架腹杆与Y1轴夹角；b和h--臂架截面宽和高，h₁--斜面桁架高，l_ω--臂架节点间距；E--弹性模量，I--惯性距。 

具体实施方式

下面结合附图和实施例对本发明作进一步的详细描述。 

如图1所示，本实施例中的利用塔式起重机运行数据获取评估点应力时域值的方法，其包括下列步骤： 

1)使用数据记录仪采集塔式起重机的起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据，并将采集到的上述实时运行参数存放到存储单元中，以获得塔式起重机较长时间的运行数据； 

2)对采集的较长时间的实时运行数据进行处理，并依据速度公式： 

求出起升、变幅和回转速度，然后根据求出的两点的速度与加速度公式： 

求出起升、变幅和回转加速度，以获得各机构的速度和加速度的参数； 

3)建立塔式起重机结构计算模型，并编制用于快速计算的计算机程序； 

4)读入结构模型计算所需参数和选择载荷组合形式的工况数据，计算某时刻塔式起重机结构杆件内力和应力，得到塔机杆件内力和应力数据，并将计算结果分类存放在数据库中； 

5)将上述第4)步骤所得的某时刻应力值加上或减去修正值，得到修正后的应力值，修正值的范围为计算值的0％～15％左右； 

6)重复第4)和第5)步骤，就可以得到塔式起重机多个时刻所对应的结构内力和应力数据，把数据整理成多个评估点的应力谱数据。 

所述选择载荷组合形式是指对起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据可全部选用或选用其中的一部分。 

所述修正值的具体取值按照杆件的截面高度与杆件长度等参数来确定。 

所述计算某时刻塔式起重机结构杆件内力和应力的步骤是： 

第一步：塔式起重机结构杆件内力计算 

塔机结构种类繁多，本实施例以桁架式的塔式起重机结构为例：臂架截面形式通常为正三角形、倒三角形；臂架为水平臂架，按照吊点分为单吊点和双吊点；塔身截面为矩形或方形，腹杆形式为六种常用腹杆形式，见附图2所示。计算载荷包括结构自重载荷、起重量载荷、风载荷、臂架回转启制动惯性载荷、吊重惯性载荷和大车行走惯性载荷等；根据需要选择其中的任何载荷组合形式。塔式起重机塔身和臂架结构为空间桁架结构，桁架节点看作是铰点，杆件为二力杆件，计算方法是力学中的力法、节点法和截面法等。最后得到各种形式塔机的每一个杆件内力计算通式，建立内力计算数据库。 

1.双吊点臂架杆件内力计算 

在起升平面内，双吊点塔式起重机臂架为一次超静定结构，采用力法计算支座反力。将BF杆断开，然后代之以多余未知力X1的基本体系如附图3所示。 

建立力法方程： 

δ_B1X₁+Δ_Bx＝0 

式中：δ_B1--表示当X1＝1单独作用时，基本结构上B点沿X1方向的相对位移；Δ_Bx--表示当臂架在载荷作用时，基本结构上B点沿X1方向的相对位移。 

1)自重载荷作用下 

(1)支反力计算 

${\mathrm{\Δ}}_{\mathrm{Bq}}=\frac{q{L}_1{L}_2(L_1^2+L_2^2+3L_1{L}_2)\sin {\mathrm{\θ}}_B}{24\mathrm{EI}}-\frac{q{L}_1{L}_2{L}_3^2(2L_1+L_2)\sin {\mathrm{\θ}}_B}{12\mathrm{EI}(L_1+L_2)}$

$+\frac{q{L}_1{L}_C{(L_1+L_2+L_3)}^2\sin {\mathrm{\θ}}_B}{2{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

${\mathrm{\δ}}_{B1}=\frac{L_l^2{L}_2^2{\sin }^2{\mathrm{\θ}}_B}{3\mathrm{EI}(L_1+L_2)}+\frac{L_B}{{\mathrm{EA}}_B}+\frac{L_1^2{L}_C{\sin }^2{\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

$X_1=-\frac{{\mathrm{\Δ}}_{\mathrm{Bq}}}{{\mathrm{\δ}}_{B1}}$

拉杆拉力： 

R_Bq＝X₁

$R_{\mathrm{Cq}}=\frac{q{(L_1+L_2+L_3)}^2-2X_1{L}_1\sin {\mathrm{\θ}}_B}{2(L_1+L_2)\sin {\mathrm{\θ}}_C}$

A点竖直方向支座反力： 

$R_{\mathrm{AqY}}=\frac{-2X_1{L}_2\sin {\mathrm{\θ}}_B+q(L_1+L_2+L_3)(L_1+L_2-L_3)}{2(L_1+L_2)}$

(2)杆件内力计算 

臂架为空间正三角形桁架结构，需分解为斜面桁架和底面水平桁架计算，自重载荷作用下，斜面桁架杆件内力计算简图如附图4所示。 

臂架自重分解在斜面桁架上的均布载荷 $q=\frac{G}{3L}(\frac{1}{\sin \mathrm{\α}}+\frac{1}{2\mathrm{con}(\mathrm{\β}/2)}),$ 分解成上、下弦杆的节点载荷 $g_m=\frac{G}{6\mathrm{mcon}(\mathrm{\β}/2)}$ 和 $g_n=\frac{G}{3n\sin \mathrm{\α}}$ (m、n分别为上、下弦杆的节点数)。 

支座反力： 

$R_{\mathrm{AY}1}=\frac{R_{\mathrm{AqY}}}{2\sin \mathrm{\α}},$ R_AX＝R_BX+R_CX

$R_{\mathrm{BY}1}=\frac{R_{\mathrm{Bq}}\sin {\mathrm{\θ}}_B}{2\sin \mathrm{\α}},$ $R_{\mathrm{BX}}=\frac{R_{\mathrm{Bq}}\mathrm{con}{\mathrm{\θ}}_B}{2}$

$R_{\mathrm{CY}1}=\frac{R_{\mathrm{Cq}}\sin {\mathrm{\θ}}_C}{2\sin \mathrm{\α}},$ $R_{\mathrm{CX}}=\frac{R_{\mathrm{Cq}}\mathrm{con}{\mathrm{\θ}}_B}{2}$

杆件内力 

当0＜X＜L1时，如附图8截面2各个杆件内力，采用截面法，取截面左端计算，分别对1和2点求矩，Y1方向合力为零得： 

上弦杆： $S_0=R_{\mathrm{AY}1}\frac{X}{h_1}-\frac{q{X}^2}{2h_1}$

下弦杆： $S_u=R_{\mathrm{AY}1}\frac{X-l_{\mathrm{\ω}}/2}{h_1}-R_{\mathrm{AX}}-\frac{\mathrm{qX}(X-l_{\mathrm{\ω}}/2)}{2h_1}$

斜腹杆： $S_d=-(R_{\mathrm{AY}1}-{\mathrm{mg}}_m-{\mathrm{ng}}_n)\frac{1}{\mathrm{con\γ}}$

L1＜X＜L2和X＞L2截面杆件内力的计算方法与上面相同。 

2)移动载荷作用下 

移动载荷F_Q位于臂架AB段距离A为X(mm) 

(1)支反力计算 

${\mathrm{\Δ}}_{\mathrm{BQ}}=\frac{-F_Q{\mathrm{XL}}_2[L_1(L_1+2L_2)-X^2]\sin {\mathrm{\θ}}_B}{6\mathrm{EI}(L_1+L_2)}+\frac{F_Q{L}_1{L}_{C}X\sin {\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

${\mathrm{\δ}}_{B1}=\frac{L_l^2{L}_2^2{\sin }^2{\mathrm{\θ}}_B}{3\mathrm{EI}(L_1+L_2)}+\frac{L_B}{{\mathrm{EA}}_B}+\frac{L_1^2{L}_C{\sin }^2{\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

$X_1=-\frac{{\mathrm{\Δ}}_{\mathrm{Bq}}}{{\mathrm{\δ}}_{B1}}$

拉杆拉力 

R_BQ1＝X₁

$R_{\mathrm{CQ}1}=\frac{F_{Q}X-X_1{L}_1\sin {\mathrm{\θ}}_B}{(L_1+L_2)\sin {\mathrm{\θ}}_C}$

A点竖直方向支反力： 

$R_{\mathrm{AQ}1}=\frac{-X_1{L}_1\sin {\mathrm{\θ}}_B+F_Q(L_1+L_2-X)}{(L_1+L_2)}$

(2)杆件内力计算 

与自重载荷作用相同，需分解为斜面桁架和底面水平桁架计算，移动载荷作用下，斜面桁架杆件内力计算简图如附图5所示。 

支座反力： 

$R_{\mathrm{AY}1}=\frac{R_{\mathrm{AQ}1}}{2\sin \mathrm{\α}},$ R_AX＝R_BX+R_CX

$R_{\mathrm{BY}1}=\frac{R_{\mathrm{BQ}1}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_B,$ $R_{\mathrm{BX}}=\frac{R_{\mathrm{BQ}1}\mathrm{con}{\mathrm{\θ}}_B}{2}$

$R_{\mathrm{CY}1}=\frac{R_{\mathrm{CQ}1}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_C,$ $R_{\mathrm{CX}}=\frac{R_{\mathrm{CQ}1}\mathrm{con}{\mathrm{\θ}}_C}{2}$

小车后轮后面桁架杆件内力，取II截面左端计算，分别对1和2点求矩，Y1方向合力为零得： 

$S_0=R_{\mathrm{AY}1}\frac{X}{h_1}$

$S_u=R_{\mathrm{AY}1}\frac{X-l_{\mathrm{\ω}}/2}{h_1}-R_{\mathrm{AX}}$

$S_d=R_{\mathrm{AY}1}\frac{1}{\mathrm{con\γ}}$

小车前轮前面到B点桁架杆件、吊点B和C之间桁架杆件和跨外桁架杆件内力计算方法同小车后轮后面桁架杆件内力计算。 

移动载荷F_Q位于臂架BC段距离A为Y(mm)和移动载荷F_Q位于臂架CD段距离A为Z(mm)时计算方法同移动载荷F_Q位于臂架AB段距离A为X(mm)时。 

3)水平桁架受移动载荷和自重载荷作用 

垂直载荷作用在水平桁架时，如附图6所示。 

竖杆A-A轴向压力： 

$S_{A-A}=\frac{G_1}{4\mathrm{tg\α}}+\frac{G}{m*3\mathrm{tg\α}}$

竖杆C-C轴向最大压力： 

$S_{c-c}=\frac{Q+G_c}{4\mathrm{tg\α}}+\frac{G}{m*3\mathrm{tg\α}}$

竖杆B-B轴向压力： 

$S_{B-B}=\mathrm{\Σ}{T}_Z-\frac{G}{m*3\mathrm{tg\α}}$ (∑T为垂直载荷在吊点B产生水平力) 

4)水平桁架受风载荷作用 

风载荷作用下的计算简图如附图7所示，臂架所受风载荷Pw(单位为N)， 按照均布的水平节点力 $g_s=\frac{p_w}{m}$ (m为节点数)作用与水平桁架一侧弦杆的各节点上。 

支反力： 

R_AZ＝mg_s $R_{\mathrm{AX}}={\mathrm{mg}}_s\frac{L}{2b}$

截面2杆件内力，取截面左端计算，分别对节点求矩得： 

下弦杆1： $S_{u1}=R_{\mathrm{AZ}}\frac{X}{b}-R_{\mathrm{AX}}-g_s\frac{X(X+l_w)}{2b{l}_w}$

下弦杆2： $S_{u2}=R_{\mathrm{AZ}}\frac{X-l_w}{b}-R_{\mathrm{AX}}-g_s\frac{X(X-l_w)}{2b{l}_w}$

水平斜腹杆内力： 

$S_d=(g_s\frac{X}{l_w}-R_{\mathrm{AZ}})/\cos \mathrm{\τ}$ (如2截面方向腹杆) 

$S_d=-(g_s\frac{X}{l_w}-R_{\mathrm{AZ}})/\cos \mathrm{\τ}$ (如1截面方向腹杆) 

5)水平桁架受臂架回转惯性载荷 

回转启制动时，臂架本身的水平惯性力引起的杆件内力计算如附图8所示。 

设臂架重量为G，微段dX的水平惯性力为： 

$\mathrm{dp}=\frac{G}{L^{\′}g}(L_i+X)*\frac{\mathrm{\ω}}{t}\mathrm{dX}$

由∑Z＝0得： 

$R_{\mathrm{AZ}}={\∫}_0^{L^{\′}}\mathrm{dp}=\frac{\mathrm{G\ω}}{2\mathrm{gt}}(L^{\′}+2l)$

由∑M_AZ＝0得： 

$R_{\mathrm{AX}}={\∫}_0^{L^{\′}}\mathrm{dp}*X=\frac{\mathrm{G\ω}{L}^{\′}}{\mathrm{Bgt}}(\frac{L^{\′}}{3}+\frac{l}{2})$

任意截面2，取右端各力对节点求矩可得： 

$S_{u2}=\frac{1}{b}{\∫}_X^{L^{\′}}\frac{G}{L^{\′}g}(X+l+u)*\frac{\mathrm{\ω}}{t}*(u+l_{\mathrm{\ω}})*d_u$

$=\frac{G*\mathrm{\ω}*(L^{\′}-X)}{b*g*L^{\′}*t}[(X+l)l_w+(l_w+X+l)\frac{L^{\′}+X}{2}+\frac{L^{\′2}+L^{\′}X+X^2}{3}]$

$S_{u1}=\frac{1}{b}{\∫}_X^{L^{\′}}\frac{G}{L^{\′}g}(X+l+u)*\frac{\mathrm{\ω}}{t}*u*d_u$

$=\frac{G*\mathrm{\ω}*(L^{\′}-X)}{b*g*L^{\′}*t}[\frac{(X+l)(L^{\′}+X)}{2}+\frac{L^{\′2}+L^{\′}X+X^2}{3}]$

取截面2左端，由∑Z＝0，计算斜腹杆内力： 

S_d＝-R_AZ/cosτ 

取截面1左，同样可得1截面方向腹杆内力： 

S_d＝R_AZ/cosτ 

2.单吊点臂架杆件内力计算 

单吊点臂架结构起升平面(即垂直平面)，作为伸臂梁计算；回转平面(即水平面)，作为悬臂梁计算。与双吊点臂架相同，先求支座反力，然后用节点法和截面法计算各个杆件内力。 

其它截面形式或腹杆形式臂架截面杆件内力计算方法与上面相同。 

3.塔身杆件内力计算 

计算的塔身截面为矩形或方形，腹杆形式为六种常用腹杆形式。具体方法是：首先将塔身空间桁架结构分解为四个平面桁架结构进行计算；接着，将塔身所受不平衡弯矩、风载荷、轴向载荷和扭矩分解到平面桁架；然后对平面桁架杆件进行计算；最后，将计算结果叠加得出塔身每个杆件的内力计算公式。 

其中主肢内力计算通式如下： 

$S_1=\frac{M_X+M_Y}{2b}-\frac{N}{2}+\frac{(g_{s1}+g_{s4})(m+1)m{l}_w}{4b}+\frac{(g_{s1}+g_{s4})(1+2Z/l_{\mathrm{\ω}}-2m)Z}{2b}$

$S_2=\frac{M_X-M_Y}{2b}+\frac{N}{2}+\frac{(g_{s1}-g_{s2})(m+1)m{l}_w}{4b}+\frac{(g_{s1}-g_{s2})(1+2Z/l_{\mathrm{\ω}}-2m)Z}{2b}$

$S_3=\frac{M_X+M_Y}{2b}-\frac{N}{2}+\frac{(g_{s2}+g_{s3})(m+1)m{l}_w}{4b}-\frac{(g_{s2}+g_{s3})(1+2Z/l_{\mathrm{\ω}}-2m)Z}{2b}$

$S_4=\frac{{-M}_X+M_Y}{2b}+\frac{N}{2}+\frac{(g_{s3}-g_{s4})(m+1)m{l}_w}{4b}-\frac{(g_{s3}-g_{s4})(1+2Z/l_{\mathrm{\ω}}-2m)Z}{2b}$

H--塔身高，单位为mm；Z--截面到支点距离，单位为mm；m--塔身单片桁架一侧节点个数；g_s、g_s1、g_s2、g_s3和g_s4--为节点风载荷，单位均为N；b为塔身截面宽，单位为mm；1标准节长，单位为mm。 

按照上面的计算过程，可以得到不同结构形式的塔机杆件的内力计算通式。 

第二步：结构杆件应力计算 

由上步计算得到的内力，经过轴心受力杆件的应力计算公式为： 

式中：N--计算的轴向力，单位为牛；A_j--杆件验算截面的净面积，单位为mm²。σ_修正为应力修正值，其具体取值按照杆件的截面高度与杆件长度等参数来确定。 

根据上述内力计算通式，以及塔机金属结构型钢参数可以得到各个杆件应力计算通式。 

第三步：编制软件 

本实施例软件编制采用VB计算机语言和Access数据库，并结合Excel工具。本系统主要分为参数输入模块、杆件应力计算模块、结果显示及输出模块三个模块。通过参数输入模块输入塔机结构参数和塔机运行数据载荷，然后将这些数据通过数据访问对象(ActiveX Data Objects，简称ADO)接口保存在Access数据库中；接着判断计算模型，然后对塔机杆件计算，当选择计算时，塔机杆件内力计算模块通过ADO接口从Access数据库中调用数据，并将计算结果保存在数据库中，并且在文本文件中备份；结果显示模块则由OLEDB接口调用数据库中的计算结果，在屏幕上显示出来，并且通过Microsoft ExcelObjects Library调用Excel绘图功能，将其结果以曲线图形显示。整个程序的系统流程如图10所示。 

(1)参数输入模块 

塔机杆件计算所需要的数据繁多，其中结构参数包括：塔机臂架、塔身、塔帽、平衡臂、拉杆等形式及尺寸；载荷参数有包括：起升重量、变幅位置、起升高度、回转位置、风速和风向以及起升机构、变幅机构和回转结构的运行速度和加速度；同时随着程序运行，数据量还在不断增加。为了避免数据重复输入，本程序系统以塔机类型和数据类型为单位进行数据管理，并且通过公用数据子集实现数据共享。在运行程序开始计算前，应该根据计算需要，填写以下各数据表格。 

塔机基本参数表、臂架参数表、拉杆参数表、塔身参数表、结构型钢和塔机运行数据载荷表。 

其程序流程如图11所示，首先选择输入方式，若选择直接输入，进行下一步输入原始数据，若选择数据库数据，从数据库中调取所需数据；接着对数据选择修改，检查是否错误，若有错，返回选择输入方式，若没有错，保存数据。 

(2)杆件计算模块 

塔机杆件计算模块是按照工况分为单工况处理和批量处理；按照计算内容分为内力计算模块和应力计算模块；按照杆件类型分为塔身杆件内力应力计算模块、臂架杆件内力应力计算模块和拉杆内力应力计算模块；按照作用载荷臂架杆件内力计算模块又分为自重载荷作用模块、移动载荷作用模块、风载荷作用模块、臂架回转惯性载荷作用模块、吊重回转惯性载荷作用模块、大车水平启制动惯性载荷模块等。在计算过程中，可以选择其中一种类型进行计算，也可以选择几种类型的组合进行计算，得到所需要的结果。该模块流程图如12所示，首先选择塔机类型及结构形式，若是单工况，读入单工况的载荷、选择组合形式并计算杆件的内力和应力；若是多工况，读入多工况的载荷、选择组合形式并计算杆件的内力和应力。 

(3)结果显示及输出模块 

计算结果的显示采用子模块的方式。与结果相对应建立塔身杆件内力模块、塔身杆件应力模块、臂架杆件内力模块、臂架杆件应力模块等。可以将塔身杆件内力、塔身杆件应力、臂架杆件内力、臂架杆件应力以文本和图形显示。 图形显示用Excel绘图功能，将其结果绘制图形并显示。显示结果按照塔机运行数据的时间历程和塔机结构杆件的编号分类。具体流程如图13所示，首先选择显示内容，接着从数据库中读取计算结果，然后选择显示方式，最后显示查看结果。 

如图9所示，所述数据记录仪包括用于监控塔式起重机参数信息的传感器采集装置1、微处理单元(STC89LE516RD)2、数据存储单元3和通讯网络接口单元4；参数信息的传感器采集装置1的信号与微处理单元2的信号输入接口连接；微处理单元2的数据输出接口连接数据存储单元3的输入接口；通讯网络接口单元4的输入端连接微处理单元2的通讯输出接口。所述参数信息的传感器采集装置1由传感器1a和模数转换芯片1b构成，传感器1a的输出端与模数转换芯片1b的输入端连接。 

下面以某型号塔式起重机为例，作进一步的说明。 

以某型号塔式起重机为例，计算塔机受任意外载荷组合下，任意杆件的内力和应力。下面是塔机臂架受自重载荷和移动载荷作用下，距离回转中心4.5m截面上弦杆、下弦杆和斜腹杆的应力计算。塔机具体结构参数如表1所示(长度为mm；重量为kg)。 

表1 

由表1可以计算得到： 

L₁＝13.125m；L₂＝25m；L₃＝11.875m；L_B＝15.42m、L_C＝39.58m；sinθ_B＝0.4086、sinθ_C＝0.1592；A_B＝B_C＝0.002375m²；h₁＝1.096m；b×h＝1.19m×0.92m；E＝2.06×10⁵MPa。 

1.距离回转中心4.5m截面杆件内力计算 

1)自重载荷作用下 

臂架自重：G＝35000N；均布载荷： $q=\frac{G}{L}=\frac{35000}{50}=7000N/m$

(1)臂架支反力计算： 

双吊点臂架为一次超静定结构，如图3所示，建立力法方程：δ_B1X₁+Δ_Bx＝0，计算支点反力。 

${\mathrm{\Δ}}_{\mathrm{Bq}}=\frac{q{L}_1{L}_2(L_1^2+L_2^2+3L_1{L}_2)\sin {\mathrm{\θ}}_B}{24\mathrm{EI}}-\frac{q{L}_1{L}_2{L}_3^2(2L_1+L_2)\sin {\mathrm{\θ}}_B}{12\mathrm{EI}(L_1+L_2)}$

$+\frac{q{L}_1{L}_C{(L_1+L_2+L_3)}^2\sin {\mathrm{\θ}}_B}{2{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

${\mathrm{\δ}}_{B1}=\frac{L_l^2{L}_2^2{\sin }^2{\mathrm{\θ}}_B}{3\mathrm{EI}(L_1+L_2)}+\frac{L_B}{{\mathrm{EA}}_B}+\frac{L_1^2{L}_C{\sin }^2{\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

$X_1=-\frac{{\mathrm{\Δ}}_{\mathrm{Bq}}}{{\mathrm{\δ}}_{B1}}$

代入上面由表1计算的参数可得：X₁＝3.35×10⁴N 

拉杆拉力： 

R_Bq＝X₁＝3.35×10⁴N 

$R_{\mathrm{Cq}}=\frac{q{(L_1+L_2+L_3)}^2-2X_1{L}_1\sin {\mathrm{\θ}}_B}{2(L_1+L_2)\sin {\mathrm{\θ}}_C}=11.17\×{10}^{4}N$

A点竖直方向支座反力： 

$R_{\mathrm{AqY}}=\frac{-2X_1{L}_2\sin {\mathrm{\θ}}_B+q(L_1+L_2+L_3)(L_1+L_2-L_3)}{2(L_1+L_2)}=0.7\×{10}^{4}N$

(2)杆件内力计算 

臂架为空间正三角形桁架结构，需分解为斜面桁架和底面水平桁架计算，自重载荷作用下，斜面桁架杆件内力计算简图如图4所示。 

臂架自重分解在斜面桁架上的均布载荷 $q=\frac{G}{3L}(\frac{1}{\sin \mathrm{\α}}+\frac{1}{2\mathrm{con}(\mathrm{\β}/2)})=417N,$ 分解成上、下弦杆的节点载荷 $g_m=\frac{G}{6\mathrm{mcon}(\mathrm{\β}/2)}=174N,$ $g_n=\frac{G}{3n\sin \mathrm{\α}}=339N$ (m、n分别为上、下弦杆的节点数)。 

斜面桁架支座反力： 

$R_{\mathrm{AY}1}=\frac{R_{\mathrm{AqY}}}{2\sin \mathrm{\α}}=4172N;$ R_AX＝R_BX+R_CX

$R_{\mathrm{BY}1}=\frac{R_{\mathrm{Bq}}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_B=8157N;$ $R_{\mathrm{BX}}=\frac{R_{\mathrm{Bq}}\mathrm{con}{\mathrm{\θ}}_B}{2}=15290N$

$R_{\mathrm{CY}1}=\frac{R_{\mathrm{Cq}}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_C=105975N;$ $R_{\mathrm{CX}}=\frac{R_{\mathrm{Cq}}\mathrm{con}{\mathrm{\θ}}_B}{2}=55800N$

距离回转中心4.5m截面杆件内力： 

上弦杆： $S_0=R_{\mathrm{AY}1}\frac{X}{h_1}-\frac{q{X}^2}{2h_1}=11600N$

下弦杆： $S_u=R_{\mathrm{AY}1}\frac{X-l_{\mathrm{\ω}}/2}{h_1}-R_{\mathrm{AX}}-\frac{\mathrm{qX}(X-l_{\mathrm{\ω}}/2)}{2h_1}=-61400N$

斜腹杆： $S_d=-(R_{\mathrm{AY}1}-{\mathrm{mg}}_m-{\mathrm{ng}}_n)\frac{1}{\mathrm{con\γ}}=-3030N$

2)移动载荷载荷作用下 

(1)臂架支反力计算 

移动载荷F_Q位于臂架AB段距离A为X(mm)时，建立力法方程：δ_B1X₁+Δ_Bx＝0，计算支点反力。 

${\mathrm{\Δ}}_{\mathrm{BQ}}=\frac{-F_Q{\mathrm{XL}}_2[L_1(L_1+2L_2)-X^2]\sin {\mathrm{\θ}}_B}{6\mathrm{EI}(L_1+L_2)}+\frac{F_Q{L}_1{L}_{C}X\sin {\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

${\mathrm{\δ}}_{B1}=\frac{L_1^2{L}_2^2{\sin }^2{\mathrm{\θ}}_B}{3\mathrm{EI}(L_1+L_2)}+\frac{L_B}{{\mathrm{EA}}_B}+\frac{L_1^2{L}_C{\sin }^2{\mathrm{\θ}}_B}{{\mathrm{EA}}_C{(L_1+L_2)}^2{\sin }^2{\mathrm{\θ}}_C}$

$X_1=-\frac{{\mathrm{\Δ}}_{\mathrm{Bq}}}{{\mathrm{\δ}}_{B1}}$

代入上面由表1计算的参数及小车变幅X＝9可得：X₁＝7.65×10⁴N 

拉杆拉力： 

R_BQ1＝X₁＝7.65×10⁴N 

$R_{\mathrm{CQ}1}=\frac{F_{Q}X-X_1{L}_1\sin {\mathrm{\θ}}_B}{(L_1+L_2)\sin {\mathrm{\θ}}_C}=5.42\×{10}^{4}N$

A点竖直方向支反力： 

$R_{\mathrm{AQ}1}=\frac{-X_1{L}_1\sin {\mathrm{\θ}}_B+F_Q(L_1+L_2-X)}{(L_1+L_2)}=1.98\×{10}^{4}N$

(2)杆件内力计算 

与自重载荷作用相同，需分解为斜面桁架和底面水平桁架计算，移动载荷作用下，斜面桁架杆件内力计算简图如图5所示。 

支座反力： 

$R_{\mathrm{AY}1}=\frac{R_{\mathrm{AQ}1}}{2\sin \mathrm{\α}}=11800N,$ R_AX＝R_BX+R_CX

$R_{\mathrm{BY}1}=\frac{R_{\mathrm{BQ}1}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_B=18630N,$ $R_{\mathrm{BX}}=\frac{R_{\mathrm{BQ}1}\mathrm{con}{\mathrm{\θ}}_B}{2}=34900N$

$R_{\mathrm{CY}1}=\frac{R_{\mathrm{CQ}1}}{2\sin \mathrm{\α}}\sin {\mathrm{\θ}}_C=13200N,$ $R_{\mathrm{CX}}=\frac{R_{\mathrm{CQ}1}\mathrm{con}{\mathrm{\θ}}_C}{2}=26700N$

距离回转中心4.5m截面杆件内力： 

$S_0=R_{\mathrm{AY}1}\frac{X}{h_1}=40370N$

$S_u=R_{\mathrm{AY}1}\frac{X-l_{\mathrm{\ω}}/2}{h_1}-R_{\mathrm{AX}}=-29000N$

$S_d=R_{\mathrm{AY}1}\frac{1}{\mathrm{con\γ}}=13590N$

3)自重载荷和移动载荷作用下杆件内力 

上弦杆：S_上＝51970N 

下弦杆：S_下＝-90400N 

斜腹杆：S_腹杆＝10560N 

2.杆件应力计算 

杆件主要是承受轴向力，为受拉或受压的轴心受力杆件，轴心受力杆件的应力计算公式为： 

式中：N--计算的轴向力，单位为牛；A_j--杆件验算截面的净面积，单位为mm²；σ_修正为应力修正值，其取值为σ_修正＝0.0Mpa。 

由杆件内力和表1中钢材型号截面积可得： 

上弦杆：S_上＝51970N 

下弦杆：S_下＝-90400N 

斜腹杆：S_腹杆＝16620N 

上面计算中载荷如表2所示，得到距离回转中心4.5m杆件应力如表3所示： 

表2 

自重载荷(kg)	起升重量(kg)	变幅位置(m)
			3500	4000	9

表3 

上弦杆(MPa)	下弦杆(MPa)	斜腹杆(MPa)
			15.5	40.0	28.1

3.批量载荷数据处理 

按照上面方法可以计算出任意载荷组合下，任意截面，任意杆件的内力和应力。对于批量载荷数据快速处理，利用上面计算过程开发的快速计算软件进行计算。如表4所示批量载荷数据。 

表4 

工况序号	起升重量   /(kg)	变幅位置   /(m)	起升高度   /(m)	回转角度   /(D)	风速   /(m/s)	风向   /(D)
							1	5	27.9	10.3	208	4.2	22
2	671	27.9	9.9	208	4.3	22
							3	2235	27.9	9.3	205	4.2	21
4	2215	25.8	8	150	4.2	23
							5	786	25.7	9.9	18	3.5	22

6	317	25.7	10	18	3.5	22
							7	71	25.7	10	18	3.6	22
8	3	22.8	7.7	214	7.5	60
							9	432	22.8	7.6	214	7.4	61
10	3168	22.8	7.4	215	7.5	63
							11	3166	16	2.2	215	7.5	59
12	3217	13	5.9	205	7.1	50
							13	3177	12.9	7.4	204	7.3	60
14	742	12.9	7.5	204	7.5	62
							15	346	12.9	7.5	204	7.2	61
16	5	48.1	15.8	186	4.3	107
							17	837	48.2	14	187	4.1	106
18	1129	48.1	13.5	186	4.2	105
							19	1213	35.6	9.6	186	4.3	107
20	1169	26.6	12.2	187	4.1	100
							21	1178	26.6	13.4	182	4.2	99
22	1235	26.5	13.6	182	4.2	108
							23	927	26.6	14.4	182	4.3	110
24	70	26.5	14.8	182	2.4	25
							25	1328	11.2	8.4	229	2.2	24
26	1347	15.8	2.7	203	2.3	23
							27	1393	24.8	8.5	172	2.2	24
28	4	25.8	8.4	172	2.5	25

按照图11所示输入表1塔机结构参数，并读入表4塔机运行载荷参数，并将其保存在Access数据库中。 

根据塔机结构参数，判读塔机结构计算模型。按照图12所示计算塔机杆 件内力和应力。 

按照图13所示输出计算结果。其中距离塔机回转中心3.8m处臂架截面杆件轴应力随工况变化如表5所示。每一杆件应力随工况变化值都可以文本或图形方式显示。 

表5 

工况序号	下弦杆   /(MPa)	上弦杆   /(MPa)	底斜杆   /(MPa)	斜腹杆   1/(MPa)	斜腹杆   2/(MPa)
						1	-41.3	-3.6	-7.1	-2.9	-1.3
2	-48.3	-2.8	-8.2	-2.1	-0.1
						3	-64	-0.9	-10.9	-0.34	2.1
4	-57.2	-1.3	-14.4	3.2	4.4
						5	-53.2	-2.4	-8.3	-0.3	-2
6	-47.3	-3.3	-7.7	-1	-2.6
						7	-43.2	-3.5	-7.3	-1.3	-2.7
8	-43	-3.7	-7.1	-2.9	-1.3
						9	-46	-3.2	-7.7	-1.4	-0.6
10	-66.6	-1	-11.3	-0.5	-2
						11	-53.7	-6.5	-9.1	-6.4	-4.4
12	-46.3	-10.6	-7.9	-10.8	-9.1
						13	-46.3	-10.6	-7.9	-10.8	-9.1
14	-45	-7.5	-7.3	-6	-5.2
						15	-44	-5	-7.2	-3.9	-3.1
16	-43	-3.7	-7.1	-2.9	-1.3
						17	-55.8	-5.5	-9.5	-4.9	-2.9
18	-59.9	-5.3	-10.2	-5.3	-3
						19	-56.2	-2.9	-9.6	-2.4	-0.3
20	-51.8	-2.4	-8.8	-1.7	0.3

21	-51.8	-2.4	-8.8	-1.7	0.3
						22	-51.8	-2.4	-8.8	-1.7	0.3
23	-49	-2.9	-8.2	-1.9	0
						24	-44	-3.5	-7.4	-2.4	-1
25	-38.9	-9	-6.7	-6	-7.4
						26	-45.1	-5.4	-7.7	-4.9	-3.2
27	-53.3	-2.2	-9.1	-1.6	0.44
						28	-41.3	-3.6	-7.1	-2.9	-1.3

Claims (4)

1.一种利用塔式起重机运行数据获取评估点应力时域值的方法，其特征是包括下列步骤：

1)使用数据记录仪采集塔式起重机的起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据，并将采集到的上述实时运行数据存放到存储单元中，以获得塔式起重机较长时间的运行数据；

2)对采集的较长时间的实时运行数据进行处理，并依据速度公式：

求出起升、变幅和回转速度，然后根据求出的两点的速度与加速度公式：

求出起升、变幅和回转加速度，以获得起升机构、变幅机构和回转机构的速度和加速度的参数；

3)建立塔式起重机结构计算模型，并编制用于快速计算的计算机程序；

4)读入结构模型计算所需参数和选择载荷组合形式的工况数据，计算某时刻塔式起重机结构杆件内力和应力，得到塔机杆件内力和应力数据，并将计算结果分类存放在数据库中；

5)将上述第4)步骤所得的某时刻应力值加上或减去修正值，得到修正后的应力值，修正值的范围为计算值的0％～15％；

6)重复第4)和第5)步骤，就可以得到塔式起重机多个时刻所对应的结构内力和应力数据，把数据整理成多个评估点的应力谱数据。

2.根据权利要求1所述的利用塔式起重机运行数据获取评估点应力时域值的方法，其特征是：所述选择载荷组合形式是指对起重量、起升高度、变幅位置、臂架回转位置、风速和风向的实时运行数据全部选用或选用其中的一部分。

3.根据权利要求1所述的利用塔式起重机运行数据获取评估点应力时域值的方法，其特征是：所述数据记录仪包括用于监控塔式起重机参数信息的传感器采集装置、微处理单元、数据存储单元和通讯网络接口单元；塔式起重机参数信息的传感器采集装置的信号与微处理单元的信号输入接口连接；微处理单元的数据输出接口连接数据存储单元的输入接口；通讯网络接口单元的输入端连接微处理单元的通讯输出接口。

4.根据权利要求3所述的利用塔式起重机运行数据获取评估点应力时域值的方法，其特征是：所述塔式起重机参数信息的传感器采集装置由传感器和模数转换芯片构成，传感器的输出端与模数转换芯片的输入端连接。

Images