贾 天 屈福政 袁任远 代一诚大连理工大学机械工程学院 大连 116024摘 要：基于相似理论分析得到核环行起重机在地震作用下的动力学相似准则，通过畸变影响缩尺模型加速度的相关物理量，使模型能较准确地反映原型在重力场下的动力学特性。根据推导的相似准则，设计了简化的核环行起重机缩尺实验模型，并通过数值模拟软件Abaqus 对模型进行模态和反应谱分析，验证了推导的相似准则的准确性，也为核环行起重机结构的抗震性能分析提供了基础。关键词：环行起重机；畸变相似；抗震性能；缩尺模型；反应谱分析中图分类号：TH113：TH215 文献标识码：A 文章编号：1001-0785（2018）10-0078-050 引言核环行起重机即核电站用环行桥式起重机（以下简称核环吊），是在核电站建造和运行阶段为反应堆厂房内的重型设备安装、维修以及反应堆换料提供吊运服务的特种重型桥式起重机。目前，发展核电是国家能源战略的重要组成部分，地震现象会威胁到核环吊的安全性。核环吊的结构庞大，受试验条件及成本等的限制，对实际结构进行抗震实验存在诸多困难，故采用缩尺模型进行核环吊的抗震性能实验研究。相似理论是近百年来一个新的科学领域，以相似理论为基础的模型实验方法，从相似三定理的建立到现在，不断地发展并日益完善，在机械工程领域得到广泛运用[1]。本文基于相似理论、依据实验设计目标，建立缩尺模型在重力场下的动力学相似准则，并通过有限元分析对基于此准则的缩尺模型及结构原型的动力学相似准则进行验证。1 核环吊缩尺实验的动力学相似理论1.1 相似理论概念相似理论主要包括相似三定理。其中，相似第一定理也称相似正定理，即“若两物理系统相似其相似准则对应的数值相等，且它们的相似指标都等于1。”相似第二定理也叫π 定理，即“当一个物理系统中有n 个物理量，且其中有k 个相互独立的基本量纲，则该系统的n 个物理量可以通过n-k 个相似准则的函数关系表示，即f （π 1，π 2，…，πn）＝ 0。”相似第三定理也叫相似逆定理，即“同一类物理现象中如果单值条件相似，且由单值条件( 主要包括系统的几何条件、边界条件、起始条件和介质条件等) 所推导出的相似准则在数值上相等，则现象相似。”其中，相似第二定理最为常用[2]。将相似定理应用于实际模型的核心是建立正确的相似准则。1.2 核环吊动力学相似准则的建立在相似准则的确定中主要采用以下两种求解方法：量纲分析法[3,4] 和方程分析法[5]。量纲分析法是求解相似准则的主要方法，其利用系统中各物理量的基本量纲建立关于量纲的齐次方程，通过因次分析和矩阵转化可求出相似准则。方程分析法是利用数学模型中的物理方程、微分方程或积分方程来求解相似准则。需要注意的是，在求解特定的动力学相似准则[6] 时，实验目标有时不能做到模型与原型完全动力相似。这时，要进一步确定并区分重要参数及可忽略参数，根据量纲分析法建立模型条件，允许有一定程度的畸变[7] 存在。如果已知原型的运动方程，则模型设计就会相对容易一些。本文结合量纲分析法和方程分析法对模型的相似比进行求解。根据核环吊的缩尺实验目的和实际工况，影响结构动态特性的主要有长度l 、密度ρ 、质量m 、时间t 、力F 、频率f 、加速度g 、弹性模量E 、截面积A 和截面积惯性矩I 等因素。建立相似准则的函数关系式采用质量纲系统[MLT] ( [M]为质量量纲，[L]为长度量纲，[T]为时间量纲) 作为基本量纲，可得到核环吊的量纲关系，α1～ α9分别代表表中各物理的幂次。变换后的指数矩阵如表2 所示。通过π 矩阵可得到相似准则π1 ～π6由量纲分析法得到实验中主要物理量的量纲为由于很难使两种材料的相关参数保持一致，模型采用与原型相同的材料，即N E ＝ 1，Nρ ＝ 1。按照量纲分析得出的相似准则结果，模型的加速度是原型的l N 倍，而原型与模型均在重力场下，为了研究小车的跳轨现象，需要保证原型与模型的加速度相似比为1，即Ng ＝ 1 [8]。因此，将时间相似比畸变为Nt ＝Nl0.5，再通过方程分析法得到畸变后的其他参数的相似比，则梁的弯曲振动的运动方程为式中：EI 为弯曲刚度，u 为垂直方向位移。对表1 中的物理量采用角标注释，其中p 表示原型，m 表示模型。模型和原型之间满足表示以下关系式同时，微分运算符满足以下关系式将上式代入式（9），整理后得到根据方程分析法，模型和原型的基本无量纲参数都一样，也就是模型和原型的方程就是完全相同的，故式（17）需满足以下关系其中NE ＝ 1，Nρ ＝ 1，Ng ＝ 1，代入式（18）中得到由于模型采用了与原型相同的材料，为了满足式（19），要将截面积惯性矩的相似比畸变为NI ＝Nl5。考虑到后续的抗震性能实验要在地震模拟系统（台面尺寸为4 m×3 m）上进行，另有运输、安装等条件的限制，将缩尺模型的几何相似比设为16:1，模型主要参数的相似准则如表3 所示。2 核环吊原型与相似模型有限元分析为了验证本文推导的核环吊抗震分析动力学相似准则的正确性，对依此相似准则设计的缩尺模型及原型实验模型进行有限元仿真分析。原型结构为320/25 t×41m 核环吊，跨度41 m，主梁长41 m，由桥架、小车、大车运行机构、起升机构等组成。2.1 有限元模型考虑结构尺寸及加工方法的限制，主梁由原箱形梁结构简化为矩形实心梁。为了保证弯曲刚度的相似，需要使梁的不同方向的截面惯性矩x y I 及I 均缩尺165 ，从而确定梁截面尺寸。依相似准则将原型简化，缩比模型的各部分参数如表4 所示。其中，大车主梁为M 1，小车承重梁为M 2，后小车架为M 3，吊重的质量为M 4。根据以上参数，在Abaqus 中建立核环吊实际结构和缩尺结构的模型，如图1、图2 所示。图1 Abaqus 原型有限元模型图2 Abaqus 缩尺模型有限元模型2.2 模态分析地震波的频率范围是0.5 ～ 10 Hz，在地震作用下，实际结构的动力学响应中低频占优，地震反应是各阶振型的迭代，而由于高振型在地震响应中的参与系数极小，一般忽略二阶以后的振型，只需保证缩尺模型与原型前2 阶频率和振型相似性[9]。模态分析[10] 是动力学分析的基础，在Abaqus 中对原型和模型进行模态分析，获得各自的固有频率和振型，如图3 所示。从图中可以看出，原型和模型的1 阶均为水平振动，2 阶均为竖直振动。核环吊原型与缩尺模型的前2 阶模态所对应的固有频率如表5 所示，从表中可以看出，缩尺模型的固有频率是原型的4 倍，也就是Nf ＝ 1/4，即Nt ＝ 4，从而验证了原型与缩尺模型的时间及频率的相似关系。(a）原型1 阶振型 （b）原型2 阶振型 （c）模型1 阶振型（d）模型2 阶振型图3 原型与模型前2 阶振型2.3 反应谱分析为了比较原型和模型在整个地震响应过程中的加速度及位移最大值，采用反应谱分析法[11] 对加速度及位移的相似性进行验证。反应谱分析是建立在振型分解反应谱理论[12] 基础上的抗震分析方法。基于振型分解理论将结构的地震响应分解为各振型分量的叠加，即对应每个振型都有一个地震作用，然后通过一定的组合方法（如SRSS、CQC、ABS 等）叠加各振型结构的地震响应得到总结构地震响应值。地震输入波的特性会影响结构的动力学特征，根据核电厂抗震设计规范[13] 中的规定，地震波选取Ⅱ类场地、用于中周期结构输入的EL Centro 1979 波。根据实验需要对地震波进行调幅处理，水平地震波幅值为0.3g ，竖向地震波幅值为0.2g ，其频谱如图4 所示。根据加速度和频率的相似准则，将频谱轴扩大4 倍，纵轴保持不变，得到图4、图5 所示最终的输入波的水平和竖直频谱。根据核电站抗震规范，依据原型结构的材料、工艺以及所处的场地类型，阻尼比选取4%。经过反应谱分析，核环吊原型和模型的加速度峰值及位移最大值均出现在主梁跨中，具体数值如表6、表7 所示。图4 原型的反应谱图5 模型的反应谱由反应谱分析结果可知，原型位移与模型位移满足16 :1 的关系。同时，原型加速度峰值与缩尺模型加速度峰值满足Ng ＝ 1，验证了加速度及位移的相似关系，说明缩尺模型在重力场下可较准确地反映原型的动力学特征。3 结论本文运用相似理论，建立了地震作用下的核环吊结构动力学相似准则，根据实验的相关限制条件和实验目的确定了影响其动力学特征的重要参数的相似比，进而得到了核环吊的缩尺模型，并通过数值仿真对比分析了模型与原型结构的动力学特征，得到以下结论：1）文中所推导的相似准则可以为今后的核环吊缩尺实验提供理论依据。同时，推导出的相似准则也可以推广到一般大型复杂的结构的振动实验模型设计中，在工程应用中具有实际应用价值。2）通过对模型各方向的截面惯性矩的缩尺比进行畸变，使原型和模型的加速度相同的方法是可靠和准确的。3）对比核环吊原型和缩尺模型的前2 阶振型、固有频率和加速度峰值可以发现，基本符合本文提出的相似准则，可以比较准确地反应原型的动态特性。参考文献[1] 胡冬奎，王平. 相似理论及其在机械工程中的应用[J].现代制造工程，2009(11)：9-12.[2] 小濱英司，安藤和彰，菅野高弘。等. 低頭型の大規模地震時挙動に関する模型振動実験及び数値解析，No.1289[R]. 横须贺市：港湾技術研究所，2014.[3] 郑哲敏. 相似理论与模化[C]. 郑哲敏文集，中国力学学会，2004.[4] 谈庆明. 量纲分析[M]. 合肥：中国科学技术大学出版社，2005.[5] DJ Schuring.Scale models in engineering fundamental and applications[M].Oxford：Pergamon Press，1977.[6] 罗忠，朱云鹏，韩清凯，等. 动力学相似理论及在结构振动分析中的应用研究评述与展望[J]. 机械工程学报，2016，52(23)：114-134.[7] Laura Diane Jacobs.Shake table experiments for the determination of the seismic response of jumbo container cranes[D]. Atlanta：Georgia Institute of Technology，2010.[8] A Otani，K Nagashima，J Suzuki.Vertical seismic response of overhead crane[J].Nuclear Engineering and Design，2002，212(1).[9] 孙广俊，李爱群. 关于结构振型参与系数和振型贡献的分析[J]. 防灾减灾工程学报，2009，29(5)：485-490.[10] 郑明刚，刘天雄，陈兆能. 基于频响函数的结构损伤检测[J]. 机械科学与技术，2001(3)：458-480，319.[11] 吴小峰，孙启国，狄杰建，等. 抗震分析反应谱法和时程分析法数值仿真比较[J]. 西北地震学报，2011，33(3)：275-278，304.[12] 李正英，李正良，范文亮. 耗能减震结构的振型分解法[J]. 重庆大学学报( 自然科学版)，2005(11)：133-137.[13] GB 50267—1997 核电厂抗震设计规范[S].举报/反馈发表评论发表作者最新文章高定位精度的第三代核环行起重机运行系统 刚柔耦合动力学仿真01-2014:40液压挖掘机转台有限元分析与疲劳强度评估01-2014:33SPMT 液压平板车车板变形有限元计算与仿真模拟01-2014:31相关文章将用5年时间打造城市副中心网红打卡地，大运河5A级景区明年展露雏形耶伦上台，弱美元趋势的终结？油价还能涨得动吗？热搜爆了！这个实验室藏不住了半月谈点赞淄博：营商土沃
本文来源于程序自动整理，版权争议与本站无关，如有侵权，请联系我们。
冀州区扬禹水工机械厂
厂家直销，支持定制，量大优惠
手机：13663289587
地址：河北省冀州区高田工业区
扫一扫加微信
