基於Matlab動力總成扭矩軸(TRA)的計算方法

【聲明①】：本公眾號轉載的文章均會在標題欄標註「轉載」；本公眾號內的原創文章則直接署名作者；其他公眾號需要轉載本公眾號的文章的，請提前私信本公眾號，謝謝！【聲明②】：該文部分信息來源於汽車之家、百度文庫等網絡，如有侵權請聯繫刪除。

「基於Matlab動力總成扭矩軸(TRA)的計算方法」

扭矩軸（TRA，Torque Roll Axis）是一根無約束的三維剛體旋轉軸。根據歐拉定理，當自由剛體收到某一特定方向的轉矩激勵時會繞著某一根軸轉動，這個軸就是扭矩軸。在動力總成懸置系統中，當施加繞曲軸轉動的扭矩時，系統不沿著任何一根主慣性軸轉動，而是繞著扭矩軸轉動。

經過小編與同事的溝通交流以及相關資料的查閱（實際上TRA的計算並不複雜，但是要理解為什麼這樣求解卻涉及到坐標軸轉換、運動方程、矩陣求解以及主慣性軸和扭矩軸區分等，奈何小編的基礎理論太差，也著實懵逼了一場~~~），基於matlab軟體給出了如下扭矩軸的計算代碼：

一，橫置動力總成的扭矩軸（MQB_TRA）

M代碼：

function [A,B,C]=TRA_MQB(Ixx,Iyy,Izz,Ixy,Ixz,Iyz)

M=[Ixx,-Ixy,-Ixz;-Ixy,Iyy,-Iyz;-Ixz,-Iyz,Izz];

N=inv(M);

p=sqrt(N(1,2)^2+N(2,2)^2+N(3,2)^2);

A=N(1,2)/p;

B=N(2,2)/p;

C=N(3,2)/p;

Cos_value=[A,B,C]

Angle_value=[acosd(A),acosd(B),acosd(C)]

命令行代碼：

Ixx=00;Iyy=00;Izz=00;Ixy=00;Ixz=00;Iyz=00;

[A,B,C]=TRA_MQB (Ixx,Iyy,Izz,Ixy,Ixz,Iyz);

二，縱置動力總成的扭矩軸（MLB_TRA）

M代碼：

function [A,B,C]=TRA_MLB(Ixx,Iyy,Izz,Ixy,Ixz,Iyz)

M=[Ixx,-Ixy,-Ixz;-Ixy,Iyy,-Iyz;-Ixz,-Iyz,Izz];

N=inv(M);

p=sqrt(N(1,1)^2+N(2,1)^2+N(3,1)^2);

A=N(1,1)/p;

B=N(2,1)/p;

C=N(3,1)/p;

Cos_value=[A,B,C]

Angle_value=[acosd(A),acosd(B),acosd(C)]

命令行代碼：

Ixx=00;Iyy=00;Izz=00;Ixy=00;Ixz=00;Iyz=00;

[A,B,C]=TRA_MLB (Ixx,Iyy,Izz,Ixy,Ixz,Iyz);

參考文獻：

郭榮，章桐.汽車動力總成懸置[M].上海：同濟大學出版社，2013.8。

     

掃碼關注