传递矩阵法分类
典型的传递矩阵计算方法有Myklestad-Prohl传递矩阵法和Riccati传递矩阵法。Myklestad-Prohl传递矩阵法有很多优点，如矩阵的维数不会随着转子系统的自由度数的增加而增加、计算效率高、程序设计简单、占用内存少等等，所以在实际工程中得到了很广泛的应用。但是，这种方法在大量应用的过程中，人们发现这种方法也存在一些问题，就是当计算的频率较高、或者结构支承的刚度很大、或者结构的自由度较多时，会出现数值不稳定的现象，从而使计算分析结果的精度大大下降[2~3,39~40]。为此，1978年Horner和Pilley提出了Riccati传递矩阵法[39]，这种方法保留了Myklestad-Prohl传递矩阵法的全部优点，且计算精度高，数值上也比较稳定。Riccati传递矩阵法在使用过程中遇到的另一个问题是在特征根的搜索过程中剩余量有许多无穷大奇点，因此可能产生增根现象，1987年王正在研究了这一现象后给出了这种奇点的消除方法[40]。
传递矩阵-matlab程序
%main_critical.m %该程序使用Riccati传递距阵法计算转子系统的临界转速及振型 %本函数中均采用国际单位制 % 第一步：设置初始条件（调用函数shaft_parameters） %初始值设置包括：轴段数N，搜索次数M %输入轴段参数：内径d,外径D，轴段长度l，支撑刚度K，单元质量mm，极转动惯量Jpp[N,M,d,D,l,K,mm,Jpp]=shaft_parameters; % 第二步：计算单元的5个特征值（调用函数shaft_pra_cal） %单元的5个特征值： %m_k：:质量 %Jp_k：极转动惯量 %Jd_k：直径转动惯量 %EI：弹性模量与截面对中性轴的惯性矩的乘积 %rr：剪切影响系数 [m_k,Jp_k,EI,rr]=shaft_pra_cal(N,D,d,l,Jpp,mm); % 第三步：计算剩余量（调用函数surplus_calculate）,并绘制剩余量图 %剩余量：D1 for i=1:1:M ptx(i)=0; pty(i)=0; end for ii=1:1:M wi=ii/1*2+50; [D1,SS,Sn]=surplus_calculate(N,wi,K,m_k,Jp_k,JD_k,l,EI,rr); D1; pty(ii)=D1; ptx(ii)=w1 end ylabel(‘剩余量’)； plot(ptx,pty) xlabel(‘角速度red/s’); grid on % 第四步：用二分法求固有频率及振型图 %固有频率：Critical_speed wi=50; for i=1:1:4 order=i [D1,SS,Sn]=surplus_calculate(N,wi,k,m_k,Jp_k,Jd_k,l,EI,rr); Step=1; D2=D1; kkk=1; while kkk<5000 if D2*D1>0 wi=wi+step;
(整理)传输矩阵法.
传输矩阵法 一、 传输矩阵法概述 1. 传输矩阵 在介绍传输矩阵的模型之前，首先引入一个简单的电路模型。如图1(a)所示， 在(a)中若已知A 点电压及电路电流，则我们只需要知道电阻R ，便可求出B 点电压。传输矩阵具有和电阻相同的模型特性。 (a) (b) 图1 传输矩阵模型及电路模拟模型 如图1(b)所示，有这样的关系式存在：E 0=M(z)E 1。M(z)即为传输矩阵，它将介质前后空间的电磁场联系起来，这和电阻将A 、B 两点的电势联系起来的实质是相似的。 图2 多层周期性交替排列介质 传输矩阵法多应用于多层周期性交替排列介质（如图2所示）， M(z)反映的介质前后空间电磁场之间的关系，而其实质是每层薄膜特征矩阵的乘积，若用 j M 表示第j 层的特征矩阵，则有： 1 2 3 4 …… j …… N
（1） 其中， （2） j δ为相位厚度，有 （3） 如公式（2）所示，j M 的表示为一个2×2的矩阵形式，其中每个矩阵元都没有任何实际物理意义，它只是一个计算结果，其推导过程将在第二部分给出。 2. 传输矩阵法 在了解了传输矩阵的基础上，下面将介绍传输矩阵法的定义： 传输矩阵法是将磁场在实空间的格点位置展开，将麦克斯韦方程组化成传输矩阵形式，变成本征值求解问题。 从其定义可以看出，传输矩阵法的实质就是将麦克斯韦方程转化为传输矩阵，也就是传输矩阵法的建模过程，具体如下：利用麦克斯韦方程组求解两个紧邻层面上的电场和磁场，从而可以得到传输矩阵，然后将单层结论推广到整个介质空间，由此即可计算出整个多层介质的透射系数和反射系数。 传输矩阵法的特点：矩阵元少（4个），运算量小，速度快；关键：求解矩阵元；适用介质：多层周期性交替排列介质。 二、 传输矩阵的基础理论——薄膜光学理论 1.麦克斯韦方程组 麦克斯韦方程组由四个场量：D 、E 、B 、H ，两个源量：J 、ρ以及反映它们之间关系的方程组成。而且由媒质方程中的参数ε、μ、σ反映介质对电磁场的影响。方程组的实质是描述电磁场的传播，即：一个变化的磁场引起邻近区域的电场变化，而此电场的变化又引起邻近磁场的变化，如此进行下去，便可抽象出电磁场的传播。如图3 所示。 ? ? ? ???==∏=D C B A M z M N j j 1)(?? ??? ?????=j j j j j j j i i M δδηδηδcos sin sin cos j j j j d N θλπ δcos 2=ε
单纯形法的计算方法
第4章 单纯形法的计算方法单纯形法求解线性规划的思路: 一般线性规划问题具有线性方程组的变量数大于方程个数, 这时有不定的解。但可以从线性方程组中找出一个个的单纯形, 每一个单纯形可以求得一组解, 然后再判断该解使目标函数值是增大还是变小, 决定下一步选择的单纯形。这就是迭代,直到目标函数实现最大值或最小值为止。 4.1 初始基可行解的确定 为了确定初始基可行解, 要首先找出初始可行基, 其方法如下。 （1）第一种情况：若线性规划问题 max z = 从Pj ( j = 1 , 2 , ? , n)中一般能直接观察到存在一个初始可行基 （2）第二种情况：对所有约束条件是“ ≤”形式的不等式, 可以利用化为标准型的方法, 在每个约束条件的左端加上一个松弛变量。经过整理, 重新对 及 ( i = 1 , 2 , ? , m; j = 1 , 2 , ? , n)进行编号, 则可得下列方程组 显然得到一个m×m单位矩阵 以B 作为可行基。将上面方程组的每个等式移项得 令由上式得 又因 ≥0, 所以得到一个初始基可行解 （3）第三种情况：对所有约束条件是“ ≥”形式的不等式及等式约
束情况, 若不存在单位矩阵时, 就采用人造基方法。即对不等式约束减去一个非负的剩余变量后, 再加上一个非负的人工变量; 对于等式约束再加上一个非负的人工变量, 总能得到一个单位矩阵。 4.2 最优性检验和解的判别 对线性规划问题的求解结果可能出现唯一最优解、无穷多最优解、无界解和无可行解四种情况, 为此需要建立对解的判别准则。一般情况下, 经过迭代后可以得到： 将上代入目标函数，整理后得 令 于是 再令 则 (1) 最优解的判别定理 若为对应于基B的一个基可行解，且对于一切 且有则 为最优解。称为检验数。 (2) 无穷多最优解的判别定理 若为一个基可行解， 且对于一切 且有 又存在某个非基变量的检验数，则线性规划问题有无穷多最优解。 (3) 无界解判别定理 若为一个基可行解，有一个> 0 ,并且对i = 1 , 2 , ?, m,有≤0 , 那么该线性规划问题具有无界解(或称无最优解)。 4.3 基变换
机床动力学建模的拓展传递矩阵法
万方数据
万方数据
万方数据
万方数据
２０１０年１１月吴文镜等：机床动力学建模的拓展传递矩阵法７３ 刀。Ｑ＝Ｆ（９）Ｑ＝Ｅ５２２＇Ｊｏ＋Ｅ６２３’１０＋Ｅ７乙１１０＋ 毛毛．１０＋岛乞Ｊ０＋Ｅｌｏｚ７＇ｊ０＋ 层Ｉｌｚ８．１０＋层１２磊．１０＋Ｅ１３ｚＦ＋Ｅ３０ｚＤ（１０） Ｆ＝Ｅ１４互．１０＋Ｅ１５乞．１０＋巨６毛＇ｌｏ＋ 巨７２５’１０＋Ｅ１８乙Ｊ０＋Ｅ９２７＇ｌｏ＋ ￡２０磊＿ｌｏ＋Ｅ２ｌｚ９．１０＋￡２２磊＋Ｅ３ｌ乞（１１） 互．ｏ＝ＥｌＺ６．１＋Ｅ２２７．Ｉ＋Ｅ３２８．１＋层４２９．１（１２） 由式（７）～（１１）得 （五ｏＥ５一Ｅ１４）互Ｚ２．ｏ＋（正ｏＥ６一Ｅ１５）五ｚ３．ｏ＋ （五ｏＥ７一Ｅ１６）五乙．ｏ＋（五ｏＥ８一Ｅ１７）毛ｚ５．ｏ＋ （墨ｏＥ９一Ｅ１８）ｒ６瓦．１＋（互ｏＥｌｏ—Ｅ１９）弓Ｚ７．１＋ （互ｏＥｌｌ—Ｅ２０）磊ｚ８，１＋（正。巨２一Ｅ２１）写ｚ９．１＋ （７ｉｏＥｌ３一Ｅ２２）ｚ－＋（７ｉｏ岛。一百３１）ＺＦ＝０（１３） 由式（６）、（１２）得 互，Ｄ（Ｅｌ乙，Ｊ＋Ｅ２２７．１＋Ｅ３２８．Ｉ＋Ｅ４毛，１）＝ｒｌ，』Ｚ１．，（１４）对于状态矢量磊＇ｌ、历＇ｌ、ｚ８＇１、而，１均为刚体１上的状态矢量，位移元素线性相关，有 易３２７．１＝Ｅ２４互，，（１５） 易３磊，ｌ＝Ｅ２５五，Ｊ（１６） ￡２３２９．１＝Ｅ２６互．，（１７）联合（１３）～（１７）将其写成矩阵的形式有 瓦ｌｚａｌＩ＝０４８×ｌ（１８）ｚａｌｌ＝（乏，ｏ召。别，。罨。烈，。 ｚ五磊。罨。ｚ０砟磊）１ 磊和Ｚｏ分别为激振点和拾振点的状态矢量，兀¨为４８×６９的高维矩阵。 ３．２结合面参数 直线进给功能部件中主要存在直线滚动导轨结合面以及电动机定子与滑板之间的螺栓结合面。对于导轨结合面模型简化为１个法向线性弹簧一阻尼系统、１个横向的线性弹簧一阻尼系统和３个转动方向的扭转弹簧一阻尼系统，以综合反映结合部各方向的微幅振动。通过锤击试验分别测定导轨法向和横向及３个扭转方向的传递函数，定义法向为Ｚ，横向为ｙ，３个坐标轴分别为Ａ、Ｂ、Ｃ。 根据单自南度系统振动方程计算出导轨各方向的接触刚度，根据半功率法计算接触阻尼。最终计算得到导轨结合面参数如表ｌ所示。电动机与滑板之问的螺栓结合面参数如表２所示。导轨结合面参数测试结果见图７。 表ｌ导轨结合部参数结果 参数数值 刚度ｋｒ／（ＭＮ?ｍ‘１２５３ 刚度ｋＪ（ＧＮ?ｍ“１２．１４ 刚度“／（ｋＮ?ｍ?ｒａｄ。。１６９３ ｝９４度ｋｓ／（ＭＮ?ｍ?ｒａｄ‘）１．７３ 刚度ｋｄ（ｋＮ?ｍ?ｒａｄ。１７２７ 阻尼ｃ；ｌ（Ｎ?ｓ?ｍ“１６４１．５ 阻尼ｃＪ（Ｎ?ｓ?ｍ’）ｌ０３４．９ 雕尼“／（Ｎ?ｍ?ｓ?ｒａｄ。）０．１４４７ 阻尼ｃ洲Ｎ?ｍ?ｓ?ｒａｄ。。）２．０１１ 阻尼Ｃｃ／（Ｎ?１ｔｌ?ｓ?ｍｄ１）０９６０２ 表２螺栓结合部参数 参数数值 刚度ｋ，／（ＧＮ?ｍ。。１ｏ．２５ 刚度ｋ，Ｊ（ＧＮ?ｍ’）０，２５ 刚度ｋｆｌ（ＧＮ?ｍ。）２．１０ 阻尼ｃ．ｒ／（Ｎ?ｓ?ｍ。）１２５ 阻尼ｅ，Ｉ（Ｎ?ｓ?ｍ。。１１２５ 阻尼ｃ∥（Ｎ?ｓ?ｍ“）２５０ （ａ）测试现场 ｛｜｜卜Ｍ以旷藩三ｈ∥ 迎卜—ｔ——专—上‘＿妻蔫ｋ套 图７导轨结合面参数测试结果 ３．３滑板有限元自由度缩减模型建－ｆｒ 创建有限元自由度缩减模型首先采用通用有限元软件得到零件的有限元法（Ｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄ，ＦＥＭ）｝－莫－型，根据零件特点选择质量集中点、 结合面连接节点、外力作用节点以及需要考察的节 万方数据
传递矩阵法在结构振动响应分析中的应用
传递矩阵法在结构振动响应分析中的应用 【摘要】传递矩阵法因其简便、快捷，已被广泛应用于机械、航空和航天等领域。本文以航空发动机低压转子临界转速分析为例，对传递矩阵法在结构振动响应分析中的应用方法和分析步骤进行了详细的介绍，并给出了某型发动机低压转子在不同支承刚度下的临界转速。 【关键词】传递矩阵；振动响应；临界转速；转子动力学 0 引言 经典传递矩阵法是20 世纪20 年代建立起来的用于研究弹性构件组成的一维线性系统振动问题的方法。经过多年的发展和完善，已经可以用于求解多圆盘轴的扭转振动问题、梁的弯曲振动模态、轴的横向振动问题、系统的静态响应和扭矩载荷响应问题、以及一维结构的振动特性分析和复合梁的振动特性等结构动力学问题。并且，由于传递矩阵法建模灵活、计算效率高等优点，已在包括光学、声学、电子学、机器人学、机械、兵器、航空、航天等诸多现代工程技术领域中得到了广泛应用[1]。 应用传递矩阵法进行分析的一般步骤为：1）结构离散化；2）建立系统传递矩阵；3）特征方程求解。 1 结构离散化 航空发动机低压转子结构简化模型见图1： 其主要组件为压气机、涡轮和低压轴。低压转子通过前、中、后3个支点与发动机转子系统相连[2]。 将该结构进行离散化处理[3-5]，并将各支点简化为线弹性体后，得到图2所示模型。 离散化处理后，整个低压转子的质量将被转换为分布式质量节点。表1给出了离散化后各质量节点的质量分布情况。 2 建立系统传递矩阵 将连续结构进行离散化处理后，实体结构将被简化成等刚性无质量梁单元及分布质量点。 3 特征方程求解 以转子转速做为变量，在不同刚度参数下对特征值进行求解。在某一给定刚
利用传递矩阵法和Riccati传递矩阵法分析转子临界转速
利用传递矩阵法和Riccati 传递矩阵法分析转子临界转速 一、 所需求解转子参数 将转子简化为如下所示： 三个盘的参数为：1232 2212322 2 1 230.0160.050.0160.0120.0250.012P P P d d d I kg m I kg m I kg m I kg m I kg m I kg m ? =?=?=???=?=?=?? 另，阶梯轴的三段轴的截面惯性矩分别为： 414243 1.73.20.9J cm J cm J cm ?=? =??=? 三段轴的单位长度轴段的质量分别为：123 2.45/ 3.063/1.587/m kg m m kg m m kg m =?? =??=? 二、 试算转轴的传递矩阵 取试算转速1200/p rad s ω== ; 则，各轴段的传递矩阵分别为： 第1段 840.061.7102.45/l m J m m kg m -=??=???=?
1 1.0006e+000 6.0007e-00 2 5.2943e-007 1.0588e-008 3.7356e-002 1.0006e+000 1.7649e-005 5.2943e-007 6.3506e+00 3 1.2701e+002 1.0006e+000 6.0007e-002 2.1170e+005 6.3506e+003 3.7356e-002 H = 1.0006e+000 ??????? 第2段 840.153.2103.063/l m J m m kg m -=??=???=? 2 1.0145e+000 1.5044e-001 1.7595e-006 8.7927e-008 3.8782e-001 1.0145e+000 2.3506e-005 1.7595e-006 4.9669e+004 2.4821e+00 3 1.0145e+000 1.5044e-001 6.6353e+005 4.9669e+00 4 3.8782e-001 H = 1.0145e+000 ??????? 第3段 840.053.2103.063/l m J m m kg m -=??=???=? 3 1.0002e+000 5.0002e-002 1.9531e-007 3.2552e-009 1.4358e-002 1.0002e+000 7.8128e-006 1.9531e-007 5.5135e+003 9.1890e+001 1.0002e+000 5.0002e-002 2.2054e+005 5.5135e+003 1.4358e-002 H = 1.0002e+000 ??????? 第4段 840.033.2103.063/l m J m m kg m -=??=???=? 4 1.0000e+000 3.0000e-002 7.0313e-008 7.0313e-010 3.1013e-003 1.0000e+000 4.6875e-006 7.0313e-008 1.9848e+003 1.9848e+001 1.0000e+000 3.0000e-002 1.3232e+00 5 1.9848e+003 3.1013e-003 H = 1.0000e+000 ??????? 第5段 840.10.9101.587/l m J m m kg m -=??=???=?
改进传递矩阵法
JOURNAL OF SOUND AND VIBRATION Journal of Sound and Vibration 289(2006)294–333 A modi?ed transfer matrix method for the coupling lateral and torsional vibrations of symmetricrotor-bearing systems Sheng-Chung Hsieh a ,Juhn-Horng Chen b ,An-Chen Lee a,? a Department of Mechanical Engineering,National Chiao Tung University,1001Ta Hsueh Road, Hsinchu 30049,Taiwan,ROC b Department of Mechanical Engineering,Chung Hua University,Taiwan,ROC Received 27January 2004;received in revised form 9August 2004;accepted 8February 2005 Available online 28April 2005 Abstract This study develops a modi?ed transfer matrix method for analyzing the coupling lateral and torsional vibrations of the symmetricrotor-bearing system with an external torque.Euler’s angles are used to describe the orientations of the shaft element and disk.Additionally,to enhance accuracy,the symmetric rotating shaft is modeled by the Timoshenko beam and considered using a continuous-system concept rather than the conventional ‘‘lumped system’’concept.Moreover,the harmonic balance method is adopted in this approach to determine the steady-state responses comprising the synchronous and superharmonic whirls.According to our analysis,when the unbalance force and the torque with n ?frequency of the rotating speed excite the system simultaneously,the en t1T?and en à1T?whirls appear along with the synchronous whirl.Finally,several numerical examples are presented to demonstrate the applicability of this approach. r 2005Elsevier Ltd.All rights reserved. 1.Introduction Rotor dynamics plays an important role in many engineering ?elds,such as gas turbine,steam turbine,reciprocating and centrifugal compressors,the spindle of machine tools,and so on.Owing to the growing demands for high power,high speed,and light weight of the rotor-bearing https://www.360docs.net/doc/f9989678.html,/locate/jsvi 0022-460X/$-see front matter r 2005Elsevier Ltd.All rights reserved.doi:10.1016/j.jsv.2005.02.004 ?Corresponding author.Tel.:+88635728513;fax:88635725372. E-mail address:aclee@https://www.360docs.net/doc/f9989678.html,.tw (An-Chen Lee).
单纯形法的计算方法
第4章 单纯形法的计算方法 单纯形法求解线性规划的思路: 一般线性规划问题具有线性方程组的变量数大于方程个数, 这时有不定的解。但可以从线性方程组中找出一个个的单纯形, 每一个单纯形可以求得一组解, 然后再判断该解使目标函数值是增大还是变小, 决定下一步选择的单纯形。这就是迭代, 直到目标函数实现最大值或最小值为止。 4.1 初始基可行解的确定 为了确定初始基可行解, 要首先找出初始可行基, 其方法如下。 （1）第一种情况：若线性规划问题 max z =n j j j=1c x ∑ 1,1,2,...,0,1,2,...n ij j i j j a x b i m x j n =?==???≥=?∑ 从Pj ( j = 1 , 2 , ? , n )中一般能直接观察到存在一个初始可行基 121(,,...,)n B P P P 0 0?? ?0 1 0 ?== ? ?0 0 1?? （2）第二种情况：对所有约束条件是“ ≤”形式的不等式, 可以利用化为标准型的方法, 在每个约束条件的左端加上一个松弛变量。经过整理, 重新对 j x 及ij a ( i = 1 , 2 , ? , m ; j = 1 , 2 , ? , n )进行编号, 则可得下列方 程组 11,1111 22,1122,1112.........,,...,0 m m n n m m n n m m m m nn n n n x a x a x b x a x a x b x a x a x b x x x +++++++++=?? +++=?? ??+++=??≥? 显然得到一个m ×m 单位矩阵