Ansys MAPDLのテスト

Geminiが出した以下のスクリプトでテストしてみた。(Goldakモデルによる熱源の移動)
一応動いてはいるらしい。何が起こっているのかはこれから見てみる。

/CLEAR
/PREP7

! --- 1. 材料特性 ---
MAT_ID = 1
MP, KXX,  MAT_ID, 15
MP, C,    MAT_ID, 500
MP, DENS, MAT_ID, 7900
ET, 1,    SOLID70

! --- 2. ジオメトリ ---
W_L = 0.100  $ W_H = 0.050 $ W_T = 0.010
BLOCK, 0, W_L, 0, W_H, 0, W_T
ESIZE, 0.004
VMESH, ALL

! --- 3. Goldakパラメータ ---
Q_EFF = 1125    ! 1500 * 0.75
AW = 0.005  $ BD = 0.004  $ CF = 0.006  $ CR = 0.018
FF = 0.6    $ FR = 1.4
PI = 3.14159265
V_W = 0.005  $ T_INC = 0.5  $ N_STP = 20

! ループ外で定数を単純化（カッコを排除）
SQ3 = 1.73205
SQPI = 1.77245
DENOM = AW * BD
DENOM = DENOM * PI
DENOM = DENOM * SQPI

C_F = 6 * SQ3
C_F = C_F * FF
C_F = C_F * Q_EFF
C_F = C_F / DENOM
C_F = C_F / CF

C_R = 6 * SQ3
C_R = C_R * FR
C_R = C_R * Q_EFF
C_R = C_R / DENOM
C_R = C_R / CR

/SOLU
ANTYPE, TRANS
TUNIF, 25
KBC, 1

! --- 5. 移動熱源ループ ---
*DO, i, 1, N_STP
    CUR_T = i * T_INC
    TIME, CUR_T
    XC = V_W * CUR_T
    YC = 0.025
    ZC = 0.010

    BFDELE, ALL, HGEN
    SFA, ALL, 1, CONV, 20, 25

    ESEL, S, LOC, X, XC-CR, XC+CF
    ESEL, R, LOC, Y, YC-AW, YC+AW
    ESEL, R, LOC, Z, ZC-BD, ZC+0.001

    *GET, E_CNT, ELEM, 0, COUNT
    *IF, E_CNT, GT, 0, THEN
        E_ID = 0
        *DO, j, 1, E_CNT
            E_ID = ELNEXT(E_ID)
            EX = CENTRX(E_ID)
            EY = CENTRY(E_ID)
            EZ = CENTRZ(E_ID)

            DX = EX - XC
            DY = EY - YC
            DZ = EZ - ZC

            ! カッコを一切使わない逐次計算
            Y_P = DY / AW
            Y_P = Y_P * Y_P
            
            Z_P = DZ / BD
            Z_P = Z_P * Z_P
            
            SUM_YZ = Y_P + Z_P

            *IF, DX, GE, 0, THEN
                X_P = DX / CF
                X_P = X_P * X_P
                DIST = X_P + SUM_YZ
                DIST = DIST * -3
                H_VAL = C_F * EXP(DIST)
            *ELSE
                X_P = DX / CR
                X_P = X_P * X_P
                DIST = X_P + SUM_YZ
                DIST = DIST * -3
                H_VAL = C_R * EXP(DIST)
            *ENDIF

            BFE, E_ID, HGEN, 1, H_VAL
        *ENDDO
    *ENDIF
    ALLSEL
    SOLVE
*ENDDO

アニメーションの方法
PlotCtrls>Animate>Deformed Results

Startをクリック

実行結果