FEAP User Forum

FEAP => Input File Issues => Topic started by: bbb on August 31, 2021, 07:03:29 AM

Title: Initial Conditions
Post by: bbb on August 31, 2021, 07:03:29 AM

Hello,

I am relatively new to using FEAP. I am trying to model a 3D poroelastic model, but I am not sure how to initialize the problem so that all the pressure in the domain, for example of a cube, is p=po.
Is there any way to initialize it without imposing the condition in INIT DISP on each node?
If not, how would be the syntax of the command?
INITial DISPlacements
1 aa 0 0 0 po
2 aa 0 0 0 po
.... until node 8?

The mesh is:
SNODES of the mesh
  1 0. 0. 0.
  2 aa 0. 0.
  3 aa aa 0.
  4 0. aa 0.
  5 0. 0. aa
  6 aa 0. aa
  7 aa aa aa
  8 0. aa aa

Thank you

Title: Re: Initial Conditions
Post by: Prof. R.L. Taylor on August 31, 2021, 09:40:28 AM

If you are using a recent version of feap you can set the initial conditions by the MATE set number.  See the user manual (Sect 11.3). 

In future it helps if you tell us which version you are using.

Title: Re: Initial Conditions
Post by: bbb on August 31, 2021, 02:34:42 PM

Thank you!