这也是官方网站上的案例(linear_hardening.mph),但是下载的DLL文件不能正确地计算雅可比矩阵,导致计算中断,需要重新编译:
gcc -shared -o elastic_plastic.dll .\elasric_plastic.c
我上传了编译的dll文件,需要的可以下载,在我的win10上能够运行成功。
前言
线弹性材料满足胡克定律,在卸载后,变形能够完全恢复,应变归零。然而,一般材料当变形超过一定程度后,会发生不可恢复变形,称之为塑性变形。线性硬化是比较简单的塑性模型。下面描述弹性-线性硬化材料模型。
模型简述
线弹性阶段:
σ = 2 G ϵ d + 3 K ϵ v o l I \bm\sigma=2G\bm\epsilon_d+3K\epsilon_{vol}\bm{I} σ=2Gϵd+3KϵvolI
其中 ϵ d \bm\epsilon_d ϵd 为应变偏量, ϵ v o l \epsilon_{vol} ϵvol 为体应变。
如果应力状态超过某个面,将发生屈服,屈服函数为:
f = σ e f − σ y s f=\sigma_{ef}-\sigma_{ys} f=σef−σys
其中 σ e f f \sigma_{eff} σeff 为等效应力,表示当前的应力状态, σ y s \sigma_{ys} σys 为屈服应力,并且该屈服面随着加载历史移动,称之为硬化函数:
σ y s =