Hi Bin,
Just recently I was experimenting with getfem and I developed a simple
example with poroelastic coupled problem (linear Biot consolidation theory).
In the attachment you can find all the files for this implementation.
If you have getfem properly installed just run "make all" from a command
line (in a bash shell for instance). This will
produce a program "poroelastic".
Then run "poroelastic pe.param". pe.param contains some description of the
problem.
Using this program you can solve both, 2D and 3D consolidation problems on
a unit cubic mesh. (the dimension is controlled with N parameter in pe.param.
But if you change it to 2D don't forget to also change the elements for the
analysis)
The program will produce .vtk files with displacements and pressure for
every increment of the consolidation analysis. View them with MayaVI. (The way
it is implemented it requires linear approximation for both pressure and the
displacement. But this is not a limitation of the analysis, but rather a
limitation of how I do the postprocessing. But you can do it in a proper way if
you want to)
Best regards,
Andriy
________________________________
Andriy Andreykiv(PhD, MSc)
Delft University of Technology
Department of Precision Engineering
Group of Fundamentals of Microsystems
Mekelweg 2
2628 CD Delft
The Netherlands
E-mail : [EMAIL PROTECTED]
Tel. : (31) 15 2786818
Fax. : (31) 15 2789475
www : http://www-tm.wbmt.tudelft.nl/~andrico
private: (31) 6 47376804
----- Original Message -----
From: Bin
To: [email protected]
Sent: Friday, May 18, 2007 7:25 AM
Subject: [Getfem-users] Coupled problem?
Hello,
Is there any example codes on coupled problem in getfem++? ( such as biot
consolidation problem in soil mechanics)
Thanks
--
Yours truly,
Bin
------------------------------------------------------------------------------
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
poroelastic.cpp
Description: Binary data
itoa.h
Description: Binary data
Makefile.dat
Description: Binary data
pe.param
Description: Binary data
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
