为了探明木薯拔起时块根最大应力的影响规律和确定块根不被拔断的允许最大拔起力,论文采用 FEM(finite element method)和 SPH(smoothed particle hydrodynamics)的耦合方法及二次回归旋转设计方法,通过构建土壤 -块根 -茎秆系统的数值模拟计算模型,进行木薯块根拔起数值模拟试验,测定各因素组合条件下的块根最大拔起力和块根最大应力,建立块根最大应力与拔起速度、块根的大小、长短和生长深度及土壤的软硬程度的多因素耦合数学模型,研究了各影响因素及交互作用对块根最大应力的影响规律,同时,通过块根最大拔起力和块根最大应力的散点图,研究了块根最大拔起力和块根最大应力的相关性,确定了块根不被拔断的允许最大拔起力,且 2014年 12月底在广西武鸣县某木薯种植地,采用随机抽样方法,进行了木薯块根最大拔起力和块根拔断率的田间试验和统计分析,对块根不被拔断的允许最大拔起力进行了验证,最大拔起力小于 0.98 kN时,块根拔断率为 2.5%。结果表明,块根最大应力与各影响因素的多因素耦合数学模型的 F检验在 0.000 1水平上显著,精度较高,可用于块根最大应力的影响分析；块根最大应力随拔起速度的增大呈先增大,后减小的变化,随生长深度、块根长度和土壤硬度的增大而增大,随块根直径的增大而减小；块根最大应力与最大拔起力相关性不强,块根允许最大拔起力约为 0.98 kN。

In order to discover the influence law of maximum stress of cassava tuber, and obtain the allowable maximumlifting force without the breakage of cassava tuber during cassava tuber lifting, the numerical simulation model of soil,cassava tuber and cassava stem system was established by adopting explicit dynamics simulation software LS DYNA.Inthe numerical simulation model, large deformations and fractures of the soils close to the tubers occurred and smalldeformations occurred in most parts of the rest.Thus, the coupling method of FEM(finite element method) and SPH(smoothed particle hydrodynamics) was applied in the numerical simulation model.The SPH method was used in inner soillayer where large deformations occurred and the FEM method was used in the outer soil layer where small deformationsoccurred.The coupling between the inner and outer soil was realized by“nodes surface” in LS DYNA.And MAT_FHWA_SOIL was selected as soil material model in the numerical simulation model, because it takes account of theinfluence of moisture content, strain softening, strain rate effect, void ratio, and pore water pressure, and obeys themodified Mohr Coulomb yield criterion.Numerical simulation tests of cassava tuber lifting were carried out by using thenumerical simulation model and the quadratic regression rotation design method.There were 5 experiment factors and 2experiment indices in the quadratic regression rotation design tests.The experiment factors were the lifting velocity, thelength, dimension and growth depth of the tuber and the soil hardness, respectively.The experiment indices were themaximum lifting force and the maximum stress of tuber which were measured under different factor combinationconditions.According to the numerical simulation tests results, the multi factor coupling mathematic models between themaximum stress of the tuber and the lifting velocity, the length, dimension and growth depth of the tuber and the soilhardness were established by using statistical software SPSS.Based on the coupling mathematic models, the differentinfluencing factors and their interactions on the maximum stress of the tuber were investigated, the relationship curveswere drawn by using mathematics software MathCAD, and the corresponding influence laws were obtained by therelationship curves.And based on the numerical simulation tests results, the scatter diagram of the maximum lifting forceand the maximum stress of tuber was plotted by using mathematics software MathCAD, and the correlation between themwas studied by the scatter diagram.According to the scatter diagram, the allowable maximum lifting force was obtainedunder which there was no breakage of tuber.Moreover, the allowable maximum lifting force under which there was nobreakage of tuber was verified by cassava tuber lifting tests, which were carried out in the cassava planting field in WumingCounty, Guangxi Zhuang Autonomous Region in the end of December 2014 by adopting the method of random sampling.Meanwhile, the maximum lifting force and the breakage rate of tuber were analyzed by statistical method.When theallowable maximum lifting force was less than 0.98 kN, the breakage rate of tuber was 2.5%.The result showed that the multi factor coupling mathematic models between the maximum stress of the tuber and the lifting velocity, the dimensionand growing depth of the tuber and the soil hardness had high precision, because the F test of the multi factor couplingmathematic models was significant at 0.000 1 level.The multi factor coupling mathematic models could be used in theeffect analysis of the maximum stress of tuber.With the increasing of lifting velocity, the maximum stress of tuberincreased at first and then decreased.With the increasing of growing depth, tuber length and soil hardness, the maximumstress of tuber increased.But the maximum stress of tuber decreased when tuber′s diameter increased.There was little correlation between the maximum stress of tuber and the allowable maximum lifting force.The allowable maximum liftingforce was 0.98 kN.

S225.7+1

国家自然科学基金项目(51365005)；国家自然科学基金项目(51065003)；广西制造系统与制造技术重点实验室课题( 13-05109S01)；广西研究生教育创新计划资助项目(YCSZ2015040)