Author Topic: user element isw  (Read 5495 times)

yw-kim

  • New Member
  • *
  • Posts: 6
user element isw
« on: January 08, 2018, 10:43:12 PM »
In my user element, if(isw.eq.0) routine didn't work.
I printed isw and isw was equal to 9.
I thought isw that changes from 0 to 26.
I don't know what is the problem.

FEAP_Admin

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 993
Re: user element isw
« Reply #1 on: January 09, 2018, 12:25:36 AM »
Your element is called with isw.eq.9 when you issue the macro command DAMPing; see Appendix D of the user manual as well as the programmers manual.

yw-kim

  • New Member
  • *
  • Posts: 6
Re: user element isw
« Reply #2 on: January 09, 2018, 12:45:12 AM »
I attached my input file.
I didn't use DAMPing.
When I printed value of isw it is equal to 9.

FEAP_Admin

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 993
Re: user element isw
« Reply #3 on: January 09, 2018, 01:06:05 AM »
There is no reason that your element should get called with isw.eq.9 with this input file.  You may have a corrupt version of the code.  You should probably download a clean version. 

The other check you can make is to change the material model to
Code: [Select]
mate
  thermal
    fourier isotropic 1
Then rerun with debugging placed in elements/thermal/therm2d.f to see what isw is appearing there.  When I do that, I see
Code: [Select]
ISW =            1
 ISW =           26
 ISW =           26
 ISW =           26
 ISW =           26
 ISW =           14
 ISW =            3
No isw.eq.9 appears.  [Test was run in version 8.5]

yw-kim

  • New Member
  • *
  • Posts: 6
Re: user element isw
« Reply #4 on: January 09, 2018, 01:17:15 AM »
I tested other user element file.
Its results are same as yours.
Thanks.