【正文】 % Y_end=T*Y_ini。F。% F=1e6*L/E1/h1% Y_end(5,1)=F。% Y_ini=zeros(6,1)。Y_end=T2*Y_part_right。Y_part_left=T1*Y_ini。 Y_ini=zeros(6,1)。 T_select=T_tol([2,6],[3,5])。 T_tol=T2*T1。T12=zeros(6,1)。%===========================================================% %簡支梁T1=exppm(G*,20)。% Y_part_right=Y_part_left+T12。% Y_ini([3,5],1)=inv(T_select)*TF_select。% TF_select=TF([2,6],1)。% TF=T2*T12。% T12(5,1)=1000。% T2=exppm(G*,20)。%======================================= G=Matrix_G(h1,L,E1,rou1,omg2)。 rou1=7800。h1=。%基梁靜力學驗證clear allclcL=。) (226)當N次循環(huán)結束后有 (227) T=exp(A)能獲得高精度計算結果的根本原因是： 數值計算的相對誤差不隨遞推過程的進行而擴散。 鐘氏精細算法中其要點是利用加法定理，取，將矩陣縮小后，保證用泰勒級數展開計算的可靠性。指數矩陣的精細算法在結構動力、優(yōu)化控制等問題中，通過變換都可以將運動（控制）微分方程寫成狀態(tài)向量形式 (220)式中，為階狀態(tài)向量；是階常數矩陣；為階載荷向量（或控制微量）。（6）驗證基層控制方程是否正確 0 1 圖7 簡支梁邊界條件分別如下： 以如圖7的簡支梁為例進行驗證，由于，故，則式可以化簡成如下的齊次方程： （其中為傳遞矩陣） 令。而式與 式左邊完全相同，由此可知： 將 式代入和，將 式代入和可得： 式和中含有未知作用力和，可利用式消去如下： 此時方程中出現(xiàn)了新的整合狀態(tài)變量，需將基梁和約束梁控制方程中出現(xiàn)的式和出現(xiàn)的式進行整合，以匹配出整合狀態(tài)向量。 基梁變形前位置 中面層 基梁變形后位置 圖6各梁的縱向位移 基梁變形前位置 基梁變形前位置 圖7各梁的橫向位移則可得：則 中面層 圖6 層合梁變形協(xié)調關系 由圖6可見，黏彈層的剪應變?yōu)椋航Y合式可得： 黏由于彈層僅計其剪切變形，故其剪切應力可表示為 式中，為諧激勵下的剪應變幅值,為粘彈層的復剪切模量。故有由于很小，所以有 由材料力學可知（其中）代入可得如圖3，又因為，整理為形式則可求得 則諧激勵下（各物理量均可表示為幅值和諧激勵因子的乘積的形式，下面方程中有上標~的各物理量均為真實物理量的幅值），其縱向振動方程和橫向振動方程化成如下形式： 對于各向同性量，其內力位移關系為： 為了處理方便，引入無量綱變量，對基梁的平衡方程和內力位移關系式進行處理，可得基梁的一階狀態(tài)向量微分控制方程組如下： 寫成矩陣的形式如下： 其中。所以，截面轉角就是y軸與繞曲線法線的夾角。桿件在軸線方向的伸長為 圖4將除以得桿件軸線方向的線應變