钢丝绳碰撞动力学模型
- 格式:pdf
- 大小:724.37 KB
- 文档页数:6
摘要吊物的摆动是影响吊车装卸效率的主要原因。
电子防摇作为一种主动防摇方式,它将减摇和运行控制结合起来考虑,不依赖于司机的操作经验,可以有效的提高吊车的装卸效率,减轻司机的工作强度,是实现港口、厂矿装卸自动化的趋势。
本文首先建立了桥式吊车运动系统的数学模型,并搭建了系统的仿真实验模型,设计了双闭环PID控制器来实现吊车系统的防摆和定位控制。
针对常规PID 控制器很难满足桥式吊车这类控制参数变化很大的复杂系统对控制精度的要求,设计了非线性PID控制器(即PID参数随误差的变化而变化),该控制方案可以消除系统静差,缩短系统响应时间,抗干扰能力较强。
运用增量型PID控制算法对PID控制策略进行了数字实现。
本系统实现了吊车到达目的地吊钩不在摆动的目的,使所吊重物在到达指定地点时,摆角为零。
从而避免了作业时,照成不必要的事故。
关键词:桥式吊车;防摆;非线性PIDAbstractGenerally speaking, the loading efficiency of crane is mainly influenced by swing of hanging objects. As an active method for avoiding swing, electronic anti-swing is not dependent on the driver’s experience and can combine swing-decreasing with movement-control to improve the efficiency of crane, and lighten the intensity of drivers. Therefore, it will be used widely for loading of port and factory.In this thesis, the mathematical model of the overhead crane motion system is established and the simulation model is also built . The two closed-loop PID controller is designed to achieve the control of anti-swing and orientation. Comparing with the controlling results of using conventional PID algorithm as a controller on different disturbance conditions, we can find that it is difficult to meet the accuracy requirements of the practical operation of crane system which has easily variable control parameters in the process of running. So we design the nonlinear PID whose variable parameters change with error. The results of simulation indicate that the steady-state error can be eliminated, and the response time of system can be shortened, in addition, the disturbance rejection ability of system can be strengthened. Then, we complish the digital realization of PID control strategy with increment PID control algorithm.This system has realized the crane hook is not swinging the purpose of arriving at our destination, make lifting heavy weights in got to the designated place, pendulum Angle is zero. Avoiding the homework, as into unnecessary accidents.Key Words:Overhead crane;Anti-swing;Nonlinear PID目录第1章绪论 (1)1.1 课题的背景及意义 (1)1.2 国内外研究现状 (2)1.3 本文研究主要内容 (4)第2章系统建模 (5)2.1 问题的提出 (5)2.2 建模机理 (5)2.3 系统模型的建立 (6)第3章方案论证 (10)3.1 系统方案的论证 (10)3.2 总体方案设计 (11)3.3 传感器的选型 (11)3.4 转换器的选型 (13)3.5 单片机的选型 (15)3.6 变频器的选型 (16)第4章硬件电路设计 (17)4.1 数据检测模块 (17)4.2 滤波电路模块 (19)4.3 转换电路 (20)4.4 单片机最小系统 (20)4.5 变频器调速模块 (21)4.6 键盘电路 (22)4.7 电源电路 (23)第5章软件设计 (24)5.1 系统流程图 (24)5.2 检测控制流程图 (25)5.3 变频器子程序 (25)5.4 报警程序 (26)5.5 PID控制的设计 (27)5.5.1 非线性PID设计 (28)5.5.2 仿真实验验证 (30)5.5.3 增量型PID (37)5.5.4 吊车系统数字PID的实现 (38)第6章总结 (40)参考文献 (41)致谢 (42)附录Ⅰ ............................................................................................... 错误!未定义书签。
基于MATLAB的臂架型起重机起升机构动力学分析摘要:本文对起重机械起升机构特点进行了分析,就臂架型起重机起升机构多质点动力学系统建立了动力学模型,确定了模型中的各个参数,利用MATLAB构造相应的仿真系统进行了动力学分析。
关键词:起重机起升机构MATLAB 机械动力学仿真1 引言起升机构是起重机械的最基本机构,目前,对臂架起重机起升机构的动力学研究还不充分,在研究中或仅考虑系统特性,或仅考虑机械系统的振动,对动力学方面的研究还不够,与实际存在一定差距,本文将起升机构简化为三质量二自由度系统,对起升过程进行了动力学分析.2 起升机构动力学模型的建立2.1 起升机构在起重机中,用以提升或下降货物的机构称为起升机构,一般称为卷扬式。
起升机构是起重机中最重要、最基本的机构。
电动机驱动的起升机构由电动机、联轴器、制动器、减速器、卷筒、钢丝绳、滑轮和吊具组成。
以轮胎式起重机伸缩臂起升机构为例,建立三质量二自由度模型,为了便于分析问题,突出主要矛盾而做如下假定:1)在起动、制动过程中时间相对较短,该过程中起重钢丝绳的长度变化,所以假定起动、制动过程中钢丝绳的长度不变;2)钢丝绳缠绕在卷筒上,假定钢丝绳在滑轮槽中不打滑;3)在起动的过程中采用预备档,假定在传动机构中各齿轮间不存在间隙冲击;4)假定动力装置运转速度只与其机械特性有关,臂架等结构振动不影响电动机运转速度;5)忽略了起升机构中各转动构件扭转对电动机机械特性的反馈影响;6)忽略地面、支腿、车架以及其他机构的挠度;7)本研究仅考虑起升平面内的振动,扭摆横向振动暂不考虑。
由1)、2)可得钢丝绳刚度为常数,在系统模型简化时由3)可将传动机构等效为一质量,从4)、5)可知起升机构等效驱动力矩由电机机械特性确定,由6)知系统各质量的运动由系统弹性件的变形所决定。
基于上述假设,将吊臂质量等效在吊臂上部的端点处、吊重及吊具质量取在吊重上,为二质量的振动系统;将驱动装置、起升卷筒及两者间的传动装置等效在卷筒轴上。
钢丝绳索具压制成形仿真分析马平;刘杰;王志勇;肖全【摘要】钢丝绳索具是广泛运用于海洋工程中起柔性连接作用的有效工具之一,具有强度高、承载能力强、韧性好等优点.由于索具中铝套壁厚的大小与索具承载力密切相关,壁厚的大小最直观的结果就是延伸率,行业中采用控制延伸率为25%~35%的方法实现索具满足承载力的要求.基于SolidWorks-ABAQUS平台,结合显示动力学分析方法对模型的压制成形过程进行仿真分析,并采用控制变量法改变索具接头中铝套的壁厚尺寸,得出索具成形所需压制力、铝套的延伸率以及钢丝绳索具疲劳寿命的变化规律.研究结果显示1X7+IWS-18mm型钢丝绳索具(对应铝套壁厚为7mm)压制成形后最大应力分布在钢丝绳自接触部分,大小为912.6MPa,成形所需压制力为1097.51kN,铝套延伸率为25.53%,疲劳次数为353次.【期刊名称】《机械设计与制造》【年(卷),期】2019(000)002【总页数】4页(P91-94)【关键词】钢丝绳索具;显示动力学;延伸率;应力;压制力【作者】马平;刘杰;王志勇;肖全【作者单位】广东工业大学机电工程学院,广东广州 510006;广东工业大学机电工程学院,广东广州 510006;建峰索具有限公司,广东广州 511356;广东工业大学机电工程学院,广东广州 510006【正文语种】中文【中图分类】TH16;TG3721 引言索具是起吊装备与重物之间连接件的统称,其具有强度高、承载能力强等优点,因此在海洋资源开采、桥梁工程、大型工程平台的起吊等领域广泛应用,目前已成为应用最广泛的挠性构件之一[1]。
压制索具采用的标准有GB/T30589―2014《钢丝绳绳端套管压制索具》[2]、GB/T16762―2009《一般用途钢丝绳吊索特性和技术条件》[3]、GB/T 6946—2008《钢丝绳铝合金压制接头》、欧盟标准:EN13411-3:2004《钢丝绳端固结法第三部分:套管和套管固结》[4]、EN13414-1:2003《钢丝绳吊索—安全第一部分:一般提升用途吊索》[5]等。
2022年第51卷第1期第20页石油矿场机械犗犐犔 犉犐犈犔犇 犈犙犝犐犘犕犈犖犜2022,51(1):20 28文章编号:1001 3482(2022)01 0020 09基于犗狉犮犪犉犾犲狓水下采油树钢丝绳下放作业窗口研究张 崇1,孟文波1,任冠龙1,唐咸弟2,赵苏文2,关清涛2,赵维青2,李新妍2(1.中海石油(中国)有限公司湛江分公司,广东湛江524057;2.中海油能源发展股份有限公司工程技术深水钻采技术公司,广东湛江524057)摘要:水下采油树的下放安装作业环境载荷非常复杂,下放作业窗口是我国南海水下采油树下放过程中亟待解决的关键问题。
针对东方1 1气田73m水下采油树钢丝缆绳下放安装方法,分析了水下采油树下放作业过程。
利用OrcaFlex建立了工程船 钢丝缆绳 水下采油树系统的耦合模型,研究了不同浪高、流速、波浪方向角等关键因素对水下采油树安装过程的钢丝绳受力和采油树运动响应的影响。
利用正交试验方法研究了各个影响因素的关键度排序,形成了水下采油树钢丝缆绳下放的作业窗口。
研究结果表明:当浪高在[5.06m,0m],流速为[0m/s,1.09m/s]的封闭区域(安全作业区域)内,适合于东方1 1气田水下采油树下放安装作业。
研究结果可为水下采油树钢丝绳下放安装过程提供一定的指导。
关键词:水下采油树;OrcaFlex;钢丝绳法;下放安装;下放作业窗口中图分类号:TE952 文献标志码:A 犱狅犻:10.3969/j.issn.1001 3482.2022.01.003犅犪狊犲犱狅狀犗狉犮犪犉犾犲狓犛狌犫狊犲犪犜狉犲犲犠犻狉犲犚狅狆犲犇狉狅狆犗狆犲狉犪狋犻狅狀犠犻狀犱狅狑犚犲狊犲犪狉犮犺ZHANGChong1,MENGWenbo1,RENGuanlong1,TANGXiandi2,ZHAOSuwen2,GUANQingtao2,ZHAOWeiqing2,LIXinyan2(1.犣犺犪狀犼犻犪狀犵犅狉犪狀犮犺,犆犖犗犗犆犆犺犻狀犪犔犻犿犻狋犲犱,犣犺犪狀犼犻犪狀犵524057,犆犺犻狀犪;2.犇犲犲狆狑犪狋犲狉犇狉犻犾犾犻狀犵犜犲犮犺狀狅犾狅犵狔犅狉犪狀犮犺,犆犖犗犗犆犈狀犲狉犵狔犜犲犮犺狀狅犾狅犵狔牔犛犲狉狏犻犮犲狊犔狋犱.,犣犺犪狀犼犻犪狀犵524057,犆犺犻狀犪)犃犫狊狋狉犪犮狋:Theenvironmentloadfortheinstallationofsubseatreesisverycomplicated,andthewindowoftheoperationisakeyissuethatneedstobesolvedurgentlyintheprocessofSubseatreesintheSouthChinaSea.Aimingattheinstallationmethodof73msubseatreewireropeinstallationinDongfang1 1gasfield,theoperationprocessofsubseatreedroppingoperationwasanalyzed,andthecouplingmodelofengineeringship wireropesubseatreesystemwasestab lishedusingOrcaFlex,andthedifferenceswerestudied.Theinfluenceofkeyfactorssuchaswaveheight,flowvelocityandwavedirectionangleontheforceofthewireropeduringtheinstallationprocessofthesubseatreeandthemovementresponseofthetreewerestudied.Theorthogonaltestmethodwasusedtostudythecriticalityrankingofeachinfluencingfactor,form ingasubseatree.Theresearchresultsshowthatwhenthewaveheightis[5.06m,0m]andtheflowvelocityis[0m/s,1.09m/s]inaclosedarea(safeoperationarea),itissuitableforthe 收稿日期:2021 07 10 基金项目:国家工信部创新专项“水下油气生产系统工程化示范应用”(2020 HGCBGXZX ZX GCFW 1780)课题2“水下采油树及其配套工具工程化技术研究”(MC_201901_S01_02) 作者简介:张 崇(1984 ),男,高级工程师,研究方向为海上油田钻完井工艺技术,E mail:zhangchong3@cnooc.com.cn。
- 18
第2章 质点动力学 一、基本要求 1.理解冲量、动量,功和能等基本概念; 2.会用微积分方法计算变力做功,理解保守力作功的特点; 3.掌握运用动量守恒定律和机械能守恒定律分析简单系统在平面内运动的力学问题的思想和方法。
二、基本内容 (一)本章重点和难点: 重点:动量守恒定律和能量守恒定律的条件审核、综合性力学问题的分析求解。 难点:微积分方法求解变力做功。 (二)知识网络结构图:
公式只有保守内力做功条件能量守恒定律公式合外力为条件动量守恒定律守恒定律动能定理动量定理基本定理能功冲量动量基本物理量)()0(
(三)容易混淆的概念: 1.动量和冲量 动量是质点的质量与速度的乘积;冲量是合外力随时间的累积效应,合外力的冲量等于动量增量。
牛二定律 动量定理 变力做功 角动量定理(不要求) 动量守恒定律 动能定理
功能原理 机械能守恒定律 - 19
2.保守力和非保守力 保守力是做功只与始末位置有关而与具体路径无关的力,沿闭合路径运动一周保守力做功为0;非保守力是做功与具体路径有关的力。
(四)主要内容: 1.动量、冲量
动量:pmv
冲量:21ttdtFI
2.动量定理: 质点动量定理:2112ttvmPPdtFI
质点系动量定理:dtPdF 3.动量守恒定律: 当系统所受合外力为零时,即0exF时,或inexFF
系统的总动量保持不变,即:niiiCvmP1 4.变力做功: drFrdFWBABAcos(为)之间夹角与rdF
直角坐标系中: )ddd( zFyFxFWzyBAx 5.动能定理:
(1)质点动能定理:k1k221222121EEmvmvW - 20
(质点所受合外力做功等于质点动能增量。) (2)质点系动能定理:niniEEWW1kio1kiinex (质点系所受外力做功和内力做功之和等于质点系动能增量。) 6.保守力、势能、功能原理: (1)保守力:做功只与始末位置有关,与经历的路径无关的力
基于MATLAB的爬行升降装置钢丝绳张力分析针对高层建筑室外便携式提升设备少见的现状,并通过利用齿轮啮合传动的工作原理,研究开发了一种新型的爬行升降装置。
运用力学的基础理论知识,研究分析钢丝绳的张力变化规律,建立了高层建筑室外新型爬行升降装置关键部件即钢丝绳的力学建模,得到了钢丝绳张力变化的时变函数。
利用MA TLAB中SIMULINK模块,进行相应数据的仿真分析,验证了新型爬行升降装置中机构设计的合理性与可行性,对高层建筑室外设备的安装与维修、清洗与维护,具有重要的意义。
标签:爬行升降装置;时变函数;动力学建模;钢丝绳;仿真引言随着我国国民经济发展和城市化进程的加快,高层建筑不断涌现,诸如:空调室外机、各类广告牌等的安装和维护工作也随之增多。
然而,用于提升的设备却少之又少,且适用性和灵活性较差,很多情况下不得不完全由人力完成,不仅劳动强度大,而且极其不安全。
故为了更好的适应社会发展,更好的保障人身安全,研究开发了一种新型爬行升降装置。
钢丝绳作为该装置的关键部件,它的性能好坏将直接关系着整个装置的安全性与可靠性。
因此,文章对该设备的钢丝绳进行受力与变形分析,并进行机构设计的可行性分析。
1 新型爬行升降装置工作原理此种爬行升降装置,主要由驅动电机、主动摩擦轮、从动摩擦轮、钢丝绳、吊装装置以及安全控制装置组成。
它利用驱动电机驱动主动摩擦轮旋转,与从动摩擦轮啮合传动,钢丝绳缠绕在主从动摩擦轮上(如图1所示),依靠钢丝绳与两个摩擦轮之间的摩擦力,实现整机设备沿钢丝绳的上升与下降。
2 爬行升降装置钢丝绳的力学模型2.1 新型爬行升降装置运动分析2.2 建立力学模型当升降装置按如图(2a)所示方向运动时(不计钢丝绳与摩擦轮之间的相对滑动),位于升降装置下面的钢丝绳长度会不断的发生变化,导致整个系统的等效质量也随之变化,从而影响系统的振动频率和振幅的变化,文章采用瑞利法与钢丝绳张力等效于刚体动张性的方法相结合,来计算钢丝绳的等效质量。
第26卷第10期 V ol.26 No.10 工 程 力 学 2009年 10 月 Oct. 2009 ENGINEERING MECHANICS197———————————————收稿日期:2008-06-16;修改日期:2008-12-09作者简介:*方子帆(1963―),男,湖北黄冈人,教授,博士,博导,副院长,从事车辆系统动力学与控制研究(E-mail: fzf@); 吴建华(1983―),男,湖北大冶人,硕士,从事机械振动与控制研究(E-mail: wujianhua83@);何孔德(1973―),男,湖北宜昌人,副教授,硕士,从事机械振动与控制研究(E-mail: hekongde@); 张明松(1965―),男,湖北荆州人,副教授,学士,从事结构设计与机械振动研究(E-mail: zms@).文章编号:1000-4750(2009)10-0197-06钢丝绳碰撞动力学模型*方子帆,吴建华,何孔德,张明松(三峡大学机械与材料学院,湖北,宜昌 443002)摘 要:以钢丝绳及其连接结构为对象,对其碰撞动力学模型进行研究。
将钢丝绳离散为单元模型,利用相对坐标关系建立其动力学模型,并将其连接结构以集中质量模型作为钢丝绳端部约束条件引入到钢丝绳动力学模型中,建立钢丝绳及其连接结构的动力学模型。
将钢丝绳的碰撞接触力引入到钢丝绳及其连接结构的动力学模型中,建立这类结构的碰撞动力学模型。
在RecurDyn 环境中建立了具有横向和垂挂空间姿态的钢丝绳及其连接结构的动力学仿真模型,并进行仿真研究。
研究结果表明这些模型可以用作刚柔混合结构的动力学分析,同时能够实现这类结构的可视化动态仿真。
关键词:钢丝绳;碰撞;动力学模型;相对坐标法;RecurDyn 中图分类号:O313; TH113.2 文献标识码:ATHE IMPACT DYNAMIC MODEL OF STEEL CABLES*FANG Zi-fan , WU Jian-hua , HE Kong-de , ZHANG Ming-song(College of Mechanical and Material Engineering, China Three Gorges University, Yichang, Hubei 443002, China)Abstract: A dynamic model of steel cable is established by a discrete elements method considering relative coordinate relationship. Its connective structures are modeled as lumped mass and incorporated into the steel cable dynamic model as end constraints. Introducing the steel cables contact-impact force into the established dynamic model of steel cable with their connective structure, the impact dynamic model of steel cables with their connective structure is established finally. An example is presented, which is steel cables consisting of a transversely placed and a vertically placed steel cable with their connective structures. The impact dynamic simulation model is established in RecurDyn. The results show that the proposed impact dynamic model can be applied in the dynamic analysis of structural systems consisting flexible bodies and rigid bodies. Key words: steel cable; impact; dynamic model; relative coordinate method; RecurDyn由于钢丝绳的材料非线性和几何非线性问题,通常采用基于Lagrange 相对坐标系模型和基于Cartesian 坐标的绝对坐标系模型建立这类结构的动力学模型。
以Song J O [1]、Simo J C [2]、Avello A [3]等为代表的学者将柔性体的大位移及弹性变形用相对惯性坐标系的单元结点坐标描述,推导出变形体的应变、位移关系。
Wu 与Haug 等[4]使用向量变分方法并结合虚功原理,采用相对坐标再叠加弹性体的模态坐标,建立了柔性多体系统的相对坐标动力学建模方法。
Chen 与Shabana [5]用绝对坐标法建立了柔性多体系统的动力学模型。
于清与洪嘉振[6]对上述两种建模方法进行了评述,认为相对坐标方法具有动力学方程广义坐标和约束方程少、计算效率高的优点,但程式化较绝对坐标方法差。
本文以198 工 程 力 学钢丝绳及其连接结构为研究对象,试图利用相对坐标法建立其碰撞动力学模型,并进行应用研究。
1 钢丝绳及其连接结构动力学模型1.1 钢丝绳动力学模型将钢丝绳离散成一系列梁单元,钢丝绳的绝对结点变形是沿着钢丝绳拓扑图的一条排序链的相对变形的累加,所以用相对位移作为广义坐标向量来描述钢丝绳的运动,建立钢丝绳的动力学模型。
图1是由两个梁单元组成的系统,在图1(b)中,结点1i −和i 分别被假定为结点i 和1i +的内侧结点。
--X Y Z 是惯性坐标系,--k k k x y z (,k i j =)是固联在结点k 上的结点参考坐标系,k r 是结点k 的位移矢量。
(1)(1)(1)--i i i i i i x y z −−−是固联在结点i 上的结点参考坐标系,第一个下脚标1i −表示第二个小脚标i 的内侧结点号。
在未变形状态时,坐标系(1)(1)(1)--i i i i i i x y z −−−和(1)(1)(1)--i i i x y z −−−的方向一致。
在惯性坐标系--X Y Z 上测量绝对结点位移,而在结点的内侧结点参考坐标系上测量相对结点位移。
(a) 两个梁单元(b) 梁单元的图形描述图1 梁单元的结点位置关系Fig.1 The relation of node position for the beam elements用(1)i i −′u 和(1)i i −′Θ分别描述一个结点的相对结点位置和方向的广义坐标。
结点i 在--X Y Z 的结点位置和方向由结点1i −和相对结点位移表示为:0(1)(1)(1)(1)()i i i i i i i r r −−−−′′=++A s u (1)1(1)(1)(1)()i i i i i i i i −−−−′=A A D ΘC (2)其中:T (1)(1)1(1)2(1)3[]i i i i i i i i θθθ−−−−′′′′=Θ (3)式(1)和式(2)中:k A (1,k i i =−)表示结点k 的参考坐标系相对惯性坐标系的转换矩阵;(1)i i −C 表示坐标系--i i i x y z 相对于(1)(1)(1)--i i i i i i x y z −−−的转换矩阵;0(1)i i −′s 表示梁未变形时相对坐标系(1)(1)(1)--i i i x y z −−−的位置矢量;(1)i i −′u 表示结点i 相对结点1i −坐标系的变形矢量;(1)i i −D 是由于坐标系(1)(1)(1)--i i i i i i x y z −−−相对结点1i −坐标系的转动转换矩阵,通过欧拉角描述:(1)1(1)12(1)23(1)3()()()i i i i i i i i θθθ−−−−′′′=D D D D (4) 式(1)的变分为:0T TT (1)(1)(1)(1)(1)(1)(1)()i i i i i i i i i i i i i r r δδδ−−−−−−−′′′′′=−++A A s u A u (5)式(5)中符号“~”表示由自身的矢量元素组成的斜对称矩阵,且(1)i i −A 满足:T (1)(1)i i i i −−=A A A (6)结点i 和1i −的虚转动关系为: T T (1)(1)(1)(1)(1)i i i i i i i i δπδπδ−−−−−′′′=+A A H Θ (7) (1)2(1)(1)1(1)1(1)2(1)1(1)1(1)210sin()0cos()sin()cos()0sin()cos()cos()i i i i i i i i i i i i i i i i θθθθθθθ−−−−−−−−′⎡⎤⎢⎥′′′=−⎢⎥⎢⎥′′′⎣⎦H (8) 合并式(5)和式(7)得到以下一对相邻单元的递归虚位移方程:(1)1(1)(1)2(1)i i i i i i i i Z B Z B q δδδ−−−−=+ (9)0T T T T T T (1)(1)(1)T (1)(1)(1)(1)1T(1)T (1)(1)2T (1)(1)[]()00000k k k i i i i i i i i i i i i i i i i i ii i i i i i Z r q I B I I B H δδδπδδδ−−−−−−−−−−−−⎧′′=⎪⎪⎡⎤′′=⎣⎦⎪⎪⎡⎤′′−+⎡⎤⎪=⎨⎢⎥⎢⎥⎪⎢⎥⎣⎦⎣⎦⎪⎡⎤⎪⎡⎤=⎢⎥⎪⎢⎥⎢⎥⎣⎦⎪⎣⎦⎩u Θs u A A A A (10) 式(10)中(1)1i i B −和(1)2i i B −仅是结点1i −和i 的相对位移(1)i i q −的函数。
整个系统在绝对结点坐标系和相对结点坐标系下的虚位移关系可以沿着系统拓扑图的一条排序链通过反复运用式(9)获得:Z B q δδ= (11)由虚功原理得到结点1i −和i 组成的梁单元i 的动力学方程为:yz xz i z i −1y i −1 x i −1r i −1r iiz (i −1)ix (i −1)iy (i −1)iy ii −1x ii − 1 i −1i +1 i ……向后排序向前排序工 程 力 学 199i i i i i qq +=M K Q (12) 式(12)中:i M 为单元i 的质量;i K 为单元i 的刚度;i Q 为单元i 上的广义力。