利用MATLAB数值计算对迭代法求解二阶频率精确度的分析结构力学小论文.docVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 利用MATLAB数值计算对迭代法求解高阶频率精确度的分析　　 XX 1103210224　 摘要：振型的求解对于了解结构的动力响应性质有着重要的意义，一般的解决方法有静力法、能量法。在实际应用中，迭代法多次逼近求解频率的方法具有形式简单，程序化的优势，因而得到了广泛的应用。迭代法求解一阶频率收敛快，精确度高，但是在计算二阶及更高阶频率时需要将之前求解的低阶频率作为初始条件代入，因此高阶频率的求解精度可能会受到前一阶频率精确度的影响，其影响可能造成收敛速度的减慢或精确度降低。同时，计算的失误也可能造成高阶频率的误差，这些误差和影响的程度如何，教材并未深入讨论，本文对此进行数值计算上的分析，得出结论：迭代法在计算高阶频率和振型时具有良好的收敛性，受以上可能因素的干扰很小，具有很好地容错性，可以满足一般的精度要求。　　 关键词：迭代法，二阶频率，MATLAB数值分析　　 计算原理　　　 由迭代法计算一阶频率的证明过程可知，不论初始迭代向量如何选取，只要迭代次数足够大，迭代结果总会收敛于一阶振型。同理，如果在假设的振型迭代向量中将一阶振型分量消除，那么经过多次迭代后，结果将收敛于二阶振型，对于更高阶振型方法相同。　　 将任意设定的振型迭代向量按振型分解：　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1)　　 利用正交性，在上式两边同时左乘，得：　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2)　 令一阶滤型矩阵为：　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3)　 得到二阶初始迭代向量：　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4)　 这样，即可利用一阶滤型矩阵在初始迭代向量中过滤掉一阶振型，保证迭代的结果收敛于二阶振型。为了与一阶振型计算方法统一，将滤型的步骤与动力矩阵合并，将滤型矩阵与动力矩阵相乘，得到二阶动力矩阵，其余计算方法与一阶振型计算完全一致。　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5)　 MATLAB 数值分析　　　 此处选取教材295页，例9-22作为分析算例。　 　　 图 SEQ 图 \* ARABIC 1 算例示意图　 首先计算一阶频率和振型，不考虑轴向变形条件下，体系有三个自由度，由静力分析图乘得到柔度矩阵和质量矩阵如下：　　　 ，　 由此得到动力矩阵为：　 　　　 设初始迭代向量为：　　 　　 为了控制迭代计算次数，定义精确度为每次迭代后归一化振型向量的第二个元素相差的比值。即：　 MACROBUTTON MTPlaceRef \* MERGEFORMAT SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 6)　 可以看出，当精确度值足够小，就可以认为迭代计算收敛。本文中，作为探究为了保证足够收敛，将精确度定为（实际工程计算不需要这样高的精确度）。　 利用MATLAB程序计算，经过四次迭代，达到精确度要求，得到：　 　　　 　　　 为了便于程序分析，将动力矩阵前的常数定为1，即：　　　 =1　　 现计算二阶振型和频率，精确度和初始振型向量设置与之前计算一致，经过四次迭代，计算得：　　 　 　　　 二阶振型振型和频率为：　 　　　 　 显然，。以下为了研究一阶误差对二阶频率和振型计算带来的影响，将此处的和作为准确值，以供对比。将一阶振型计算的误差表