有限矩形源表面变形：Okada 1985 年断层错位位移、倾斜及应变解的计算（Matlab开发）

简介：
Okada [1985]模型分析了弹性半空间中因剪切和拉伸断层引起的表面变形的解析解。此模型被广泛应用以模拟由构造断层（如地震）或火山岩脉（如岩浆侵入）等局部扰动所引发的地表变形现象。基于给定矩形断层几何参数（长度、宽度、深度、走向及倾角）和三个分量位错幅度（前倾度、滑动量和张开量），该模型可用于计算自由表面的位移、倾斜程度及应变值。建议的Matlab脚本是基于冈田方程的文字转录，仅需将坐标系进行旋转即可应用——其中断层参数由相对于北的方向角定义，并采用rake、slip和opening（而非U1、U2、U3）参数，遵循Aki & Richards [1980]的定义标准。所有计算均基于断层质心为基准点展开。拉梅常数λ和μ则以泊松比ν替代（适用于各向同性介质，默认值取0.25），因为该公式与弹性参数无直接关联。这些方程则可应用于分析(x,y)平面上除倾角之外的所有输入参数及输出结果。