Author Topic: Thermal properties in truss elements  (Read 3533 times)

Paolo Di Re

  • Jr. Member
  • **
  • Posts: 14
Thermal properties in truss elements
« on: March 28, 2017, 01:54:51 AM »
Dear All,
I was running the attached file in FEAP 8.4, which is a simple linear elastic truss with a uniform temperature variation DT = -600.

I expect the truss to shorten because I'm assigning T0 = 600 under MATErial, which is what I want, but something goes wrong. The truss shorten too much. If I check the displacement of node 2, it seems the model is using a value of the coefficient of linear thermal expansion alpha = 1e-5*E, where E is the Young's modulus.

What am I doing wrong?

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2647
Re: Thermal properties in truss elements
« Reply #1 on: March 28, 2017, 10:38:16 AM »
Yes there is a bug!  The following changes can be made:

1. In ./elements/material/inmate.f -  where d(3) is set change to read:
             if(oned) then
               d(3) = alp(1)
             else
               d(3) = alp(1) + alp(2) + alp(3)
             endif

2. In .elements/material/small/elas1d.f:  Change d(47) back to d(3).

Hopefully this fixes your problem.

Paolo Di Re

  • Jr. Member
  • **
  • Posts: 14
Re: Thermal properties in truss elements
« Reply #2 on: March 29, 2017, 01:33:37 AM »
Thank you very much. Everything works great now!