Hi Robert,

So, I've implemented the elasticity part of the problem using Sfepy, see,

https://pages.nist.gov/pfhub/benchmarks/benchmark4.ipynb/

for the problem definition. My implementation is contained in

https://github.com/wd15/sfepy-bm4/blob/73df8d6a929c9ee70019785aef770cef57938...

and

https://github.com/wd15/sfepy-bm4/blob/73df8d6a929c9ee70019785aef770cef57938...

Basically, I'm using both the dw_lin_elastic and dw_lin_prestress to implement the two terms in the stress / strain equation. Does what I've done make sense to you with respect to how those terms are being used?

I'm not currently solving the phase field equation for $\eta$. I'm simply imposing an $\eta$ field. Initially, I'm going to solve the phase field equation using FiPy and then use Sfepy for the elasticity. Eventually, I'd like to learn how to do it all in Sfepy.

Any advice would be greatly appreciated.

Cheers,

Daniel

On Thu, Sep 6, 2018 at 7:58 AM Robert Cimrman cimrman3@ntc.zcu.cz wrote:

On 09/05/2018 06:50 PM, Daniel Wheeler wrote:

Thanks for the references. I'll let you know how it goes. I might solve the elasticity part with Sfepy and the CH part with FiPy for now.

What is the best example in Sfepy to understand how to solve multiple PDEs?

Try looking at the examples in examples/multi_physics/ - it is usually elasticity + something.

r.