関数get_axisynmec_comp マニュアル

(The documentation of function get_axisynmec_comp)
Last Update: 2025/8/7


◆機能・用途(Purpose)

モーメントテンソルの3つの固有値のうちの2つが共通となるソースメカニズムについて、 モーメントテンソルの1つの成分の値を求める。
Calculate a moment tensor component of a source mechanism for which two of the three eigenvalues are equal to each other.

◆形式(Format)

#include <sourcemec.h>
inline double get_axisynmec_comp
(const char ∗component,const double M1,const double M2,
double theta,double phi)


◆引数(Arguments)

component 計算するモーメントテンソル成分名を表す文字列。 "Mxx", "Myy", "Mzz", "Mxy", "Myz", "Mzx"のうちのいずれか。
A string representing the moment tensor component name to calculate, which must be either "Mxx", "Myy", "Mzz", "Mxy", "Myz", or "Mzx".
M1 縮退していない方の固有値。
The eigen value which is not degenerated.
M2 縮退している方の固有値。
The eigen value which is degenerated.
theta \(M_1\)に対応する固有ベクトルと\(z\)軸とのなす角。 単位は°で、0以上90以下の値を指定する。
Angle between the eigenvector corresponding to \(M_1\) and the \(z\)-axis, in degrees greater than or equal to 0 and less than or equal to 90.
phi \(M_1\)に対応する固有ベクトルの\(xy\)平面への投影と \(x\)軸とのなす角を反時計回りに測った角度。 単位は°で、0以上360未満の値を指定する。
Angle between the eigenvector corresponding to \(M_1\) projected on \(x\)-\(y\) plane and the \(x\)-axis measured counterclockwise, in degrees greater than or equal to 0 and less than 360.


◆戻り値(Return value)

引数M1,M2,theta,phiの値をそれぞれ \(M_1\), \(M_2\), \(\theta\), \(\phi\)として、 引数componentの値に応じて以下のようになる。
The return value depends on the argument component and is given by the table below, where \(M_1\), \(M_2\), \(\theta\), and \(\phi\) are the values given by the arguments M1, M2, theta, and phi, respectively.
引数componentの値
Value of argument component 		戻り値
Return value
Mxx \((M_1-M_2)\sin^2\theta\cos^2\phi+M_2\)
Myy \((M_1-M_2)\sin^2\theta\sin^2\phi+M_2\)
Mzz \((M_1-M_2)\cos^2\theta+M_2\)
Mxy \((M_1-M_2)\sin^2\theta\sin\phi\cos\phi\)
Myz \((M_1-M_2)\sin\theta\cos\theta\sin\phi\)
Mzx \((M_1-M_2)\sin\theta\cos\theta\cos\phi\)


◆使用例(Example)

double Mxx=get_axisynmec_comp("Mxx",3.0e+12,1.0e+12,30.0,60.0);

この例では\(\theta=30^{\circ}\), \(\phi=60^{\circ}\), 地震モーメントが\(10^{12}\) N mの 開口クラックメカニズムの\(M_{xx}\)成分の計算が行われる。
In this example, an \(M_{xx}\) component of a tensile crack source mechanism with \(theta=30^{\circ}\), \(\phi=60^{\circ}\), and a seismic moment of \(10^{12}\) N m is calculated.