在剛體運(yùn)動學(xué)、飛機(jī)飛行、衛(wèi)星姿態(tài)等領(lǐng)域，歐拉角是一個非常重要的概念和控制參數(shù)。

通俗講，一個直角坐標(biāo)系Sys-1，繞自己的X、Y、Z坐標(biāo)軸依次轉(zhuǎn)動一個角度，變?yōu)樽鴺?biāo)系Sys-2，則轉(zhuǎn)動的角度為歐拉角。歐拉角的值與轉(zhuǎn)動順序相關(guān)，例如從Sys-1轉(zhuǎn)到Sys-2，轉(zhuǎn)序X→Y→Z與轉(zhuǎn)序Z→X→Y的歐拉角通常不同。因此，講歐拉角時必須明確轉(zhuǎn)序。除了繞3個軸轉(zhuǎn)的歐拉角，還有繞2個軸轉(zhuǎn)的歐拉角，例如轉(zhuǎn)序X→Y→X。

除了運(yùn)動學(xué)，精密設(shè)備的指向變化，有時也用歐拉角表示。例如，某精密指向、定位設(shè)備安裝面為X-Y平面，指向?yàn)閆軸。由于外載荷（力、熱等）作用，引起該設(shè)備安裝面、指向發(fā)生變化，這時，可以用指定轉(zhuǎn)序下的歐拉角表示設(shè)備指向狀態(tài)的變化量。此處指定轉(zhuǎn)序下的歐拉角稱為絕對歐拉角。

對于有一定角度關(guān)系的兩臺精密設(shè)備，當(dāng)關(guān)注兩臺設(shè)備相對指向變化時，引入相對歐拉角描述。

下面以簡單例子介紹平面變形、指向歐拉角的計(jì)算，包括絕對歐拉角、相對歐拉角。

1 簡化模型

下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側(cè)3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。

2 計(jì)算要求

計(jì)算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。

3 數(shù)據(jù)處理

使用平面節(jié)點(diǎn)坐標(biāo)、位移數(shù)據(jù)計(jì)算平面變形歐拉角。可以使用Python腳本輸出平面節(jié)點(diǎn)編號、節(jié)點(diǎn)坐標(biāo)（X、Y、Z）、節(jié)點(diǎn)位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數(shù)據(jù)文件，打開的文本文件中7列數(shù)據(jù)為節(jié)點(diǎn)編號、坐標(biāo)、位移。

三個平面10個工況的節(jié)點(diǎn)數(shù)據(jù)文件如下圖所示。每個文件中包含一個工況一個平面的節(jié)點(diǎn)編號、坐標(biāo)、位移數(shù)據(jù)。

4 絕對歐拉角計(jì)算

使用PyQt+Python開發(fā)了一個簡單的小軟件，計(jì)算絕對歐拉角、相對歐拉角。

首先計(jì)算各平面的絕對歐拉角。

計(jì)算平面1的10個工況的絕對歐拉角。

平面1變形的絕對歐拉角計(jì)算結(jié)果如下圖所示。

伴隨絕對歐拉角計(jì)算結(jié)果，軟件同時寫出了平面變形前后的坐標(biāo)系數(shù)據(jù)，如下圖。每行18個數(shù)據(jù)，每3個數(shù)據(jù)為一個坐標(biāo)軸向量，變形前后2個坐標(biāo)系，6個坐標(biāo)軸，18個數(shù)據(jù)。

5 相對歐拉角計(jì)算

利用計(jì)算絕對歐拉角時得到的坐標(biāo)系文件，計(jì)算平面變形相對歐拉角，如下圖所示，計(jì)算平面2相當(dāng)于平面1、平面3相對與平面1的相對歐拉角。

計(jì)算結(jié)果如下圖所示。

6 小結(jié)

上述軟件用的算法申請了發(fā)明專利，軟件申請了軟著。CAE工程師，也可以自制軟件工具，解決重復(fù)性、復(fù)雜性數(shù)據(jù)處理等工作痛點(diǎn)。

2021年8月24日于西昌衛(wèi)星發(fā)射中心