• 7.26 MB
  • 71页

塔式起重机臂架结构安全评定方法研究

  • 71页
  • 当前文档由用户上传发布,收益归属用户
  1. 1、本文档共5页,可阅读全部内容。
  2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,可选择认领,认领后既往收益都归您。
  3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形,可联系本站下载客服投诉处理。
  4. 文档侵权举报电话:19940600175。
'国内图书分类号:TH213.5密级:公开国际图书分类号:621西南交通大学研究生学位论文塔式起重机臂架结构安全评定方法研究年级二〇一五级姓名陈富申请学位级别工稈硕士专业机械工程指导老师吴晓副教授王彬教授二?一八年四月 ClassifiedIndex:TH215U.D.C:621SouthwestJiaotonUniversitgyMasterDereeThesisgRESEARCHONSAFETYASSESSMENTMETHODFORJIBSTRUCTUREOFTOWERCRANEGrade:2015Candidate:ChenFuAcademicDegreeAppliedfor:MasterDegreeSpeciality:MechanicalEngineeringSuervisor:AssociateProf.XiaoWUpProf.BinWanggAril.2018p 西南交通大学学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定,同意学校保留并向国家有关部门或机构送交论文的复印件和电子版。本人授,允许论文被査阅和借阅权西南交通大学可以将本论文的全部或部分内容编入有关数据库进行检索,可以采用影印。、缩印或扫描等复印手段保存和汇编本学位论文本学位论文属于1.保密□,在年解密后适用本授权书;。2.不保密<使用本授权书“”(请在以上方框内打V)学位论文作者签名:指导老师签名:?日期S/f曰期: 西南交通大学硕士学位论文主要工作(贡献)声明本人在学位论文中所做的主要工作或贡献如下:1.对塔式起重机臂架结构的安全评定流程和方法进行研究。从断裂评定和疲劳分析两方面着手,结合塔式起重机臂架结构的结构特点和工作环境等因素,研究用于塔式起重机臂架结构安全评定的流程和方法。2,.通过建立不含缺陷的塔式起重机臂架结构模型,确定含缺陷部位应力分布情况利用数值计算的方法确定各评定参数,进行结构安全评定分析,同时得出结构的失效6评定图,从而完成对塔式起重机臂架结构断裂安全评定的分析,对结果进行分析3.结合有限元分析结果推导出该结构在随机载荷下的临界裂纹尺寸,同时对Paris一aris曲线裂纹扩展曲线进行优化,分别釆用阶段和二阶段的P,计算结构在不同裂纹初始尺寸下的疲劳扩展寿命,并对结果进行对比分析,得出相关参数对于塔式起重机臂架结构安全评定结果的影响。4.以含单边裂纹的有限宽度平板为例,对相互作用积分法和扩展有限元法两种裂纹分析方法的裂纹仿真结果进行了对比分析,并结合有限元分析软件ABAQUS对塔式起重机臂架结构裂纹结构进行了仿真模拟,得到裂纹区域结构的断裂参数以及裂纹单元的扩展方向及状态。本人郑重声明:所呈交的学位论文,是在导师指导下独立进行研究工作所得的成果。除文中己经注明引用的内容外,本论文不包含任何其它个人或集体己经发表或撰写过的研宄成果。对本文的研究做出贡献的个人和集体,均已在文中作了明确说明。本人完全了解违反上述声明所引起的一切法律责任将由本人承担。学位论文作者签名:日期: 西南交通大学硕士研究生学位论文第I页摘要在塔式起重机中,臂架结构为主要的承载部位,长期在交变载荷下工作。结合近几年塔式起重机安全事故,臂架结构为塔式起重机断裂破坏发生的高危构件。塔式起重机的臂架结构主要是焊接桁架结构,存在初始焊接缺陷,并且经常在高强度的随机载荷下工作,很容易发生疲劳破坏,从而使得臂架结构在工作使用中存在安全隐患。因此,针对塔式起重机臂架结构安全性的研究是确保塔式起重机安全运行的保障。本文以某建筑施工工地服役的TC6020塔式起重机为研究对象,从结构完整性评定理论出发、,以断裂力学理论可靠性理论、损伤演化准则以及疲劳损伤理论等为基础,结合有限元分析结果和数值分析等技术手段,展开对塔式起重机臂架结构安全评定方法的研宄。论文的主要内容及结论有:(1)基于结构完整性理论来对塔式起重机臂架结构的安全评定流程和方法进行研宄。从断裂评定和疲劳分析两方面着手,结合塔式起重机臂架结构的结构特点和工作环境等因素,研宄用于塔式起重机臂架结构安全评定的流程和方法。(2)对塔式起重机臂架结构的断裂安全评定方法进行研宄。以某建筑工地使用的TC6020塔式起重机臂架结构为例,利用有限元分析软件ABAQUS建立该塔式起重机1臂架结构的有限元模型;同时,以英国标准BS790三级断裂评定准则为参考,结合臂架结构裂纹缺陷的参数信息以及有限元应力分析结果,对裂纹缺陷进行断裂安全评定,从而得出该塔式起重机臂架结构的断裂安全评定结果。(3)对塔式起重机臂架结构的疲劳安全评定方法进行研究。同样以TC6020塔式起重机臂架结构为研宄对象,以名义应力疲劳寿命分析法和疲劳裂纹扩展评定法为基础,结合塔式起重机臂架结构实际的疲劳载荷工况和裂纹缺陷特点,对比分析两种方法的优缺点,结合结构的有限元分析结果,对疲劳评定参数进行计算,预测臂架结构的疲劳寿命,从而得出该塔式起重机臂架结构的疲劳安全评定结果,并讨论各参数对于疲劳评定结果的影响。(4)塔式起重机臂架结构的裂纹扩展有限元仿真。首先,对相互作用积分法和扩展有限元法进行对比分析一,以有限大含缺陷平板为例,分析两种方法在裂纹仿真中的优劣势;其次,结合塔式起重机臂架结构的有限元分析结果,提取裂纹附近的边界条件,建立裂纹所在位置的局部有限元模型,采用扩展有限元法对含裂纹缺陷臂架结构进行裂纹扩展仿真分析,得到含裂纹缺陷臂架结构的断裂参数以及裂纹单元的扩展方向及状态。关键词:塔式起重机;臂架结构;结构安全性;断裂评定;疲劳评定;裂纹扩展 西南交通大学硕士研究生学位论文第II页AbstractInthetowercranetheibstructureisthemainartofthebearinanditwithstands,jpg,-Ittlontermalternatinloads.ncombinationwihhesafetyaccidentsofthetowercraneinggrecen-tearstheibstructureisthehihriskareawheretowercranesbreakdown.Theiby,jgjstructureofthetowercraneismainlyaweldedtrussstructuresothereareinitialweldin,gdefec-ts.Inadditiontheworkisoftenerformedinahihintensitcclicrandomload,pgyy,whichisronetofatiuefailuresothereisasafethazardintheworkoftheibstructure.pg,yjThereforetheresearchonthestnxcturesafetoftheibstructureisarticularlimortant.,yyjppTakingtheTC6020towercraneservingonabuildingconstructionsiteastheresearchobectaccordintothetheorofstructuralinteritassessmentandbasedonfracturej,gygy,mechanicstheory,reliabilitytheory,damageevolutioncriterionandfatiguedamagetheory,aresearchonthemethodofevaluatinstructuralInterityofibstructurehasbeendoneggjcombinedwithfiniteelementanalysisresultsandnumericalanalysisandothertechnicalmeans.Themaincontentoftheresearchinthisaerincludesthefollowinasects:ppgp(1)Basedonthetheoryofstnxcturalintegrity,thesafetyassessmentprocessandmethodsfortheibstructurearestudied.Fromtwoasectsoffractureassessmentandjpfatiguedamageasetofrocessesandmethodsthataremoresuitableforstructuralsafet,pyassessmentofthejibstructureareformedcombinedwiththestructuralcharacteristicsandworkinenvironmentoftheibstnxcture.gj(2)Researchonthefracturesafetyassessmentofthejibstructure.TakingthejibstructureofaTC6020towercraneasanexamlethefiniteelementanalsisiscarriedout.p,ytake-AtthesametimetheBritishstandardBS7910threelevelfractureassessmentcriteriaas,areferencecombininthefiniteelementanalsisresultoftheibstructureandtheosition,gyjparameterofthecrackdefectsthefracturesafetassessmentiscarriedout.Thustheresultp,y,ofthesafetassessmentoftheibstructureisobtained.yj3Researchonthefatiuesafetyassessmentoftheibstructure.Takintheib()gjgjstructureofaTC6020towercraneasanexamleanalzintheadvantaesandp,yggdisadvantagesofnominalstressfatiguelifeanalysismethodandfatiguecrackpropagationassessmentmethodcombininwiththeactualfatiueloadinconditionsandcrackdefectgggcharacteristicsoftheibstructure.Thefatiuelifeoftheibstructureisredictedcombinedjgjpwiththeresultsofthefiniteelementanalysisofthestructure.Andtheinfluenceofvariousparametersonfatigueassessmentresultsisalsodiscussed. 西南交通大学硕士研究生学位论文第III页4Researchonthecrackroaationanalsisoftheibstructure.Takinafinitelate()ppgyjgpwithalaredefectasanexamletheinteractioninteralmethodandtheextendedfinitegp,gelementmethodareusedtosolvetheforceintensityfactorrespectively.Andtheaccuracyofthetwomethodsforcracksimulationiscomaredandanalyzed.Combinedwiththefinitepelementanalysisresultsofthetowercraneboomstructure,thelocalfiniteelementmodelofthecrackregionisestablished.Thecrackroaationsimulationanalsisiscarriedoutonppgythedefectibstructureandthefracturearametersofthedefectibstructureandthej,pjpropagationdirectionandstateofthecrackcellareobtained.Keywords:Towercrane;Jibstructure;Structuralsafety;Fractureassessment;FatigueassessmentCrackroaation;ppg 西南交通大学硕士研究生学位论文第IV页目录-1m#懸1U研宄背景1.2国内外研究现状21212..疲劳强度的发展及研宄现状1.2.2塔式起重机结构安全评定现状314.3结构完整性理论1.3.1结构完整性概念41.3.2结构完整性评定原理41.3.3结构完整性评定的主要技术51.4课题研宄的目的和意义61.5论文研宄的主要内容6第二章塔式起重机安全评定理论与方法821T.C6020塔式起重机基本参数82.1.1臂架结构参数82.1.2臂架材料力学性能92.1.3起重特性曲线102.2断裂评定分析理论与方法102.2.1线弹性断裂力学理论与方法102221..弹塑性断裂力学理论与方法32.3疲劳评定分析理论与方法152.3.1名义应力疲劳寿命分析法162.3.2疲劳裂纹扩展评定法162.4本章小结18第三章塔式起重机臂架结构的断裂评定方法193.1断裂评定流程193.2裂纹位置及参数193.3塔式起重机臂架结构有限元分析213.3.1臂架结构三维模型213.3.2边界条件213.3.3载荷分析223.3.4有限元分析结果24 西南交通大学硕士研究生学位论文第V页3.4断裂评定准则26—3.4.1级简化评定263.4.2二级常规评定273.4.3三级延性撕裂评定293.5断裂评定计算实例303.5.1断裂评定参数计算303.5.2臂架结构的断裂评定323.6本章小结33第四章塔式起重机臂架结构的疲劳评定方法344.1疲劳评定流程344.2裂纹尺寸及载荷分析3514.2.臂架结构裂纹尺寸354.2.2疲劳载荷分析364.3疲劳评定准则374.4疲劳评定计算实例3814.4.确定疲劳评定参数384.4.2臂架结构的疲劳评定394.5本章小结40第五章臂架结构疲劳裂纹的扩展有限元仿真425_1裂纹扩展仿真方法425.1.1相互作用积分法42125..扩展有限元法435_2裂纹仿真分析实例445.2.1裂纹模型445.2.2裂纹模型有限元分析455.2.3计算结果对比分析485.3臂架结构裂纹模型495.3.1局部实体模型495.3.2裂纹模型5015.3.2边界条件55.4断裂失效准则525.5裂纹扩展结果535.6本章小结53结论与展望55 西南交通大学硕士研究生学位论文第VI页l¥Ctf5758攻读学位期间发表的学术论文及科研成果61 西南交通大学硕士研究生学位论文第1页第一章绪论1.1研究背景近年来,高层建筑物不断增加,塔式起重机得到了,随着房地产行业的大力发展越来越广泛的应用。基于工作效率高、使用范围广、操作方便等诸多特点,塔式起重机在各种工程项目中,尤其是现代工业与民用建筑中的应用越来越广泛,逐渐成为施1[]工建筑中的主要大型机械设备。与此同时,随着投入使用的塔式起重机的逐渐增多,由塔式起重机结构破坏所引起的案例频频发生,造成了巨大的经济损失和人员伤亡,令人触目惊心。据相关统计数据显示,在近年来相关塔式起重机金属结构断裂事故中,超过三分之二是由于长期2[]承受复杂的交变应力而产生的疲劳破坏的原因造成的。塔式起重机疲劳破坏发生的主要位置为:塔身标准节结构、臂架结构、拉杆结构等,这是由于在施工过程中这些构件主要承受随机载荷和循环交变载荷的作用造成的3[]1-1。其中,臂架结构的疲劳破坏时有发生,图为两起在施工现场发生的塔式起重机。结构断裂的事故,造成了巨大的经济损失和人身伤害,同时社会影响极为恶劣1鼷_图-工现场塔式起重机断裂事故11施一目前,大多投入工程使用的大型机械中,具有定缺陷的结构普遍存在。随着服役期限的不断增加,加上,长时间的承受交变载荷使结构的性能退化和承载能力降低初始缺陷的存在,从而使得塔式起重机在继续使用中存在安全隐患,无法保证结构继续安全可靠地工作。目前,但由,质监部门已有塔式起重机结构安全评定的相关规范于不确定因素太多,计算分析困难,导致评定难度大大増加。臂架结构作为塔式起重机中最重要的受力构件,长时间承受交变载荷,是塔式起重机疲劳断裂破坏的重灾区。 西南交通大学硕士研究生学位论文第2页因此,对于臂架结构疲劳破坏的研宄是确保塔式起重机安全运行的保障。结构的疲劳破坏过程根据其发展过程分为疲劳裂纹的萌生阶段、疲劳裂纹的扩展5[]阶段和失稳断裂阶段。塔式起重机臂架结构疲劳破坏起初只是很小的裂纹,随着服役期的逐渐增加一,在交变载荷的不断作用下,这些微小裂纹不断扩展,逐渐延伸达到定程度从而导致结构的疲劳破坏,才会引发安全事故,如经常检查,及时发现和处理。结构裂纹缺陷,并进行相关的结构安全评定,可大大减少此类事故的发生1.2国内外研究现状1.2.1疲劳强度的发展及研究现状疲劳是指在循环载荷的作用下,材料在某点或某些点逐渐产生局部的永久性损伤,6一一继而在定的循环次数后逐渐萌生出裂纹或进步扩展直至发生断裂的现象[]。疲劳强一已开始引起人们的重视度问题,作为结构安全评定中的个重点问题,在很早便。第一次工业革命后,以蒸汽机为动力的机动运载工具和机械设备逐渐得到广泛的使用,“”人类社会进入蒸汽时代。随着越来越多的金属结构投入大规模的使用,结构的破坏现象也时常发生,这些破坏事故也始终困扰着当时的工程师们。1829年,德国的工程师Albert在对矿山卷扬机焊接链条的疲劳试验中,研宄发现随着外加载荷次数的逐渐增多,链条的结构性能也会随之减弱,这首次解释了裂纹破坏的原因,从此人们开“”一始关注这种特殊的现象。直到十九世纪三十年代末,疲劳这专业术语才由法国工程师J.V.Poncelent首次提出,用来描述材料在外加交变载荷下发生结构变形直至最一终发生断裂破坏的过程,自此这词汇便沿用至今。直到1850年德国工程师Wohler一台疲劳试验机-研究设计出了第沃勒疲劳试验机,首次在金属结构上进行了疲劳试验,得出了金属结构疲劳寿命与载荷应力的关系,进而获得了材料的曲线,并首次提出了疲劳极限的概念,奠定了金属疲劳理论的基础。20世纪前,对于疲劳的研究主要集中在宏观的结构变形上,进入20世纪后,科学家们开始通过金相显微镜来研究疲劳破坏的过程。1910年,A.A.Grlffith以玻璃为对象。1iner来研究脆断的理论和实验,这也成为了断裂力学的基础945年,M.A.M提出了-MPin线性累积损伤理论的理论公式,形成了almagramer线性累积损伤法则。到了1963年,美国的科学家C.Paris将断裂力学的方法应用到疲劳裂纹的计算中,提出了估算裂纹扩展规律的Paris公式和损伤容限设计方法,使断裂力学和疲劳逐渐结合了起来。20世纪60年代,科学家们开始将统计学的方法应用于疲劳寿命计算。此后,可靠性理论和损伤容限法也逐渐在结构的疲劳设计中得到应用。相比于国外疲劳研究工作的快速发展,国内的疲劳强度研宄起步的相对较晚,直到20世纪50年代,国内对于疲劳方面的理论以及实验研究才逐渐得到发展。目前, 西南交通大学硕士研究生学位论文第3页经过多年的研宄发展,国内对于疲劳强度的研宄逐渐取得了较大的成绩,有些研宄成7[]Pari果甚至步入世界领先水平。吴晓等在对桥式起重机金属结构疲劳寿命估算中,将s公式与线性累积理论相结合一,提出了种简便的疲劳研宄模型来估算起重机在变幅和8随机载荷下的疲劳寿命[]。范小宁等人采用人工神经网络的技术结合雨流计数法和线性9一[]累积理论,提出了种估算疲劳剩余寿命的新方法。徐永春等人应用断裂力学的思想,建立含裂纹缺陷的结构模型,再结合有限元分析的方法估算该结构的疲劳裂纹扩展寿1〇一[]命。董赞等人运用了种生死单元网格的技术来模拟真实裂纹的扩展过程,为疲劳裂纹的寿命评估提供了新方法一。这些研宄成果在定程度上填补了国内对于疲劳强度方面研究的空白,对金属结构的稳定安全及事故预防运行提供了保障。1.2.2塔式起重机结构安全评定现状塔式起重机的发展已有六十多年了,随着建筑行业的大力兴起,对于塔式起重机的需求越来越多。作为大型高空作业起重设备,随着高层建筑的不断增加,塔式起重机的安全运行尤为重要。为了保证塔式起重机的安全运行,需要对结构进行专业的结构疲劳安全评定计算,而这需要投入大量的人力物力。为了对塔式起重机结构进行安全评定,同时又能简化工作量,提髙工作效率,相关安全评定规范亟需制定和完善。对于金属结构安全评定规范,在国外已经有了比较成熟的发展,目前在压力容器及石油化工设备等领域已得到广泛的应用一些发达国家先后已制定出了多种相关,在欧美标准,目前使用比较广泛的如美国API579、欧洲SINTAP标准和英国的BS7910标准等。同时,我国在2004年发布了适用于含缺陷压力容器安全评定的标准-2004标准GB/T19624,填补了国内对于结构安全评定标准方面的空白。塔式起重机作为较为常见的大型机械设备一,其结构安全问题直都是国内外结构安全领域和质检部门深度研究的对象。目前,国内外对于塔式起重机结构安全方面的研宄主要围绕疲劳寿命预测和稳定性等方面,具体分为基于基础实验的数据分析和数值计算两个方面。赵威威等人基于现场测得的实验数据得到塔式起重机典型工况下12的载荷-时间历程[],再结合疲劳分析软件估算在役塔式起重机的疲劳寿命。张朋等人利用塔式起重机相关设计规范一,提出了种基于反向推算法来估算塔式起重机剩余寿13[]命的方法。郑夕健等人通过采用神经网络的方法搭建塔式起重机应力识别网络,提一出了种通过人工智能获取金属结构载荷谱的方法,进而求得塔式起重机的疲劳寿命。14[]工作数据王进等人通过检测得到塔式起重机的,利用有限元分析的方法得到结构的15疲劳载荷谱[],进而预测塔式起重机结构的疲劳寿命。张宏伟等人通过有限元分析软件ANSYS来模拟塔式起重机臂架结构节点载荷,结合断裂力学及Paris扩展曲线,提出了一种能够用来快速预测出塔式起重机臂架结构疲劳寿命的数值分析方法。总的来说,随着可靠性理论、断裂力学理论、人工智能和有限元分析方法在疲劳 西南交通大学硕士研究生学位论文第4页安全评定中的不断应用一,目前国内在塔式起重机结构安全评定方面的研宄已取定了定的成果。但是,结合国内具体情况,这与真正应用于建筑施工中塔式起重机的结构安全评定还存在一定的距离,主要是由于随着服役时间的增加,国内很多服役塔式起重机缺少完整的的相关工作记录,同时对于疲劳试验数据也不是很充分,缺乏相关的理论基础。,造成较大的计算误差因此,开展对于适用于服役塔式起重机结构安全评定的方法的研宄具有重要的意义和工程应用价值。1.3结构完整性理论1.3.1结构完整性概念随着现代工业的快速发展,对于工程结构可靠性和安全性的要求越来越重要,作为现在工程结构的两大关键特征,可靠性和安全性与工程结构的安全运行息息相关。一结构的可靠性与安全性相互关联,同时也存在着定的区别,研宄可靠性的目的是为了减少结构故障的发生,它的研宄对象是故障,而安全性的研宄对象是危险,以达到降低事故发生的目的。一种对含缺陷结构能否继续正常工作的定量评价结构完整性是,体现了结构抵抗4一[]二种是破坏的能力。目前,现代工程上存在四类标准:第种是产品的设计标准,第建造标准或“质量控制标准,第三种是运行过程安全监控标准,以及第四种建立在适”“”“”用性或者合于使用基础上的标准。结构完整性评定准则正是以合于使用为原则而建立的一种新的方法体系,目前己在航空发动机、压力容器以及石油管道的结构安全评定中得到了广泛的应用。因此,在现代工程结构中,安全性、可靠性和完整性三者相互关联又彼此约束,往往结构的完整性很大程度上决定了结构的安全性和可5[]靠性。1.3.2结构完整性评定原理结构完整性是以确保结构的安全运行为先决条件,来评定结构能否达到规定功能的一种度量指标一。结构完整性评定过程是个多学科交叉的应用,共涉及到工程结构、无损检测、材料科学等多个学科:,具体的评定步骤如下(1)对结构进行无损检测,确定缺陷的种类、形状和大小;(2)结合结构所承受的外部载荷情况,对结构进行应力分析,进而通过计算分析得到含缺陷部位所受的应力;(3)查阅相关资料,确定含缺陷部位所受残余应力的大小;(4)根据相关实验数据或者查阅文献数据来确定材料的力学性能,得出结构的性能参数,如断裂强度,屈服强度〇、断裂軔性尤、以及材料的疲劳扩展速率等;%;/e通过以上分析步骤,提取结构的应力大小、材料性能以及缺陷情况等评定参数, 西南交通大学硕士研究生学位论文第5页对含缺陷结构进行结构完整性评定分析计算,从而判别该含缺陷结构的安全性和可靠性。1.3.3结构完整性评定的主要技术11.3.3.失效评定图技术失效评定图技术又名FAD技术,是综合考虑了结构的脆性断裂和塑性失稳两种失效形式,是目前对含缺陷结构进行安全评定方面应用最广泛的方法。目前,国际上建“”“”4立在[]适用性或合于使用基础上的标准,大多数采用的都是失效评定图技术。失效评定图技术如今已在多种缺陷评定标准得到应用-,如英国推出的BSIPD649391、-REPORTAPV-瑞典的SA/FOU91/01等,我国也有推出相关的缺陷评定标准,如S95规范等。失效评定图技术的核心就是得出结构的失效评定图(FAD),目前比较通用的结1-2所构失效评定图如图示。在图中,两坐标轴分别表示结构脆性断裂和塑性失稳两种失效形式,横坐标&代表了结构抵抗塑性失稳的能力,纵坐标&代表了结构抵抗脆性断裂的能力。通过对含缺陷结构的分析结果,将评定点绘制于结构的失效评定图内,若评定点处于两坐标轴与评定曲线所包围的曲线内时,表示结构的缺陷是可以接收的,结构可以继续安全稳定的工作,;否则结构为不安全的,会发生破坏。牛非安全区^7''K'l'■失效评定曲线-安全区\、IIIIIIILr图-21通用失效评定图1.3.3.2概率安全评定方法随着含缺陷金属结构安全评定方法和技术的不断发展,对于结构安全评定的精度要求也越来越高。由于实际工程下机械结构工况和外部环境的不确定性,导致相关评定参数的选取往往具有随机性,概率安全评定方法正是充分考虑了实际工程情况中外部工作环境和评定参数的不确定性,将随机理论与安全评定相结合,对安全评定中具有统计不确定性的随机参数建立随机干涉模型,从而求解结构的失效概率。概率安全 西南交通大学硕士研究生学位论文第6页“”评定的首次提出是在R6新补充的附录概率断裂评定意见和附录中,随后,其他含缺陷结构评估标准也相继添加了概率安全评定的内容。对含缺陷构件的概率安全评定过程中:,主要需要考虑以下两方面的内容(1)考虑在结构加工和制造过程中以及使用中所造成的非人为的不确定性;(2)考虑在结构安全评定中由于人的经验判断和相关的限定条件所造成的参数的不确定性;研究以上两方面的不确定性,对安全评定的结果进行校正,以确保评定结果的可靠度。将通过计算所得到的失效概率与所要求的失效概率相比较,判别式为P<P-11〇()f式中,4为实际情况下通过计算得到的结构失效概率;巧为安全评定标准所要求的结构失效概率-1);当计算结果符合判别式(1的要求时,认为结构是安全的,可以继续安全可靠的工作,否则,则不安全。1.4课题研究的目的和意义随着建筑行业的快速发展,高空建筑物逐渐增多,塔式起重机的投入使用越来越广泛。随着塔式起重机服役年限的逐渐增加,由塔式起重机结构断裂引起的安全事故也越来越多,造成巨大的经济损失的同时也造成了严重的人身伤亡。因此,开展对于塔式起重机结构安全评定方法的研宄势在必行。在塔式起重机中,臂架结构为承载的主要构件,长期在交变载荷下工作。结合近几年塔式起重机安全事故,臂架结构为塔式起重机断裂破坏发生的高危构件。塔式起重机的臂架结构主要是焊接桁架结构,存在初始焊接缺陷,并且经常在高强度的循环随机载荷的工作,很容易发生疲劳破坏,从而使臂架结构在工作中存在安全隐患。因此,针对塔式起重机臂架结构安全性的研宄是确保塔式起重机安全运行的保障。以建筑工地使用的服役塔式起重机为研宄对象,在前人对于塔式起重机疲劳安全研宄的基础上,结合断裂力学、疲劳损伤理论和可靠性理论,借助于有限元分析软件ABAQUS,通过数值计算和有限元分析的手段,研宄用于塔式起重机臂架结构安全评定的方法,对塔式起重机臂架结构的疲劳寿命进行预测,使其在服役期限内能够安全的运行,节约资源的同时为社会带来更多的经济价值。1.5论文研究的主要内容本文以某建筑施工工地服役的TC6020塔式起重机为研宄对象,从结构完整性评定理论出发,以断裂力学理论、可靠性理论、损伤演化准则以及疲劳损伤理论等为基础,结合有限元分析结果和数值分析等技术手段,展开对塔式起重机臂架结构安全评定方法的研究。论文的主要研究内容如下: 西南交通大学硕士研究生学位论文第7页(1)基于结构完整性理论来展开塔式起重机臂架结构的安全评定流程和方法方面的研究。从断裂评定和疲劳分析两方面着手,结合塔式起重机臂架结构的结构特点和工作环境等因素,研宄用于塔式起重机臂架结构安全评定的流程和方法。(2)展开塔式起重机臂架结构断裂安全评定的研究。以英国标准BS7910三级断裂评定准则为参考,结合塔式起重机臂架结构的结构特点和裂纹缺陷的位置参数,提出适用于塔式起重机臂架结构断裂评定的流程和方法。以某建筑工地使用的TC6020塔式起重机臂架结构为例,对该臂架结构进行有限元分析计算,提取断裂评定参数,从而得出该塔式起重机臂架结构的安全评定结果。(3)展开塔式起重机臂架结构的疲劳安全评定的研究。结合相关的疲劳理论,提出适用于臂架结构的疲劳分析流程和方法。同样以TC6020塔式起重机臂架结构为研宄对象,结合结构的有限元分析结果,估算臂架结构的疲劳寿命,从而得出该塔式起重机臂架结构的疲劳评定结果,并讨论各参数对于疲劳评定结果的影响。(4)塔式起重机臂架结构的裂纹扩展仿真分析。对相互作用积分法和扩展有限元法进行对比分析一,以有限大含缺陷平板为例,分析两种方法在裂纹仿真中的优劣势;结合塔式起重机臂架结构的有限元分析结果,建立裂纹区域的局部有限元模型,对含缺陷臂架结构进行裂纹扩展仿真分析,得到含裂纹臂架结构的断裂参数以及裂纹单元的扩展方向及状态。 西南交通大学硕士研究生学位论文第8页第二章塔式起重机安全评定理论与方法2.1TC6020塔式起重机基本参数塔式起重机主要由塔身、塔帽、起重臂架、平衡臂架、回转支承架、底架、台车12[]。架等结构组成,对于不同的塔式起重机,由于构造上的差异,部件也有所增减本文所研究的塔式起重机为某建筑工地使用的TC6020塔式起重机-,结构简图如图21所-示,其主要性能参数见表21。.mm-?W^p1_身§戀ifiII-(m-(y图2-TC16020塔式起重机表2-TC16020塔式起重机主要性能参数独立高度变幅速率回转速度/最大起重额定起重小车质量/吊钩重量/臂架重量//m/m/minr/min量/1量/1kgkgkg-53.813/38/5800.710232030067752.1.1臂架结构参数TC6020塔式起重-m。机臂架结构由八节桁架结构组成,总长为60如图22所示每,节主要是由上弦杆、下弦杆和腹杆组成。上弦杆采用圆钢,下弦杆采用方钢,空间斜腹杆与水平斜腹杆均为圆钢-2所示,水平直腹杆为角钢,具体参数如表2,结构简图 西南交通大学硕士研究生学位论文第9页-3所示如图2。第一节第二节第三节第西节第五节第六节第七节第八节、八八八/\八八八八八八八八八八7\八八八7^7^八八八八八/v\A八八八八八八八八八八八八0001QQ0010000500500QCQQ1Q0QQ50005001___t^l60000^^2-2图TC6020塔式起重机臂架结构表2-2TC6020塔式起重机臂架结构参数名称型号材料---12节36节78节上弦杆070¥70055Q345B下弦杆U108X10d98x8d68x6Q345B空间斜腹杆057x3.5057x3.5032^3Q235B051x3¥51x3051x3水平斜腹杆042x3?&42x3042x3Q235B水平直腹杆槽钢63x6槽钢50x5槽钢5〇x5Q235B长拉杆0)60Q345B短拉杆0>55Q345B肢上弦靡-4\\\\\\\\/\/\/\/\/\/\/\/\/臂架上弦图2-3臂架结构简图2.1.2臂架材料力学性能表2-3力学性能参数材料参数弦杆许用应力腹杆许用应力材料密度P弹性模量泊松比■MP3〇/aa/MPa/(k/mm)EfMPaV[][]g-65参数值2571757.85xl〇2.1xl〇0.3参照GB/T-113752992《塔式起重机设计规范》结合塔式起重机的结构构造和工,作环境,本文选用载荷组合B对该TC6020塔式起重机臂架结构进行分析计算,其中 西南交通大学硕士研究生学位论文第10页安全系数为1.34。已知所选用的塔式起重机臂架结构所有弦杆材料为Q345B,所有腹235B2-杆材料为Q,材料的各力学特性参数如表3所示。2.1_3起重特性曲线起重特性曲线作为塔式起重机作业性能的一个重要指标,表示该塔式起重机最大起重量随工作幅度变化的关系。该TC6020塔式起重机采用的是机械式力矩限制器,可通过计算起升载荷对臂架结构上铰点的力矩来求得起重特性曲线。通过计算分析,该塔式起重机起重特性曲线可表示为"10000R<R,{=2-1Q\1551564()-620R>R,'1-Li?0.78。式中,A为额定起重量所能达到的最大幅度将该表达式在函数图象中表示出来,便得到了该塔式起重机起重特性曲线。求得的TC6020塔式起重机起重特性曲线如图2-4所示。1200010000???#??8000J\?6000J.'刪j\4000j2000J,0- ̄!111\1ir010203040SO6070工作幅度兑'm图2-4TC6020塔式起重机起重特性曲线2.2断裂评定分析理论与方法断裂评定分析的基本原理是断裂力学理论一。断裂力学是门用应用力学的理论及一18[]方法来研宄含缺陷材料和结构的破坏问题的门学科。随着理论技术和实验条件的不断发展,应用力学在工程上的应用越来越广泛,可用来计算各种条件下疲劳裂纹的扩展问题等。断裂力学是从线弹性力学和弹塑性力学两方面来分析,把裂纹作为断裂破坏的边界条件,计算裂纹尖端的应力应变场和位移场,通过建立这些场与断裂控制参数之间的物理模型来模拟裂纹的扩展过程。2.2.1线弹性断裂力学理论与方法一线弹性断裂力学理论是断裂力学里面发展比较早的个分支,主要是在线弹性理 西南交通大学硕士研究生学位论文第11页论的基础上对含裂纹结构进行力学分析。在工程实际中,结构所受到的外部载荷不同,导致裂纹的扩展也会发生变化,根据裂纹体的受力情况和裂纹的扩展方式,通常把裂19[]I:型(张开型)裂纹、II型(滑开型)裂纹以及III型(撕开型)裂纹纹分为三类,一如图2-5所示。线弹性断裂力学的研究方法主要有两种,种是从能量的角度分析,一基于能量的G判据准则;另种是从应力应变场的角度来描述裂纹的扩展过程,基于应力应变的K判据准则。TTf(i(滑开型型(张开型)裂纹II型)裂纹m型(撕开型)裂纹图2-5三种常见裂纹类型2211...G判据准则裂纹的形成可以看做一个能量转化的过程,在裂纹的形成过程中,将系统释放的能量转化成新裂纹界面的表面能一,这过程符合能量守恒定律。可以通过能量转化的过程来描述这一过程,表达式为G=G2-2jIC()式中,为材料的应变能释放率;G/e为材料常数,描述材料抵抗裂纹扩展的能力,可通过实验获得。该准则用能量释放率来描述裂纹的扩展过程,即为材料脆性断裂的G判据准则。2.2丄2:J判据准则应力强度因子尤作为断裂力学中一个重要的控制参数,通常用来描述材料抵抗脆2()性断裂的能力[]。判别裂纹发生失稳断裂的临界条件为K=K2-3:IC()式中,尤,是材料的应力强度因子;^/e为裂纹临界条件下的应力强度因子,即材料的断裂初度,描述材料抵抗脆性断裂的能力,可通过材料的拉伸实验获得。这就是材料的脆性断裂足判据准则,用应力强度因子来描述材料的断裂性能以及裂纹扩展情况。线弹性断裂力学的两种方法是可以相互转化的,G和^之间符合下式的关系, 西南交通大学硕士研究生学位论文第12页平面应变状态下EG=2-37()S2,平面应力状态下IE式中,E为材料的弹性模型;v为泊松比。同时,在弹性条件下,Gy与尺之间也符合下式的关系式G=^2-5ICf()由此可见,在对结构进行线弹性断裂力学分析时,采用两种方法完全是等效的,可通过具体可获得的参数来选取合适的方法进行线弹性断裂分析。2.2.1.3应力强度因子在实际生产中,在材料的内部难免的会有缺陷的存在,如夹杂、气孔等,由于这些缺陷的存在,使得材料内部存在不连续性,这种不连续性可等效为材料裂纹在裂纹尖端产生的应力集中,形成裂纹尖端的应力应变场,应力强度因子正是断裂力学中用来表示这个应力应变场的参数。对于三维裂纹,参照弹性力学理论,裂纹尖端的应力应变场的表达式为9.9.3^1■〇SW[^co---—---_slsmsinsin2+coscos^()()〇222222f^19.e.W.99W0n,co--—--—^sl+smsmsmcoscos^()222222〇21uIU〇----<><2>><>=<2/sin26f//cos+〒/+()f^xl7ir22丌ry4qc〇se.ewen.d.3tr^--—co---—Lcossmcoss1sinsinV()222222e一〇〇也了—00-JV.>■A:A!A:式中,a〇〇为裂纹尖端的应力分量分别为I型裂纹、II型;/,//,///;,,;,,裂纹、III型裂纹对应的应力强度因子;为裂纹尖端极坐标分量;p为材料的泊松比。-r趋向于零时由式25可以看出,对于实际裂纹,裂纹尖端的应力存在奇异性,当,理论应力将趋向于无穷大,这显然与实际情况不相符,因此用应力强度因子的值来表示裂纹尖端的应力场的强弱。在线弹性力学理论中,对于应力强度因子的计算主要有三种方法:解析法、数值模拟法和实验法。(1)解析法对于形状规则并且边界条件简单的常见裂纹,往往可通过解析法求得裂纹的应力强度因子。解析法求解应力强度因子简单易操作,为日常中最常见的求解应力强度因子的方法。针对于不同的裂纹类型,有不同的应力强度因子求解公式。查阅应力强度 西南交通大学硕士研究生学位论文第13页21一[]般情况下因子手册,,对于常见的I型张开型裂纹,应力强度因子的通用解析式为=ra-K7cr^/^27;()式中,7为修正系数,与裂纹的几何形状、位置和载荷情况有关;C7为结构不含裂纹情况下裂纹尖端所受的名义应力;a为裂纹尺寸。(2)数值模拟法随着数值分析技术和有限元技术的不断发展,这些学科在断裂力学上面的应用也越来越广泛。对于不规则的结构和复杂的加载条件下,用解析法往往无法准确的得到裂纹结构的应力强度因子,数值分析的方法应用而生,可以用来模拟裂纹尖端的应力应变场,从而得到裂纹尖端的应力强度因子。目前,应用数值分析的方法有很多,Chan2223一SK[]a[]等人提出了基于唯方法的位移外推法,BrS〇um等人提出了1/4节点位移法,2425[][]311也等人提出了修正位移法,刘明尧等人分别采用1/4节点位移法、虚拟裂纹闭合法以及位移外推法对含穿透裂纹的平板有限元模型进行断裂分析,求解应力强度因子,并在此基础上研宄了载荷情况、裂纹尺寸和构件尺寸对于应力强度因子的影响。对于上述求解应力强度因子的数值方法可根据求解对象的不同分为直接求解法和间接求解法。直接法计算过程简单易操作,直接利用节点的位移和应力公式,从而推算应力强度因子;而间接法则通过求解单元的应变能,再经过转化间接求解得到应力强度因子。(3)实验法实验法是通过实验的方法得到能够提供应力强度因子的准确值,精度最高,但受限于实验条件,往往只是用来验证其他方法的准确性。近些年,随着解析法和数值法不断地优化发展,但实验法的地位还是无法代替。目前,计算应力强度因子最可行有效的方法是将解析法、数值模拟法和实验法相互结合,通过数值模拟法来模拟裂纹的扩展情况,得到结构的应力强度因子,大大提高工作效率,简化计算流程,再通过实验法提供的实验数据对数值模拟法的结果进行验证,增加结果的可行度,并通过解析法来总结推导出计算应力强度因子的经验公式。2.2.2弹塑性断裂力学理论与方法对于线弹性断裂力学理论,当且仅当主要表现为弹性和脆性断裂的材料来说才完26[]全适用。然而,在实际工程中,完全弹性体是很少存在的,结构在承载过程中,裂纹尖端会出现大范围的塑性变形的情况。这时,采用线弹性断裂力学的分析方法对结构进行断裂分析,会造成结果误差很大,得到的结果不可信。近年来,随着断裂力学的不断发展,为了更准确的得到表征裂纹塑性变形过程的一参量,弹塑性断裂力学作为断裂力学的个重要分支逐渐取得卓越的成效。目前,弹塑性断裂力学理论己经比较成熟,在实际工程上逐渐得到了普遍的应用,其中使用最 西南交通大学硕士研究生学位论文第14页27[]多的两种方法为COD(裂纹张开位移)法和J积分方法。2221D...CO弹塑性断裂力学的COD法用裂纹尖端张开位移即COD来表征裂纹的弹塑性断裂过程,当结构在外部载荷作用下,裂纹尖端的张开位移3逐渐增大到其极限值&时,。结构便发生失稳断裂,判别依据可表示为S<S2-8C()式中,先为裂纹扩展至开裂时裂纹尖端张开量的临界值,可通过对材料的拉伸试验获得,为材料的断裂軔性参数。对于裂纹尖端的张开位移值5的计算方法在很早之前便有人己经提出过。1960年,27[]■Duda-M模型gle等人提出了著名的D,在该模型中,与裂纹面垂直的均布拉力〇作用在无限大的平板上,穿透裂纹的长度为仏,裂纹尖端存在长度为i?的塑性区,塑性2-区作用均布拉应力h,如图6所示。MMa2c「DT图2-6-DM模型示意图该模型认为裂纹尖端为塑性区,外部区域为弹性区,由此推倒的无限大平板上穿透裂纹尖端张开位移COD可表示为^=2-7ln(sec)(9)^^nE2cjs式中,3即为裂纹尖端张开位移。此模型用裂纹尖端张开位移即COD来表征裂纹的弹塑性断裂过程。2222/积分法._./积分方法是与尤判据方法比较相似,也是通过裂纹尖端的应力应变场来描述裂纹的扩展过程。然而,用弹塑性断裂力学的理论很难直接获得裂纹尖端应力应变场,并且在裂纹尖端存在应力奇异性,在断裂分析时需要避开这个区域,这在弹塑性断裂一一力学领域直是个棘手的问题。直到?/积分的提出,利用积分的线路无关性很好地解决了这个问题-d为。?/积分定义如图27所示C为围绕裂纹尖端的积分路线,,线路内部面积。 西南交通大学硕士研究生学位论文第15页X丨■裂纹尖端;,c2-图7?/积分定义示意图?/积分可以从能量的角度来描述,即沿着线路C外部区域对内部区域做功的速率,4一[]不小于储存在X区域中内能的变化率和不可逆的损耗能量率之和。通过进步研宄,得出的J积分的表达式为J=wdy-—-(Tds210)()jc式中,w为应变能密度;m为位移矢量;r为线路c上的应力矢量;办为沿着积分线路c的线路元素。通过该公式所求的的/积分可用来表示裂纹尖端的应力场。当裂纹尖端J积分的值小于或等于结构开裂的临界值;时裂纹结构为安全的,不,一旦当J积分达到这个临界值会发生开裂,,则裂纹会发生开裂,这就是/积分裂纹结构弹塑性分析方法?/积分法可表示为。J<J2-11()C。式中,义为/积分的临界值,可通过实验获得对于大多数金属材料,必须对试件尺寸严格要求下,才能测出准确的;的值。对于脆性材料,可通过线弹性断裂力学理论实现J积分与应力强度因子的转化,转化公式为=^-J212ICf()I■。式中,为型裂纹的/积分临界值对于非脆性材料,此公式不在适用,对于不易通过实验测得J积分临界值的金属材料,可通过阻力曲线来描述?/积分的临界值。2.3疲劳评定分析理论与方法金属结构的疲劳破坏按照缺陷的发展过程通常分为三个阶段:疲劳裂纹萌生阶段、疲劳裂纹扩展阶段和失稳断裂阶段,疲劳寿命主要取决于疲劳裂纹的扩展阶段。目前,对于塔式起重机疲劳裂纹的计算方法主要有基于材料的W曲线的名义应力疲劳寿命分析法和基于断裂力学的疲劳裂纹扩展评定法。 西南交通大学硕士研究生学位论文第16页2.3.1名义应力疲劳寿命分析法名义应力疲劳寿命分析法是将名义应力作为计算评定参数,参照结构材料的W曲线-i,考虑工作环境和加载情况等因素的影响,得出结构构件的SV曲线,结合疲劳累积损伤理论,从而得出结构的疲劳寿命。2.3.1.1材料的曲线一一材料的5WV曲线可由标准小试样在疲劳试验机上实验测得。将组标准试件(一 ̄般812个试件),在定的平均应力下,施加不同的应力幅,从而测出试件断裂时的循环数iV。然后在坐标系中画出这些点,从而绘制出拟合曲线即为材料的基本曲-。线,如图28所示曲线^1 ̄r--- ̄ ̄■- ̄IIrriIII''rri>f?■''"n||疲劳寿命N图2-8材料曲线2.3.1.2疲劳累积损伤理论疲劳累积损伤理论从提出到现在已经经过了多年的研究,相继有数十种损伤假设-M的问世,但在工程上能够得到广泛应用的并不多。目前,应用最广泛的为Palmgeniner法则,可表示为=2-#13()=*1■式中,iV为构件总寿命;%为某应力水平下的循环次数;况为该应力作用下结构发生疲劳破坏时对应的循环次数。当结构累积损伤等于1时,零件发生破坏,可表示为=-&/^.1(214)=il11若总损伤小于,则构件是安全的,则,;若总损伤大于构件不安全结构将会发生疲劳破坏。2.3.2疲劳裂纹扩展评定法一疲劳裂纹扩展评定法是将断裂力学与疲劳两门学科结合起来的种疲劳评定方法, 西南交通大学硕士研究生学位论文第17页该方法假定结构存在初始裂纹或缺陷,基于断裂力学理论,结合裂纹缺陷的尺寸等参数,根据对裂纹扩展速率公式的积分,从而对结构的疲劳寿命进行预测。2.3.2.1疲劳裂纹扩展速率疲劳裂纹扩展速率表示疲劳裂纹随着结构的应力强度的变化而继续扩展的速率,与结构应力强度因子具有一定的函数关系。该函数关系曲线可分为不扩展区、条纹扩-。展区和快速扩展区三个区域,如图29所示ma=dldNC<M)1:I/■■不扩駆I快速扩展区M-g?:条纹扩展区f:/m■/ ̄ ̄- ̄ ̄.u,.,nv^':(MPa应力强度因子幅值M(mm)图2-9裂纹扩展曲线在恒幅载荷下,目前国际上采用最为广泛的为著名的Paris公式,表示为mda/dN=AAK2-15()()式中,而/狀为疲劳裂纹扩展速率;J、m为材料常数,可通过实验测得;为应力强度因子范围,可由下式求得=-YAcr^m216()式中,Act为结构应力幅;7为几何修正参数,可结合塔式起重机臂架结构的具体结构特点求得;a为裂纹尺寸。2.3.2.2疲劳寿命估算对于含缺陷构件,,可通过无损检测技术测得裂纹的初始尺寸你数值计算得到相应的临界裂纹尺寸ae,由载荷信息计算得到相应的应力强度因子值,结合材料的疲劳裂纹扩展速率表达式即可计算出构件的疲劳寿命-1s。将式(26)代入Pari公式,并对裂纹尺寸进行积分便可得到恒幅载荷下构件的疲劳裂纹扩展寿命,即°= ̄da7V—.fmja〇AYA<jra(47)、a(M〇=U2?^AYAcr()2-7=<,(1)”—,A]—m关2士,2-mAYAcrV^()() 西南交通大学硕士研究生学位论文第18页当所求得的裂纹扩展寿命大于所要求的疲劳寿命,同时裂纹还没有达到极限尺寸时,则构件是安全的;反之,则构件不安全,结构将会发生疲劳破坏。2.4本章小结首先,本章对TC6020塔式起重机的结构和基本参数进行了介绍。其次,结合臂架结构的结构特点,从断裂评定和疲劳评定两方面详细阐述了对塔式起重机臂架结构安全评定的具体理论和方法,并对各种分析方法的优劣情况和适用条件进行了分析比较,从而为后续塔式起重机臂架结构的安全评定奠定了理论基础。 西南交通大学硕士研究生学位论文第19页第三章塔式起重机臂架结构的断裂评定方法本章参照英国标准BS7910三级断裂准则,提出适用于塔式起重机臂架结构断裂评定的方法和流程。以TC6020塔式起重机为研宄对象,运用ABAQUS建立该塔式起重机臂架结构有限元模型,同时结合断裂力学理论,对塔机上出现的平面裂纹进行断裂安全评定,计算评定参数,绘制失效评定图(FAD),从而得出该臂架结构的断裂安全评定结果。3.1断裂评定流程对于含缺陷架结构的断裂评定方面的研宄一系列对应的,国内外相关规范已给出评定流程。本文参考BS7910标准中对于金属材料结构安全评定过程,结合塔式起重机。臂架结构的实际情况,提出适用于塔式起重机臂架结构断裂评定的方法和流程对于-1所示的流程进行含裂纹塔式起重机臂架结构的断裂安全评定分析,可按如图3。裂纹位置及参数确定建立臂架结构有限元模型裂纹位置应力分析选择评定准则^^一级简化评定二级常规评定三级延性撕裂评定^确定评定参数r^安全评定计算(Lr,Kr)^得出评定结论图3-1塔式起重机臂架结构断裂评定流程图3.2裂纹位置及参数对塔式起重机臂架结构进行结构安全评定,首先得采用相关技术手段对塔式起重。机臂架结构进行裂纹缺陷检测,目前常见的裂纹检测方法有无损检测和有损检测常见的裂纹无损检测技术有:超声波检测、渗透检测、射线检测、磁粉检测、涡流检测、 西南交通大学硕士研究生学位论文第20页28[]漏磁检测、声发射等。通过对塔式起重机臂架结构缺陷信息统计分析,塔式起重机293()[],-2所臂架结构疲劳损伤的危险位置主要集中于下弦杆与腹杆的焊接处,如图3示。裂纹缺陷在臂架结构上的位置如图3-3所示。图3-2臂架结构下弦杆与腹杆焊接处裂纹'S0、^\微图3-3臂架结构裂纹位置经无损检测得到的实际裂纹缺陷形状复杂多样,不便于进行计算分析,为了便于对裂纹缺陷的规范化,需要对实际测得的裂纹缺陷进行规则化处理。依照BS7910,可将-常见的结构缺陷表征为规则的裂纹状表面裂纹、埋藏裂纹或穿透裂纹,如图34所示。T?,...,(a)穿透裂纹(b)埋藏裂纹(c)表面裂纹图3-4常见的结构缺陷类型本文以TC6020塔式起重机臂架结构为研究对象进行断裂评定,根据臂架结构的特点,先初步确定缺陷可能出现的位置,再结合具体位置的缺陷类型,选用合适有效的无损检测法来进行探测,针对焊接部位可采用超声波检测或射线检测,针对结构表面 西南交通大学硕士研究生学位论文第21页的裂纹缺陷可采用渗透检测的方法。对于检测得到的缺陷信息,需要进行规则化处理,参照裂纹缺陷的规范化标准将裂纹转化为规则的形状,从而得到臂架结构的裂纹参数来进行安全评定计算。在臂架结构下弦杆与腹板焊接处,裂纹缺陷规则化处理后的裂-1所示纹尺寸如表3。表3-1臂架结构裂纹表征尺寸裂纹编号裂纹高度裂纹长度almm2c!mm1362383483.3塔式起重机臂架结构有限元分析3.3.1臂架结构三维模型塔式起重机臂架结构主要是由桁架结构组成,其中,主弦杆受力后会产生拉、压、弯、扭,而腹杆只受轴向拉、压力。针对塔机起重臂的结构特点,基于ABAQUS在非线性分析方面的强大功能,本文使用ABAQUS软件对TC6020塔式起重机起重臂进行有限元仿真分析。定义主弦杆为可承受扭矩的B32梁单元,腹杆和拉杆均采用T3D2三维模型图如图-。杆单元,35所示L图3-5TC6020塔机起重臂模型3.3_2边界条件在塔式起重机的实际使用中,臂架结构的始端通过销轴与塔身相连接,可以绕着水平轴在竖直平面内转动,上端拉杆与塔顶固接。因此,需要对拉杆上端节点进行全约束,对六个方向自由度进行完全约束,具;对臂架结构的始端节点约束五个自由度-体的约束情况如图36所示。 西南交通大学硕士研究生学位论文第22页L图3-6臂架结构边界条件3.3.3载荷分析本文选用载荷组合B对该TC6020塔式起重机臂架结构进行分析计算。在塔式起重机工作状况下,需考虑的的载荷主要有起升载荷、自重载荷、风载荷和水平惯性载34[]荷。3.3.3.1起升载荷起升载荷指塔机起升的最大起重量、起升装置以及起升过程中其他设备的重力之和。起升载荷可按下式进行计算:P=-<pG+G+G312(0)c()eg式中,为各工况的起重量;为小车自重;为吊钩重;为起升动载荷系-112008数,根据起重机设计规范GB38,考虑塔式起重机工况下的载荷,取起升动载系。数妁为1.253.3_3.2白重载荷从地面提升起升质量会对塔式起重机臂架结构的自重产生冲击作用,因此,在计算35[]自重载荷时需要将臂架结构自重乘以起升冲击系数%。结合该塔式起重机工况载荷,取起升冲击系数&为1.1。则自重载荷可按下式进行计算:P=3-2G9fi()式中,G为塔式起重机臂架结构总重量。同时,由于臂架结构上安装有其他设备,但是在有限元建模过程中,为了简化模型,提高计算效率,会对结构模型进行必要的简化和忽略处理,因而使得所建的臂架结构有限元模型自重小于实际实体自重。在此,引入质量修正系数X,通过增大密度来模拟自重载荷,修正后的密度为'p=Xp3-3()式中,Z为实际求得的自重载荷与有限元模型的自重之比。3.3_3_3风载荷塔式起重机一般是在露天环境下使用,因而需要考虑风载荷对臂架结构的影响,并把风载荷看作是作用在臂架结构上沿任意方向的水平力。本文主要研宄塔式起重机 西南交通大学硕士研究生学位论文第23页处于工作状态下的受力情况,故只需要计算在正常工作情况下结构所承受的风载荷,表达式为P=CA3-4wq()式中,C为风力系数;g为工作状态下的计算风压;J为臂架结构的迎风面积,通2TC6020=过计算分析16.08m。查询,得到该塔式起重机臂架结构迎风面积J起重机=17工作状态计算风压设计手册,结合该塔式起重机工作环境,取风力系数C.,=Pa-g250,带入式34中,求得该TC6020塔式起重机臂架结构所受到的风载荷为,=1x250x=^.716.086833^实际上,风载荷以均布载荷的方式作用在臂架结构的迎风面上,在ABAQUS有限元模型中,可等效成作用在迎风面所有节点上的集中载荷。通过对该臂架结构的有限元模型分析,得到该臂架结构迎风面共含有728个节点,则每个节点上的风载荷为=^^^=翌=94#F'w迎风面节点数7283.3.3.4水平惯性载荷在塔式起重机启制动时,会对回转中心产生惯性载荷作用,所产生的回转惯性载荷为=-P(pmsr36H5()=。式中,内为弹性振动载荷系数,1.0<<2.0,本文取1.5m为塔式起重机%外旋转部分的质量,为臂架结构与拉杆的质量总和;r为臂架结构重心到回转中心的距离;s为回转角加速度,可通过下式求得,£=—-37()60,式中,《为臂架结构的回转角速度查询该TC6020塔式起重机工作参数,t为该塔式起重机的制动时间,查询TC6020塔式起重机性能参数表,得到制动时间为10s。通过计算分析,求得水平惯性载荷为^^-=P(pmrH5eot2;rxQ=15x8795x-x25.6.60x10=2475.557V对于臂架结构上的水平惯性载荷,可等效成所有节点上集中载荷。通过对该臂架结构的有限元模型分析,得到该塔式起重机臂架结构共有1136个节点,则分布在每个节点上的水平惯性载荷为4247555,=FH='总结点数1136 西南交通大学硕士研究生学位论文第24页3.3.4有限元分析结果对TC6020塔式起重机臂架结构进行三维建模,将不同工况下的载荷情况加载与模型上进行计算,可得出在不同工况下臂架结构的有限元结果。挑选该塔式起重机四种-典型工况进行有限元分析,四种典型工况数据如表33所示。表3-2TC6020塔式起重机典型工况起重量g/t工作幅度i?/m工况一260工况二2469.75.工况三4.6930工况四1013.9根据上述塔式起重机典型工况的相关数据,在ABAQUS有限元分析软件中对TC6020塔式起重机臂架结构进行静力分析,得到的各工况的载荷情况和应力云图如图3-7至3-10所示。S.HsaMfn.!S|fe5etiKipo5■個論^4l2L2feMG-msimrntl-4387^itZ-lEfiMI-KBLrt?fila6,2-H03?!)?3iaaIHIK)331BsJ?W3?3I]C(l.(Sl;:+:^f4-.vssa:saci^A?-?Xn=IrareBt1:StepmeI,KBWta^'iSrS.ifees你?UDefwmdjnSc*Faaur:?一一(a)工况载荷示意图(b)工况应力云图一图3-7工况有限元分析结果S.lfeesMiiibfcssdKfpaEsJ41.?s?82f■H237H!E-葙組纏量l-47,e*?i.4Kfl5.167■:ta'^A-L?-?=Xlrfl?nerEl:StspTme1.IB3MpyVSrS,蛛swBaarsmiBfic&a:ciRfrDfifVl)DflkSFa+l,376HHa工况二载荷示意图b工况二应力云图()()图3-8工况二有限元分析结果 西南交通大学硕士研究生学位论文第25页S.H5ESMibpfesaflaipMB挪75S)■■nsfeseH-4-!25felSe-+yi,iStse1■■招3誠;Its■-233HB+1Ztemfcs^6m—涵改必繼MSI:K_js/?(.t*■*LSW丨*fina^WrS.lfes'rasimieniKcrcor-Dfc^vJDfatS*Fttl,771MBa工况三载荷示意图b工况三应力云图()()图3-9工况三有限元分析结果S.H5BM|feHtiK〇pwE-n^mmsamnssam.?S-].s7wa4^4124fe?-BtniOtTfl:f97UM■+.!.招32%姐■).W?+5336■I--6393HI'.'*'?、-、、,|t2.77MJi;S-l-,vaairstci^Atn=HiIrari*rel:StEpTmeI,WORTayifrS.UseDefemeiVs:UCStereKnSdeFalx+154&tKa工况四载荷示意图b工况四应力云图()()图3-工况四有限元分析结果10对于以上四种典型工况的应力云图进行分析,可知该塔式起重机臂架结构的最大应力值均小于材料的许用应力,满足结构的设计和使用要求。在第三种工况下,即在两拉杆中间处加载时,臂架结构在下弦杆与腹杆焊接处出现应力集中,所受应力最大为167.8MPa。对于研宄臂架结构下弦杆与腹杆的焊接处的裂纹,此工况为最危险工况,故将此工况作为后续评定分析的计算工况。将缺陷部位的主应力作为断裂评定中的计算参数。在进行断裂评定计算时,假定。臂架结构不存在缺陷,提取缺陷对应位置的应力分布根据应力的作用区域和性质,一次应力P和二次应力0可将其划分为,可通过线弹性计算方法计算得到。在对TC6020塔式起重机有限元模型四种典型工况的分析中,发现工况三为其最危==4i?.6免工作幅度30m。险工况,即起重量0,依据有限元分析结果,结构的最大等■=效应力〇167。.8A/Pfl_一=o=对于塔式起重机臂架结构,弯曲应力可忽略不计,即p,〇则a;次应力6==16二〇7。主要为裂纹缺陷位置的应力,为结构的最大应力,即恳.8^^对于次_应力==,主要为焊接应力,依照BS7910标准,0crs235MPa。m 西南交通大学硕士研究生学位论文第26页3.4断裂评定准则“”目前,对于含缺陷结构的安全评定分析,基于合于使用的原则建立的结构安全评定方法以及与之对应的安全评定标准越来越走向成熟,国际上比较通用的有欧洲工业安全评定准则SINTAP标准和英国标准委员会推出的BS7910标准等,国内也有GB/T19624来规范对含缺陷压力容器进行安全评定的流程。由于含缺陷结构的材料性能、载荷以及评定精度要求等的区别,因此在对不同结构进行安全评定时,所选取的。失效评定曲线也会有所不同,进而所选用的评定等级亦不相同在常用安全评定标准一般将平面缺陷的安全评定方法分为三个等级中,。本文结合最新的英国标准BS7910《金属结构中缺陷验收评定方法准则》,将含缺陷结构的安全评定方法分为一二级简化评定、级常规评定以及三级延性撕裂评定,可31[]根据具体的评定要求选择最合适的评定等级。随着评定等级的增髙,所需要的评定参数也会增多,相应的评定精度也会随之增高。评定等级越低,结果越趋向于保守,当用不同评定等级所得到的评定结构发生矛盾时,以运用更高级评定等级所得到的结果为准。3—.4.1级简化评定一对含缺陷结构的级简化评定方法是最简单的评定方法,对结构的性能参数等要求最低,同时结果也更趋于保守。当能获取的对含缺陷结构的信息比较少时,可以选用一级简化评定对结构的安全进行初步判定一。按照采用的方法不同,级简化评定可分为1A评定和1B评定两种。3.4.1.11A级失效评定1A失效评定方法采用的是简化失效评定图来对结构进行安全评定,其简化评定图如图3-11所示。?*_评样I非安全区0,07A-:。_5米安全区00--:|110-80.00,51.0图3-111A级失效评定图评定曲线由两条分别平行于坐标轴的直线组成,两条直线的表达式为 西南交通大学硕士研究生学位论文第27页°3'8:()f〇s{r式中,#为材料的断裂比,由材料的性能决定;5;为材料的载荷比;对于结构的载荷比可按下式求得,S=^ ̄-r(39)式中,0为材料的流变应力;C7为材料的参考应力,根据结构的载荷情况求得。>re/结构的断裂比^的表达式为'灰语(31〇)式中,5为裂纹尖端张开位移(CTOD)值;为CTOD表示的材料的断裂軔度。通过对含缺陷结构的评定参数进行分析计算,,将所得的结果表示在失效评定图中当评定点位于评定曲线与坐标轴所围成的矩形区域内时,则认为结构为安全的;否则,不能确定结构的安全性,需增加结构的评定参数,选用更高级的评定方法对结构进行安全评定。3.4.1.21B级失效评定1B失效评定没有采用失效评定图技术,是通过计算含缺陷结构的最大容许等效裂纹尺寸,再将裂纹进行当量处理为等效裂纹,将等效裂纹尺寸与最大容许等小裂纹尺寸进行比较,根据比较结果来对含缺陷结构进行安全评定。其中,最大容许等效裂纹尺寸为E8—^<,当1时3-nHL()^'—- ̄-—>cr>a+a,当时2s(zlZ2)式中,五为弹性模量;%为屈服应力;&为膨胀效应系数;为结构拉伸应力。当所评定缺陷的等效裂纹尺寸小于结构的最大容许等效裂纹尺寸时,则认为结构可继续安全运行,选用更;否则,不能确定结构的安全性,需增加结构的评定参数高级的评定方法对结构进行安全评定。3.4.2二级常规评定二级常规。二评定是现在应用最广泛的安全评定方法,也最为大家所接受的级常规评定也分为2A级失效评定和2B级失效评定两种,两种方法均采用的失效评定图技 西南交通大学硕士研究生学位论文第28页术,主要是评定曲线的区别。3.4.2.12A级失效评定2A级失效评定1A采用失效评定图技术,和级简化评定相比,对评定曲线进行了优化-12,评定结果更加精确。该评定曲线如图3所示。其中,评定曲线由曲线部分和垂直于横坐标的截断线两部分组成,曲线部分的表达式为26--^=-:10.14〇.3+0.7ex0.65Z312r(4)p(r)(){}式中,为断裂比;总为载荷比。竖直截断线的表达式为=-313L()rLr^式中,为载荷比h的极限值,取决于结构的材料特性,对于奥氏体不锈钢==材料,可取1.8;对于无屈服的低碳钢及奥氏体不锈钢材料,可取&1.25;_对于无屈服的低合金钢材料=,可取1.15对于其他不能直接按钢材类别确定;的材料,可按下式进行计算,L3-14()。式中,为材料的抗拉强度12-",平定曲线1.〇:0--8\非安全区K-\r0-.6\-安全区\0-.4\垂直线':厂〇,0 ̄.0IIiiIiI0.00.5,DL1.52,01maxLr图3-122A级失效评定图与1A级失效评定的判别方法类似,结合结构的实际情况,对含缺陷结构的评定参数进行计算,将所得到的评定参数以坐标的形式表示在失效评定图上,当评定点落在、评定曲线截断线以及两坐标轴所围成的区域内时,则认为结构是安全的,不;否则能确定结构的安全性。3A2.22B级失效评定2B级失效评定也是采用的失效评定图技术,且评估精度比2A级失效评定的精度高,所得到的评定结果更加准确。但是,相比于2A级失效评定比较普遍的应用,2B 西南交通大学硕士研究生学位论文第29页级失效评定主要应用于焊接结构的各种母材及其焊缝金属材料,且不适用于热影响区2B级失-应变数据域。效评定需要获得材料在适当温度下的应力,经过计算从而得到结构的失效评定曲线。该评定曲线的表达式为°'5+鸟以j-15(3)〇>时,当44max式中,为材料在单轴拉伸状态下的真实应变。2B级失效评定中的评定参数与2A级失效评定中的相同,但需要提供材料真实状态下的应力应变数据。材料的真实应力应变曲线及其相对应的2B级失效评定图如图3-。13所示U^评定曲线■评定点,^||,,r厂二1麵.非安全区\zf0*,5安全区 ̄^0.0|III0/■「max£c.000.350.701.051.404(a)材料应力应变曲线(b)2B级失效评定图图3-132B级失效评定图3.4.3三级延性撕裂评定在二级常规评定中,当评定点位于失效评定曲线上方且A<&时,不能判断结胃构是否安全,这就需要采用三级延性撕裂评定方法。该评定方法也是采用的失效评定图技术,但只是用于具有延性撕裂特征的材料,比如奥氏体和铁素体不锈钢等。三级延性撕裂评定按评定曲线的区别,分为3A评定、3B评定以及3C评定三种安全评定方法。3A3.13A级失效评定该评定方法是在2A级评定方法基础上的延伸-,同样不需要材料的应力应变数据,-1所采用的评定曲线与2A级评定曲线也很类似,如图34所示。当评定点位于失效评定曲线、垂直截断线和两坐标轴围成的区域内时,认为结构是安全的;否则,则不安全。相比于2A级评定-应变曲线的初,这种评定方法考虑了误差的影响,但对于应力始硬化率较高的材料时,评定结果会过于保守,误差较大。 西南交通大学硕士研究生学位论文第30页15娜卜獅'U丨,"\;—\非安全区-0.6-0.4安全区\断线y截-I-101■?■■01110?00.51-0L「max1_52.0L图3-143A级失效评定图3.4.3.23B级失效评定3B-应变曲线数据级失效评定图于2B级别的评定图相类似,同样需要材料的应力,如图3-15所2B级示。相对于安全评定,该评定方法不再受约束,可适用于所有金属材料的安全评定过程。评定曲线■评定点||费季鄉.\■0,5安全区-〇■〇—I111r0.000.350.701.05/_「max1.40K3-B图153级失效评定图3.4_3.33C级失效评定3C级失效评定方法与前面所述的评定方法均有所不同,该评定曲线需要引入材料的弹塑性断裂参数J积分,J积分的取值可以通过对结构进行弹塑性分析求得,进而绘制相应的失效评定曲线。该失效评定曲线的表达式为°.5-C/e/J,当々以):匪尤{哼n-l6)()叫〇,当4>4_时式中一,人和/的数值对应的是在同载荷下的取值,是载荷比总的函数,;因此该评定曲线还是断裂比与载荷比的函数。3.5断裂评定计算实例3.5.1断裂评定参数计算参照BS7910三级断裂准则,其中得到普遍应用的为采用常规评定中的2A级失效 西南交通大学硕士研究生学位论文第31页评定,该评定方法采用通用失效评定图(FAD)的方法进行,相比于1A级简化评定,该失效评定图中对于失效评定曲线做了一定的优化处理,结果更加精确,同时,2A级失效评定法不需要提供材料的应力-应变曲线,这大大简化了评定流程,降低了对于评定结构材料属性以及作用应力等信息的要求,同时也大大降低了计算的繁琐程度。对于本文所研究的某建筑工地用的TC6020塔式起重机,针对臂架结构上的裂纹缺陷,选一用2A级失效评定方法来进行断裂安全评定,评定曲线由曲线部分和条垂直于横坐标轴的竖直线组成。3.5.1.1载荷比计算载荷比反映了结构的塑性失稳程度,用来描述裂纹尖端的塑性量,可通过结构含缺陷位置的参考应力与材料的屈服强度的比求得。在对载荷比进行计算的过程中,可31一[]以忽略考虑二次应力的影响,只需要得到结构的次应力。臂架结构裂纹缺陷部位的载荷比可通过下式进行计算,讀2,俄2L.二17'(3)—31幻(式中,%为材料的屈服强度,<可通过下式求得,2aC=-々18(3)Bc+B()式中,a为裂纹高度;c为裂纹长度;3为裂纹所在处结构的厚度。结合臂架结构上的裂纹缺陷信息,求得的TC6020塔机起重臂结构的载荷比及相关-参数值见表33。表3-3载荷比及相关参数裂纹编号系数屈服强度载荷比fo^/MPaLr100692350849..20.0862350.864301142350892..3.5.1.2断裂比计算断裂比反映了结构的断裂失效程度,用来描述裂纹尖端扩展至断裂的能力。通过对结构裂纹缺陷部位进行计算分析,可求得裂纹尖端的应力强度因子的值,与材料断裂初度的比值即为结构的断裂比。臂架结构裂纹缺陷的断裂比可按下式进行计算,=^-K+319P(),Ykmat=0式中,为塑性修正因子,经简化计算得,取p;,为材料的断裂韧度,& 西南交通大学硕士研究生学位论文第32页为裂纹尖端的应力强度因子,可通过如下解析式求得,K=FY-c7+<y4^a(320i()p()^)[]一式中r二1,Te咖和T〇次应力与次应力,可按照BS790标准((>分别表示修正后的中附录M进行计算。通过安全评定计算,得到TC6020塔式起重机起重臂结构的断裂比及相关参数见表-34所示。表3-4断裂比及其相关参数裂纹编号修正一次应力修正二次应力应力强度因子断裂軔度断裂比050'5.ToPmm-T咖/MPa^s/MaX//MPaKmaMPamm((t/149141.0156.745939.0616052.60.1552205.728216.21295.3106052.60.2143152.0485160.0351106.3056052.60.1833.5.2臂架结构的断裂评定根据BS7910标准中2A级失效评定准则,结合臂架结构的材料属性,绘制失效评定图。对于该TC6020塔式起重机臂架结构,所有弦杆材料为Q345B,所有腹杆材料235B=。三处裂纹的为Q,均为低碳钢,则取垂直线方程为1.25将以上分析得到载荷比Zr和断裂比心所构成的评定点(0.849,0.155)、(0.864,0.214)、(0.892,0.183)在失效评定图中表述出来,即可得到含缺陷臂架结的2A级失效评定结果图,-如图316所示。L2□1j评定曲线a2…Lz0O3'°-8\非安全区Kr-\-0.6\-安全區\-垂直线0.4\△.〇,: ̄0.0iiiiiir0.00.51.0L腿1.52.0Lr-图316失效评定图(FAD)通过对臂架结构的2A级失效评定结果图进行分析,可以得到该TC6020塔式起重机臂架结构的断裂评定结果。观察发现,三处裂纹缺陷的评定点均处在安全区域内,即三处裂纹缺陷均是安全可以接受的。其次,评定点的位置接近评定曲线,在该TC6020 西南交通大学硕士研究生学位论文第33页塔式起重机的后续工作中,需要对塔式起重机工作状态进行实时的监控与信息的反馈,确保臂架结构能够安全稳定的运行。3.6本章小结本章以TC6020塔式起重机为研究对象,对该塔式起重机的臂架结构进行断裂评定分析,参照英国标准BS7910,采用二级常规断裂评定方法中的2A级安全评定方法,对臂架结构上的裂纹缺陷进行了断裂安全评定研宄,得出的结论如下:(1)对TC6020塔式起重机臂架结构的2A级安全评定结果显示,该结构上的三处裂纹缺陷均处在可安全接受的范围内,该臂架结构仍可继续正常运行,并且在后续工作中,需要对塔式起重机工作状态进行实时的监控与信息的反馈,确保臂架结构能够安全稳定的运行。(2)结合塔式起重机臂架结构的结构特点和裂纹缺陷信息,使用ABAQUS对塔式起重机臂架结构进行建模分析,同时将结构有限元分析方法应用于结构的断裂安全评定过程中,在保证评定结果精确性的同时,大大提高了评定工作的效率。 西南交通大学硕士研究生学位论文第34页第四章塔式起重机臂架结构的疲劳评定方法本章以基于材料的曲线的名义应力疲劳寿命分析法和基于断裂力学的疲劳裂纹扩展评定法两种方法为基础,结合TC6020塔式起重机臂架结构进行分析,得出该臂架结构疲劳寿命的计算结果,并进行疲劳评定分析,从而研究用于塔式起重机臂架结构的疲劳评定方法。4.1疲劳评定流程本文第二章给出了两种常用结构疲劳评定方法计算流程,名义应力疲劳寿命分析法需要先得到材料的曲线,核心内容是累积损伤理论,是假定材料内没有缺陷和裂纹,而疲劳裂纹扩展评定法是假定结构存在初始裂纹。对于塔式起重机臂架结构,由于焊接缺陷的存在,使臂架结构在投入使用时已经存在着缺陷和裂纹,从而导致利用名义力疲劳寿命分析法计算得到的结果存在误差;但累积损伤理论的分析方法更接近塔式起重机臂架结构实际工作中的破坏情况。本文参考两种方法的优缺点,结合塔式起重机臂架结构实际工况,建立针对塔式起重机臂架结构疲劳评定的流程图,如图-41所示。裂纹位置确定计算分士裂纹尺、ajc疲劳载荷分析确定危截面应^选择评定准则..1一Par阶is曲线二阶Paris曲线^确定评定参数A,MAO(,AK〇,AKUJ计算疲劳寿命#得出评定结果图4-重机臂架结构疲劳评定流程图1塔式起 西南交通大学硕士研究生学位论文第35页4.2裂纹尺寸及载荷分析4.2.1臂架结构裂纹尺寸4.2.1.1初始裂纹尺寸%对于塔式起重机臂架结构,由于焊接缺陷及自身材料缺陷的存在,在工程使用中不可避免的会存在初始裂纹。对于塔式起重机臂架结构初始裂纹尺寸fl,因结构尺寸、。工艺等因素的不同一材料以及焊接,目前还未形成统的标准。根据相关规范和经验,一35[]?2mm0。采用II级焊缝质量标准的初始裂纹般在.5本文为了研宄初始裂纹尺寸对疲劳寿命结果的影响,分别取初始裂纹尺寸为0.5mm,1mm,1.5mm,2mm。4.2.2.1临界裂纹尺寸在第三章中已经知道,对于在塔式起重机臂架结构中产生的裂纹,其应力强度因子可通过解析式求得,为-K-=Ya^m41()式中,7为形状修正因子;a为结构所受应力;a为裂纹长度。基于线弹性断裂理论,当含裂纹构件应力强度因子达到材料的断裂軔度^,时的裂纹尺寸即为构件的临-界裂纹尺寸。代入式41可得结构的临界裂纹尺寸+为a= ̄4-2c()式中,为臂架结构材料的断裂軔度;F为结构的形状修正因子,结合臂架结构=112的结构特点和工作环境.。,取形状修正参数7由此,可得到结构临界裂纹-2所示尺寸+和与结构最大等效应力〇的关系图,如图4。^40-35-\E30-,\基25-m\眛2〇-—15140160180200220240260最大等效应力cw、(MPa)图4-2临界裂纹尺寸与最大等效应力的关系图 西南交通大学硕士研究生学位论文第36页■BS有限元计算结果〇=16-结合第三章中A7,42中查得QU得到.8」_代入图,胃=该TC6020塔式起重机臂架结构的临界裂纹尺寸&31.6_。4.2.2疲劳载荷分析一疲劳载荷分析是疲劳评定计算中个关键的部分,结构的载荷谱的确定对疲劳寿?1。命影响很大,如果载荷谱误差为0%,将会导致寿命误差达到23倍结构的载荷谱包含载荷的大小和出现的频数等信息,是作用在构件上外载荷的统计表示。对于在役的塔式起重机,其变幅小车的运动大多为随机的,因此结构所承受的真实载荷往往是不规则的。对于承受随机载荷的结构进行疲劳寿命估算时,必须先进行载荷谱的分析。载荷谱的分析步骤为:(1)通过对服役的塔式起重机进行现场观察,对其工作情况进行统计分析,得到该塔式起重机臂架结构的载荷-时间历程;(2)在计算时必须对这些不规则的载荷谱进行处理,通过各类循环计数法将载荷一-时间历程简化为系列的全循环或半循环的过程,从而得到设计计算所需要的载荷谱,如峰值法、蒙特卡洛法或雨流法等。(3)结合该塔式起重机臂架结构的有限元模型,进行加载分析,得到各载荷情况一下的应力情况,绘制应力谱,通过系列简化处理该臂架结构危险部位的应力幅及其对应的频率。通过以上方法对本文研究的TC6020塔式起重机进行分析,发现其危险截面的应力幅值服从正态分布,采用八级阶梯分布应力幅,来模拟该塔式起重机臂架结构的工作-应力谱,如图43所示。对于变幅载荷和随机载荷下工作的结构,结构的载荷分析过程复杂繁琐,不便于直接用于对结构的疲劳评定计算。而恒幅载荷是预测结构疲劳寿命最简单实用的载荷类型,可采用Miner法则对塔式起重机臂架结构的随机载荷进行等效处理,可得到TC6020塔式起重机臂架结构应力幅为1=ik3^=-Acr43()n-式中.3,Vo为各级对应的应力幅; ̄为各级应力幅对应的频数;代入图4各项数;20MPa据可求得结构应力幅为60.。 西南交通大学硕士研究生学位论文第37页丰1000—mS00—I.I^g600_蠢.彡§<y494-杳多観yyw.400_1MK-3多多?^p^i^200-44^i4,M隱II?-H9:iaama^^010.6228_7643.8760.3270_69S3J290_29.1012应力巾iA<?MPa4-3TC图6020塔式起重机臂架结构的应力幅谱4.3疲劳评定准则目前,对于疲劳评定的相关准则已有多年的研宄,目前最为大家所接受的是基于Paris曲线的疲劳评定准则。Paris公式是由美国的C.Paris在断裂力学方法的基础上提出的估算裂纹扩展规律的著名关系式一。自提出以来,由于其形式简单的特点,直得到广泛的应用,它能够较好的描述裂纹扩展的第二阶段,通常作为预测结构疲劳裂纹扩展寿命方法的依据。对于塔式起重机臂架结构,其主要为焊接结构,由于在制造工艺和使用环境等因素的影响,结构中可能已经存在初始裂纹和缺陷。对于含缺陷的结构,Paris公式能够很好地描述结构中裂纹的扩展过程,这对估算结构的疲劳寿命具有重大意义。tK=^m〇(CI3人^od^dN=4^r/(=!hK^KtKLo〇ga—ari()阶段Ps曲线 西南交通大学硕士研究生学位论文第38页iT1(:过渡点l/?^ml=2/dadN.4A^為J(Ax.?bK=LKLogKQ(b二阶段Paris曲线)图4-4疲劳裂纹扩展规律在用对数表示的疲劳裂纹扩展速率与应力强度因子范围的函数中,Paris公式通常一一-可表示为两条垂直线和条斜线的阶段曲线来表示,如图44(a)所示。该曲线将裂纹的扩展过程用规则的直线表示出来,大大简化了计算过程,这也降低了计算的精度。结合结构载荷历程、平均应力、温度和腐蚀介质等因素的影响,对Paris曲线进行一-了适当的修正,种双斜线段的二阶Paris曲线也逐渐得到应用,如图44(b)所示。本文结合塔式起重机的工作环境和使用情况,分别使用两种曲线来进行裂纹寿命计算,并对结果进行分析和比较。4.4疲劳评定计算实例4.4.1确定疲劳评定参数Paris公式中的常数X和w为材料属性,需要通过实验确定,对于TC6020塔式起235B一重机臂架结构,其主要材料为Q345B和Q,查阅文献,得出阶段与二阶段Paris-曲线相关评定参数如表41所示。表4-1Paris曲线相关参数门槛值过渡点过渡点对应裂常数常数裂纹扩展准则AK〇厶幻纹尺寸Am3232(N/mm)(N/mm)cij/mm_13一阶段x2.61l〇3.0063//17I2.1〇xl〇5.10两阶段631441.45u_12II1.29xl〇2.88 西南交通大学硕士研究生学位论文第39页4_4.2臂架结构的疲劳评定通过上述疲劳评定计算,结合TC6020塔式起重机臂架结构实际工作情况,将所得到的临界裂纹尺寸、臂架结构载荷应力幅以及Paris曲线相关参数带入疲劳裂纹扩展寿一命估算公式2-17中,依次代入不同的初始裂纹尺寸,分别采用阶段和两阶段的裂纹扩展规律计算,得到不同初始裂纹尺寸和采用不同裂纹扩展规律条件下的疲劳裂纹扩-。展预测寿命,如表42所示表4-2臂架结构疲劳寿命预测结果初始裂纹尺寸疲劳裂纹扩展预测寿命AT(cycle)一a〇(mm)阶段Paris曲线二阶段Paris曲线660.52.74x102.97xl〇66xx11.81l02.0〇l0661.51.39xl0l.l〇xl〇6621.15xl00.92x10结合上述疲劳评定结果,可较准确的估算出TC6020塔式起重机臂架结构的疲劳裂纹寿命。依照疲劳评定准则,当结构的疲劳裂纹扩展预测寿命大于设计所要求疲劳寿命,同时裂纹尺寸又没有达到临界尺寸时,则裂纹缺陷是安全可以接受的,该臂架结构仍可正常运行,则裂纹是不,需要;反之可接受的对臂架结构裂纹部位进行维护和优化后才能继续投入使用。为了更准确的研究不同初始裂纹尺寸对疲劳裂纹扩展寿命的影响一,可采用阶段一裂纹扩展规律,分别得到不同初始裂纹尺寸阶段Paris曲线下的疲劳裂纹的扩展曲线,如图4-5所示。UJJ05COOOO10000W)15000002JXJOOOO2500000300WSX)裂纹t广屑寿命办:图4-5不同初始裂纹尺寸下裂纹的扩展曲线同时一,取初始裂纹为0.5_,分别得到采用阶段和两阶段疲劳裂纹扩展准则条-件下裂纹的扩展曲线,如图46所示。 西南交通大学硕士研究生学位论文第40页一阶段裂绞扩展曲线二35阶觸绞扩展曲线15震II顏II10參如m^0050000010000CX3150(XMfi200000025000K)30500003SOOOOO■广展寿命欣图4-6不同扩展准则下裂纹的扩展曲线结合表4-2中评定4-5中疲结果和图劳裂纹扩展曲线可知,初始裂纹尺寸的取值对于疲劳裂纹的扩展寿命估算结果影响很大,且当初始裂纹取值很小时,对于结果的影响尤为明显,即在相同的裂纹扩展量下,初始裂纹尺寸越小其循环次数越多,同时,随着初始裂纹尺寸的增加而疲劳寿命急剧下降-。从图46中可看出,采用二阶疲劳裂一阶疲劳扩展准则的疲劳寿命纹扩展准则所得到的疲劳寿命要小于采用,可结合实际情况选用更合适的方法来估算疲劳裂纹的剩余寿命。同时可从图中看出,随着裂纹尺寸的增加,裂纹的扩展速率逐渐增大,且对于不同的初始裂纹尺寸以及不同的裂纹扩一展规律下,裂纹的扩展曲线大致相同。当裂纹增加到定值时,扩展曲线会出现急剧攀升,裂纹扩展速率急剧增加。因此在对塔式起重机的例行检查中,当探测到臂架结构上的裂纹尺寸接近或者大于20mm时,必须适当减少该塔式起重机工作中的负载,并对塔式起重机工作状态进行严密监控、实时跟踪和及时维护。4.5本章小结本章详细介绍了两种常用的疲劳寿命估算方法,在两种方法的基础上,结合TC6020塔式起重机臂架结构的分析实例,估算出了该臂架结构的疲劳寿命,并进行了疲劳评定分析:,得出的结论如下(1)针对塔式起重机臂架结构的疲劳裂纹,将累积损伤理论与断裂力学理论应用于结构的疲劳评定过程中,对裂纹缺陷进行了疲劳裂纹分析,提高了计算结果的准确性和可靠性,同时该评定流程和方法简单易操作,对于塔式起重机臂架结构的疲劳安全评定具有重要的工程价值。(2)在塔式起重机臂架结构的疲劳评定中,初始裂纹尺寸的选择对评定结果的影响很大,应结合实际情况谨慎选择。同时,随着初始裂纹尺寸的增加而疲劳寿命逐渐20_时下降,在裂纹尺寸接近,结构的疲劳寿命随裂纹尺寸的变化速率急剧增加。因此,在实际工程工作中,当探测到裂纹尺寸接近或者大于20mm时,为了阻止裂纹 西南交通大学硕士研究生学位论文第41页的继续扩展,可对裂纹缺陷位置进行加筋处理,或者进行阻断裂纹工艺避免应力集中,或者在严格的工艺下进行补焊,加强塔式起重机工作状态进行严密监控、实时跟踪和及时维护。(3)结合塔式起重机的工作环境等因素,对传统的Paris曲线进行了适当的修正,ari一运用双斜线段的二阶段Ps曲线进行疲劳寿命的估算,并与传统阶段的Paris曲线所得的结果进行对比分析,采用二阶疲劳裂纹扩展准则所得到的疲劳寿命要小于采用一。阶疲劳扩展准则的疲劳寿命,评定结果更安全 西南交通大学硕士研究生学位论文第42页第五章臂架结构疲劳裂纹的扩展有限元仿真目前,对于进行结构裂纹扩展仿真的方法多种多样,最为大家所接受的主要有相一互作用积分法和扩展有限元法。应力强度因子作为结构安全评定中个重要的物理量,一往往被拿来作为评定种分析方法准确性的参数。本文主要针对两种常见的裂纹仿真一含单边裂纹的有限宽度平板为例方法,以,比较两种分析方法在裂纹扩展仿真中的优劣势,最终选择扩展有限元法对TC6020塔式起重机臂架结构裂纹局部模型进行裂纹扩展分析,得到臂架结构裂纹的扩展路径,以便做出相应的裂纹扩展防御措施,提高塔式起重机臂架结构的安全性。5.1裂纹扩展仿真方法5.1.1相互作用积分法相互作用积分法的基本原理是以?/积分的理论思想为基础,通过在裂纹尖端建立辅助场的方法来代替实际裂纹尖端的真实应力场和位移场,从而通过对辅助场的分析53[]计算来模拟实际裂纹的扩展过程。在有限元计算分析中,将虚拟裂纹尖端的应力场和位移场作为辅助场,实际裂纹尖端的应力场和位移场作为裂纹尖端的真实场。一相互作用积分法的表达式与?/积分的表达式较类似,是个在将裂纹尖端包含在内的积分回路上的能量积分,可表示为2))麵—麵-1(5)14,417JJUJdxx-o式中,为结构真实的应力应变场的参数变量;为附加f,#,#()的辅助应力-应变场的参数变量。当积分回路逐渐靠近裂纹尖端区域时,辅助场与附加场之间的关系可表述为1212))))/=i:-r+(52)r(;J44)^E22))^式中,为真实场下的应力强度因子;为辅助场下的结构应力(<,4)强度因子一;矿为与材料参数相关的常数量,与材料的弹性模量和泊松比有关,关系式可表述为'五,平面应力*E=E5-3\_(k);平面应变2l-vl 西南交通大学硕士研究生学位论文第43页2)5-2可以看出:==通过式,当辅助场的应力强度因子满足条件10时,通过<,<?相互作用积分法可以直接求得实际裂纹的真实场的I型应力强度因子<\同理,当辅足=0尤=1II助场的应力强度因子满足时,即可求得实际裂纹的真实场的型应力强产,^度因子X。i5.1.2扩展有限元法扩展有限元法是以传统的有限单元法为基础,对裂纹面或者裂纹尖端附近的单元54[]节点通过裂纹近场位移解来进行増强,从而模拟裂纹出现的过程。扩展有限元法改进了有限单元法的形状函数,使得在单元内部的形状函数可存在不连续的现象,从而可以更好地分析不连续的问题。应用扩展有限元法分析含缺陷结构断裂问题时,可先不考虑裂纹面的位置,直接,然后再引入单元分解法的思想对结构进行网格划分,将反应局部特性的附加函数带59[]入到离散位移表达式。对于分析裂纹缺陷问题,可在扩展有限元法分析中引入Heav-iside函1数来进行加强,从而反映出裂纹面的位移不连续,如图5所示(空心点表示节点)-2所,用裂纹渐进位移逼近场函数来加强含裂纹尖端的单元,如图5示(实心点表示节点),进而反映裂纹尖端区域的局部特性。?M;丨Itt][il——[十一士°.i+-—二二广14士^ ̄一—I—丨卜」丨OcymX一一--—HiI勹,十,———U ̄ ̄O—Q———_-图51单元内任意位置裂纹节点图图52裂纹尖端单元加强节点图对于分析裂纹不连续的问题时,扩展有限元法采用节点扩展的方法,分别对被裂纹完全穿过的单元节点和裂纹尖端附近的单元节点采用不同的扩充形函数,裂纹节点模型如图5-3所示。一裂纹□波裂纹完全穿过的单元节点—————O波裂纹完全穿-------^[M3E3E]E3EJ43E^^^(3过的单元节点£—?OO图5-3裂纹节点模型图 西南交通大学硕士研究生学位论文第44页其中,对于被裂纹完全穿过的单元,裂纹面附近的位移场会发生跳跃,扩充形函数可按下式计算,=-cx(54)p人xNxH(J)、人、Xy式中,表达式,为常用的节点位移函数;丑为沿裂纹面间断跳跃函数为1X ̄°,=-//(5)jc5(){/〇c为结构的水平集函数,表达式为)'-=—--—xvo^nVxxVsinnxx(56)f()gi())+式中,《为裂纹线上的单位法向量。/(jc)表示节点到裂纹面的最短距离,规定当+一致时取正值节点所在位置与《的方向,否则取负值。对于裂纹尖端周围节点,裂纹结构的扩充形函数可按下式计算<x=xx-PAN<(57))j{)p{)式中,#00为各向同性材料的裂纹尖端的渐进函数,可由下式函数的线性组合求得,-x=s-cs-sincos-(8)<j)()Vrin,Vro,Vrsin^sin,^^5^^^^式中,r和0为节点在裂纹尖端极坐标系中的位置参数。一通过引入扩充形函数的概念,把,可以在对裂纹结构进行计算分析时些已知解的信息通过形函数的形式添加到有限元中,从而可以大大减少计算时间,同时提高计算精度。使用两种扩充形函数,可表示出含裂纹结构的位移场函数,表达式如下,h=NNHtN-uxU+xxa+x<xbt(59){)^jJ()(f())J()K()^)K()^^^j===ilil/1式中,为常用的节点位移函数;M为常规形状函数的节点自由度;,&和&为节点的附加自由度。5.2裂纹仿真分析实例5.2.1裂纹模型ABAS裂纹分析的方法一本文为了研究QU,以含单边裂纹的有限宽度平板为例,借助ABAQUS软件,分别采用相互作用积分法和扩展有限元法法对含裂纹平板进行裂纹分析,计算不同裂纹长度下的裂纹尖端的应力强度因子值,并与解析法求得的理论5-4-1值进行比较分析。该含裂纹平板几何模型如图所示,结构的各项参数如表5所示。Z===100mm#20mm52mw。,宽度,厚度 西南交通大学硕士研究生学位论文第45页wLl““““t图5-4裂纹板几何模型表5-1含裂纹平板各项参数物理量参数值L=mm=B=100W20mm2mm尺寸参数,,5材料性能参数五=x二2.1l〇M%v0.3,=载荷参数cr5QMPa=4mm1裂纹长度a2mm6mm8mm0mm,,,,在对裂纹模型进行计算分析的过程中一,应力强度因子是至关重要的个参数量,往往对于应力强度求解的准确性决定了该种裂纹分析方法结果的可信度。对于求解应力强度的方法,在第二章中已给出了最常见的三种方法,其中解析法能够准确得到应力强度的理论精确值,往往将求得的应力强度因子的计算结果与解析法算得的理论值进行对比,来验证该计算方法的准确性。对于应力强度因子的理论值,可按照解析式-)(27来计算。其中,对于修正系数y的值,通过查阅英国标准BS7910可知,对于有限大含裂纹平板的修正系数可按下式进行计算,234=-------71.12.210217(503+.6,+30.410)()()()()BBBB式中,a为裂纹长度;万为平板的厚度。5.2.2裂纹模型有限元分析5.2.2.1相互作用积分法计算借助ABAQUS软件建模模块可以快速建立不含裂纹平板的几何模型,再对裂纹部位进行分割,模拟含裂纹结构。通过相互作用积分法对裂纹模型划分网格时,裂纹必须在网格边界上,否则会造成计算结果不收敛。为了使网格划分规则,需要对裂纹尖端区域进行特殊处理一,对模型整体采用六面体单元C3D8R,对裂纹尖端区域(即第3D6-层围线区域)采用楔形单元C,得到的有限元模型如图55所示。 西南交通大学硕士研究生学位论文第46页■ilil圏_士1|||翁3_麵__lifffll(a)裂纹板有限元模型(b)裂纹尖端局部放大图图5-5相互作用积分法有限元模型=以裂纹长度a2mm为例对模型施加载荷和裂纹,通过相互作用积分法求解计算,得到的应力云图和裂纹区域局部放大图如图5-6所示。通过观察发现,相互作用积分法所得的形变模型发生裂纹断裂区域未发生单元的破坏现象,这些区域的网格只是发生了挤压变形。sesS,MiA75%(vg:)+2,678e+02+245702H.e++2.235+026-+2P^,014e+02—-+179302.e++1.571+02e■-e+1.350+02-十1■■.1296+02+0-+9.0741ee+6.860+01+4.647e+01+2.434+01e-eH2.206+00,--ODB:Jobxh2yf.odbAbaqusStandard&.121Fri:LGMT+08:002018—j/Ist-vep:'=It1StTi1,000ncremen;epmePiVSMirmaryar:,sesDeformedVar:eformaoncaleacor:6UDtiSFt+7.915+02(a)裂纹板应力云图邏(b)裂纹张开区域图5-6相互作用积分法计算结果图 西南交通大学硕士研究生学位论文第47页一Cmck-1在分析步模块下创建个新的历史输出,将输出变量设为裂纹裂纹尖端的应力强度因子,便能在输出文件中查得用相互作用积分法求得的裂纹尖端的应力强度-因子,当裂纹长度为2mm时,输出的应力强度因子结果如图57所示。'.KFA:CT0:MATES:RESTI:G)CRACKCRACKFRONTG;,NT0;.URSNAMEHODESET1t:A4H-- ̄0UTPUT2^CRACK1-5-:1535154.5154.5Kl,153,S''---%H..K2;1.iQ344g87EDS57204£:03aSE03S232E''----;--:K3.:.:4.99S2E035.2:757E03l9C26E〇366385E;0SMTSDIRECTION(DEG):lOOtfd.OOd;0.000(tOOCSJfromKs:t.1D21tl.1026H.10350.1D35图5-7相互作用积分法计算应力强度因子结果5.2.2.2扩展有限元法计算采用扩展有限元法进行裂纹分析,裂纹模型与平板模型被看做两个相互独立的个立模型一体,需要分别建,划分网格,然后再通过设置接触属性将两者相互作用到起。在对平板进行网格划分时,不需要考虑裂纹位置的影响,直接选用把节点六面体减缩59[]单元C3D8R对模型进行网格划分。网格的密度越大,计算结果的精度越髙。本文主50x250要为了研宄两种方法计算结果的比较,,为了减少网格密度对结果的影响采用5-8所示的网格密度对模型进行网格划分,所得到的有限元模型如图。■___1111腦鄉卿&___|111111K,:1_麵|1图5-8扩展有限元法有限元模型=此处也取裂纹长度a2mm为例对模型施加载荷和裂纹,采用扩展有限元的方法-对模型求解计算,得到的应力云图和裂纹区域局部放大图如图59所示。通过观察发现,这是,扩展有限元法所得到的形变模型中,在裂纹开裂位置处的单元发生了破坏由于在建立有限元模型时裂纹即穿过这些单元内部,裂纹模型与平板模型相互独立造成的。 西南交通大学硕士研究生学位论文第48页SMisesjv:(Ag75%)+2265+02■.e-207702+.6++.6+189002j-+十1.702e02-,十1.514e+02-「+.e+:■132702-+.e+■113902-■+9.51日e+01-+.e+g764201U-+5.766e+01?!+3:.B90e+01+.e+H201401+.+1387e00Y-.,-+ODB:JobxfemodbAbausStandard6121FriMarGMT08:002018q/Ste-p:Step1Increment1:StepTime=1.0C0ZXrmarar:sesPiyVS,MiDeformedVar:UDeformationScaleFactor:+8.094e+02(a)裂纹板应力云图7mb裂纹张开区域()图5-9扩展有限元法计算结果图当裂纹为2mm时ck-,用扩展有限元法求得的裂纹Cm1裂纹尖端的应力强度因子-1结果如图50所示。KFACTORESTIMATESCRACKCRACKFRONTCONTOURSNAMENODESET1234CRACK-1XFEM_1Kl:195.8161.016L4153.7K2--:08954305021341297....K3:3,20612-791L931,318MTSD--IRECTION(DEG):0.52392.1691.5159.516fromKs:,16620.11340.J0—113801031图5-10扩展有限元法计算应力强度因子结果5.2.3计算结果对比分析2mm4mm为了研究不同裂纹长度对于计算结果的影响,取裂纹长度a分别为,,6mm,8mm,10mm,分别采用相互作用积分法和扩展有限元法求得裂纹尖端的应力强度因子K的值5-2,再与解析法求得的应力强度因子理论值进行比较分析,结果如表所示。 西南交通大学硕士研究生学位论文第49页表5-2裂纹板应力强度因子计算值裂纹长度fl理论值相互作用积分法扩展有限元法0'5mmiPa-mm/CM模拟值误差模拟值误差0505-mrne-^/MPa7/%^/MPamme2/%214844147020961520624.....4243.37239.11.73248.972.363611310370482.555.37.7..58529.76520.981.65545.182.910796.61779.852.1817.362.6从表中数据可以得出以下结论:(1)使用相互作用积分法和扩展有限元法求得的应力强度因子值与理论值的误差均小于5%的允许误差范围之内,所以这两种方法求得的应力强度因子值均为可信的,可以根据实际情况选用合适的方法求解应力强度因子。(2)随着裂纹尺寸的增加,两种计算方法求得的结果的误差值并没有明显的变化。(3)相互作用积分法求解应力强度因子的值要比扩展有限元法求得的值的误差更小,结果更为准确,如果对应力强度因子的精度要求比较高的情况下,可采用相互作用积分法求得。(4)扩展有限元法求得的应力强度因子值要略大于理论值,而相互作用积分法所求的的应力强度因子模拟值要略小于理论值,如果要用来进行断裂评定可选用扩展有限元法模拟值用于计算,会使评定结果更为保守,结构安全性更为可靠。5.3臂架结构裂纹模型5.3.1局部实体模型为了更详细的反映臂架结构裂纹局部扩展情况,并节省计算时间,只需建立裂纹区域(即下弦杆与水平腹杆焊接处)的局部模型,从而研宄该处裂纹的扩展规律。由于本文第三章对臂架结构进行整体分析时,采用的是线单元建立的臂架结构的整体模型,故不能直接采用子模型的技术进行分析,需要借助三维分析软件CATIA重新建立-11所示下弦杆与水平腹杆焊接处的局部实体模型,如图5。将局部实体模型导入ABAQUS中进行网格划分,由于采用扩展有限元法计算,所以裂纹与实体模型是相互独立的,故不需要对裂纹尖端区域进行网格加密处理,采用六面体单元C3D8R将模型划分为15793个单元,臂架结构下弦杆与水平腹杆连接处局-部结构网格模型如图512所示。 西南交通大学硕士研究生学位论文第50页图5-11局部实体模型zi7r5-图12局部结构网格模型5.3.2裂纹模型对于研宄裂纹的扩展问题一,研宄对象是裂纹扩展这动态的过程。对于相互作用一一积分法,通过上述分析裂纹的过程可看出,需要把实体模型离散成组有限的且按定方式相互连接的单元组合体,并且在单元划分过程中必须严格要求单元的界面和裂一纹边界等几何界面保持致,为了保证计算结果收敛,需要对裂纹尖端进行局部的网一-所示格加密,如图55(b)。当裂纹进行扩展时,裂纹尖端单元会发生变形,每扩展次,需要对裂纹尖端重新划分网格,这大大增加了工作量;而在扩展有限元法中,裂-8中所示且并不纹模型与实体模型是相互独立的,如图5,裂纹边界可在单元内部,影响计算精度,故在裂纹的扩展中不需要再对实体模型重新进行网格划分,这在保证计算精度的同时也大大提髙了效率。同时,对于裂纹的扩展方向,在传统有限元法中,在加载裂纹时必须预先设定裂纹的扩展方向,裂纹只能沿着预设路径的单元边界进行扩展,而在实际情况中,裂纹的扩展路径是任意的,我们往往无法预知裂纹的真实扩展方向,;而利用扩展有限元法不需要预设裂纹扩展方向可以很方便的模拟沿任意方向扩展的裂纹模型,这与实际情况更为贴合。因此,本文选用扩展有限元法来对臂架结构进行裂纹扩展分析。根据第四章中分析结果-a,取预埋裂纹长度为10_,裂纹形状如图513()所示,5-13将裂纹装配在局部模型的下弦杆与水平腹杆连接处,如图(b)所示,裂纹类型选用XFEM,允许裂纹生长。设置分析步,打开非线性,设置场输出变量。 西南交通大学硕士研究生学位论文第51页<.>(a)预埋裂纹形状(b)裂纹位置5-图13裂纹模型5.3.3边界条件本文第三章己给出臂架结构的有限元分析结果,从结果应力云图可以看出,在工况三中,最大应力出现在短拉杆下端右侧的下弦杆与水平腹杆处,将此处作为裂纹局部161、163、模型研宄的区域,从臂架结构的整体模型中提取在局部模型边界上(节点287-1、289)的位移结果,设置局部模型的边界条件,如图54所示。EditBdaCnditiXounryoon|Name--iBC161Name;BC163Type:Displacement/RotationType:DispiacemenVRotationSS-ilS-ilteprtep1Statc,GeneraSteite1Statc,Genera()pp()--Region:Set7Region:Set6CSYSGlblCSYSGlbl:(oa)X:(oa)JLDisribuion:Uniformp权x》Disribuion:Uniform^权ttttjj|0U1--:55.86?80U1:64.3554||]0U2:2分.13340U2:33.58570U3:2.660620U3:2.80358j|radians□radians||UR1:UR1;丨|radians□radians1IUR2;UR2:丨|didiURraansURraans□3:□3:AlidRm\7AliudeRmptue:(ap)Fmpt;(amp)|卜卜|Note-:ThdisltlueillbNote:Thdltlilbepacemenvaweeispacemenvaueweitadibtt.tdibtt.maninensusequensepsmsmalnensusequensepsOKCancelOKCancel||| 西南交通大学硕士研究生学位论文第52页NameBC-287NBC-289:ame:Tttttttype:Displacemen/RoaionType:Displacemen/Roaion-Genera-Step;Step1Static,lSte;Ste1Static,General()pp()R----iSt4RiSt5egon:eegon:eCSYSrGlobalCSYS:GlobalJ()JL()Ltt"妖x^f^x)D-isribuion:Uniform}Distribution:Uniformp|--U..1:429257U1:579520|]0|0U2:24.77940U2:33.4217110U3:6.062970U3:6.9S6Q51]d国URidi1:raansUR1:raans1I||□|g□radiansradiansUR2:QUR2:[ ̄didi||UR3:raans□UR3:raans||Amplt\/Altd\yiude:RampFmiue:RampF()卜p()QNote:ThedisplacementvaluewillbeNote:Thedisplacementvaluewilbemainainedinsubseuenses.mainainedinsubseuenses.tqttptqttpOKCancelCancel[||OK|jj|-图514局部模型边界条件5.4断裂失效准则采用扩展有限元法来模拟裂纹的扩展过程中,需要定义两个重要的参数:材料的断裂失效准则和材料在破坏损伤中的相关参数。材料的失效准则作为材料单元发生破坏的起始判据一,对于宏观材料,断裂失效准则般可分为四类:基于能量的断裂失效准则、基于经验公式的断裂失效准则、基于应力应变的断裂失效准则以及基于损伤演61_63化的断裂失效准则[]。本文中,对于局部裂纹模型的裂纹扩展分析,采用最大主应力失效准则来作为材料损伤的起始判据,最大主应力设为84.4MPa。同时,采用最大主应力的失效准则需要设置材料的损伤演化规律和损伤稳定性系数,来保证计算结果的收敛性。对于损伤演化的类型选用基于能量的演化类型,线性软化、极大值退化、幂次法则指数型混合模式的损伤演化规律,法向方向和两个切向方向的断裂能为:===W-(422wm1。^仏仏.i,具体的损伤演化参数设置如图55所示损伤稳定性系_5数设为5.〇xl〇。*SubotionEditorpDamatgeEvoluionTpr.Enerygy日Softening;LirwafQdti:?DeraaonMaxtmufngj|MiddbhtPxemoeeavonowerLawrao:rrModemixtfEnegyjjSwiPo*f1Usetem-dedenttperafairependaa*brffiv:0Numeodd〇ft*bl*sDtaah?rod.SMNomlMod*ra__FfactureEnergyFractureEnergy"FirstDirectionSecondOirectiorOKCancel|j]图5-15损伤演化参数设置图 西南交通大学硕士研究生学位论文第53页5.5裂纹扩展结果通过ABAQUS扩展有限元方法计算分析,得到裂纹局部模型的断裂参数PHILSM-16-1变量的计算结果如图5所示,裂纹扩展单元状态如图57所示。'SSSS's'FWLSM刑-湖sIvIIJI丨_ICH4TN41不|、-KStep;Step1f\II7—?%/ncremen:eme=.lWki;TS=1It24StpTi1000/2*fi¥245tTCO?FTimarVar:PHUSMv'y』IIKDd"eformeVanUDefotiSFt100C00rmaoficaleacor:+.6+7X汗_」\@丨^_■图5-PHILSM16断裂参数计算结果和75S_::CC<3H<Mt■-i15-^1r■磅-,:3Si:2-■+lHT.:C(Ejl2-^St3f<:it?E3Ji1!ri331t^SfcpTKi!33fW-1?---'-1ri../r-T:jyT-.-,XD?^ris-iVaUCrf^n=i(SateFacir+:XD:?i3图5-17裂纹扩展单元状态-16和图5-17可以看出从图5,下弦杆与水平腹杆连接处裂纹主要受到拉力作用,为一I型张开型裂纹,与实际情况相致。同时,从裂纹的扩展路径可看出,裂纹逐渐向结合部位深处扩展一,扩展方向与裂纹并非沿着裂纹方向,并与裂纹方向形成定夹角,可在裂纹扩展路径上进行特殊处理来阻止裂纹继续扩展。5.6本章小结本章借助非线性分析软件ABAQUS对塔式起重机臂架结构裂纹扩展过程进行仿:真模拟,具体的结论如下(1)以含单边裂纹的有限宽度平板为例,分别运用相互作用积分法和扩展有限元法对该含裂纹结构的应力强度因子进行分析计算,并将计算结果与通过解析法得到的理论值进行对比分析,验证了两种方法所求结果的准确性。 西南交通大学硕士研究生学位论文第54页(2)通过对两种方法的对比分析,最终选用扩展有限元法对臂架结构上裂纹进行裂纹扩展分析。通过建立裂纹局部区域模型,得到断裂参数PHILSM变量的计算结果和裂纹扩展单元的状态,并验证了裂纹扩展方向的实际性。 西南交通大学硕士研究生学位论文第55页结论与展望随着建筑行业的蓬勃发展,作为建筑施工现场最重要的起重运输机械,塔式起重机在建筑作业中扮演着越来越重要的角色。随着服役塔式起重机年限的逐渐增多,由塔式起重机断裂引起的安全事故也越来越多,造成巨大的经济损失的同时也造成了严重的人身伤亡。在塔式起重机中,臂架结构为承载的主要部位,长期承受交变载荷,。结合近几年塔式起重机安全事故,臂架结构为塔式起重机断裂破坏发生的高危构件本文以某建筑施工工地服役的TC6020塔式起重机为研究对象,从结构完整性评定理论出发、,以断裂力学理论可靠性理论、损伤演化准则以及疲劳损伤理论等为基础,结合有限元分析结果和数值分析等技术手段,展开对塔式起重机臂架结构安全评定方法的研宄。论文的主要内容及结论有:(1)塔式起重机臂架结构断裂评定流程与方法。以某建筑工地使用TC6020塔式起重机臂架结构为例,对该塔式起重机臂架结构进行有限元分析;同时,以英国标准BS7910三级断裂评定准则为参考,结合臂架结构有限元分析结果和裂纹缺陷的相关参数,计算断裂评定参数,对臂架结构进行断裂安全评定,通过对绘制的失效评定图分析,该结构上裂纹缺陷均处在可接受的范围内,该臂架结构仍可正常运行。(2)塔式起重机臂架结构疲劳评定流程和方法。同样以TC6020塔式起重机臂架结构为研究对象,结合塔式起重机臂架结构实际的疲劳载荷工况和裂纹缺陷特点,采用线性累积Miner法则对塔式起重机臂架结构的随机载荷进行等效处理,结合结构的有限元分析结果,采用基于Paris曲线的疲劳评定准则对该臂架结构的疲劳寿命进行预测,并分析不同裂纹初始尺寸和评定曲线对于评定结果的影响。一(3)塔式起重机臂架结构的裂纹扩展分析。以含单边裂纹有限大含缺陷平板为例,分别采用相互作用积分法和扩展有限元法对应力强度因子求解,对比分析两种方法对于裂纹模拟的准确性。结合塔式起重机臂架结构的有限元分析结果,建立裂纹区域的局部有限元模型,采用扩展有限元法对含缺陷臂架结构进行裂纹扩展仿真分析,得到含裂纹臂架结构的断裂参数以及裂纹单元的扩展方向及状态,并验证了裂纹扩展方向的实际性。一目前,对于塔式起重机臂架结构的安全评定方法还未具有统的执行标准,本文在对TC6020塔式起重机臂架结构的安全评定的基础上研究适用于塔式起重机臂架结一构的安全评定方法和流程,有助于进步对于塔式起重机安全评定标准或规范的制定。但由于实际生产、使用具体情况的复杂多变,对于塔式起重机臂架结构的安全评定仍一是个复杂的问题,受限于作者自身水平与工程实践的经验,本文的研究还存在不足与需要完善的部分。具体后续的研究可从以下几个方面开展: 西南交通大学硕士研究生学位论文第56页(1)对于臂架结构的疲劳与断裂特性真实数据,可通过具体的疲劳断裂实验来获得,通过对于实际模型同比缩小,结合实际加载情况,可获得真实的疲劳与断裂数据,为本文的理论结果提供实验数据支持。(2)本文所采用的评定参数都使用的是常数值,但在实际工程使用中,由于材料的不均匀性、随机载荷以及在制造使用中人为因素的影响,评定参数往往是不确定的并服从一定的分布规律,可引入概率学以及可靠性理论,来研宄这些不确定性的规律以及对于评定结果的影响,从而更加真实地反映结构的真实情况。 西南交通大学硕士研究生学位论文第57页致谢一研宄生的时光转眼就要接近尾声,在此,衷心地感谢在这里遇到的每个人,是你们见证了我的成长。首先,在这里要感谢我的导师吴晓副教授,在论文选题、研究思路、课题实施、论文撰写及修改过程中,我得到了吴晓副教授的悉心指导。吴晓老师知识渊博、治学严谨、具有高度的责任心,论文的顺利完成无不倾注着他的心血。借此机会,真诚表达我对恩师的无比尊敬和由衷感激!非常感谢工程机械系可敬的老师们包括王金诺教授、程文明教授、管会生教授、王少华教授、于兰峰教授、张则强教授、黄松和副教授、王海波副教授和许志沛副教授等等,谢谢你们在研宄生三年里对我的帮助和指导。感谢西南交大,感谢交大机械学院及群体教职员工为我们提供的良好学习平台和生活环境,同时感谢熊刚、罗涛、罗友红、丁正男等师兄的指导帮助,感谢同门冯晓华、杨海艳和杨成洪、何更旺、高春丽、汪赢、袁雨阳、王昕、吴青科等师弟师妹的鼎力支持,感谢在研宄生期间鼓励和支持我的同学和朋友。一最后,由衷的感谢在背后直默默支持我的家人,是你们默默地辛勤付出,是你们毫无保留的爱一,是你们在我的每个人生节点处的鼓励与支持,才让我在人生道路上鼓足勇气、顺利前行!陈富谨致二零一八年四月 西南交通大学硕士研究生学位论文第58页参考文献1于兰峰.塔式起重机结构刚性及动态优化研究D.成都2006.[]:西南交通大学,[]2徐永春林云-王者静.塔式起重机结构件安全寿命评估J.建筑机械200311:5456.,,,[][]()32006-李联光.塔式起重机疲劳裂纹的无损检测[J].广东建材10:128129.[],()4陈国华.结构完整性评估M.北京:科学出版社,2002.[][]20068339-45董月香髙增梁.疲劳寿命预测方法综述[J].大型铸锻件:1.[],,,()6李舜酩.机械疲劳与可靠性设计M.北京:科学出版社,2006.[][]7吴晓罗薇刘璐等.在役桥(门)式起重机金属结构疲劳寿命预测分析[J].中国安全科学学报[],,,201020295-99:.,()8范小宁,徐格宁,王爱红.基于人工神经网络获取起重机当量载荷谱的疲劳剩余寿命估算方法[]J工程学报-?机械20114720:6974.[],,()9-.塔机结构件的疲劳性能分析和质量控制J.200527477.[]郑荣跃,徐永春[]建筑机械,():[10]董赞,蔡敢为,郑战光等.基于损伤力学的疲劳寿命和裂纹扩展的数值分析[J].中国机械工20-程,10,21(20:24122415.)11赵威威.基于试验的在役塔式起重机疲劳寿命分析D.大连:大连理工大学,2012.[][]12张朋李梅山孟宪颐等J2000542-44[].塔式起重机钢结构疲劳寿命估算.建筑机械:.,,[],()13郑夕健何劝云费烨.RBF神经网络和均匀设计在塔式起重机安全状态模式识别中的应用J.,,[][]机械设计与制造2005-l:7880.,()[14]王进.基于有限元法的塔式起重机钢结构疲劳寿命研宄[D].重庆:重庆大学,2005.15张宏伟江西建材-.塔机起重臂疲劳寿命的数值分析J.201615:284289.[][],()[16]张彦华.焊接力学与结构完整性原理[M].北京:北京航空航天大学出版社,2007.17MilneI,AinsworthRA,DowlinAR,etal.R6:AssessmentoftheInteritofStructures[]ggyContainingDefects,Revision3J.InternationalJournalofPressureVessels&Piin,1998,[]pg32-1:3104.()1J工程学院学报-8康颖安.断裂力学的发展与研究现状.湖南200616l:3942.,,[][]()19丁遂栋M工业出版社,199[].断裂力学[].机械7.20IrwinGR.AnalsisofstressandstrainsneartheendofacracktraversinalateJ.J.al.mech”[]ygp[]pp-195724:361364.,PI]中国航空研究院.应力强度因子手册[M].科学出版社,1993.22ChanSKsTubaIS,WilsonWK.OnthefiniteelementinlinearfracturemechanicsJ.Enineerin[][]ggFraM-ctureechanics197021:117.,,()23BarsoumRS.AlicationofuadraticisoarametricfiniteelementsinlinearfracturemechanicsJ.[]ppqp[] 西南交通大学硕士研究生学位论文第59页Inrnanaurnar114-tetiolJolofFacture9740:603605.,,()[24]ShihCF,LorenziHGD,GermanMD.Crackextensionmodelingwithsingularquadraticnnaurnau-isoarametricelementsJ.IternatiolJolofFractre1976124:647651.p,,[]()@5]刘明尧,柯孟龙,周祖德等.裂纹尖端应力强度因子的有限元计算方法分析m.武汉理工大学20-学报,11(33):116121.P6]V.库默,M.D.杰曼,C.F.施.弹塑性断裂分析工程方法[M].北京:国防工业出版社,1985.M出版社200607]程勒.断裂力学.北京:科学.,赵树山[],28周志新工程-?机械裂纹无损检测方法综述J.机电20173410:11391143.[][],,()92014R]孙环.塔式起重机起重臂疲劳裂纹扩展及寿命研宄[D].西安:西安建筑科技大学.,00]林华,赵毅.在用塔式起重机安全性鉴定中常见裂纹形式及处理方法[J].建筑安全,2046-48117:.()B-methodn31S79102005GuidetosforassessitheacceptabilitofflawsinmetallicstructuresS.[],gy[]Londonrndardnun11-466:BitishStasIstittio203:.,32赵章焰,吕运冰,孙国正.J积分法测量低碳钢Q235的断裂軔性KICJ.武汉理工大学学[][]20-报,12,244:111112.()[33]周凯旋,朱守寨,吴志纯.基于ABAQUS的动态塔机起重臂有限元分析[J].建设机械技术与管理20-143:105107.,()[34]王金诺,于兰峰.起重运输机金属结构[M].北京:中国铁道出版社,2002.[35]范俊祥,陆念力.塔式起重机[M].北京:中国建材工业出版社,2004.7.36GB/T-137521992.塔式起重机设计规范[S.[]]37KlinerC.FailuresofcranesduetowindinducedvibrationsJ.EnineerinFailureAnalysis,2014,[]g[]gg43-198220:.38何君儒.基于断裂力学的塔式起重机疲劳性能分析[D].南昌:南昌航空大学2013.[],39熊刚晓涛等.桥门式起重机箱形梁结构安全评定方法研宄J.机械[],吴,罗()[]强-度:2016386:12531257,()40王进.基于有限元法的塔式起重机钢结构疲劳寿命研究D.重庆:重庆大学,2005.[][]D541朱惠华.TC2040型塔机结构有限元分析.武汉:武汉大学.200.[][]42马灿,刘树林,温磊等.基于ABAS的塔式起重机结构分析及实验验证[J.机械制造,[]QU]20-1755634:8699.,()[43]ZaiJ,CaoJ,BellAJ.Thefatiguestrengthofboxgirdersinoverheadtravellingcranes[J].Sne440-405tructuralEnier1997223:1.g,,()[44]RettenmeierP,RoosE,WeiheS,etal.Assessmentofmixedmodecrackpropagationofcranerunwaubnnuren-yirderssectedtocyclicloadinJ,EieerinFractMechaics,2016,153:1124.gjg[]gg45Han,WanY,YinY,etal.DeterminationofstressintensitfactorformodeIfatiuecrackbased[]Qgyg 西南交通大学硕士研究生学位论文第60页onnnracureM201511-12finiteelementanalysisJ.EieerinFtechanics:86.[]gg,46KlinerC.FailuresofcranesduetowindinducedvibrationsJ.EnineerinFailureAnalsis,2014,[]g[]ggy43-:198220.[47]熊刚.基于结构完整性理论的桥(门;)式起重机箱形梁结构安全评定方法研究[D].成都:西南交通大学,2016.[48]罗涛.含缺陷起重机箱形梁结构的弹塑性断裂参数模拟及安全评定[[D].成都:西南交通大学,2016.49罗友红.桥(门)式起重机箱形梁概率安全评定与疲劳寿命估算D.成都:西南交通大学,2017.[][][50]T.Belytschko,T.Black.Elasticcrackgrowthinfiniteelementswithminimalremeshing.IntetilJlf45-onaournaorNuicalMethodsinEnineerin.1999.5601620.rnamer():gg51J.Dolbow,T.Belytschko.Afiniteelementmethodforcrackgrowthwithoutremeshin.[]g-InternationalJournalforNumericalMethodsinEnineerin.1999461:131150.gg,()52nSSYantTddtress-intitfactorfinfuntillWaJFCorenHDosRH.Senssorsurfacecrackscona[]gs,,yyradedmaterialsundermodeltnermomenchanicalloadinJ,InternationalJournalofSolidsandgg[]Structures,2004.53宫经全张少钦李禾等.基于相互作用积分法的应力强度因子计算J.南[],,[]昌航空大学学报,20-15291:4248.,()[54]庄茁,柳占立,成斌斌等.扩展有限单元法[M].北京:清华大学出版社,2012.55DauxC,MoesN,DolbowJ,etal.Arbitrarybranchedandintersectincrackswiththeextended[]gfiniteelementmethodJ.InternationalJournalforNumericalMethodsinEnineerin2000,,[]gg481211-1:74760.()[56]余天堂.扩展有限单元法[M].北京:科学出版社,2014.57BabuskaI,MelenkJM.ThePartitionofunitmethodJ.InternationalJournalforNumerical[]y[]thin404-MeodsEnineerin1997727758:.gg,,()58董玉文任青文.基于XFEM的混凝土开裂数值模拟研宄[J].重庆交通大学学报:自然科学版[],,200928-401:36,,()余天堂-59董玉文.直接计算应力强度因子的扩展有限元J.计算力学学报2008251:7277.,,,()[][][60]SMohammadi.Extendedfiniteelementmethod[M].London:BlackwellPublishingLtd,2008.61李彩霞.基于扩展有限元法的裂纹扩展分析研宄D.成都:西南交通大学,2015.[][][62]杨志锋,周昌玉,代巧.基于扩展有限元法的弹塑性裂纹扩展研宄[J].南京工业大学学报:自然2014364-科学版:5057.,,()63余天堂-.模拟三维裂纹问题的扩展有限元法J.岩土力学20103110:32803285.,,[][]()[64]张彦华.焊接力学与结构完整性原理[M].北京:北京航空航天大学出版社,2007.65GB/T-137521992.塔式起重机设计规范[S.[]] 西南交通大学硕士研究生学位论文第61页攻读学位期间发表的学术论文及科研成果1.发表学术论文11.陈富吴晓杨海艳冯晓华.基于BS790的塔式起重机臂架结构安全评定方法[],,,研究[J].起重运输机械(已录用).2.参与科研项目“1.四川省国际合作计划项目,服役起重机金属结构强度退化机理与安全性能评估[]”研宄,2014HH0022。'