Hi Ronghai,

On 11/05/2014 09:42 PM, Ronghai Wu wrote:

Dear Sfepy developers and users,

I am wondering if it is possible to solve a non-constant prestress problem with sfepy, as show

where all fields exist in the whole domain(no subdomain is considered). is a constant stress field, however is a non-constant scale field(having different values at different nodes, but values are known). So, therefore prestress is also non-constant. I know the "linear_elasticity/prestress_fibres" example in Sfepy documentation. But the prestress in that example is constant(having same values at different nodes). So, I am wondering if it is possible to implement non-constant prestress problem. Thanks.

It should work out of the box. Look at [1] to see how to define a non-constant material parameter (the prestress here) by using functions.

r.

[1] http://sfepy.org/doc-devel/users_guide.html#defining-material-parameters