详情

原文传递一种基于剪应变的土体破裂面确定方法和系统

专利名称：	一种基于剪应变的土体破裂面确定方法和系统
摘要：	本发明涉及一种基于剪应变的土体破裂面确定方法和系统，涉及土壤勘测领域。包括以下步骤：S1：获取一个装有土体作为模型箱，并对所述土体进行变形，将所述土体的变形的时间划分为至少2个时间段，观测所述土体的变形并得到与每个所述时间段对应的土体示踪点的位移场；S2：获取一个所述时间段作为计算时段，并在所述计算时段的所述模型箱内构建三节点应变计算单元；S3：根据所述土体示踪点的位移场得到示踪点位移场插值，根据所述示踪点位移场插值、所述三节点应变计算单元计算所述三节点应变计算单元的节点位移坐标。本方案解决了如何准确的确定土体破裂面的形态的技术问题，适用于确定土体破裂面的形态。
专利类型：	发明专利
国家地区组织代码：	湖南;43
申请人：	湖南理工学院
发明人：	刘晓平;姜豪;徐川;刘王苗;刘晓红;刘正夫;张家强;曾永庆;陆仲阳;张政;张建伟;彭锟
专利状态：	有效
申请日期：	2019-07-04T00:00:00+0800
发布日期：	2019-09-17T00:00:00+0800
申请号：	CN201910598198.0
公开号：	CN110243666A
代理机构：	重庆市信立达专利代理事务所(普通合伙)
代理人：	陈炳萍
分类号：	G01N3/00(2006.01);G;G01;G01N;G01N3
申请人地址：	414006 湖南省岳阳市岳阳楼区奇家岭学院路439号
主权项：	1.一种基于剪应变的土体破裂面确定方法，其特征在于，包括以下步骤： S1：获取一个装有土体作为模型箱，并对所述土体进行变形，将所述土体的变形的时间划分为至少3个时间段，观测所述土体的变形并得到与每个所述时间段对应的土体示踪点的位移场； S2：获取一个所述时间段作为计算时段，并在所述计算时段的所述模型箱内构建三节点应变计算单元，所述三节点应变计算单元包括第一节点、第二节点和第三节点； S3：根据所述土体示踪点的位移场得到示踪点位移场插值，根据所述示踪点位移场插值、所述第一节点、所述第二节点和所述第三节点计算所述三节点应变计算单元的节点位移坐标，所述节点位移坐标包括与所述第一节点对应的第一位移坐标，与所述第二节点对应的第二位移坐标，与所述第三节点对应的第三位移坐标； S4：根据所述第一位移坐标、所述第二位移坐标和所述第三位移坐标计算所述三节点应变计算单元的正应变和剪应变； S5：根据所述正应变和所述剪应变得到主应变和最大剪应变； S6：获取所述计算时段前一个时间段的最大剪应变作为原始剪应变，根据所述最大剪应变和所述原始剪应变得到最大剪应变增量； S7：重复步骤S2-步骤S6直至遍历所有所述时间段； S8：获取一个所述时间段作为目标时段，将所述目标时段内的所述模型箱均匀分成至少2个竖向长条区域； S9：获取所述竖向长条内的所有所述最大剪应变增量中的最大值，并将所述最大值的位置坐标记录为最大坐标； S10：根据最小二乘法和所述最大坐标拟合成曲线，并记录为所述目标时段内土体塑性区的几何形态； S11：重复步骤S8-步骤S10直至遍历所有所述时间段； S12：按照所述时间段的先后顺序，记录第一个贯通所述模型箱的所述几何形态为土体破裂面几何形态。 2.根据权利要求1所述的基于剪应变的土体破裂面确定方法，其特征在于：步骤S1具体为，获取一个装有土体作为模型箱，并对所述土体进行变形，将所述土体的变形的时间划分为至少2个时间段，基于PIV技术观测所述土体的变形并得到与每个所述时间段对应的土体示踪点的位移场。 3.根据权利要求1所述的基于剪应变的土体破裂面确定方法，其特征在于，步骤S4具体为： S41：获取所述三节点应变单元的三个节点的坐标，所述三个节点的坐标分别为第一节点坐标(xi,yi)、第二节点坐标(xj,yj)和第三节点坐标(xm,ym)，并根据下列公式计算出与所述第一节点坐标对应的第一节点位移坐标(ui,vi)、与所述第二节点坐标对应的第二节点位移坐标(ui,vi)和与所述第三节点坐标对应的第三节点位移坐标(ui,vi): 其中，A为3个节点所围成的面积： 4.根据权利要求3所述的基于剪应变的土体破裂面确定方法，其特征在于，步骤S5具体为：根据下列公式、所述第一位移坐标、所述第二位移坐标和所述第三位移坐标计算所述三节点应变计算单元的正应变εx,εy和剪应变γmax：其中表示u对x求偏导数，表示v对y求偏导数，表示u对y求偏导数，表示v对x求偏导数。 5.根据权利要求4所述的基于剪应变的土体破裂面确定方法，其特征在于，步骤S6具体为：根据下列公式、所述正应变εx,εy和所述剪应变γxy得到主应变ε1,ε2和最大剪应变γmax： 6.一种基于剪应变的土体破裂面确定系统，其特征在于：包括观测模块、计算模块和几何形态模块，所述观测模块用于获取一个装有土体作为模型箱，并对所述土体进行变形，将所述土体的变形的时间划分为至少2个时间段，观测所述土体的变形并得到与每个所述时间段对应的土体示踪点的位移场；所述计算模块用于获取一个所述时间段作为计算时段，并在所述计算时段的所述模型箱内构建三节点应变计算单元，所述三节点应变计算单元包括第一节点、第二节点和第三节点；所述计算模块还用于根据所述土体示踪点的位移场得到示踪点位移场插值，根据所述示踪点位移场插值、所述第一节点、所述第二节点和所述第三节点计算所述三节点应变计算单元的节点位移坐标，所述节点位移坐标包括与所述第一节点对应的第一位移坐标，与所述第二节点对应的第二位移坐标，与所述第三节点对应的第三位移坐标；所述计算模块还用于根据所述第一位移坐标、所述第二位移坐标和所述第三位移坐标计算所述三节点应变计算单元的正应变和剪应变；所述计算模块还用于根据所述正应变和所述剪应变得到主应变和最大剪应变；所述计算模块还用于获取所述计算时段前一个时间段的最大剪应变作为原始剪应变，根据所述最大剪应变和所述原始剪应变得到最大剪应变增量；所述计算模块还用于重复获取所述时间段直至遍历所有所述时间段并得到与每个所述时间段对应的最大剪应变增量；所述几何形态模块用于获取一个所述时间段作为目标时段，将所述目标时段内的所述模型箱均匀分成至少2个竖向长条区域；所述几何形态模块还用于获取所述竖向长条内的所有所述最大剪应变增量中的最大值，并将所述最大值的位置坐标记录为最大坐标；所述几何形态模块还用于根据最小二乘法和所述最大坐标拟合成曲线，并记录为所述目标时段内土体塑性区的几何形态；所述几何形态模块还用于重复获取所述时间段直至遍历所有所述时间段并得到与每个所述时间段对应的土体塑性区的几何形态；所述几何形态模块还用于按照所述时间段的先后顺序，记录第一个贯通所述模型箱的所述几何形态为土体破裂面几何形态。 7.根据权利要求6所述的基于剪应变的土体破裂面确定系统，其特征在于：所述观测模块用于获取一个装有土体作为模型箱，并对所述土体进行变形，将所述土体的变形的时间划分为至少2个时间段，基于PIV技术观测所述土体的变形并得到与每个所述时间段对应的土体示踪点的位移场。 8.根据权利要求6所述的基于剪应变的土体破裂面确定系统，其特征在于：所述计算模块用于获取所述三节点应变单元的三个节点的坐标，所述三个节点的坐标分别为第一节点坐标(xi,yi)、第二节点坐标(xj,yj)和第三节点坐标(xm,ym)，并根据下列公式计算出与所述第一节点坐标对应的第一节点位移坐标(ui,vi)、与所述第二节点坐标对应的第二节点位移坐标(ui,vi)和与所述第三节点坐标对应的第三节点位移坐标(ui,vi): 其中，A为3个节点所围成的面积： 9.根据权利要求8所述的基于剪应变的土体破裂面确定系统，其特征在于：所述计算模块还用于根据下列公式、所述第一位移坐标、所述第二位移坐标和所述第三位移坐标计算所述三节点应变计算单元的正应变εx,εy和剪应变γmax：其中表示u对x求偏导数，表示v对y求偏导数，表示u对y求偏导数，表示v对x求偏导数。 10.根据权利要求9所述的基于剪应变的土体破裂面确定系统，其特征在于：所述计算模块还用于根据下列公式、所述正应变εx,εy和所述剪应变γxy得到主应变ε1,ε2和最大剪应变γmax：
所属类别：	发明专利