Apply a changing Dirichlet boundary condition

[Jie Wu]

Dear All,

I would like to apply a Dirichlet boundary condition like {u = [ye^x, 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

[Prof. S. Govindjee]

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.

[Prof. R.L. Taylor]

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. 

[Jie Wu]

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