硕士学位论文-大型造纸机基础动力设计的理论研究与应用.pdf

湖南大学硕士学位论文大型造纸机基础动力设计的理论研究与应用姓名：吴昊申请学位级别：硕士专业：工程力学指导教师：刘腾喜20090418硕十学位论文摘要造纸工业是国民经济重要基础原材料产业之一。近年来，我国引进了若干国外先进的大型造纸成套设备，其中大型造纸机是最关键的设备。造纸机运行中要受到很大的动力荷载的作用，因此对造纸机基础的动力设计研究随之开展。本文以实际工程项目为背景，重点介绍了造纸机基础的动力分析方法和步骤，对工程实际有重要的应用价值。本文首先利用有限元方法计算了单桩和群桩的动力刚度，并和动力机器基础设计规范中公式计算得到的结果进行了比较。有限元结果验证了地基刚度和阻尼的频率相关性，并证明按动规方法求得的地基刚度和阻尼可用于工程实际。然后采用子结构法，通过动规方法求得桩基刚度和阻尼，对大型造纸机的桩基机器系统进行了动力分析。重点介绍了动力分析的基本步骤，并对振动系统的地基的模拟和荷载的处理做了详细描述。通过对振动系统的模态和简谐分析和得到的速度响应谱曲线确定了满足纸机制造商动力设计要求的纸机基础设计。最后讨论了在造纸机动力分析中的相关问题，包括基础模型的简化、有限元中阻尼矩阵的确定以及对于处理造纸机某特殊部分动力分析的方法。这些问题的处理方法已用于多项实际项目的动力设计中，经过了实践的考证，可以作为同类项目设计的参考。关键词：造纸机；桩基础；模态分析；简谐分析；子结构法；阻尼矩阵犬型造纸机基础动力设计的理论研究与应用AbstractPaper-makingindustryisoneoftheimportantbasicmaterialindustriesinnationaleconomyInrecentyears，severalgiantpaperequipmentshavebeenimportedfromabroad，inwhichthegiantpapermachinesarethekeycomponentsThepapermachineissufferedtogreatdynamicloadsinoperation，SOthedynamicdesignandresearchofthefoundationofpapermachineareveryimportanttasksBasedonapracticalsubject，thedynamicanalysisiscarriedoutprimarilyinthispaperwhichisofgreatbenefittotheengineeringapplicationUsingthefiniteelement，thedynamicstiffnessofthesinglepileandgrouppilesarecalculated，andthecomparisonofthepresentresultswithonesintheDynamicMachineBaseDesignSpecification)ismadeSincetheinDynamicMachineBaseDesignSpecificationisbasedonthemassspringdampingmodelTheresultsfromFEMshowsthatthestiffnessandthedampingoffoundationswiththecharacteristicofthefrequencydependenceandthestiffnessandthedampingcalculatedaccordingtoDynamicMachineBaseDesignSpecification)canbeusedinengineeringpracticesInthispaperthestiffnessofthepilefoundationanddampingisfirstlycalculated，andaccordinglythedynamicanalysisofthepilefoundationmachinesystemofthegiantpapermachineisinvestigatedbyadoptingthesubstructureThenthebasicstepofthedynamicanalysisandthewaytomodelthefoundationofthevibrationsystemandtheloadareprimarilyandelaboratelyintroducedBasedonthespectrumresponsecurveobtainedthemodalandharmonicanalysisofthevibrationsystem，thebasicpaper-makingmachinemodelmeetingthemanufacturingcompanySdynamicrequirementonthemachineisreasonablydesignedAtlastsomeproblemsinthedynamicanalysisarediscussed，whichincludethesimplificationofthebasicmodel，thedeterminationofthedampingmatrixofthefiniteelementmodelandthedynamicapproachofsomespecificcomponentsofpapermachinesThetreatmentsofsuchproblemshavebeenutilizedinthedynamicdesignsofseveralpracticalsubjectsTheengineeringpracticeshowsthatthedynamicanalysisschemeinthispapercanbeusedforthedynamicdesignsoffoundationsofpapermachinesII硕十学位论文Keywords：giantpapermachine；pilefoundation；modalanalysis；theharmonicanalysis；substructure；dampingmatrixIII湖南大学学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名：美景日期：汐呷年f月加日作者签名：云灵日期：汐呷年s月加日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权湖南大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。本学位论文属于1保密口，在年解密后适用本授权书。2不保密囹。(请在以上相应方框内打“)作者签名：导师签名：日期：日期：厂月汐日岁月阳日硕lj学位论文11研究背景及意义第1章绪论造纸工业是国民经济重要基础原材料产业之一，纸及纸板的消费水平已成为衡量一个国家现代化水平和文明程度的一个重要标志。在经济发达国家，纸及纸板消费量增长速度与其国内生产总值增长速度同步。造纸工业在现代经济中所发挥的作用已越来越多地引起世人瞩目，被国际上公认为“永不衰竭的工业，在美国、加拿大、日本、芬兰、瑞典等经济发达国家，造纸工业已成为其国民经济十大支柱制造业之一。造纸工业是一个与国民经济发展息息相关的重要产业，它涉及林业、农业、机械制造、化工、电气自动化、交通运输、环保等多个产业，对上下游产业的经济有一定拉动作用。同时，造纸工业又是一个技术密集、工艺复杂、资源消耗量大、产生污染物多的工业。20世纪下半叶以来，中国造纸工业迅速崛起，取得了举世瞩目的成就。目前，我国纸及纸板的生产量和消费量均居世界第2位，仅次于美国，产品种类齐全，质量已达到或接近世界先进水平，在世界造纸工业竞争格局中已具有相当的影响力。但是，我们在充分肯定成绩的同时，必须看到我国造纸行业对环境污染的日益严重性，这已经引起了全社会的关注。随着世界经济的发展，纸产品消费的增长，今后全球纸业的发展将面临资源约束、环境压力等问题，我国造纸工业面临的这些问题尤显突出。优化结构，整合资源，引进先进技术与装备是缓解上述矛盾的有效办法。采用现代化大型造纸设备具有生产率高、纸张质量好、对环境污染小等优点。近年来，我国引进了并仍将引进若干套国外先进的大型造纸成套设备，其中大型造纸机是最关键的设备n。81。大型造纸机通常采用框架式基础，包括地面以上的框架，地面以下的基础梁或板(即群桩的承台)以及桩。沿机器长度方向(下面称机器方向)，通常隔20-30米设置沉降缝，沉降缝将基础分为若干部分，各部分的机器、基础(包括框架及桩)以及地基构成了独立的振动系统。造纸机机器、基础以及地基(包括桩和土)构成的振动系统，属于典型的土一结构动力相互作用问题。根据动力设计的要求，应对于上述的每一部分(振动系统)建立包括机器、基础以及地基的分析模型。目前，国内引进了相当数量的先进大型造纸机，而机器的基础是不能引进的，需根据场地条件进行设计。基础设计中最主要的内容是动力设计。如果设计不周，轻则不能保证纸张质量，重则极易造成机器的损坏。因此，国外机器制造商对基础的动力设计提出了严格的要求。如果基础设计由国内设计人员承担，大型造纸机基础动力设计的理论研究0应用外商将对动力设计成果进行严格审查，这也从一个方面证明了基础动力设计的重要性。基础动力设计的实质是系统的动力分析，初步设计一般很难满足机器制造商提出的动力要求，通常需进行多次动力分析，才能完成动力设计的工作。111土一结构动力相互作用研究意义大型动力机器基础、核电站结构物、桥梁、海上采油平台以及高层建筑等结构要受到诸如机器扰力、地震、风和波浪等动力作用。早在三十年代后期，人们就认识到在动力作用下，上部结构和地基是相互影响的关联系统，应作为整体系统来研究其动力行为。以动力机器基础为例，动力机器产生的动荷载(例如不平行惯性扰力)使机器及基础产生振动，振动能量一部分通过地基传递，使地基变形相当于弹簧阻尼系统，于是机器、基础和地基组成了复杂的振动系统。动力机器的基础设计则建立在地基一基础一机器系统的动力分析基础之上1。以地震作用为例，一方面震源产生的地震波，经过多层介质的反射、折射和滤波作用，将振动能量传给结构物，引起结构物的振动；另一方面，结构振动时，又将一部分能量回输到地基中。这样，上部结构物和地基形成了一个复杂的动力系统，结构在地震作用下的行为是这个复杂动力系统的综合响应。大量的理论分析和震害调查都表明，土一结构动力相互作用对结构的地震响应具有重要的影响。因此，土一结构动力相互作用的研究受到工程界高度重视，一方面其研究成果可直接为工程建设服务，具有重要的应用价值；另一方面其研究思想和方法丰富了工程力学的内容，具有重要的理论意义n0。1。112土一结构动力相互作用及效应对于最一般的情况(例如地震问题)，土一结构动力相互作用问题可以按以下三个方面来描述n2l。(1)若建筑物场地为覆盖在基岩上的土层，土的放大和滤波作用将改变由基岩输入的地震波幅值和频率成份。这种未建建筑物以及未经开挖的场地运动称为自由场运动，由基岩输入运动来确定自由场地面运动的过程称为自由场反应分析。(2)建筑物的存在将改变上述自由场运动。以箱形基础为例，基础相对于土而言可以认为是刚性的。土与基础交界面上各点的运动应近似满足刚体运动的关系，或者说基础相对于上述自由场而言是一种位移边界条件。显然，这种刚体运动关系或位移边界条件将改变自由场的运动，这种作用一般为运动相互作用。(3)施加于结构的惯性荷载通过基础与土发生作用，使土体变形，从而构成了土一结构动力系统。在惯性荷载作用下，该动力系统的动力相互作用一般称为惯性相互作用。2硕士学位论文土一结构相互作用的效应体现在以下三个方面：(1)建筑物的存在将改变自由场运动，即运动相互作用效应。在实际分析中，考虑到地震输入的很多不确定因素，一般不考虑运动相互作用效应，而直接采用自由场运动作为土一结构惯性相互作用分析的地震输入。但如果结构的存在明显地改变了自由场运动，则这种近似分析可能导致较大的误差。(2)相对于基础固定的情况(刚性地基假设)，地基上的变形使土一结构系统固有频率和振型都发生变化，软弱地基上这种变化更加显著。显然，动力系统这些重要动力特性参数的改变将对系统的动力反应有重要的影响。(3)振动能量中相当一部分能量通过地基土介质的阻尼而消散，对于近似于弹性半空间的场地而言，这种阻尼作用将使结构的动力反应减小很多。如果为动力机器基础问题，则仅有上述的惯性相互作用。对于软弱地基上采用群桩基础的框架式基础，一般桩基承台埋置于土中，因此必须考虑土、桩和承台的相互作用。因此，土一桩一承台一框架一机器组成了复杂的动力系统，尽管只有惯性相互作用其分析也是十分复杂的。12国内外相关领域的研究现状121单桩一土动力相互作用研究现状在桩基动力学中，研究的最早、较为成熟的问题是单桩一土系统在动力荷载作用下的相互作用。根据桩周土计算模型和模拟方式的不同，现有的单桩一土相互作用分析理论和研究方法可归纳为解析方法、离散模型方法、数值方法和实验方法四类n31。(1)解析方法解析方法通常以线弹性或粘弹性均匀连续介质中的三维波传播理论为基础确定地基反力。这类方法能正确地表示几何阻尼和土层的共振现象，但无法反映桩一土界面上的几何非线性行为。Novakn4叫副等得到了桩在竖直方向和水平方向简谐荷载作用下桩身的动力响应。竹宫宏和与山田善n61在Novak方法的基础上针对层状土提出了一种有效的近似解析方法。刘宗贤和李玉亭n71给出了桩基础的横向自振特性以及在横向动荷载与地震荷载作用下强迫振动的解析解。王慧、杨光辉和张鸿儒引建立了桩一土运动相互作用连续介质力学模型，获得了桩周土的位移和应力解析表达式。何玉敖n引针对水平振动时可能出现的桩一土局部分离现象，提出了具有第三维变形特征的无限平面内圆孔辐射应力混合边值问题的计算理论。(2)离散模型法离散模型法起源于Winkler地基梁模型，这种方法将桩视为埋置于土介质中3大型造纸机基础动力设计的理论研究j应用的梁，忽略土的连续性而将桩周土的阻抗效应用分布的相互独立的弹簧和阻尼器代替，因而可以考虑土性沿深度的变化以及材料的非线性性质。这类方法简便实用，物理概念清楚，计算工作量小，因此，在桩基理论中一直受到广泛重视。同时它也存在着一定的局限性，如弹簧和阻尼器系数的取值以及参振土体体积或质量的确定等问题尚未得到很好解决，不能详细描述土中应力波的传播、土的屈服破坏的发展过程和桩一土界面上的破坏等复杂的物理现象。Novak”们利用Winkler模型求出了桩的动刚度和阻尼，并应用于摩擦桩的垂直振动分析中他。后来，Novak和AboulElla乜21对此进行了改进，改进后的模型不仅考虑了地基的几何辐射阻尼，而且考虑了材料滞回阻尼的影响。为了研究桩周土体的非线性、桩与土的相对滑移和分离等复杂的物理现象，Novak和Sheta陋31将桩周土体定义为弱化区，并将弱化区内的剪切模量和阻尼比分别加以适当的降低和增大。近年来，Naggar和Novak基于Kondner他41建立的非线性应力应变关系提出了非线性桩一土相互作用模型，并进行了时域计算凹引。Gazetas、Fan和Kaynia心引进行了全面的变动参数研究和桩一土相互作用动态分析。Galzetas恤73提出了动力相互作用因子的概念，对均质和非均质土层中的桩一土动力相互作用问题进行了深入探讨。陈云敏、陈仁朋和朱斌心81给出了单桩在桩顶偏心竖向动力荷载作用下的动力响应解。蒯行成、沈蒲生和陈军心刚提出了层状土中单桩动力阻抗的计算方法。范敏、解明雨和邬瑞锋凹叫分析了土的非线性性质对桩一土一结构相互作用体系地震响应的影响。(3)数值方法在桩一土动力相互作用分析中，广泛使用的数值方法主要包括有限元法和边界元法及其耦合方法。其中有限元法相对比较成熟，得到了广泛应用，目前已由简单的线性分析发展至复杂的非线性分析，由平面分析发展至三维分析。在有限元分析中，一般将桩基础附近较大范围的地基土与桩一起视为整体，在基底面输入地震动作用，直接计算桩基和地基土的地震响应。这种方法能够用相容单元分别模拟地基土和桩体，可以真实地考虑桩基和土体的力学性质、复杂的几何形状以及任意荷载作用，使许多复杂的实际问题得以求解，但是所需计算机容量较大，计算工作量大。而边界元方法只需对桩一土界面进行离散，同时将桩离散为梁柱单元，因此，计算量相对较小。但这类方法一般仅适用于线性分析，迄今未能应用于非线性分析。Blaney口利用轴对称有限元方法对桩的动力特性进行了三维分析。Berger翱利用有限元程序ALUSH对地震荷载作用下的桩一土相互作用进行了数值分析。考虑桩的弯曲特性，Kagawa和Kraft凹33对有限元程序ALUSH进行了改进，并发展了新的有限元程序PLLUSH。Angelides和Roesset41利用发展的非线性有限元程序对桩一土体系的动力相互作用特性进行了研究。赵振东和傅铁铭口51对在施4硕十学位论文加于桩顶的侧向脉冲动荷载作用下桩的非线性动力特性进行了分析。雷国辉、赵维炳和施建勇阳63利用有限元方法对单桩的锤击贯入过程进行了数值模拟。张崇文、赵剑明和张社荣71提出了一种求解桩一土相互作用体系动力非线性反应的动力层元分析模型。姜忻良、黄艳和丁学成81用样条有限元法分析了承受建筑物荷载的桩基和土体，用半解析无限元模拟桩基周围半无限土体，分析了相邻建筑物的相互作用。肖晓春、迟世春和林皋口钉对横向荷载作用下的单桩和群桩分别用协调单元和非协调单元进行了有限元模拟计算。(4)实验方法由于桩一土动力相互作用问题的复杂性，以往的研究工作主要局限于理论分析和数值模拟，缺乏实验验证。近年来，动力相互作用的实验研究日益受到重视。桩一土动力相互作用问题中最有意义的实验研究是现场的强震观测和现场原型实验，但由于地震的不可预见性，有针对性的强震观测难以人为控制实验进度，等待周期往往很长。而现场原型振动实验成本很高，至今得到的有价值的实验数据非常有限。为了更深入地认识地震荷载作用下桩一土动力相互作用的机理，必须了解桩一土动力相互作用的规律，进行桩一土动力相互作用体系的振动台模型实验研究成为必要，这将为验证桩一土动力相互作用分析理论和计算方法的合理性与可靠性提供参考依据。Kobayashi、Yao和Yoshida们通过实验探讨了桩一土间的相对滑移和相对分离以及桩周土的非线性变形特性等复杂的物理过程。NomuraH对可液化土中的桩一土相互作用效应进行了振动台模型实验研究。楼梦麟和王文剑H列等通过振动台模型试验，探讨了桩一土相互作用对结构动力特性和地震反应的影响。122群桩一土动力相互作用研究现状在实际工程中，桩基往往是以群桩的形式出现的。当群桩桩间距很小时，群桩效应会对单桩的刚度和阻尼产生很大的影响。对于静载下的桩群，如果某一根桩位于其相临桩位移场范围内，则该桩位移会由于相邻桩位移场影响而增加，因此，群桩的总位移总是大于单桩在平均荷载作用下的位移，这里所说的平均荷载是指群桩荷载除以群桩的桩数。在动力荷载下的群桩与在静力下的类似，只是群桩的动力阻抗与频率相关。目前，关于群桩动力特性的分析方法可分为两大类，即直接分析方法和近似分析方法。直接分析方法，将整个桩群作为一个整体进行分析，同时考虑桩群中各桩的相互作用效应，该方法具有物理概念清楚的特点。基于叠加原理的近似分析方法，首先利用相互作用系数的概念考虑群桩中每两根桩之间的相互作用效应，然后将群桩中每两根桩之间的相互作用效应进行叠加，得到整个桩群的动力响应。此法简单，计算工作量较小，便于工程实用而受到广泛重视。Polous和DavisH3。443首先大型造纸机基础动力设计的理论研究j应用引入相互作用因子的概念来考虑两桩之间的静力相互作用，并证明群桩静力效应可以通过叠加两桩的相互作用效应而求得，对群桩的相互作用问题做了开创性的工作。自Poulos提出静力相互作用系数、群桩效应后，有人借助于边界积分方程、有限单元法和简化近似方法求得了桩顶各自由度的相互作用系数，并得出在静荷载下，群桩效率系数总是小于或等于1。但Kaynia和Kausel等(1981)一些学者的研究表明，动力群桩系数是随频率波动的，可以是正的或负的，可以大于1或小于1。因此，在激振频率较高的情况下，静力相互作用系数对群桩动力稳态反应不能提供有用的信息。为此，Kaynia和Kausel引提出了动力相互作用系数的概念，认为群桩动力相互作用效应也可以通过叠加桩一桩两两动力相互作用效应来估计。Dobry与GazetasH鲥提出了一种计算均质土中动力相互作用系数的简化算法，其动力相互作用系数的定义和Poulos的静力相互作用系数的定义式是一致的，仅仅是用动位移取代了静位移。Gazatas和Markis7叫81分别对轴向和侧向荷载作用下的动力相互作用分析方法进行了进一步的完善，并将其推广到非均质地基中。这些研究均是在线性分析理论基础上建立的，不能直接应用于动力非线性分析，因而具有较大的局限性。桩的动力研究虽然是基于静力分析成果而开展起来的，但是由于两者所涉及问题领域不同，以及它们研究目的不同，从而导致所采用分析手段也不尽相同，所以两者发展并不互为依附关系。迄今为止，桩基动力特性研究也渐趋成熟，由最初定性概念逐步向量化计算分析过渡，群桩的动力相互作用要比单桩的动力分析复杂的多，因此直到十九世纪七十年代才用弹性波动理论对群桩的动力效应进行研究。一个比较实用的分析方法就是考虑有效土体质量，在文克尔地基桩土系统模型中，通过有效土体质量来考虑土体惯性力作用，但是这种方法的缺点是没有考虑桩土相互作用的辐射阻尼，而辐射阻尼是桩土相互作用动力分析的一个重要参数。正是群桩中桩一土一桩相互作用，使得群桩的动力分析极为复杂，近年来国内外对群桩阻抗的研究受到了重视，并做过一些群桩在动力荷载下的实验。Novak等人H叫通过等效群桩静力因子(考虑相互作用的频率因子)的方法来反映群桩动力效应，这种调整方法使得数值结果与精确解有了较为一致的逼近，但仍不能合理地解释桩土间的波动特征过程。Nogamit哺们在利用等效单桩的概念求解两桩间动力反应时提出了动力作用效应因子的雏形。Mylonakis及其合作者喵系统地将这种采用计算方法应用到桩群竖向、水平及摇摆各项阻抗分析中，对均匀和非均匀介质中桩一土一桩动力相互作用问题进行了较全面而广泛的研究，取得了共识性的成果。蒯行成等瞄剐进一步地将上述思想推广到了层状介质，这种以简化动力文克尔模型计算桩土体系的动力相互作用因子法是比较适用于工程实际的。另外对群桩也作过一些试验，Novak和Grigg(1976)、ElSharonouby和Novak(1984)、6硕一f二学位论文ElMarsafawi(1992)分别测试了群桩在动力简谐荷载下的位移，实验结果表明，位移与激振频率极为相关。Hasini和Woods(1989)测试了不同桩间距群桩的反应，实验表明，桩间距4d在1014以下时，群桩效应是不能忽略的。JamesPBurry和MichaelJPender33等人，在不同的地基条件下，改变桩间距，得到一些对工程实践有指导意义的结果，对于硬土地基当桩间距s朋12时，群桩效应可以忽略，对于软土地基当桩间距4d16，群桩效应可以忽略。123地基一桩基一上部结构动力相互作用研究现状土一结构动力相互作用问题的研究内容十分广泛，文献非常多，动力阻抗的计算方面，关于圆形、矩形、任意地面形状、明置、埋置等基础的动力阻抗己研究得相当成熟了。地基的变形将影响土一结构系统的动力特性和动力反应，地基土越软弱，这种影响越大，而软土地基广泛采用桩基础。因此，地基桩基上部结构动力相互作用更值得重视，此外，桩一土一桩的动力相互作用比其它类型的基础(如箱基、筏基)与土的相互作用要复杂得多。土一结构动力相互作用分析方法可以分为两大类，即直接法和子结构法。直接法是将结构、基础以及周围有限范围内的土体作为整体分析对象，将地震波作用在虚拟边界上，这种直接法通常采用有限元对整个研究对象进行离散，为了模拟波向无穷远处辐射的半无限域特征，需在土体虚拟边界上设置人工边界。一般情况下，对地基一桩基一上部结构动力相互作用直接方法分析，是将土(半无界域)、桩、承台、上部结构统统离散为有限元，这种方法无须作更多的简化，是一种能处理较为复杂系统的方法，但由于地基为三位半无限域，其计算工作量十分庞大，因此这种全数值方法并不能取代子结构方法。工程中广泛采用子结构法，将土一桩一承台作为一个子结构，该子结构为一半无界域系统，上部结构作为另一个子结构，土桩系统对上部结构的作用可等效为弹簧阻尼器系统，弹簧等效刚度和阻尼可通过地基动力阻抗的概念求得，一旦求得了地基动力阻抗，惯性相互作用就类似于通常的有限域振动问题。因此，采用子结构方法，如何求得地基系统的动力阻抗成为关键问题，大多数的相互作用研究也集中在这一方面。一般情况下，由于运动相互作用对自由场运动的改变并不显著，因此，地基一桩基一上部结构动力相互作用研究工作大部分集中在惯性相互作用方面，且其中大部分文献是关于动力阻抗计算方面的内容。确定桩基动力阻抗的方法可分为两大类，即较精确的数值方法和各种简化方法。数值方法包括有限元和边界元等方法，这类方法的特点是所作的假定比较少，因而从简历分析模型的角度看，数值方法被认为是精确的方法。数值方法分析桩一土一桩相互作用的优点是能较方便地考虑土的非均匀性，分层等因素，特别是有限元方法较便于考虑非线性因素，数值方法由于所做的假定较少，其结构被认为是较精确的(理论意义上)。但由于7人型造纸机荩础动力设计的理论研究j应用其计算量十分庞大，输入数据多，因此通常用于理论研究。简化方法的种类较多，其主要特征是在建立分析模型时作了或多或少的各种假定，不同方法之间的区别主要是所作假定不同。一般而言，各种简化方法最后也要进行数值计算，目前简化方法中广泛采用动力文克尔地基梁模型(BDWF)和动力相互作用因子(DynamicInteractionFactor)来计算单桩和群桩的动力阻抗。采用数值方法分析桩一土一桩相互作用的工作相当多，由于边界元特别适用于半无限域边界，因此80年代以来，采用边界元(包括数值边界积分法)研究这类问题的工作得到重视和发展。这里列出了一些有代表性的近年的工作，文献54利用层状介质中的格林函数，通过积分变换技术和数值边界积分，建立了求解层状介质中桩一土一桩动力相互作用的分析体系并编制了程序，该程序既可用于动力阻抗的计算，可以进行运动相互作用分析，并常被用来计算动力相互作用因子。文献55采用边界单元法对均匀介质中的群桩进行分析。文献56采用边界单元法研究土性随深度作线性变化的非均匀介质中群桩的动力阻抗。文献57提出了一种准三维有限单元来分析桩基础，这种准三维单元可以大大减少计算时间，使其能用于实际计算。文献E58将这种方法推广到桩基的非线性分析，分析中考虑土剪切模量和阻尼随应变水平改变、屈服等非线性效应，对一试验单桩和22群桩进行了数值计算，并与试验结果进行了对比，结果是令人满意的。计算桩动力阻抗的另一类重要的方法为各种简化方法，这类方法在深入分析桩一土一桩相互作用机理的基础上，做出某些简化，提出较简单但在物理上合理的模型，计算简单，结果良好，便于实际应用。数值方法和简化方法是研究相互作用问题的两条途径，都值得重视和继续深入研究。在单桩分析中，广泛采用动力文克尔地基梁模型哺引，用动弹簧系数和阻尼系数来模拟桩周土的作用，不同的作者采用不同的方法来确定弹簧和阻尼系数阳0。611，用解析方法和其他的近似方法计算单桩的动力阻抗。对于群桩，简化方法中采用动力相互作用因子来考虑群桩效应，动力相互作用因子一般利用数值方法计算，其性质与静力相互作用因子完全不同。针对目前实际中仍广泛采用静力相互作用因子来计算群桩动力阻抗的现状，文献62认为群桩中桩一土一桩之间的作用是强烈的，采用静力相互作用因子的结果还不如不考虑群桩效应的结果(即认为群桩的刚度等于各单桩刚度之和)。而动力相互作用因子法是目前可用于工程实际的唯一值得推荐的方法。动力相互作用因子一般采用数值方法计算，但文献中可查到的数值结果十分有限，而要编制可适合于各种桩数和桩间距等情况的计算动力相互作用因子的图表资料是不现实的，因此研究计算动力相互作用因子的简化方法近年来受到重视阳3。871。但目前这些方法只能用于以下的情形，即水平振动只能用于均匀介质，轴向振动只能用于均匀或土性随深度为线性或幂次变化的非均匀介质。简化模型也被用来进行桩一土的非线性动力相互作用分析，文献68-69对均匀介质中的单桩，引入桩周8硕士学位论文软化层的概念来考虑土的非线性以计算单桩的动力阻抗，在此基础上通过引入等效的线性相互作用因子来考虑群桩效应以计算群桩的动力阻抗。文献70利用上述的动力阻抗对一座近海钻塔的模态特性进行分析，这是采用简化方法研究考虑非线性效应的土一结构动力相互作用比较系统的工作，是比较前沿的课题。国内研究桩基动力阻抗的工作要少得多。文献71-74采用Novak等人的分析方法对单桩及群桩进行了分析，并将分析结果与若干试验结果进行了对比。由于桩基试验工作并不多，故这些试验结果是很宝贵的。文献75对软土地基上加桩的块式基础(用于动力设备)进行了研究，计算了动刚度，分析中考虑了回填土的影响，但没有考虑群桩效应。对于低承台群桩基础，承台埋置在土中，在计算地基系统的动力阻抗中，除了土一桩相互作用外，土与承台的相互作用必须考虑。目前关于土一桩一承台动力相互作用的研究并不多见，事实上，土一桩一承台静力相互作用的问题就已经很复杂了，文献76-84研究了群桩一土一承台的共同相互作用，并利用承台对群桩刚度的贡献对一些桩基事故进行了处理。对于群桩一土一承台动力相互作用研究比较复杂。而这一问题在软土地区采用群桩基础的动力机器基础设计中必须进行研究，这方面报道十分罕见。124大型造纸机动力设计国内外研究现状近年来，我国的大型造纸设备基本上由国外引进。为了保证机器的正常运行以维护机器厂家自身的声誉和经济利益，外商一般要求由他们进行机器基础的设计，而设计费相当高昂。如果由国内自行设计，则要求设计单位提供详尽的动力分析结果，这些动力分析内容必须经外商审核认可。在此背景下，近年来国内的有些单位为引进的机器进行基础设计，国内在这方面的研究也随之于展开。但是国内在这一方面的研究还处于初步阶段，相关的理论研究与分析手段还有待发展。特别是基础的动力分析，目前国内的设计单位无法完成。为了满足国内造纸业与相关领域生产设备更新时对配套辅助设施提出新的要求，目前国内已开始这一方面的研究，也取得了一些积极的成果，在这一方面的研究已接近国外先进水平。13本论文的研究内容本文以实际工程项目为背景，研究土一桩一承台一机器系统的动力相互作用。本文首先介绍了桩基础的理论和动力阻抗计算方法以及动力分析中的模态分析理论和简谐响应分析理论。然后比较有限元法和动力机器基础设计规范阳51(以下简称动规)公式计算得到的地基刚度结果，用实例验证了动力分析中地基刚度与阻尼的频率相关性。本文重点介绍了根据纸机制造商提供的造纸机基础资料利用有限元软件ANSYS对其进行建模并进行动力响应分析，用实例详细介绍9大型造纸机基础动力设计的理论研究j戍用了关于动力机器振动系统的分析方法和步骤，此方法对实际工程应用有重要价值。最后本文讨论了从造纸机的动力计算中发现的一些值得研究的问题，可以作为在应用ANSYS建模以及关于动力机器基础动力设计同类项目的一些参考。本文分5章，具体内容如下：第一章：介绍了土一桩基一机器一上部结构动力相互作用和单桩、群桩基础以及大型造纸机基础动力分析的一些理论和研究现状，并简要说明了本文章的研究意义和研究内容。第二章：详细介绍了桩基础和桩基动力阻抗的理论知识以及动力分析中的模态分析理论和简谐响应分析理论。第三章：根据桩基动力阻抗理论，应用有限元法对桩基础的阻抗进行了计算，用实例验证了动力分析中刚度与阻尼的频率相关性，并和由动规公式计算得到的刚度与阻尼进行了对比。第四章：对大型造纸机某部分振动系统进行动力分析，结果满足造纸机制造商的动力设计要求，提出了桩基设置在该项目中的必要性。第五章：通过实际模型考察了ANSYS动力分析中不同阻尼矩阵取值的结果，讨论了采用空间杆件模型简化三位实体模型的方法和优点，最后就造纸机某特殊部分的动力分析方法进行了讨论。10硕卜学位论文第2章基本理论21桩基础的基本理论211基本概念桩基础是最常用的一种深基础，当地基浅层土质不良，采用浅基础无法满足结构物对地基强度、变形及稳定性方面的要求时，往往需要采用桩基础。承受竖向荷载的桩如图21所示，通过桩侧摩阻力以及桩端阻力将上部荷载传递到深部土(岩)层，而承受横向荷载的桩如图22所示，则由桩身材料和桩侧土(岩)的弹性力来抵抗。桩基础具有承载力高、稳定性好、沉降量小而均匀、抗震性能好等优点。而且桩基础能以不同的桩型和施工方法适应不同的水文地质条件、荷载性质和上部结构特征，因此具有良好的适用性。在港口、桥梁、高耸塔形建筑，海上采油平台，支撑建筑及核电站等大型结构物中，桩基础得到了广泛的应用。桩基础由基桩和联接于桩顶的承台共同组成。若桩身全部埋于土中，承台底面与土体接触，则称为低承台桩基；若桩身上部露出地面而承台底位于地面以上，则称为高承台桩基。建筑桩基通常为低承台桩基础。高层建筑中，桩基础应用广泛。图21竖向荷载下的桩HIl，I，。IIIII土层IIII上图22横向荷载下的桩212单桩动力阻抗定义采用子结构法或动规方法分析土一结构相互作用问题时，桩一土系统等效为弹簧一阻尼器系统，此时需计算桩基的动力阻抗。土一结构动力相互作用问题从大型造纸机幕础动力设计的理论研究i应用早期labm随引对弹性半空间表面上作用集中振动力的位移积分解和reissner阳73所创的弹性半空间刚性圆形基础振动的计算模式，都是以“阻抗的形式描述的。目前，对单桩动力阻抗的计算可以采用数值方法(有限元法、边界元法等)和简化方法。有两种基本的计算模式：一是弹性板空间模式；二是质弹阻计算模式，前者相对于后者来说比较复杂，质弹阻计算模式是简化方法中的一种被广泛采用的模式。桩一土系统动力模型为动力文克尔地基梁模型(beamondynamicwinkler-foundation)，可简写为bdwf模型瞄31。各方面的研究表明桩一土的动力相互作用使桩的刚度和阻尼随着频率而变化。设桩头受到简谐力eoe泐作用，单桩动力阻抗计算就是求桩头力和桩头位移之间的比值，设桩头位移为瞩e栅，则单桩动力阻抗可定义为：艮鲁(21)巩、将式(21)写成复刚度的形式，即K=KI+吗上式中，实部K，代表桩一土系统的动力刚度，记为Kl=K；虚部墨代表桩一土系统的阻尼，可写成如下的形式：如=caC(22)因此，单桩动力阻抗也可写为如下的形式：K=K+i(oC(23)213群桩的阻抗支持上部结构的桩基一般是群桩基础，而不是单桩。在基础设计时不能不考虑的就是群桩效应，群桩效应就是指群桩基础受竖向荷载后，由于承台、桩、土的相互作用，使其桩侧阻力、桩端阻力、沉降等性状发生变化而与单桩明显不同，承载力往往不等于各单桩承载力之和。设群桩的桩数为咒，各桩尺寸和材料相同，一般可以认为各桩受到承台的刚性约束，图23为群桩示意图，图中B为无质量的刚性板，矿是群桩的水平运动位移，形是群桩的竖向运动位移。12硕Jj学位论文I图23群桩刚度计算简图B现简述一种计算群桩阻抗的简化方法，来阐述群桩相互作用原理和阻抗计算。在群桩上的竖向简谐荷载为F，各桩头轴向力为E，承台可看作一个刚性板对各桩起到一个刚性约束的作用，使得各单桩有着相同的位移。因此，群桩的竖向位移等于各单桩的竖向位移彬o，即=彬o(汪1，玎)(24)第i根桩的位移等于该桩在自身荷载的作用下产生的位移加上其余以一1根桩对该桩的作用产生的位移，即彬o=彬+陟+陟爹+陟翟(25)其中，眩表示桩f在自身荷载下的位移；昭表示由于第J根桩的振动(在作用下)引起的第i根桩的位移。引入单桩的动力阻抗和动力相互作用因子来表示群桩的竖向位移，可以得到甩个方程如下：彤=W0去(E+乃+巳)(扣1，刀，V=1，刀)(26)式中，K为单桩的轴向动力阻抗；岛为桩的桩头荷载；(J=1，1)为桩i与其余各桩的动力相互作用因子，且口：=1。群桩的竖向荷载F，由各桩分担，即E=毛(27)则群桩的竖向动力阻抗为：砰2老(2名)同理，可以定义群桩的水平动力阻抗表示如下所示：群=争(29)13大型造纸机基础动力设计的理论研究+j应用22有限元中的阻尼矩阵在所有的振动系统中都具有阻尼，阻尼的主要效应是从系统中消去能量。阻尼系统的能量损失会使振动系统的振动幅值衰减。各种工程结构的阻尼力及其产生的机理是非常复杂的。从宏观看，阻尼有两种主要形态。一种是由结构周围粘性介质产生的阻尼，称为粘性阻尼。粘性阻尼的阻尼力一般近似地认为与运动速度成正比。另一种是由结构材料内部分子间摩擦产生的阻尼，称为结构阻尼或材料阻尼或滞变阻尼。结构阻尼的阻尼力一般近似地认为与弹性体的应变速率成正比。在实际中，要精确地计算各种单元的阻尼矩阵是很困难的，通常是从结构的整体上考虑阻尼的影响，近似地估计阻尼力做功消耗的能量。结构动力分析的有限元方法中，线性振动系统的运动方程为：Mii(t)+Cu(t)+砌(f)=P(t)(210)式中，M为系统的质量矩阵；C为系统的阻尼矩阵；五为系统的刚度矩阵；4、办和五分别为系统的位移向量、速度向量和加速度向量；p为系统的荷载向量。方程(21O)为采用线性粘性阻尼的运动方程的形式。阻尼矩阵C可通过不同的阻尼理论确定，常用的阻尼理论主要有瑞利阻尼理论、复阻尼理论等。阻尼矩阵的确定将在第5章中详细说明。23大型造纸机动力分析方法造纸机机器、基础以及地基(包括桩和土)构成的振动系统，属于典型的土一结构动力相互作用问题。分析这样的振动系统通常有两类方法p831：一类方法是将上部结构及地基全部按有限元离散伸4l。为了模拟地基的半无限特性，尽管可以选取合理的人工边界，但地基的计算区域仍需取得相当大，这样的三维有限元模型自由度太多，实际动力设计中工作量太大，不适于工程设计。另一类为子结构方法。子结构分析方法把系统分成结构和无限地基两部分考虑(这里的结构指广义结构，包括实际结构及其附近有时存在的部分不规则地基部分)，先计算无限地基与结构相接触边界的阻抗和散射特性，然后把这些特性作为结构部分的边界条件，通过二者接触界面上的力和位移连续条件加到结构部分的动力分析中去。结构部分的动力分析可以用有限单元法很容易地完成，然而无限地基边界的阻抗和散射特性的分析却不那么容易，它成为子结构分析方法的关键。国内外学者们围绕这一问题，提出许多解决的办法，使得子结构方法被认为是目前土一结相互作用分析的最有效工具。子结构法在对模型的处理上又可分为两种方法，一种是取上部结构及桩建立有限元模型，承台下及桩间土则模拟为弹簧和阻尼器，这种方法在对桩的处理上似乎更精细，但将土模拟为弹簧和阻尼器并不是一件容14硕上学位论丈易的工作(主要困难是不易确定弹簧的刚度和阻尼参数)，此外也增加了模型的自由度。子结构对模型处理的另一种方法是仅对上部结构建立有限元模型，而将地基(包括桩一土)模拟为弹簧一阻尼器系统。这种模型的自由度最少，同时弹簧一阻尼器的参数可依据动规得到。造纸机动力分析的具体方法和步骤将在第4章中详细说明。24模态分析理论结构的自振特性(固有频率和振型)是结构动力计算中的主要内容。在运动方程(21o)中，令系统的荷载向量p(f)=0，便可得到结构的自由振动运动方程，即刎(f)+翻(f)+砌(f)=0(211)通过求解方程(211)，可以确定结构的固有频率和振型。在工程实际中，结构的阻尼对结构的固有频率和振型的影响不大，因此，可以进一步忽略阻尼的影响，得到结果无阻尼自由振动的运动方程为：刎(f)+胁(f)=0(212)上述方程的解具有如下的形式，即脚(f)=COS国t(213)式中，为位移“(f)的振幅列向量，它与时间f无关，称为振型；彩为与结构的固有频率。将式(213)代入方程(212)，约去cos(or，可得方程(K一国2M)uo=0(214)因此，求解方程(212)就是寻找满足方程(214)的功2值和非零向量，这种问题称为广义特征值问题。记t=国2，则彳和分别称为广义特征值和广义特征向量。方程(214)可写成：(置一tM)uo=0对于上述齐次线性方程组，有非零解的充分必要条件是方程组下的系数行列式等于零，即I置一五Ml=0(215)方程(215)称为结构的特征方程，其中刚度矩阵置和膨都是n阶的方阵，其中n等于自由度的数目，所以方程是关于co：的n次代数方程，由此可解出结构的玎个固有频率。1525简谐响应分析设机器上作用有n组简谐荷载，每组简谐荷载中各简谐力的圆频率和幅值相同，但各组荷载之间的频率互不相同。任一简谐力可写成：P(f)-me02e删(216)对于具有圆频率为国的一组简谐荷载，方程式(210)的右边项可表示为P(f)=删(217)P=mere2I(218)式中，P是简谐荷载的幅值向量；，是单位列向量。由方程式(218)可以看出，当荷载频率在感兴趣的范围内扫描时，幅值向量P随缈变化，但在ANSYS的简谐分析中，P必须为常数向量，为此令P=eo=屁(219)在eo=屁删作用下，设方程式(21O)的解为UiUo-Uoe删(220)将式(219)和式(220)代人式(210)得到眉Vo=I(221)式中，K=IK一彩2M+icoC)(222)求解复数方程式(221)，可得到“。于是，根据叠加原理，系统对该组荷载P=mec02I的简谐响应为U=忱=mec02Uoe泐(223)求得每组荷载的简谐响应后，将这些简谐响应迭加，即可得到系统对n组简谐荷载的总响应。设基础上某量测点对圆频率为q的第k组荷载的响应为=4sin(cod-孽ok)(224)那么，该点对n组简谐荷载的位移和速度响应分别为“(f)=4sin(cod-Ck)(225)k=l1，(f)=q4cos(cokt一婊)(226)k=i一般而言，方程式(225)SFll方程式(226)；14再是简谐运动，因此需根据方程式(220)和方程式(221)画出响应的时间历程曲线。然后，就可由这些曲线确定