精确时程分析在管道水(汽)锤计算中的应用探讨
李琳*山东电力工程咨询院有限公司济南250013
曲顺利山东省化工规划设计院济南250013
摘要在基于应力计算软件如国际通用的CAESAR和水利计算软件如PIPENET、AFT等被普遍应
用的基础上,介绍对管道进行水(汽)锤分析的一种更为精确合理的方法-时程分析。通过实例,详细介绍时程分析的过程,并与现在普遍应用的响应谱分析方法结果进行比较。
关键词管道应力分析时程分析响应谱分析
工业管道在运行中,由于阀门突然关闭或开启、泵等设备突然停止等引发的管道振动、发生啸声,流体力学研究中称这种现象为水(汽)锤现象。水(汽)锤是由于管道内部介质流动状态突然改变,以流体中的声速沿管道轴向传播形成压力波,压力波在传播中,由于管系中的压力不平衡将在管道轴向产生冲击力。这类冲击力的瞬时值可能达到惊人的数值,作用于管系、支架和设备接口,易形成管系共振、支架损坏等后果,影响系统的安全运行,所以对主要管道(如大口径水管道、系统中有快关或开阀管系)进行水(汽)锤计算是非常必要的。
对管系进行水(汽)锤计算分两步: 需要计算管系中各管段的动态力数值及这个力和时间的曲线关系;!知道动态力的大小和变化关系,计算动态过程对管系应力、支吊架荷载等的影响。
式中,F为最大不平衡力;C为压力波在管道中的传播速度;L为直管段长度;dw/dt为阀门动作过程中流量的变化率。
以上公式基于如下假定: 整个过程温度恒定;!蒸汽视为理想气体;#流动无摩擦;∃流体为单一流动;%由压力波产生的蒸汽特性变化是等熵的;&在阀门、弯头、三通、大小头处的动力损失忽略;∋不考虑波的干扰。
鉴于手算不精确,引入PIPENET软件计算,避免不准确性。可得到较准确的最大动态力数值及动态力和时间的变化关系曲线。
有了准确的输入数据,就可以把数据输入应力分析软件CAESAR,分析各管段上的动态力对管系造成的影响。
2模拟动态力对管系的破坏
实际上管道的水(汽)锤计算不外乎两个目的: 计算管系的应力是否超过规范规定的许用值,通常为防止管道的破坏,在管道的适当位置要加限位装置和液压阻尼器,以提高管系的刚度,防止瞬间破坏;!计算限位装置和阻尼器所要承受的荷载,选择合适的阻尼器和限位装置。
CAESAR软件使用动态运动方程的数值积分法来模拟荷载作用期间的系统响应。
用动态运动方程来求解分析整个管系的响应如下:
Mx((t)+Cx)(t)+Kx(t)=F(t)
1流体瞬态分析
由于水(汽)锤是由于管道内介质流动状态突然改变而形成压力波,并在管道产生轴向冲击力。动态分析外在表现是管道受冲击,而实际上是由于内部流体产生的力。所以对流体进行精确的动态计算也就相当于计算管道受到的冲击力。在没有引进精确瞬态计算软件如PIPENETTRANSIENT前,计算管系的不平衡力多用手算,手算公式如下
[1]
:
F=C∀(L/C)∀dw/dt=L∀dw/dt
式中,M为系统的质量矩阵;x((t)为系统加速
*李琳:工程师。1999年毕业于东南大学热能与动力工程专业。一直从事热能专业勘测设计工作。联系电话:85182720,Emai:llilin@sdleccom。
(0531)
2010,20(6)李琳等精确时程分析在管道水(汽)锤计算中的应用探讨
19
度矩阵;C为系统阻尼矩阵;x)(t)为速度矢量;K为系统刚度矩阵;x(t)为位移矢量;F(t)为作用荷载矢量。
(1)解这个方程要用数值积分的方法求得系统的速度、加速度、位移、反力和应力。通过反复迭代来求解管道动态运动方程。
荷载数据被汇编成动力分析过程中的几个时间层,然后被模拟系统动态响应的计算机程序读取,求解每个时间段的系统动力方程。存入每个节点和振动模态的动力参数(加速度、速度、位移等等),然后积分时间段的这些值,这就是时程分析。时程分析通过检查所有时间段的解来确定有荷载作用情况下,任何时间,管道上任意位置所发生的最大位移、力、应力等等
[2]
3例题
以某工程的高温蒸汽管道(汽轮机再热热段管道)为例,介绍用CAESAR进行精确时程分析的步骤及一些关键参数设置,并与传统响应谱分析方法结果比较。
在完成管道静力计算的基础上,要对管道进行动态分析。假设中压联合气门的关闭时间是015s,管道布置及支吊架设置见图1。
。
如果能借助计算机,时程分析是动态分析最精确的方法。但这种时程分析的方法在管道水(汽)锤计算中未得到广泛的应用,原因是这种
方法需要迭代,如果手算工作量将非常大。
(2)现在沿用的多是响应谱分析方法。响应谱分析法是把一个动态力用响应因子和频率之间的关系来描述(每个力生成一个DLF文件),管系的每个振型都和一个频率上的响应因子对应,多个动态力的响应叠加起来就是整个系统的响应。这个方法优点为不需要迭代,运算量相对小,这是响应谱计算被普遍应用原因之一;缺点是这种方法在原理上不如数值积分法(时程分析)方法精确。
以上两种方法都已经嵌套进CAESAR软件,并与PIPENET软件有接口,可以直接将PIPENET生成的每个力转化成响应谱曲线(DLF文件),这也是现在普遍用响应谱分析方法的原因之一。
随着计算机技术的发展,原来时程分析的弱点如计算速度慢、需要消耗更多资源等已经克服。但由于CAESAR软件没法直接读入PIPENET软件生成的力和时间的文件,且这种曲线关系靠人工输入工作量相当大,也使的这个方法的应用受到了。
鉴于此,笔者编制了一个将PIPENET数据转化成CAESAR时程分析输入数据的软件,解决数据输入的瓶颈,使此类问题采用精确的时程分析成为可能。图1热段管道立体图
热段管道参数:设计压力为522MPa;设计温度为574∗;再热蒸汽流量为1744t/h;主管为914+32;支管为660+22。
图1中圆圈代表承受动态力的阻尼器和限位支架编号,方框为弯头节点号,即各管段不平衡动态力的作用点。
4时程分析步骤
动态计算首先要算出各管段的不平衡动态力和时间的关系,这是由PIPENET软件的瞬态分析模块完成的,软件生成了每个管段不平衡力和时间关系的文件,即*FRC文件。下面以生成的*FRC文件为基础,导入CAESAR软件进行计算。如何生成*FRC文件本文不再赘述。
(1)CAESAR软件做时程分析需要输入每个荷载作用点的力和时间的数据,为了方便读取20CHEMICALENGINEERINGDESIGN化工设计2010,20(6)
PIPENET生成的*FRC文件,笔者编制了将*FRC文件离散化的程序,将*FRC文件离散成一个个单独点的力和时间的文件。转化程序界面见图2。
PIPENET里的力生成的方向是和管道箭头方向(建模方向)相同的,没有坐标系,主要目的是要把这个力转到CAESAR应力计算坐标系中。
实际上,汽锤力是通过弯头体现出来的,对弯头来讲,第一个波汽锤力的作用方向是顺着汽流方向,所以要把PIPENET计算力的方向转化到CAESAR应力计算坐标系中。可以通过Force、direction一栏,设置在CAESAR分析时力的作用方向,CAESAR动态分析输入设置力的方向界面见图5。
图2PIPENET数据转化程序
界面使用简单,离散化的数据文件生成在指定目录下面(本例为E:
散化结果界面见图3。
\\),PIPENET数据离
图5CAESAR动态分析输入设置力的方向
(4)CAESAR动态分析输入设置阻尼器刚度界面见图6。
图3PIPENET数据离散化结果
(2)在CAESAR动态分析界面建立时间和力作用曲线,CAESAR动态分析输入各时程力界面见图4。
图6CAESAR动态分析输入设置阻尼器刚度
(5)CAESAR动态分析输入参数设置界面
见图7。
图4CAESAR动态分析输入各时程力
(3)要把PIPENET里的力转化成应力计算
坐标下的力,涉及到一个方向转化的问题。因为图7CAESAR动态分析输入参数设置2010,20(6)李琳等精确时程分析在管道水(汽)锤计算中的应用探讨
21
完成以上设置参数后,即可进行时程分析计算,时程分析结果界面见图8。
少等于作用荷载的最大持续时间加上第一个求解振型的时间。对热段管道来讲,汽机至锅炉的距离大约100m,压强波从阀门行至锅炉大约只用
01~02s,管系最低阶固有频率0396Hz,则持续荷载时间应为1/0369+02=291s,取3s。54阻尼系数
CAESARII默认临界阻尼比为003。这个值在管道系统中可以在001~003之间修改。阻尼系数取为002意味着计算更保守。当将数值
图8时程分析结果
由003改为002时,如果数据发生了显著的变化,那么有必要对整个系统进行微调,不止是临界阻尼比,还包括系统参数和载荷参数的调整。USNuclearRegulationGuide161中推荐对直径大于12(的管道,阻尼系数取002。该工程阻尼系数取002。55时程输出工况
如果持续荷载时间取3s,要求输出3个工况,将得到在1s、2s、3s时的动态输出结果和一个最危险工况的输出结果。输出工况的个数至少是1个。可根据需要设置输出的工况,输出工况个数对结果没有影响。56质量模型
选连续质量模型比离散质量模型更为精确,且考虑了挠度和转动效应对质量分配的影响,缺点是计算时间比不连续质量模型长,且耗用计算机内存量大。为结果精确起见,计算选用连续质量模型。笔者对比连续质量模型与离散计算模型的计算速度,连续计算模型约是离散计算模型计算时间的3倍,计算耗时很长,约20min。
[4]
计算结果包括纯动态工况、冷态和动态工况组合、热态和动态工况组合。并把每个工况的最
大响应自动计算出来。最大响应的结果包括热位移、应力和动态力等。
从以上设置步骤来看,时程分析的设置要比响应谱分析设置简单得多,工况组合简单,动态只有一个工况,动静组合也只有两个工况,与响应谱多则几十个工况相比,大大节省了管道动态分析的输入时间,且不容易出错。
5时程分析参数设置
由4(5)节中得知,时程分析需要进行参数设置。参数设置的合理与否对时程分析结果影响很大。本节详细介绍时程分析方法的各参数设置
[3]
。
51截止频率
CAESAR技术参考手册上推荐的优化截止DLF为105,基本上包括了共振范围内所有的振型
[1]
。DLF文件可由CAESAR自带响应谱生
成器生成,查看每个力的DLF文件可确定合理的截止频率。从工程每个作用力的DLF曲线来看,时程分析的截至频率拟取至60Hz。52时程时间步长
一般来讲,时间步长越小,计算结果越精确。推荐的时间步长应该这样选择:时间步长X截止频率(Hz)小于01。该工程截止频率为60Hz,最大的时间步长应为01/60=000167s。53荷载持续时间
CAESAR技术参考手册上推荐持续时间至6时程分析和响应谱分析结果比较
时程分析和响应谱分析是两种不同的理念,所以有时结果上存在较大差异是可能的。从计算
原理上,时程分析更符合实际。笔者对热段管道也进行了响应谱分析,和时程分析是基于一样的PIPENET结果文件,为了保证两种方法输入的最大力是相同的,响应谱分析取力的组合系数为10。将两种分析方法的结果比较,供参考。需要说明该比较结果仅限于该管道,可能不具有普遍性。22CHEMICALENGINEERINGDESIGN化工设计2010,20(6)
表1管系应力结果对比
61管系应力
该管道最大应力比较,时程分析最大应力水平为616%,而响应谱分析为662%。
选部分三通点与弯头点作比较,管系应力结果对比见表1。
由表1可见,用响应谱分析与时程分析计算管系的应力基本在一个数量级上,总体来讲响应谱分析方法计算的管道应力水平偏大,结果趋于保守。势必会选择更多的阻尼器和设置更保守的弹簧刚度。62汽锤力
汽锤力计算结果对比见表2。
节点2401203041991341461611691792743492944
1584314741239416394040941292771436044226844324523435395249797
响应S1时程分析S2应力最大发相对误差%
(kPa)(kPa)生的时刻(ms)(S1S2)/S1*100
14671110043619436788355382333975450207421024221849251338521
12133344546113173714354484321901205842686
740255912231165978742705-15140394751118485256
表2汽锤力计算结果对比
支吊架点3070阻尼87阻尼3025阻尼107阻尼139阻尼198阻尼321阻尼281阻尼6阻尼125刚性吊架128X向限位155刚性吊架271Z向限位657X向带间隙限位
作用方向
ZZXXXYXXXYXYZX
响应谱分析动态力
F1(绝对值)
(N)
82492850871020908638019066915706210519881738104458173510529435427030717161458
时程分析动态力F2(矢量)
(N)
5928755204-26921-27755114916-107314-461885271863-175442-43377-304421248500-130971
时间(ms)387392584355307213165242328315354234176533
相对误差%(|F1|-|F2|)/|F1|*100
28133512736367873973316756104063-7892-0588014071910-11311
由表2看出,两种方法计算出的动态荷载差别很大,大部分支架处响应谱方法计算的力要大,但在某些支架处又不够保守。在动静结合工况力的计算公式:
响应谱分析力=动态力绝对值+静态力矢量时程分析力=动态力矢量+静态力矢量
工况可以单向受力。但如果动态力+静态力要改变原静态力的方向,单向限位应该改为双向限位。如125号为刚性吊架,静态力是-44803,动态合力是-175442,合力未改变方向。但如果合力改变方向,就应该改成,Y,支吊架形式由刚性吊架改为限位拉杆。这一点是响应谱分析得不到的。63动态位移
从公式上看,时程分析显然更合理。对限位装置来讲,有些为单向限位,在静态
在这里只截取阻尼器的热态位移做比较,热态工况动态位移对比见表3。
2010,20(6)李琳等精确时程分析在管道水(汽)锤计算中的应用探讨
表3热态工况动态位移对比
23
支吊架点30708730251071391983212816
响应谱分析位移D1(绝对值)(mm)X163245210711131156505520112066951418032681155696587264309
Y14540096346921452145700152146020730418543736309475
Z15672231585675180116517471593915091452117380071290859468226
X
时程分析位移D2(矢量)(N)
Y1149029635335534436171-49112-2142466-736039-554861-314019
Z1563583158098717566141708373-357027-145108-372572288415-448857
-16143521054748-113933550623111958751456503675112693427-267930
从表3的结果可以看出,响应谱分析得到的动态最大位移值是没有方向的,严格来说结果不满足要求。时程分析带有时间参数,可以准确得到最大力作用时刻,并且可动态观测汽锤发生时管系动作情况,比响应谱分析更符合实际情况。
响应谱分析计算的所有数据,而且结果更合理。(5)用CAESAR进行时程分析已不是问题,这时PIPENET计算结果的合理性就显得非常重要。关键是力和时间的关系要准。例题仅将输入数据的时间段取在0~800ms(包括第一个波的作用周期)。
(6)理论上讲,对已知力随时间变化关系的荷载,都可应用时程分析方法。它不适用非预知荷载,如地震等。
参考文献
1陈国宇.蒸汽管道的汽锤暂态响应及其荷载的简化分析原理[J].热机技术,2001,
4(2).
2北京艾思弗计算机软件技术有限责任公司.TechnicalReferenceManual技术参考手册.
3北京艾思弗计算机软件技术有限责任公司.户手册(中文版).4
EstablishingTimeHistoryInput?AnExample.COADEMechanicalEngineeringNews,Vol18(6/94).
(修改回稿2010-06-17)
CaesarII用
7结语
(1)在基于水利计算软件如PIPENET对各管段进行不平衡动态汽锤力计算的基础上,本文
详细介绍了用时程分析方法对管系进行分析,校核管系动态应力、得出动态位移和各刚性架、阻尼器承受动态荷载等。比传统的手算时程分析更精确,考虑因素更全面。
(2)精确时程分析的原始数据处理及输入已不再是瓶颈,如果PIPENET结果调整,也可很方便地进行修改计算,便于调试。
(3)从原理上讲,时程分析是汽锤力分析的标准方法,比响应谱分析方法更精确。从数据输入上来讲,也比响应谱分析更为简便。
(4)从计算结果来看,时程分析可计算出
(上接第34页)
6结语
(1)我国正面临产业结构调整大好时机,精细化工产业的发展前景广阔,其相关设计防火技术标准的研制正当其时,国家标准委已将精细化工列入标准研制的三年计划中。
(2)最近获悉上海市消防协会为适应上海
市精细化工产业发展的需要,正筹备先行编制地方标准−精细化工企业设计防火技术规程.,这是令人高兴的上下合力之举,盼望标准尽早出台,填补空白,促进我国精细化工产业的健康持续发展。
(收稿日期2010-06-05)
2010,20(6)ABSTRACTSOFCHEMICALENGINEERINGDESIGN
Keywords
pipe
1
stressanalysisTimeHistoryAnalysis
ABSTRACTS
SummaryofDesignTechnologyfor800kt/aDomesticmade
LargescaleSulfurBasedSulfuricAcidPlant
YangShaoxing,etal
(YunnanChemicalIndustryDesignInstituteCo.,Ltd.,Kunming650041)
Brieflyintroducetheproductionmethod,
technologicalprocess,
mainequipmenttypeselectionandthelevelofautomaticcontroltechnologyforthedomesticmadelargescalesulfurbasedsulfuricacidplant.Compareandsummarythedesigntechnologywiththeresultsoftheproductionperformancetest.Keywordssummary
ApplicationofEnergySavingTechnologyinPotassium
NitrateProduction
ZhangGang
(YueyangInstituteofPotassiumSaltScience,Yueyang414000)Themajorenergysavingtechnologiesadoptedinpotassiumnitrateproductionincludethepotassiumchloridedissolutionbysteamdirectheating,crystallizationbyaircooling,dissolvingammoniumnitratefollowedbyevaporation,twoeffectevaporation,thesteamfedinseriestotheairheatersofdrier,
etc.
Theenergyconsumptionand
productioncostaresignificantlyreducedbyusingtheenergysavingtechnologiesmentionedabove.Thesteamconsumptionofpertonpotassiumnitratecanbereducedby2622kg.Keywordspotassiumnitrateenergysavingaircoolingtwoeffectevaporation
ProbeonTypeSelectionandDesignEmphasisofPiperack
inChemicalIndustrialPlant
PanFeng
(WuhuanEngineeringCo.Ltd.,Wuhan430223)
Probeonthetypeselectionanditsdesignemphasisofthepiperackinchemicalindustrialplant.Putforwardthesuggestionofreasonableselectionforstructuredesignscheme,optimaldesign.Keywords
piperackforoverallplant
structuresystem
piperackoptimaldesign
ProbeonApplicationofAccurateTimeHistoryAnalysis
inWater(Steam)HammerCalculationforPipe
LiLin,etal
(ShandongElectricPowerEngineeringConsultingInstituteCo.
soastoachievetheaimof
crystallizationby
largescale
sulfuricacid
plantdesign
technical
ResponseSpectrumAnalysis
SafeFactortoBeNoticedinFlareSystemDesign
WangXiaoxia,etal
(PetroChinaNortheastRefining&ChemicalEngineeringCo.Ltd.JiLinDesignInstitute,JiLin132002)
Flaresareusedtotreatthecombustiblegaseswhichisunabletocollectandreprocess,andaretheimportantfacilitytoensurethesafe
productionandreducetheenvironmentpollution.Thisarticledescribesthesafefactorstobenoticedinflaresystemdesign.Keywordssafefactor
flaredischargepipingnet
firerain
safetyvalve
completecombustionbackfire
StudyonConstructionMethodofDoubleFloatingPan
withoutAssemblyJigforLargescaleStorageTank
RenHongl,ietal
(CPECCAbuDhabiCrudeOilPipelineProject,Beijing100120)Theconstructionmethodofdoublefloatingpanwithoutassemblyjigisanalyzedandstudiedthoroughlyinthispaperbasedontheexample
of160000m3floatingrooftankconstructioninAbuDhabicrudeoilpipelineproject.Mainlydiscusstheinstallationandconstructionmethodandthebuoyantforcecalculationforthedoublefloatingpan.Theverificationcalculationresultsagreewellwiththeactualvalue.Byusingthisconstructionmethodforlargescalefloatingrooftank,thefabricationofassemblyjigcanbeomitted,theerectionofoverheadoperationandthedifficultyforconstructiontechniquecanbereduced.Theinstallationqualityandconstructionefficiencyforthedoublefloatingpancanbeimproved.Theconstructioncostcanbesaved.Thesafetyofconstructioncanbeenhanced.
Keywordsfloatingrooftankdoublefloatingpanassemblyjiginstallationmethod
ConsiderationofCharacteristicsandFirePreventionDistance
forFineChemicalIndustry
WuJun,etal
(ShanghaiChemEngineeringIncorporation,Shanghai200235)AnalyzethecharacteristicsforthefinechemicalindustryenterprisesbuiltinShanghaiChemicalIndustryPark.Putforwardthesuggestionthatthedesignstandardofthefireprotectionforthefinechemicalindustryenterpriseshouldbelaiddownassoonaspossible,especiallyforthefirepreventiondistances.Keywords
finechemicalindustry
characteristics
design
standardfirepreventiondistance
CurrentResearchandApplicationProspectofBiodiesel
LvJun,etal(CNOOCZhejiangNingboLNGCo.,Ltd.,Ningbo315800)
Biodieselisakindofhighqualityandenvironmentalfriendlyregeneratedbiomassnewenergy.Theresearchanddevelopmentoftheregeneratedenergyhasbecomeanenthusiasmintheworldwidescale.Althoughthebiodieselhasbeenthepracticalproductcommonly
acceptedbythemarketintheindustrializedcountriesinEuropeandAmerica,theregularregionalmarkethasnotbeenformedinourcountryandtheindustryisrelativelybackward.ThisarticleintroducesandanalyzesthebiodieseltechnologyandthecurrentstatusandtrendLtd.,Jinan250013)
BasedonthestresscomputationalsoftwaresuchastheuniversalCAESARandthewaterconservancycomputationalsoftwaresuchasthePIPENET,AFTandsoonbeingwidelyused,accurateandreasonablemethod.ie.
introduceamore
theTimeHistoryAnalysisto
analyzethewater(steam)hammerforpipe.Introduceindetai,lbyexamples,theprocedureoftheTimeHistoryAnalysis,andcomparewiththeresultsoftheResponseSpectrumAnalysiswhichisgenerallyapplicableatpresent.
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- huatuo0.cn 版权所有 湘ICP备2023017654号-2
违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务