转子动力学计算分析与优化

Ansys转子动力学

基于ANSYS的转子动力学分析 1、题目描述 如图1-1所示,利用有限原原理计算转子临界转速以及不平衡响应。 图 1-1 转子示意图及尺寸 2、题目分析 采用商业软件ANSYS进行分析,转子建模时用beam188三维梁单元,该单元基于Timoshenko梁理论,考虑转动惯量与剪切变形的影响。每个节点有6个(三个平动,三个转动)或7各自由度(第七个自由度为翘曲,可选)。 轴承用combine214单元模拟。该单元可以模拟交叉刚度和阻尼。只能模拟拉压刚度,不能模拟弯曲或扭转刚度。该单元如图2-1所示,其有两个节点组成,一个节点在转子上,另一个节点在基础上。

图 2-1 combine214单元 对于质量圆盘,可以用mass21单元模拟,该单元有6个自由度,可以模拟X,Y,Z 三个方向的平动质量以及转动惯性。 3、计算与结果分析 3.1 转子有限元模型 建模时,采用钢的参数,密度取37800/kg m ,弹性模量取112.1110pa ,泊松比取0.3。轴承刚度与阻尼如表1所示,不考虑交叉刚度与阻尼,且为各项同性。 Kxx Kyy Cxx Cyy 4e7N/m 4e7N/m 4e5N.s/m 4e5N.s/m 将转子划分为93个节点共92个单元。有限元模型如图3-1所示。

图 3-1 转子有限元模型 施加约束时,由于不考虑纵向振动与扭转振动,故约束每一节点的纵向与扭转自由度,同时约束轴承的基础节点。施加约束后的模型如3-2所示。 图 3-2 施加约束后的有限元模型 3.1 转子临界转速计算 在ANSYS中可以很方便的考虑陀螺力矩的影响。考虑陀螺力矩时,由于陀螺矩阵是反对称矩阵,所以求取特征值时要用特殊的方法。本文考虑陀螺力矩的影响,分析了在陀螺力矩的影响下,转子涡动频率随工作转速的变化趋势,其Campell图如图3-3所示。同时给出了转子的前四阶正进动涡动频率与反进动涡动频率以及固有频率。如表3-2所示。

哈工大结构动力学大作业2012春

