I am doing transient heat tranfer analysis, where specific heat is a function that needed to be define as temperature dependent. Any help instead of defining the specific heat in tabular form in Abaqus Cae.
It is a common practice in Chemical Engineering to define the constant pressure specific heat (Cp) as a polynomial function of temperature (T),
Cp=a+bT+cT2+dT3
Constants a, b, c, and d can be estimated with linear regression from tabular data.
Polynomial model are easy to use in transient heat transfer by numerical solution of Fourier equation.
Thank You, Miguel Angel Garcia-Alvarado for your help. How can I implement it in Abaqus, using subroutine?
Excuse me Ankit I have not any experience in Abaqus. I have solved many transient heat transfer equations in FORTRAN or in MatLab. I think that you can implement the Cp as function of Temperature in Abaqus in the differential equation. As example, the Fourier law is,
d(RoCpT)/dt=kdiv(T)
where “Ro” is the medium density, “Cp” the medium constant pressure specific heat, “T” temperature, “k” heat conductivity and “div” is the divergence (you must define the appropriate coordinate system). The left side derivative must be partial.
Like Cp is a polynomial temperature function Cp=a+bT+cT2+dT3 the Fourier law is re-written as,
aRodT/dt + 2bTRodT/dt+3cT2RodT/dt+4dT3RodT/dt=kdiv(T) or
dT/dt=kdiv(T)/(aRo+2bTRo+3cT2Ro+4dT3Ro)
And you can follow the Abaqus protocol for differential equations.
- Badges
- Science topic