计算边界条件弹性地基梁结构表面在x和y方向上的弯矩值，利用Matlab数值拟合软件制作了三维表面内力分布的等值线图

对于边界条件弹性地基梁结构，在x和y方向上的弯矩值可以通过以下公式计算：

Mx = -EI(d^2w/dx^2)
My = -EI(d^2w/dy^2)

其中，Mx和My分别表示在x和y方向上的弯矩值，E表示弹性模量，I表示截面惯性矩，w表示梁的挠度，d^2w/dx^2和d^2w/dy^2表示挠度在x和y方向上的二阶偏导数。

为了制作三维表面内力分布的等值线图，需要进行以下步骤：

1. 对梁的挠度进行离散化，将其分割为若干个小段。

2. 计算每个小段在x和y方向上的弯矩值Mx和My。

3. 将每个小段的弯矩值Mx和My转化为等效应力值，利用等效应力值制作出三维表面内力分布的等值线图。

4. 利用Matlab数值拟合软件对等值线图进行拟合，得到一个函数关系式，用于预测任意一点的内力分布情况。

需要注意的是，在计算每个小段的弯矩值时，需要考虑边界条件和梁的受力情况。具体的计算方法可以参考结构力学相关的教材和资料。