急求流固耦合单元阵的MATLAB程序

zhy

单元流固耦合矩阵公式为： [Re]n*3n=~~[Np]n*1{n}[Nq]3*3ndS
[Np]为流体中声压的插值函数,大小为n*1；   [Np]  为结构的位移插值函数，小为3*3n
{n}为交界面出的方向向量，设为{0，0，1};~~为两单元交界面处的二重积分
n为单元节点数，流体和结构单元均采用六面体八节点单元，即n=8。
S为两单元的交界面，最后转化到自然坐标系下应为边长为2的正方形。为一平面，z为常数。
本人认为此耦合矩阵程序将涉及两个单元（流体+结构）以及12个节点（交界面处公用4个节点），节点坐标值已知，怎样在自然坐标系下编程？
由于公式复制不上去，给大家带来不便 ，请谅解
请各位指教，其他语言也可或推荐相关书籍，不胜感激！