FEAP User Forum

FEAP => Input File Issues => Topic started by: Jie Wu on July 29, 2016, 02:20:19 PM

Title: Apply a changing Dirichlet boundary condition
Post by: Jie Wu on July 29, 2016, 02:20:19 PM
Dear All,

I would like to apply a Dirichlet boundary condition like {u = [yex, sin(y)]} to a 1*1 square. It looks the command 'CSURface' may be helpful, tut I don't know to to implement it in FEAP.

Any suggestions would be deeply appreciated!

Best,
Jie
Title: Re: Apply a changing Dirichlet boundary condition
Post by: Prof. S. Govindjee on July 29, 2016, 04:03:21 PM
Have a look at this:  http://www.ce.berkeley.edu/~sanjay/FEAP/umesh1.html

Just beware that this was written when pointers into hr( ) and mr( ) were 32bit so you will have to the updated memory and pointer management system.
Title: Re: Apply a changing Dirichlet boundary condition
Post by: Prof. R.L. Taylor on July 30, 2016, 07:27:57 AM
A version of the user mesh routine umesh1.f that is up to date is attached.  It shows how pointer arithmetic can be avoided and be a little clearer on what is being done.  It also may correct for some compiler warnings.  Note that the calling argument for "umesh*" routines has changed and includes another parameter "tx" (which is not being used in the current module). This would  permit having commands like "TORS DISPs".  Otherwise it is the same formulation as the previous post. 
Title: Re: Apply a changing Dirichlet boundary condition
Post by: Jie Wu on August 02, 2016, 07:29:08 AM
Dear Prof. Govindjee and Prof. Taylor,

Thank you very much for your directions and especially the recommended sample umesh1.f file. They helped me a lot! Now I can apply the changing Dirichlet boundary conditions.

All the best,
Jie