(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111533524.3 (22)申请日 2021.12.15 (71)申请人 河北建筑工程学院 地址 075000 河北省张家口市经开区朝阳 西大街13号 (72)发明人 谢军　庞京京　张斌　阎杰　 黄久鹏　郭飞　王贺贺　张华帅　 梁金晓　王政　罗岩　赖镇泉　 陆超　王啸天　邢国斌　王丽君　 于旭涛　 (74)专利代理 机构 沈阳之华 益专利事务所有限 公司 21218 代理人 韩凌宇 (51)Int.Cl. G06F 30/13(2020.01)G06F 119/14(2020.01) (54)发明名称 一种修正Daviden kov本构模型的构建方法 (57)摘要 本发明公布了一种修正Davidenkov本构模 型的构建方法， 包括以下步骤： S1： 将 “阻尼比退 化系数”的概念引入到用分段函数表示的 Davidenkov本构模型当中； S2： 推导出修正 Davidenkov本构 模型的动应力 ‑应变骨架曲线和 滞回曲线公式； S3： 基于修正Dav idenkov应力 ‑应 变滞回曲线推导出修正Davidenkov本构模型的 阻尼比退化系数、 阻尼比和动剪切模量公式。 本 发明在用分段函数表示Davidenkov本构模型的 基础上引入了 “阻尼比退化系数 ”的概念， 不仅继 承了用分段函数表示的Davidenkov本构模型的 优点： 应力 ‑应变骨架曲线剪应力不再随剪应变 无限增大， 还改善了 由Davidenkov应力 ‑应变滞 回曲线计算的阻尼比大于土体实际阻尼比的现 象， 更加真实的体现出土体的耗能情况， 为后续 研究提供了一种合理可 行的本构模型。 权利要求书3页 说明书11页 附图8页 CN 114357573 A 2022.04.15 CN 114357573 A 1.一种修正Davidenkov本构模型的构建方法， 所述Davidenkov本构模型用分段函数表 示， 所述Daviden kov本构模型的动剪切模量比G/Gmax的表达式如下： G/Gmax＝1‑H(γ)    (1) 式中， G为剪切模 量； Gmax为初始最大剪切模量； γ为剪应变； A、 B和γ0为与土性有关的拟 合参数； 所述用分段函数表示的Daviden kov本构模型的动应力 ‑应变骨架曲线表示 为： 式中， τ(γ)为剪应力； τc为应力‑应变滞回曲线加卸载转折点处对应的剪应力； τult为破 坏剪应力； 其特征在于， 包括以下步骤： S1： 将“阻尼比退化系数 ”的概念引入到用分段函数表示的Daviden kov本构模型当中； S2： 推导出修 正Daviden kov本构模型的动应力 ‑应变骨架曲线和滞回曲线公式； S3： 基于修正Davidenkov本构模型的应力 ‑应变滞回曲线推导出修正Davidenkov本构 模型的阻尼比退化系数、 阻尼比和动剪切模量公式。 2.根据权利要求1所述的一种修正Davidenkov本构模型的构 建方法， 其特征在于： 步骤 S1所述的“阻尼比退化系数 ”是指试验阻尼比与由滞回曲线计算阻尼比的比值： K(γ)＝DT/ DM    (4) 式中， K(γ)为阻尼比退化系数； DT为通过试验得到的阻尼比； DM为由滞回曲线得到的阻 尼比。 3.根据权利要求1所述的一种修正Davidenkov本构模型的构 建方法， 其特征在于： 步骤 S2所述修正Daviden kov本构模型的动应力 ‑应变骨架曲线表示 为： 式中， γ'＝| ±γm‑γc|/2， γm为应力‑应变滞回曲线所能达 到的最大剪应 变。 4.根据权利要求1所述的一种修正Davidenkov本构模型的构 建方法， 其特征在于： 步骤 S2所述修正Daviden kov本构模型的动应力 ‑应变滞回曲线表示 为：权　利　要　求　书 1/3 页 2 CN 114357573 A 2式中， τm为应力‑应变滞回曲线所能达到的最大剪应力； τc为为应力 ‑应变滞回曲线加卸 载转折点处对应的剪应力。 5.根据权利要求1所述的一种修正的Davidenkov本构模型的构 建方法， 其特征在于： 步 骤S3所述 修正Daviden kov本构模型的阻尼比退化系数表示 为： 式中， DD为分段函数表示的Daviden kov本构模型的计算阻尼比。 6.根据权利要求1所述的一种修正的Davidenkov本构模型的构 建方法， 其特征在于： 步 骤S3所述的修 正Daviden kov本构模型的阻尼比表示 为： 式中， DMD为修正Davidenkov本构模型的 阻尼比； SLOOP为一个滞回圈的面积； SE表示三角 形的面积。 7.根据权利要求1所述的一种修正的Davidenkov本构模型的构 建方法， 其特征在于： 步 骤S3中所述推导动剪切模量公式包括以下步骤： (1)： 对(5)式 中的剪应 变γ求导， 得到初始加载 段动剪切模量的计算公式(| τ|≤ τult)：权　利　要　求　书 2/3 页 3 CN 114357573 A 3