等值刚度法管道静力计算程序及其使用

等值刚度法管道静力计算程序及其使用编制：蔡立颖校对：施杰峰审核：时伟杭州锅炉集团有限公司（杭州锅炉厂）二零零三年一月1前言：等值刚度法属力法的计算范畴，与之对应的还有位移法（即有限元法）。根据电力部标准《火力发电厂汽水管道应力计算技术规定》（SDGJ6-78）编制，只能计算树状分支；如遇环形，用刚性支点分割。原理是将两个分支等效为一个分支，依此类推将多分支等效为单支结构。（一）计算程序的编制本计算程序是根据“管道应力计算技术规定”的要求编制的，计算内容与计算原则如下：（1）管道在工作状态下，由内压和持续外载（包括自重）作用产生的一次应力验算，并给出工作状态下各支点的荷重。（2）管道由热胀、冷缩和其它位移受约束而产生的二次应力计算，按计算管道沿坐标轴的全补偿值和钢材在20oC时的弹性模数计算，并考虑弹簧的附加力的影响。（3）管道在工作状态下对设备（或端点）的推力和力矩计算。考虑热胀、端点附加位移、有效冷紧和自重，按钢材在计算温度下的弹性模数计算。（4）管道在冷状态下对设备（或端点）的推力和力矩计算。考虑冷紧、自重和弹簧附加力，按钢材在20oC时的弹性模数计算。并必须比较运行初期冷状态下的力和力矩同应变自均衡后的自拉力和力矩，取其大者（绝对值）作为管道在冷状态下对设备（或端点）的推力和力矩值。（5）管道由冷状态到工作状态的热位移的计算。按计算管道沿坐标的全补偿值和刚才在20oC时的弹性模数计算，并考虑弹簧附加力的影响。根据各支吊点的热位移和工作载荷，选择弹簧并给出弹簧的予压缩值。（6）管道由冷紧和弹簧附加力作用的位移值计算，以作为管道支吊架设计和安装调整的一个依据。按冷紧值、自重和钢材在20oC时的弹性模数计算，并考虑冷紧口位置和弹簧附加力的影响。根据计算内容和等值刚度法的特点，计算程序分为四个工况：工况1，简称“自重工况。进行自重分配，计算弹簧支吊架的工作载荷，同时计算出管道分支的力，以作为一次应力验算的依据。工况2，简称“初热”工况。计算管道运行初期在工作状态下对设备（或端点）的推力，同时算出刚性支吊架的荷重。工况3，简称“由热自冷”工况。进行热胀二次应力验算，并根据工况1算出的自重作用力和力矩，进行一次应力验算，以及一次加二次应力验算；这个工况还计算出弹簧吊架的安装荷重；根据应变自均衡的情况，求出应变自均衡后冷状态管道对设备（或端点）的推力和力矩，以及这个状态下刚性支吊架的荷重。在这个工况中，还给出各支吊点的热位移并选择支吊架弹簧。工况4，简称“初冷”工况。计算管道安装状态下对设备（或端点）的推力和力矩以及刚性支吊架的安装荷重，同时给出弹簧支吊架的冷紧位移值。程序总框图如图1-26所示。在总框图中：SR框——输入原始数据并打印原始数据。SSJB框——形成GXB（管种计算系数表）、GDB（管端表）、GDWYB（管端位移表）、LJB（连接表）、LLB（冷拉表）、YJB（元件表）、SJTJB（数据统计表），有给定荷重时，把值送入PCB（弹簧荷重表）。SFZYJ框——形成THB（弹簧表）、BAB（坐标旋转矩阵表）、BAZB（Ba指示表）、MZBB（未坐标表）、PWSB（自重等值表）、DTB（热膨胀位移表）、DWB（自重位移表）、FB（分支集中荷载柔度矩阵表）。SEFZ框——算出每一个分支的始端力P和未端位移D。XTH框——选弹簧，形成YSB（约束表）。SYL框——计算一次应力、二次应力和一次加二次应力，形成YLB（约束类型表），并打印应力表。STL框——算推力，形成PGB（分支始端力表）、PCB（弹吊荷重表）、PMTB（对管端热态推力表）、PMOB（对管端冷态推力表）、PMTSB（对管端热松弛推力表）、PMOSB（对管端冷松弛推力表）、PTB（热推力表）、POB（冷推力表）、POSB（冷松弛推力表）、LWYB（冷位移表）、RWYB（热位移表）。DY框——打印支吊架表、推力表和冷拉表。2KS条件图1-26SRSSJBSFZYJ1GkGk4SGkKZ算下一道题SEFZSPMDIGk=SkKZGJYSS0XTH合格SYLSTLGk=4DYGk+1--GkGk=3SZ3（二）计算程序的使用范围用等值刚度法编制的这个多分支管道静力计算程序，可用于TQ—16型电子计算机，计算“树形”多分支管道。同时，必须注意管道的一个分支节点只能直接连接三个分支。当一个分支节点上连接四个及以上分支时，其连接点处可错开一个微小直元件。由于TQ—16型电子计算机储存容量的限制，管道算题的容量应控制在下列范围内。（1）原始数据数（SJS）1500；（2）管道分支数（FZS）95；（3）管子品种数（GZS）30；（4）需校验管端坐标的管种数（GDS）20；（5）有附加位移的管端数（GDMY）20；（6）元件数（YJS）300；（7）约束分支数（支吊架数——YSS）80；（8）冷紧口数（LLS）20；（9）元件改变走向次数（Ba）50；（10）环内分支数（DFZS）20；（11）环节点数（JDS）10。（三）原始数据的填写方法1．原始数据填写前的准备工作1）绘制计算管道的单线立体图。2）选取基准坐标系。基准坐标系用右旋直角坐标系（OXYZ—即四指由X--Y，拇指方向为Z；四指由Y--Z，拇指方向为X；四指由Z--X，拇指方向为Y）。坐标系的原点O可以选在任意点上（一般宜选在计算管道的某一端点），规定Y轴正方向为垂直向上，一般取X轴为主厂房纵向的水平轴，Z轴为主厂房横向的水平轴。3）划分管道分支。管道分支一般以分支节点来划分，分支点可以是管道分叉点、支吊点、冷紧点、端点、不同管材管子的连接点，或者管道上指定的任一验算点，凡是分支点与节点之间的管道或支吊架作为一个分支。分支划分以后，对分支节点编号，节点号可取0至999中的任一正整数，但不应重号。每一个分支，不能以弧元件开始或结束。如遇此种情况，则要在其前或后加一个微小直元件，譬如0.00001，这样并不影响计算结果的正确性。同样，两个弧元件相连时，中间亦须加一个微小元件，其方向为两弧的公切线方向。三通件不能直接与刚性元件相连，末梢分支的元件不能全部是刚性元件，如遇这两种情况，亦可加一个微小元件。4）元件的划分。计算管道中的元件，按其性质、形式、形状和管子规格来划分：直管、弯管、（包括热压弯管、焊接弯管）、热挤压三同、各类焊制三通件为柔性元件；阀门、对焊法兰、铸钢三通、铸钢弯头等刚度大的元件作为刚性元件（铸钢弯头用两个相交的刚性直元件代替）。元件划分以后，标出每个元件在各坐标轴上的投影长度（或角度）。对于大小头，一般可按大头管子计算。2．原始数据填写的格式原始数据填写分四个数据段，顺序按下列四张表格填写1）第一数据段表1-1名称符号单位备注设备通道号LP逻辑设备，随意，一般填4指打印机等计算年份YY四位整型数计算月、日M．D月占两位整数，日占两位整数42）第二数据段表1-2符号备注KZ（1）填1表示取消冷紧，要冷紧填0KZ（2）填1表示输出GXB，每个分支始、末端以及末端位移KZ（3）填1印FiEKiEMDiDiDWiDtiKZ（4）填1印YJBLJBKZ（5）填1印YZBZZBMZBBαFWJFJPWJKZ（6）填1不选弹簧KZ（7）填1打印选弹簧过程KZ（8）填1打印中间点力和应力KZ（9）任填，暂空KZ（10）填1表示弹簧归和选取第一版本填2表示弹簧归和选取第二版本填3表示弹簧归和选取第三版本，一般用此版。表1-2说明：此表主要功能有二，一是控制二级输出，KZ（1）--KZ（8）填1，如填0则表示一级正常输出；二是控制弹簧选用版本，程序中根据KZ（10）取1、2、3可分别采取三个版本之一，三种版本的弹簧规格见表2-1，2-2，2-3。表2-1第一版本弹簧规格表（第II类max=70mm）NoPmax（公斤）自由高度Ha（cm）柔度FT2NoPmax（公斤）自由高度Ha（cm）柔度FT2119712.60.072218205022.50.003421219715.80.035519242026.80.002891329216.80.024020336027.60.002081451420.60.013621435028.10.001611581518.40.0085922560029.00.0012516115521.60.0060623649027.20.0010817156220.30.0044824717529.30.00098表2-2第二版本弹簧规格表（第II类max=70mm）NoPmax（公斤）自由高度Ha（cm）柔度FT2NoPmax（公斤）自由高度Ha（cm）柔度FT21011431370.849510922022380.03181022521480.227811028342560.02471033981660.175911133602830.02081045661750.123711239232800.01701057661840.091411347522980.01471069971910.070211455623160.012010712962180.054011564693380.010810816372200.042811679733580.0088表2-3第三版本弹簧规格表NoPmax（公斤）自由高度Ha（cm）柔度FT2NoPmax（公斤）自由高度Ha（cm）柔度FT2101501400.1511115752300.00476102771500.097411221002500.003571031181600.063611328002500.002681041821600.041211437502600.002001052821600.026811550003200.001501063751800.020911666603300.001131075001800.015011788803700.0008451086661900.0113118118004300.0006361098882000.00845119157504400.00040611011802000.00636120210004600.00035753）第三数据段表1-3名称符号单位备注题目号TMH用正整数的阿拉伯数字表示，现没有此行。数据数SJS填以下数据的总个数，现可任意。计算安装温度Taz℃最大允许弹簧串联数MaxCLS见说明①，上限为3。荷重变化系数Kp一般取0.35另一种为0.25适用范围小。管件标记GJ见表1-3说明②。填0表示要输入一种管子的下面12个参数。计算压力Pjkg/cm2计算温度Tj℃计算温度下的线膨胀系数αj10-4cm/m℃(=10-6mm/mm.℃)钢材在20℃时的弹性模数E20106kg/cm2=(KN/mm2/0.098)x10-3钢材在计算温度下的弹性模数Ej106kg/cm2凡GJ=0的管件，均应填满该13个参数。所有GJ=0的管件，应先填；然后才能填GJ=1、2、…的管件管子外径Dwmm管子壁厚Smm弯管弯曲半径Rm管道单位重量（含水、汽、保温）WGkg/m钢材在计算温度下基本许用应力[σ]jkg/mm2钢材在20℃时的基本许用应力[σ]20oCkg/mm2环向焊缝系数h管件标记GJ1见表1-3说明②管种号（GH）取其自然顺序号。从GJ=0算起，所有GJ=0~9管件的编号，第一个管件为1，第二个管件为2，依此类推，GJ=0~9混编。以便在后边的FZBT中指定管段的归属情况。凡GJ=1的管种应填满该6个参数，没有则不填。主管所属管种号GH管子外径Dwmm管子壁厚Smm弯管弯曲半径Rm管道单位重量WGkg/m管件标记GJ2见表1-3说明②同上表凡GJ=2的管种应填满该3个参数，没有则不填。主管所属管种号GH弯管弯曲半径Rm管件标记GJ3~8见表1-3说明②同上表凡2＜GJ＜9的管种应填满该6个参数，没有则不填。主管所属管种号GH1支管所属管种号GH2壁厚S1mm热挤压三通过渡区平均半径r1m管道单位重量WGkg/m管件标记GJ凡GJ=9