✍️ 朱邓达  •  🗓️ 2026-05-19 (创建时间)

计算动态解(面波解)

计算 动态面波解 的主要计算流程如下(目前仅支持 C 模块)。其中值得注意的是, 在经过 modsum 模块计算得到不同震源的面波格林函数后, 后续的合成位移/应力等流程完全和 计算动态解(全波解) 是一样的,共享计算模块。

关于面波计算的基本理论详见 姚振兴和谢小碧 (2026) ,这里不再做过多介绍。

        flowchart TB

    VV(["eigenv"])
    FN(["eigenfn"])
    MS(["modsum"])
    SS(["syn"])
    EE(["strain"])
    RR(["rotation"])
    TT(["stress"])

    V["Compute Dispersion Curves"]
    F["Compute Eigenfunctions, Energy Integrals,
       Group Velocity and Dispersion Sensitivity"]
    M["Compute Surface-wave Green's Functions (and its Spatial Derivatives)"]

    VV --> V
    V --> FN --> F
    V --> MS --> M
    M --> SS
    SS --> EE
    SS --> RR
    SS --> TT

    classDef cmdcls1 fill:#A1E3F9,stroke:#006BFF,stroke-width:2px,color:#333;
    classDef cmdcls2 fill:#f9f2d9,stroke:#e8d174,stroke-width:2px,color:#333;
    class VV,FN,MS cmdcls1
    class SS,EE,RR,TT cmdcls2