ANSYS中的速率状态依赖摩擦

亲爱的同事

一个快速的问题:-我正在玩接触的连续模型(赫兹，韦斯特加德等)和摩擦定律。但是，我更喜欢让我的学生使用ANSYS有限元代码。在Ansys中有标准库仑，速率依赖，但不是速率状态(除非有人编码一些用户子程序)。看到

https://www.sharcnet.ca/Software/Fluent14/help/ans_ctec/Hlp_ctec_realkey..。

Dieterich(1979)和Ruina(1983)(另见Dieterich(2007))提出了一个速率和状态依赖的摩擦定律来解释大量岩石摩擦的实验室数据。其他正则化的摩擦定律也存在，特别是Cochard和Rice(2000)用实验得到的Prakash和Clifton(1993)摩擦定律的简化版本，证实了在不稳定范围内数值解不收敛通过网格尺寸减小。最后，Ben-Zion和Huang(2002)发现，对于足够大的传播距离，自锐化和发散行为也存在正则化定律，并且“正则化摩擦的参数必须经过微调，才能在给定的传播距离内产生明显的稳定性”。

你知道谁可以为经典速率状态定律实现这些用户子程序吗?

对您的帮助表示诚挚的问候和感谢。

迈克

参考文献

迪特里希，j.h.(1979)。岩石摩擦建模实验结果和本构方程。j .地球物理学。第84号，2161-2168号决议
陈志伟(2007)基于摩擦速率和状态的断层滑动模型的应用，地球物理学报，2007，第4章，第109 - 109页
普拉卡什V.，克利夫顿r.j. (1993)“压力-剪切中的时间分辨动态摩擦测量，可变形固体动力学的实验技术”，应用力学Div，卷165,33 -48。
Cochard A.， Rice J.R. (2000)“不同物质之间的断层破裂:病态、正则化和滑动脉冲响应”，地球物理学报，105,25 891。