温度依赖线弹性及热胀冷缩考虑的UMAT源代码分享

温度依赖线弹性

1 本构理论

1.1 率形式

本构方程为：

2 UMAT代码

umat代码用C++编写。其中材料的性质与温度相关，根据插值得到，相关代码如下：

// compute elastic properties by interpolating
double interpolate_elastic_properties(double x0, double x1, double y0, double y1, double x)
{
    double fac0 = 1.0 - (x - x0) / (x1 - x0);
    double fac1 = (x - x0) / (x1 - x0);
    if (x < x0)
    {
        return y0;
    }
    else if (x > x1)
    {
        return y1;
    }
    else
    {
        return fac0 * y0 + fac1 * y1;
    }
}

3 测试

3.1 带孔板拉伸

Abaqus的计算结果为：

利用umat计算的结果为：

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删