摘要：从弹性力学基础理论出发,采用刚度矩阵法,推导了应用于直角坐标系下的三维多层弹性层状体系静力学数值解法。引入二维傅里叶变换及高斯积分求解法,基于MATLAB数学软件平台编制计算程序,实现三维多层弹性层状体系理论计算方法的数值求解。针对典型有砟轨道轨下基础结构,采用提出的计算方法和编制的相应计算程序对其进行静力学分析,并将所获得的计算结果与采用通用有限元程序ABAQUS的计算结果进行对比。分析结果表明：采用提出的计算方法和通用有限元计算方法获得有砟轨道轨下基础最大竖向位移分别为1.50、1.95 mm,最大竖向应力分别为0.34、0.21 MPa,计算结果较为接近,计算反映出来的各状态分量变换规律基本一致,提出的计算方法及其相应计算程序可应用于多层弹性层状体系的静力学计算。

关键词：路基工程 弹性层状体系 刚度矩阵法 数值解法 

单位：西南交通大学土木工程学院; 四川成都610031; 西南交通大学道路工程四川省重点实验室; 四川成都610031