室内土工格栅拉拔试验物理模型边界及加载装置通常是刚性的，但土工格栅工程应用的边界条件往往不能与室内拉拔试验物理模型边界相对应。为了研究边界条件对拉拔试验结果的影响，在室内三轴和拉拔试验的基础上，运用三维离散元方法，对土工格栅加筋风积沙的拉拔试验展开数值模拟，研究加载顶面刚、柔性加载方式和模型前壁刚、柔性边界4种组合条件即刚性顶面刚性前壁（RR）、柔性顶面刚性前壁（FR）、刚性顶面柔性前壁（RF）、柔性顶面柔性前壁（FF）对筋土界面宏细观特性的影响。分析了不同组合条件下拉拔力与拉拔位移的关系、界面剪切强度、试样内部应力链、孔隙率分布和颗粒旋转规律，分析了刚性加载板和柔性边界颗粒位移、FF组合下土工格栅的变形以及拉拔过程中剪切带的演化规律。研究表明，模型前壁的刚、柔性边界对拉拔力和拉拔位移曲线模式以及界面摩擦角的影响很大。刚性前壁边界条件下，拉拔力和拉拔位移曲线呈加工软化型；柔性前壁边界条件下，拉拔力和拉拔位移曲线近似呈双折线型。当模型前壁由刚性变为柔性时，界面摩擦角减小了7º～8º。当法向应力小于90 kPa时，建议FR组合、RF组合和FF组合的拉拔力峰值折减系数分别取0.9、0.6和0.7。4种组合在拉拔试验过程中试样体积增加，呈剪胀特性。RR和FR组合剪切带厚度分布为颗粒平均粒径的5.38倍和10.79倍；RF与FF组合剪切带的分布范围更广。该研究成果有助于深入揭示刚、柔加载方式以及刚、柔前壁边界下土工合成材料和风积沙的相互作用机制