结构动力学大作业 对于如下结构,是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。 1 以上为一个简支梁结构。集中质量块放于梁上,质量块距简支梁的左端点距离为L. 将该简支梁简化为欧拉伯努利梁,并离散为N 个单元。每个单元有两个节点,四个自由度。 单元的节点位移可表示为: ]1122,,,e v v δθθ?=? 则单元内一点的挠度可计作: 带入边界条件: 1 3 32210)(x a x a x a a x v +++=0 1)0(a v x v ===3 322102)(L a L a L a a v L x v +++===1 10 d d a x v x ===θ2 321232d d L a L a a x v L x ++===θ1 0v a =

[]12 3 4N N N N N = 建立了单元位移模式后,其动能势能均可用节点位移表示。单元的动能为: 00111()222 l l T T T ke e e e e y E dx q N Ndxq q mq t ρρ?===??? 其中m 为单元质量阵,并有: l T m N Ndx ρ=? 带入公式后积分可得: 222215622541322413354 1315622420133224l l l l l l l m l l l l l l ρ-?? ??-??= ?? -?? ---? ? 单元势能可表示为 22 200 11()()22 2 T l l T T e pe e e e q y E EI dx EI N N dxq q Kq x ?''''== =??? 其中K 为单元刚度矩阵,并有 ()l T K EI N N dx ''''=? 2 23 2212 612664621261266264l l l l l l EI k l l l l l l l -????-??=??---??-?? 以上为单元类型矩阵,通过定义全局位移矩阵,可以得到系统刚度矩阵和系统质量矩 1 1θ=a )2(1)(3211222θθ+--=L v v L a )(1)(22122133θθ++-= L v v L a 1232133222231)(θ???? ??+-+???? ??+-=L x L x x v L x L x x v 2 2232332223θ??? ? ??-+???? ??-+L x L x v L x L x 2 4231211)()()()()(θθx N v x N x N v x N x v +++=

转子动力学知识

2转子动力学主要研究那些问题 答:转子动力学是研究所有不旋转机械转子及其部件和结构有关的动力学特性,包括动态响应、振动、强度、疲劳、稳定性、可靠性、状态监测、故障诊断和控制的学科。这门学科研究的主要范围包括:转子系统的动力学建模与分析计算方法;转子系统的临界转速、振型不平衡响应;支承转子的各类轴承的动力学特性;转子系统的稳定性分析;转子平衡技术;转子系统的故障机理、动态特性、监测方法和诊断技术;密封动力学;转子系统的非线性振动、分叉与混沌;转子系统的电磁激励与机电耦联振动;转子系统动态响应测试与分析技术;转子系统振动与稳定性控制技术;转子系统的线性与非线性设计技术与方法。 3转子动力学发展过程中的主要转折是什么 答:第一篇有记载的有关转子动力学的文章是1869年Rankine发表的题为“论旋转轴的离心力”一文,这篇文章得出的“转轴只能在一阶临界转速以下稳定运转”的结论使转子的转速一直限制在一阶临界以下。最简单的转子模型是由一根两端刚支的无质量的轴和在其中部的圆盘组成的,这一今天仍在使用的被称作Jeffcott转子的模型最早是由Foppl在1895年提出的,之所以被称作“Jeffcott”转子是由于Jeffcott教授在1919年首先解释了这一模型的转子动力学特性。他指出在超临界运行时,转子会产生自动定心现象,因而可以稳定工作。这一结论使得旋转机械的功率和使用范围大大提高了,许多工作转速超过临界的涡轮机、压缩机和泵等对工业革命起了很大的作用。但是随之而来的一系列事故使人们发现转子在超临界运行达到某一转速时会出现强烈的自激振动并造成失稳。这种不稳定现象首先被Newkirk发现是油膜轴承造成的,仍而确定了稳定性在转子动力学分析中的重要地位。有关油膜轴承稳定性的两篇重要的总结是由Newkirk和Lund写出的,他们两人也是转子动力学研究的里程碑人物。 4石化企业主要有哪些旋转机械,其基本工作原理是什么 汽轮机:将蒸汽的热能转换成机械能的涡轮式机械。工作原理:在汽轮机中,蒸汽在喷嘴中发生膨胀,压力降低,速度增加,热能转变为动能。作用与功能:主要用作发电用的原动机,也可直接驱动各种泵、风机、压缩机和船舶螺旋桨等。还可以利用汽轮机的排汽或中间抽汽满足生产和生活的供热需要。 燃气轮机:是一种以空气及燃气为介质,靠连续燃烧做功的旋转式热力发动机。主要结构由三部分:压气机,燃烧室,透平(动力涡轮)。作用与功能:以

转子动力学

转子动力学是固体力学的一个分支。本文主要研究转子支承系统在旋转状态下的振动,平衡和稳定性,特别是在接近或超过临界转速的情况下转子的横向振动。转子是涡轮机,电动机和其他旋转机械的主要旋转部件。 200多年来,工程和科学界一直关注转子振动。w.j.m. 1869年英格兰的兰金(Rankin)和1889年法国的拉瓦尔(c.g.p.de Laval)对挠性轴的测试是研究此问题的先驱。随着现代工业的发展,高速细长转子逐渐出现。由于它们通常在柔性状态下工作,因此它们的振动和稳定性变得越来越重要。转子动力学的主要研究内容如下: ①临界速度 由于制造误差,转子每个微小部分的质心与旋转轴略有偏离。当转子旋转时,由上述偏差引起的离心力将使转子产生横向振动。在某些速度(称为临界速度)下,这种振动似乎非常强烈。为了确保机器不会在工作速度范围内产生共振,临界速度应适当偏离工作速度,例如大于10%。临界速度与转子的弹性和质量分布有关。对于具有有限集总质量的离散旋转系统,临界速度的数量等于集总质量的数量;对于具有连续质量分布的弹性旋转系统,临界速

度是无限的。传递矩阵法是计算大型转子支撑系统临界转速的最常用数值方法。要点是:首先,将转子分成几个部分,每个部分左右两端的四个部分参数(挠度,挠度角,弯矩和剪切力)之间的关系可以通过传递来描述。该部分的矩阵。以此方式,可以获得系统的左端和右端的横截面参数之间的总传递矩阵。然后,根据边界条件和自然振动中非零解的条件,通过试错法求出各阶的临界速度,得到相应的振动模式。 ②通过临界速度的状态 通常,转子以可变速度通过临界速度,因此通过临界速度的状态是不稳定的。与以临界速度旋转时的静止状态不同,有两个方面:一是振幅的最大值小于静止状态的振幅,速度越大,振幅的最大值越小。另一个是振幅的最大值不会在像静止状态那样的临界速度下出现。在不稳定状态下,频率转换干扰力作用在转子上,这使分析变得困难。为了解决这种问题,在数值计算或非线性振动理论中必须使用渐近法或级数展开法。 ③动态响应

结构动力学大作业

结构动力学作业 姓名: 学号:

目录 1.力插值法 (1) 1.1分段常数插值法 (1) 1.2分段线性插值法 (4) 2.加速度插值法 (7) 2.1常加速度法 (7) 2.2线加速度法 (9) 附录 (12) 分段常数插值法源程序 (12) 分段线性插值法源程序 (12) 常加速度法源程序 (13) 线加速度法源程序 (13)

1.力插值法 力插值法对结构的外荷载进行插值,分为分段常数插值法和分段线性插值法,这两种方法均适用于线性结构的动力反应计算。 1.1分段常数插值法 图1-1为一个单自由度无阻尼系统,结构的刚度为k ,质量为m ,位移为y (t ),施加的外力为P (t )。图1-2为矩形脉冲荷载的示意图,图中t d 表示作用的时间,P 0表示脉冲荷载的大小。 图1-1 单自由度无阻尼系统示意图 图1-2 矩形脉冲荷载示意图 对于一个满足静止初始条件的无阻尼单自由度体系来说,当施加一个t d 时间的矩形脉冲荷载,此时结构在t d 时间内的位移反应可以用杜哈梅积分得到: 0()sin ()2 (1cos )(1cos ) (0) t st st d P y t t d m t y t y t t T ωττω πω=-=-=-≤≤? (1-1) 如果结构本身有初始的位移和速度,那么叠加上结构自由振动的部分,结构的位移反应为: 02()cos sin (1cos ) (0 )st d y t y t y t t y t t T πωωω =+ +-≤≤ (1-2)

图1-3 分段常数插值法微段示意图 对于施加于结构任意大小的力,将其划分为Δt 的微段,每一段的荷载都为一个常数(每段相当于一个矩形的脉冲荷载),如图1-3所示,则将每一段的位移和速度写成增量的形式为: 1cos t sin t (1cos t)i i i i y P y y k ωωωω +=?+ ?+-? (1-3) i+1/sin t cos t sin t i i i y P y y k ωωωωω =-?+ ?+ ? (1-4) 程序流程图如下

高等结构动力学大作业

Advanced Structural Dynamics Project The dynamic response and stability analysis of the beam under vertical excitation Instructor:Dr. Li Wei Name: Student ID:

1.Problem description and thepurpose of the project 1.1 calculation model An Eular beam subjected to an axial force. Please build thedifferential equation of motion and use a proper difference method to solve this differentialequation. Study the dynamic stability of the beam related to the frequency and amplitude of the force. As shown in the Fig 1.1. Fig1.1 1.2 purpose and process arrangement a.learninghow to create mathematical model of thecontinuous system and select proper calculation method to solve it. b.learning how to build beam vibration equation and solve Mathieu equation. https://www.360docs.net/doc/b33992000.html,ing Floquet theory to judgevibration system’s stability and analyze the relationship among the frequency and amplitude of the force and dynamic response. This project will introduce the establishment of the mathematical model of the continuous system in section 2, the movement equation and the numerical solution of using MATLAB in section 3,Applying Floquent theory to study the dynamic stability of the beam related to the frequency and amplitude of the force in section 4. In the last of the project, we get some conclusions in section 5.

转子动力学知识

转子动力学知识 2转子动力学主要研究那些问题? 答:转子动力学是研究所有不旋转机械转子及其部件和结构有关的动力学特性,包括动态响应、振动、强度、疲劳、稳定性、可靠性、状态监测、故障诊断和控制的学科。这门学科研究的主要范围包括:转子系统的动力学建模与分析计算方法;转子系统的临界转速、振型不平衡响应;支承转子的各类轴承的动力学特性;转子系统的稳定性分析;转子平衡技术;转子系统的故障机理、动态特性、监测方法和诊断技术;密封动力学;转子系统的非线性振动、分叉与混沌;转子系统的电磁激励与机电耦联振动;转子系统动态响应测试与分析技术;转子系统振动与稳定性控制技术;转子系统的线性与非线性设计技术与方法。 3转子动力学发展过程中的主要转折是什么? 答:第一篇有记载的有关转子动力学的文章是1869年Rankine发表的题为“论旋转轴的离心力”一文,这篇文章得出的“转轴只能在一阶临界转速以下稳定运转”的结论使转子的转速一直限制在一阶临界以下。最简单的转子模型是由一根两端刚支的无质量的轴和在其中部的圆盘组成的,这一今天仍在使用的被称作Jeffcott转子的模型最早是由Foppl在1895年提出的,之所以被称作“Jeffcott”转子是由于Jeffcott教授在1919年首先解释了这一模型的转子动力学特性。他指出在超临界运行时,转子会产生自动定心现象,因而可以稳定工作。这一结论使得旋转机械的功率和使用范围大大提高了,许多工作转速超过临界的涡轮机、压缩机和泵等对工业革命起了很大的作用。但是随之而来的一系列事故使人们发现转子在超临界运行达到某一转速时会出现强烈的自激振动并造成失稳。这种不稳定现象首先被Newkirk发现是油膜轴承造成的,仍而确定了稳定性在转子动力学分析中的重要地位。有关油膜轴承稳定性的两篇重要的总结是由Newkirk和Lund写出的,他们两人也是转子动力学研究的里程碑人物。 4石化企业主要有哪些旋转机械,其基本工作原理是什么? 汽轮机:将蒸汽的热能转换成机械能的涡轮式机械。工作原理:在汽轮机中,蒸汽在喷嘴中发生膨胀,压力降低,速度增加,热能转变为动能。作用与功能:主要用作发电用的原动机,也可直接驱动各种泵、风机、压缩机和船舶螺旋桨等。还可以利用汽轮机的排汽或中间抽汽满足生产和生活的供热需要。

转子动力学

固体力学的分支。本文主要研究转子轴承系统在旋转状态下的振动,平衡和稳定性,特别是在接近或超过临界速度的运行状态下转子的横向振动。转子是旋转机械(例如涡轮机和电动机)中的主要旋转部件。 工程和科学界一直关注转子振动已有200多年的历史了。1869年英国W.J.M Rankin撰写的有关离心力的论文以及法国C.G.P.de Laval于1889年对挠性轴进行的测试是研究此问题的先驱者。随着现代工业的发展,高速细长转子逐渐出现。由于它们通常在柔性状态下工作,因此振动和稳定性问题变得越来越重要。转子动力学的主要研究内容如下: 由于制造误差,转子的每个微段的质心通常会略微偏离旋转轴。当转子旋转时,由上述偏差引起的离心力将导致转子横向振动。在某些转速(称为临界转速)下,这种振动非常强烈。为了确保机器不会在工作速度范围内产生共振,临界速度应偏离工作速度超过10%。临界速度与转子的弹性和质量分布有关。对于具有有限集中质量的离散旋转系统,临界速度的数量等于集中质量的数量。 对于具有连续质量分布的弹性旋转系统,存在无限的临界速度。用于计算大型转子支撑系统的临界转速的最常用数值方法是传递矩阵法。要点如下:首先,将转子分成几个部分,每个部分左右两端的四个部分参数(挠度,挠度角,弯矩和剪切力)之间的关系可以用下式描述:

本节的转移矩阵。以此方式,可以获得系统的左端面和右端面的截面参数之间的总传递矩阵。然后,根据边界条件和自然振动中存在非零解的条件,通过试错法求出各阶的临界速度,然后得到相应的振动模式。 由于Jeffcott转子的特殊性,唯一的轮盘位于两个刚性支撑之间,因此可以忽略陀螺力矩对临界转速的影响。Jeffcott转子在无阻尼状态下的临界速度可以看作是其固有频率,但是对于其他类型的转子,陀螺力矩对临界速度的影响是不能忽略的,这是与结构动力学的差异之一。和振动力学。 就转子动力学而言,在存在外部阻尼的情况下,Jeffcott转子的临界速度高于其在非阻尼状态下的固有频率,该结论也适用于其他类型的转子。 对于结构动力学和振动力学,具有外部阻尼的结构系统的固有频率低于没有阻尼的结构系统的固有频率。

ANSYS中的转子动力学计算

ANSYS 中的转子动力学计算 安世亚太 许明财 转子动力学是固体力学的一个重要分支,它主要研究旋转机械的“转子-支承”系统在旋转状态下的振动、平衡和稳定性问题,其主要研究内容有几个方面:临界转速、动力响应、稳定性、动平衡技术和支承设计。在旋转机械研究设计中,转子动力学的性能分析是极其重要的一个方面。 旋转机械广泛应用于以下领域: y 涡轮机械 y 能量站 y 机械工具 y 汽车 y 家用器械 y 航空领域 y 船舶推进系统 y 医疗器械 y 计算机设备 传统的转子动力学分析采用传递矩阵方法进行,由于将大量的结构信息简化为极为简单的集中质量—梁模型,不能确保模型的完整性和分析的准确度;而有限元在处理转子动力学问题时,可以很好地兼顾模型的完整性和计算的效率,但多年来转子的“陀螺效应”一直是 制约转子动力学有限元分析的“瓶颈”问题。 ANSYS 很好地解决了动力特性分析中“陀螺效应”影响的问题,而且陀螺效应的考虑不受计算模型上的限制,使得转子动力学有限元分析变得简单高效。 本文对ANSYS 的转子动力学计算功能进行简要介绍。 1 ANSYS 转子动力学的理论基础 ANSYS 转子动力学分析中,两种参考坐标系可供选择:静止坐标系和旋转坐标系。空间点P 在静止坐标系(其原点在O′)下的位置矢量为r′,在旋转坐标系(其原点在O)下的位置矢量为r。 在静止坐标系下转子的动力方程为: [][][]{}F {u}K }u ]){gyr [C C (}u {M =+++&&& 式中:为陀螺效应矩阵 ]gyr [C

在旋转坐标系下转子的动力方程为: [][][]{}F }r ]){u spin [K K (}r u ]){cor [C C (}r u {M =?+++&&& 式中:为哥氏效应矩阵, 为旋转软化效应刚度矩阵 ]cor [C ]spin [K 2 ANSYS 转子动力学的计算功能和新技术 ANSYS 转子动力学计算包含如下功能: y 无阻尼临界转速分析 y 不平衡响应分析 y 阻尼特征值分析 y 涡动和稳定性预测 典型的应用包括: y 轴的弯曲变形 y 扭转振动 y 转子轴未对准 y 旋转部分的平衡 y 流制振动 为了分析时计入哥氏效应、陀螺效应和支承的影响,ANSYS 发展了下列新技术单元: SHELL181 4节点有限应变壳单元 PLANE182 二维4节点结构实体单元 PLANE183 二维8节点结构实体单元 SOLID185 三维8节点结构实体单元 SOLID186 三维20节点结构实体单元 SOLID187 三维10节点四面体结构实体单元 BEAM188 三维一次有限应变梁单元 BEAM189 三维二次有限应变梁单元 SOLSH190 三维8节点层合实体-壳单元 COMBIN214 二维轴承单元(可变刚度和阻尼) ANSYS 考虑陀螺效应时没有计算模型上的限制,故可选择一维(梁、管)、二维(轴对称)和三维复杂计算模型进行分析。同时,ANSYS 还提供了一系列功能以完善转子动力学的计算,包括: y CMOMEGA 可以通过组(component CM_NAME)对多个转子指定不同的转速 y CORIOLIS 可以考虑哥氏效应在不同参考坐标系下的影响 y PLCAMP 可绘制坎贝尔图,为临界转速的确定提供了方便 y PRCAMP 可打印固有频率和临界转速 y CAMPB 可绘制预应力结构的坎贝尔图。 y PRORB 可打印出转子涡动幅值 y PLORB 可绘制转子不同截面的涡动轨迹 y SYNCHRO 可以指定激励频率是否与结构速度同步 y 多轴转子 y 正向/反向旋转和不稳定性 3 实例 实例1:梁单元与实体单元比较

结构力学大作业

结构力学大作业——五层三跨框架结构内力计算 专业班级:土木工程XXXX班 姓名 XXXXX 学号:XXXXX 指导教师:XX

目录 一、题目 (3) 二、任务 (5) 三、结构的基本数据 (5) 1.构件尺寸: (5) 2.荷载: (5) 3.材料性质: (5) 四、水平荷载作用下的计算 (5) 1.反弯点法 (6) 2.D值法 (8) 3.求解器法 (12) 五、竖直荷载作用下的计算 (15) 1.分层法 (16) 2.求解器法 (21) 六、感想 (24)

二、题目 结构(一) 1、计算简图如图1所示。 4 . 2 m 3 . 6 m 3 . 6 m 3 . 6 m 3 . 6 m 图1

’ 图2 q’ 图3

二、任务 1、计算多层多跨框架结构在荷载作用下的内力,画出内力图。 2、计算方法: (1) 水平荷载: D 值法、反弯点法、求解器,计算水平荷载作用下的框架 弯矩; (2) 竖向荷载:迭代法、分层法、求解器,计算竖向荷载作用下框架弯矩。 3、对各种方法的计算结果进行对比,分析近似法的误差。 4、把计算过程写成计算书的形式。 三、结构的基本数据 E h =3.0×107kN/m 2 柱尺寸:400×400,梁尺寸(边梁):250×600,(中间梁)300×400 竖向荷载:q '=17kN/m 水平荷载:F P '=15kN 构件线刚度:)12 (,3 bh I l EI i == 柱子:43-3 10133.212 400400m I ?=?= 柱 第一层:m kN i ?=???= -152382.410133.2100.33 71 第二--五层:m kN i ?=???= -177786.310133.2100.33 72 梁: 边梁:43-3105.412 600250m I ?=?=边梁 m kN i ?=???=-225006105.4100.3373 中间梁:43-3106.112 400300m I ?=?=中间梁 m kN i ?=???=-228571 .2106.1100.3374 四、水平荷载作用下的计算 水平荷载: F P =16kN ,F p '=15kN

2018西南大学[0729]《结构力学》大作业答案

1、结构的刚度是指 1. C. 结构抵抗变形的能力 2、 图7中图A~图所示结构均可作为图7(a)所示结构的力法基本结构,使得力法计算最为简便的 C 3、图5所示梁受外力偶作用,其正确的弯矩图形状应为()C 4、对结构进行强度计算的目的,是为了保证结构 1. A. 既经济又安全 5、改变荷载值的大小,三铰拱的合理拱轴线不变。 1. A.√ 6、多余约束是体系中不需要的约束。 1. B.×

7、结构发生了变形必然会引起位移,结构有位移必然有变形发生。 1. B.× 8、如果梁的截面刚度是截面位置的函数,则它的位移不能用图乘法计算。 1. A.√ 9、一根连杆相当于一个约束。 1. A.√ 10、单铰是联接两个刚片的铰。 1. A.√ 11、虚功原理中的力状态和位移状态都是虚设的。 1. B.× 12、带拉杆三铰拱中拉杆的拉力等于无拉杆三铰拱的水平推力。 1. A.√ 13、瞬变体系在很小的荷载作用下会产生很大的内力,所以不能作为结构使用。 1. A.√ 14、虚位移原理中的虚功方程等价于静力平衡方程,虚力原理中虚功方程等价于变形协调方程。 1. A.√ 15、体系的多余约束对体系的计算自由度、自由度及受力状态都没有影响,故称多余约束。 1. B.× 16、力矩分配中的传递系数等于传递弯矩与分配弯矩之比,它与外因无关。 1. A.√ 17、当上部体系只用不交于一点也不全平行的三根链杆与大地相连时,只需分析上部体系的几何组成,就能确1. A.√ 18、用力法计算超静定结构时,其基本未知量是未知结点位移。

B.× 19、静定结构在非荷载外因(支座移动、温度改变、制造误差)作用下,不产生内力,但产生位移。 1. A.√ 20、力法和位移法既能用于求超静定结构的内力,又能用于求静定结构的内力。() 1. B.× 21、静定结构在非荷载外因(支座移动、温度改变、制造误差)作用下,不产生内力,但产生位移。()1. A.√ 22、位移法和力矩分配法只能用于求超静定结构的内力,不能用于求静定结构的内力。( ) 1. B.× 23、 图2所示体系是一个静定结构。() 1. B.× 24、力矩分配法中的分配系数、传递系数与外来因素(荷载、温度变化等)有关。 1. B.× 25、三铰拱的水平推力不仅与三铰的位置有关,还与拱轴线的形状有关。 1. B.× 26、三铰拱的主要受力特点是:在竖向荷载作用下产生水平反力。 1. A.√ 27、两根链杆的约束作用相当于一个单铰。 B.× 28、不能用图乘法求三铰拱的位移。

结构动力学大作业

目录 一、结构特性矩阵 1.1框架设计 (2) 1.2截面尺寸 (2) 1.3动力自由度 (2) 1.4结构的一致质量矩阵 (3) 1.5结构的一致刚度矩阵 (13) 二、频率与振型 2.1简化的质量矩阵 (25) 2.2简化的刚度矩阵 (25) 2.3行列式法求频率与振型 (27) 2.4Stodola法求频率与振型 (27) 三、时程分析 3.1框架资料 (31) 3.2地震波波形图 (31) 3.2瑞利阻尼 (32) 3.4操作步骤 (33) 3.5各楼层位移时程反应图 (37)

一、结构特性矩阵 1.1框架设计 框架平面图如图1所示,跨度均为6.0m,层高均为3.6m,混凝土采用C30。 图1 框架平面图 1.2截面尺寸 梁均为300mm600mm ? ?,柱均为500mm500mm 1.3动力自由度 框架结构可以理想化为在节点处相互连接的梁柱单元的集合。设梁、柱的轴向变形均忽略不计,只考虑横向平面位移,则该框架有3个平动自由度和12个角自由度,共15个自由度,并对梁柱单元分别编号,如图2所示: 图2 单元编号及自由度

将结构分成在有限个节点处相互连接的○1~○21个离散单元体系,通过计算各个单元的一致质量矩阵、一致刚度矩阵,并将相关的单元叠加求得整个单元结构的一致质量矩阵、一致刚度矩阵。 1.4结构的一致质量矩阵 在节点位移作用下框架梁和柱上所引起的变形形状采用三次Hermite 多项式,因此均布质量梁的一致质量矩阵为: ??? ???? ???????4 3 2 1 I I I I f f f f =420L m ?? ? ?? ???????------222 2432213341322221315654132254156 L L L L L L L L L L L L ???? ????????? ????? (4) .. 3 2 1 v v v v 梁:m =250060.030.0??=450kg/m, L=6m;

结构动力学大作业

结构动力学大作业 问题描述 《建筑结构抗震设计》高振世P247 该建筑为一幢六层现浇钢筋混泥土框架房屋,屋顶有局部突出的楼梯间和水箱间。混泥土强度等级:梁为C20,柱为C25。混凝土密度为2500kg/m3 本题目将对该梁柱结构的框架房屋进行模态分析,求解出该结构的前8阶固有频率及其对应的模态振型。框架的平、剖面见图1,图2。构件尺寸参见表1、表2。 其材料为混凝土,相关参数为:杨氏模量C20为2.55e10N/m2,C25为2.8e10 N/m2。 图 1 平面视图 图 2 剖面图 表 1 梁的几何尺寸 部位断面 b×h (m×m) 跨度L (m) 屋顶梁0.25×0.60 5.7 楼层梁0.25×0.65 5.7 走道凉0.25×0.40 2.1 表 2 柱的几何尺寸 层次柱高(m)断面(m×m) 1 4 0.50×0.50 2,3,4,5,6 3.6 0.45×0.45

建模 利用有限元商业软件ANSYS8.0建模和有限元分析。 1单元类型(Element Type ): Beam4 图 3 单元BEAM4的特征 2材料模型: 本题目中所有材料都假定是各向同性线弹性体。 表 3 MARERIAL MODEL No. 梁/柱混凝土标号弹性模量EX N/m2 泊松比PXY 密度DENS kg/ m3 1 柱C25 2.8e10 0. 2 2500 2 梁C20 2.55e10 0.2 2500

3单元实常数: 表 4 REAL CONSTANT No. 部位 横截面面积 AREA(m2) Z轴惯性矩 IZZ(m4) Y轴惯性矩 IYY(m4) 高度 TKZ(m) 宽度 TKY(m) 1 底层柱0.25 5.208e-3 5.208e-3 0.5 0.5 2 其余柱0.2025 3.417e- 3 3.417e-3 0.45 0.45 3 走道梁0.1 1.333e-3 5.208e- 4 0.4 0.25 4 楼层梁0.162 5 5.721e-3 8.464e-4 0.65 0.25 5 屋顶梁0.15 0.45e-3 7.8125e-4 0. 6 0.25 4几何模型,网格划分,施加边界条件: 柱划为个6单元,走道梁划为3个单元,楼层梁划为5个单元,合计3029个单元 图 4

非线性转子 动力学

航空发动机非线性转子碰磨研究 XXX (XXXX 机械工程上海200072) 摘要:综述了国内外非线性转子动力学的研究现状,讨论了非线性转子动力学研究中的7个主要问题,并引述了大量相应的国内外文献,包括:非线性转子动力学研究的一般方法;求解非线性转子动力学问题的数值积分方法;大型转子-轴承系统高维非线性动力学问题的降维求解;基于微分流形的动力系统理论方法;转子非线性动力学行为的机理研究和实验研究;高速转子-轴承系统的非线性动力学设计,最后讨论了非线性转子动力学研究中存在的问题及展望。 关键词:非线性;高速转子;数值积分法 The research for Aeroengine nonlinear rotor WANG Qing-long (Shanghai university mechainal engineering 20072 shanghai) Abstract: Reviewed the research status of nonlinear rotor dynamics both at home and abroad, discusses the seven main in the study of nonlinear rotor dynamics. To questions, and cited a large number of relevant literature both at home and abroad, include: common methods of nonlinear rotor dynamics; To solve the non-linear. Rotor dynamics problems of numerical integral method; Rotor - bearing system of large dimension reduction solution for high dimensional nonlinear dynamics; In the theory of differential dynamic system of the manifold method; Rotor nonlinear dynamics behavior of mechanism research and experiment research; High speed rotor shaft. Bearing system of the nonlinear dynamics design, and finally discusses the problems of nonlinear rotor dynamics research and prospects. Key words: nonlinear; High speed rotor; The numerical integral method. 由于旋转机械系统中各种异常振动的存在,常常引发灾难性的事故。过去研究转子-轴承-基础系统大多采用基于线性转子动力学理论。例如传统转子动力学对转子-轴承系统稳定性问题的研究,一般采用8个线性化的刚度与阻尼特性系数的油膜力模型。对于大型旋转机械中存在的油膜力、密封力、不均匀蒸汽间隙力等严重的非线性激励源,由于数学模型不够完善,以致系统中存在的许多由非线性因素引起的多种复杂动力学行为尚没有彻底搞清,不能满足现代工程设计的需要,迫切需要建立转子-轴承系统的非线性动力学理论,揭示系统存在的各种非线性动力学行为,提出转子-轴承系统的非线性动力学设计方法,研究旋转机械中存在的各种实际问题,这对提高旋转机械运行的稳定性、安全性、可靠性具有重要的现实意义和实际工程背景。 随着非线性动力学理论的发展,非线性转子动力学理论和方法也受到了关注,大量的研究成果使转子动力学面貌一新。但现有的非线性动力学理论和方法在解决高维动力系统方面还存在困难,而工程实际中的转子-轴承-基础系统是一个复杂的高维系统,从而吸引了更多的研究者从事这方面的研究,特别是现代非线性动力学理论在转子动力学中的应用,已成为当今国

结构力学大作业(华科)

一、任务 1.求解多层多跨框架结构在竖向荷载作用下的弯矩以及水平荷载作用下的弯矩和 各层的侧移。 2.计算方法: (1)用近似法计算:水平荷载作用用反弯点法计算,竖向荷载作用采用分层法和二次力矩分配法计算。 (2)用电算(结构力学求解器)进行复算。 3. 就最大相对误差处,说明近似法产生误差的来源。 4. 将手算结果写成计算书形式。 二、结构形式及各种资料 1. 计算简图:如图1所示。 2. 基本计算参数 底层柱bXh(mm) 其它层bXh(mm) 边梁bXh(mm) 中间梁bXh(mm) 500X500 450X450 250X450 250X450 材料弹性模量: 72 3.210/ h E kN m =? 竖向荷载: 2 1 =23/ g kN m,2 2 =20/ g kN m 水平荷载: =32 p F kN 1,2 =18 P F kN 3. 荷载分组: (1)计算水平荷载(见图2);(2)计算竖向恒载(见图3); L1L2H1 H2 H2 H2 H2 F F F F F 图1 计算简图图2 水平荷载作用

g2 g1 g1 g1 g1 q2 q1 图3 竖向荷载作用 三、计算内容 ?水平荷载 1、反弯点法 (1)求柱的剪力 由所给数据可得各层梁柱的线刚度(单位:kN·m)如下表: i底柱i其它柱i左梁i右梁 34792363331270825417 第五层柱;F Q14 = F Q25 = F Q36 = 18/3kN = 6kN 第四层柱;F Q47 = F Q58 = F Q69 = 50/3kN 第三层柱;F Q710 = F Q811 = F Q912 = 82/3kN 第二层柱;F Q1013 = F Q1114 = F Q1215 = 114/3kN 第一层柱;F Q1316 = F Q1417 = F Q1518 = 146/3kN (2)求柱的弯矩 第五层柱;M 14 = M 41 = M 25 = M 52 = M 36 = M 63 = 6×3/2 = 9kN·m 第四层柱;M 47 = M 74 = M 58 = M 85 = M 69 = M 96 = 50/3×3/2 = 25kN·m 第三层柱;M 710 = M 107 = M 811 = M 118 = M 912 = M 129 = 82/3×3/2 = 41kN·m 第二层柱;M 1013 = M 1310 = M 1114 = M 1411 = M 1215 = M 1512 = 114/3×3/2 = 57kN·m 第一层柱;M 1316 = M 1417 = M 1518 = 146/3×4.8/3 = 77.87kN·m M 1613 = M 1714 = M 1815 = 146/3×2×4.8/3 = 155.74kN·m (3)求梁的弯矩 分别取结点1、2为隔离体 1 M12 ∑M1=0 M12=M14=9kN·m M14

转子动力学

转子动力学 工程界和科学界关心转子振动的历史已有二百多年,1869年英国的W.J.M.兰金关于离心力的论文和1889年法国的C.G.P.de拉瓦尔关于挠性轴的试验是研究这一问题的先导。随着近代工业的发展,逐渐出现了高速细长转子。由于它们常在挠性状态下工作,所以其振动和稳定性问题就越发重要。转子动力学的研究内容主要有以下5个: ①临界转速 由于制造中的误差,转子各微段的质心一般对回转轴线有微小偏离。转子旋转时,由上述偏离造成的离心力会使转子产生横向振动。这种振动在某些转速上显得异常强烈,这些转速称为临界转速。为确保机器在工作转速范围内不致发生共振,临界转速应适当偏离工作转速例如10%以上。临界转速同转子的弹性和质量分布等因素有关。对于具有有限个集中质量的离散转动系统,临界转速的数目等于集中质量的个数;对于质量连续分布的弹性转动系统,临界转速有无穷多个。计算大型转子支承系统临界转速最常用的数值方法为传递矩阵法。其要点是:先把转子分成若干段,每段左右端4个截面参数(挠度、挠角、弯矩、剪力)之间的关系可用该段的传递矩阵描述。如此递推,可得系统左右两端面的截面参数间的总传递矩阵。再由边界条件和固有振动时有非零解的条件,籍试凑法求得各阶临界转速,并随后求得相应的振型。

②通过临界转速的状态 一般转子都是变速通过临界转速的,故通过临界转速的状态为不平稳状态。它主要在两个方面不同于固定在临界转速上旋转时的平稳状态:一是振幅的极大值比平稳状态的小,且转速变得愈快,振幅的极大值愈小;二是振幅的极大值不像平稳状态那样发生在临界转速上。在不平稳状态下,转子上作用着变频干扰力,给分析带来困难。求解这类问题须用数值计算或非线性振动理论中的渐近方法或用级数展开法。 ③动力响应 在转子的设计和运行中,常需知道在工作转速范围内,不平衡和其他激发因素引起的振动有多大,并把它作为转子工作状态优劣的一种度量。计算这个问题多采用从临界转速算法引伸出来的算法。

ANSYS_中的转子动力学计算

ANSYS 中的转子动力学计算 雷先华 安世亚太 转子动力学是固体力学的一个重要分支,它主要研究旋转机械的“转子-支承”系统在旋转状态下的振动、平衡和稳定性问题,其主要研究内容有几个方面:临界转速、动力响应、稳定性、动平衡技术和支承设计。在旋转机械研究设计中,转子动力学的性能分析是极其重要的一个方面。 传统的转子动力学分析采用传递矩阵方法进行,由于将大量的结构信息简化为极为简单的集中质量—梁模型,不能确保模型的完整性和分析的准确度;而有限元在处理转子动力学问题时,可以很好地兼顾模型的完整性和计算的效率,但多年来转子的“陀螺效应”一直是制约转子动力学有限元分析的“瓶颈”问题。ANSYS 很好地解决了动力特性分析中“陀螺效应”影响的问题,而且陀螺效应的考虑不受计算模型上的限制,使得转子动力学有限元分析变得简单高效。 本文对ANSYS 的转子动力学计算功能进行简要介绍。 1 ANSYS 转子动力学的理论基础 ANSYS 转子动力学分析中,两种参考坐标系可供选择:静止坐标系和旋转坐标系。空间点P 在静止坐标系(其原点在O ′)下的位置矢量为r ′,在旋转坐标系(其原点在O )下的位置矢量为r 。 在静止坐标系下转子的动力方程为: [][][]{}{}([]){}{}M u C C u K u F gyr +++=&&& []C gyr 为陀螺效应矩阵。 在旋转坐标系下转子的动力方程为: [][][]{}{}([]){}([]){}M u C C u K K u F cor spin r r r +++?=&&& []C cor 为哥氏效应矩阵,[]K spin 为旋转软化效应刚度矩阵。 2 ANSYS 转子动力学的计算功能和新技术 ANSYS 转子动力学计算包含如下功能: – 无阻尼临界转速分析 – 不平衡响应分析 – 阻尼特征值分析

重庆大学_结构动力学大作业

研究生课程考核试卷 (适用于课程论文、提交报告) 科目:结构动力学大作业教师:刘纲姓名:学号: 专业:土木工程类别:学术上课时间:2016 年9 月至2016 年12 月考生成绩: 阅卷评语: 阅卷教师(签名) 重庆大学研究生院制

土木工程学院2016级硕士研究生考试试题 科目名称:结构动力学考试日期:2016年12月总分:20分 1、按规定设计一个2跨3层钢筋混凝土平面框架结构(部分要求如附件名单所示;未作规定部分自定)。根据所设计的结构参数,求该结构的一致质量矩阵、一致刚度矩阵; 2、至少采用两种方法求该框架结构的频率和振型; 3、输入地震波(地震波要求如附件名单所示),采用时程分析法,利用有限元软件或自编程序求出该框架结构各层的线性位移时程反应。 要求给出: (1)框架结构图,并给出一致质量矩阵和一致刚度矩阵; (2)写出两种方法名称及对应的频率和振型; (3)输入地震波的波形图,计算所得各楼层位移反应时程图。 第 1 页共1页

1框架概况 1.1框架截面尺寸 框架立面图如图 1.1所示,各跨跨度为13900L mm =,各层建筑层高均为 33700L mm =,对应的梁截面分别为2200400mm ?,柱截面均为2300300mm ?。设楼层进深为24200L mm =,板厚为100mm ,顶层板厚取最小值120mm 。 图1.1框架立面图 1.2动力自由度 框架结构可以理想化为在节点处相互连接的单元(梁和柱)的集合。设 梁、柱的轴向变形均忽略不计,只考虑横向平面位移,则该框架有3平动自由度

和9角自由度,共12自由度。自由度编号及梁柱单元编号如图1.2所示。 图1.2自由度编号和梁柱单元编号图 1.3力学参数计算 梁柱混凝土都采用C30,其弹性模量为42 3.010/ E N mm =? 2 200400mm ?梁截面惯性矩: 33 4 1 200400 1212 bh I mm ? == 2 300300mm ?柱截面惯性矩: 33 4 2 300300 1212 bh I mm ? == 2 200400mm ?梁刚度: 3 442 1 200400 3.010 3.210 12 EI kN m ? =??=?? 2 300300mm ?柱刚度: 3 442 2 300300 3.010 2.02510 12 EI kN m ? =??=?? 2 200400mm ?梁的单位长度质量(按照计算重力荷载代表值的方法计算):

相关文档
最新文档