DPC子程序密度变化实现方法？

请教各位大神，在使用vusdfld过程中，进行相对密度的变化写入计算，因为本构方程各个参数直接与相对密度有关，因此要计算单元的相对密度，需要计算单元体积，可是真的特别困难，孩子急疯了，求求 总结:如何在显示计算过程中直接引用单元体积，或者进行单元体积计算？