用語集 — 記号・変数名の対応表
理論編(02/04)の数式記号と、実装編(03/05)のVBA変数名の対応をまとめます。コードを読んでいて変数名の意味に迷ったら、まずここを確認してください。
空力解析(02_空力理論編.md/03_空力実装編.md)
| 数式記号 |
VBA変数名 |
意味 |
単位 |
| θ |
geo.Th(i) |
コサイン分布の角度パラメータ(0〜π) |
rad |
| y |
geo.Y(i) |
スパン方向座標(翼根=0、正が翼端方向) |
m |
| b |
geo.BSpan |
全スパン(翼端〜翼端) |
m |
| c(y) |
geo.Cord(i) |
局所翼弦長 |
m |
| Aₙ, Bₙ |
Bvec(q, 1) |
循環のフーリエ係数(奇数次のみ) |
— |
| n(フーリエ次数) |
XX(q)(= 2q-1) |
フーリエ級数の次数(奇数) |
— |
| K |
geo.K(= n\2) |
フーリエ係数の個数、制御点の個数 |
— |
| μ(θ) |
(minfP * cordP / (4*BSpan)として展開) |
モノプレーン方程式の無次元パラメータ |
— |
| a₀ |
minfJohen(i) |
局所揚力傾斜(2次元翼型データ由来) |
1/rad |
| −α_L0 |
zeroAJohen(i) |
ゼロ揚力角の符号反転値(2次元翼型データ由来)。注意: 標準的な意味でのα_L0そのものではなく、その符号を反転した値がここに格納されている(キャンバー翼型ではα_L0は負なので、zeroAJohenは正の値になる) |
deg |
| α(幾何迎角) |
geo.aoa(i) |
取付角+ねじり下げ角 |
deg |
| (α − α_L0) |
aba(i)(= geo.aoa(i) + zeroAJohen(i)) |
絶対迎角(ラジアン)。zeroAJohenが−α_L0を格納しているため、コード上は「足し算」になる |
rad |
| α_i(誘導迎角) |
ai(i) |
吹き下ろしによる迎角の減少分 |
deg |
| αe(有効迎角) |
Ae(i) |
幾何迎角 − 誘導迎角 |
deg |
| w/V(吹き下ろし比) |
bn(i), W(i) |
無次元吹き下ろし速度、WはV1を掛けた実速度 |
— , m/s |
| CL |
result.CL |
全機揚力係数 |
— |
| δ(スパン損失係数) |
del_ |
楕円分布からのズレ(高次フーリエ項の寄与) |
— |
| e(スパン効率) |
result.ew(= 1/(1+δ)) |
誘導抗力の効率、1が理想 |
— |
| Cl(y) |
Cl_dist(i) |
局所揚力係数分布 |
— |
| Cd(y) |
CdDisJohen(i) |
局所プロファイル抗力係数分布 |
— |
| Cm(y) |
CmDisJohen(i) |
局所モーメント係数分布(1/4翼弦基準) |
— |
| CDi(y) |
CdInd(i) |
局所誘導抗力係数分布 |
— |
| Cp(y) |
Cp(i) |
局所圧力中心位置(翼弦割合、0〜1)。「空力中心(aerodynamic center)」とは別概念なので注意 |
— |
| Re |
reArr(i) |
局所Reynolds数(局所翼弦・V1・動粘性係数から) |
— |
| V1 |
result.V1 |
つり合い飛行速度(トリム後) |
m/s |
| Weight |
Weight |
機体重量(構造+パイロット) |
kgf |
| Wtot |
Wtot |
現在のV1での総揚力(kgf換算) |
kgf |
| Cmg |
result.Cmg |
主翼全体の縦揺れモーメント係数 |
— |
| ReqPower |
result.ReqPower |
必要空力パワー(=全抗力×V1) |
W |
構造解析(04_構造理論編.md/05_構造実装編.md)
| 数式記号 |
VBA変数名 |
意味 |
単位 |
| Ex, Ey |
yex, yey |
材料主軸のヤング率(繊維方向・直角方向) |
Pa |
| Gxy |
yes_ |
材料主軸の面内せん断弾性率 |
Pa |
| νxy |
rmux |
材料主軸のポアソン比 |
— |
| U1〜U5 |
u1〜u5 |
Q行列の材料不変量(角度に依存しない) |
Pa |
| Q̄11, Q̄12, Q̄16, Q̄22, Q̄26, Q̄66 |
q11(j)〜q66(j) |
層jの、繊維角度で回転した縮約剛性 |
Pa |
| θ(積層角) |
thetas(j) |
層jの繊維角度 |
deg |
| A11〜A66 |
a11〜a66 |
面内剛性行列(厚み比重み付き平均) |
Pa |
| Ex_eff, Gxy_eff |
exEff, gxyEff(Mod11) / exEff, eyEff(Mod17) |
積層板全体の等価物性(A行列の逆行列成分由来) |
Pa |
| EA |
sec.EA(i) |
軸剛性 |
N |
| GK |
sec.GK(i) |
ねじり剛性 |
N・m² |
| EIy, EIz |
sec.EIy(i), sec.EIz(i) |
曲げ剛性(v-θz平面、w-θy平面) |
N・m² |
| Iy, Iz |
sec.Iy(i), sec.Iz(i) |
断面二次モーメント(応力計算用の物理値) |
m⁴ |
| Q16(連成剛性) |
sec.Q16(i) |
曲げ-ねじり連成剛性 |
N・m² |
| z(層境界半径) |
z2(j) |
プライ境界の半径方向位置(内径基準) |
m |
| σ1, σ2, σ6 |
sig1Plate, sig2Plate, sig6Plate |
板座標系での応力(1=軸方向、2=周方向、6=面内せん断) |
Pa |
| σx, σy, σs |
sigmaX, sigmaY, sigmaS |
材料主軸へ回転した後の応力 |
Pa |
| stx, stxd, sty, styd, sts |
stx(j)〜sts(j) |
強度allowable(繊維/直角×引張/圧縮、面内せん断) |
Pa |
| R(強度比) |
TsaiWuRatioの戻り値、ratio |
Tsai-Wu複合則の破壊安全率 |
— |
| ε1, ε6 |
eps1Beam, eps6Beam |
梁理論から求めた材料主軸方向のひずみ |
— |
| u, v, w, θx, θy, θz |
Ux/Uy/Uz/THx/THy/THz(clsSparDeformation) |
節点の変位・回転(6自由度) |
m, m, m, rad, rad, rad |
| BMD |
deformation.BMD(i) |
曲げモーメント分布(v-θz平面) |
N・m |
| SFD |
deformation.SFD(i) |
せん断力分布(v-θz平面) |
N |
| Fx, Fy, Fz |
deformation.Fx/Fy/Fz(i) |
要素の断面内力(軸力・せん断力2方向) |
N |
| MxInt, MyInt, MzInt |
deformation.MxInt/MyInt/MzInt(i) |
要素の断面内力(ねじり・曲げモーメント2方向) |
N・m |
| zakutu1〜4 |
zakutu1〜zakutu4 |
座屈応力(4方式) |
Pa |
| b1〜b4 |
b1〜b4(Mod16) |
要素剛性行列の6×6ブロック(節点i-i, i-i+1, i+1-i, i+1-i+1) |
— |
共通の添字・接頭辞の慣習
- 変数名の末尾
Arr・(): 要素または節点ごとの配列であることを示す
〜Johen(ジョーヘン=如変、「ブレンド後」の意味の造語): 2つの翼型/材料を補間比でブレンドした後の値〜Fem: 空力側のコサイン分布グリッドから、構造側の等間隔グリッドへスプライン変換した後の値(Mod15_桁荷重)sec: clsSparSection(断面剛性)のインスタンス変数名として頻出geo: clsWingGeometry(翼形状)のインスタンス変数名として頻出