On 06/03/2015 03:41 PM, Kid Guo wrote:
at first thank you for your help, I have solved the problem.I have seen this example,sinbc()function.and I add a new function in biot.py. but it something wrong. 1.I added a new facet named Right and select its region.Like "'Right' : ('vertices in (x > 0.999999)', 'facet')," 2.a new function added in biot.py. 3.I load a stress to the left and right facet. if the stress is a constant,it worked.but when I added a new function on the left facet and another new function on the right facet,it didn't work.
What exactly does not work? I am getting (with biot.py from your next message) the linear pressure distribution in z direction on the left and right facets, as the functions for Dirichlet conditions on pressure enforce.
BTW. what you apply is the fluid pressure in the pores. To apply a stress load on the boundary, you would need to use the 'dw_surface_ltr' term (as in [1]), not the Dirichlet BC on pressure.
r. [1] http://sfepy.org/doc-devel/examples/linear_elasticity/linear_elastic_tractio...