Hi Robert,

Now I encounter a more complicated case. What if the stiffness tensor is not constant, but a function of eta, as shown in PDF "non-constant elastic tensor.pdf". Is it possible to implement it based on the modified code "Ronghai.py"?

Regards
Ronghai

在 2014年11月15日星期六UTC+1下午7时02分31秒,Robert Cimrman写道:
On 11/14/2014 06:02 PM, Ronghai Wu wrote:
> Hi Robert,
>
> I have some explanations of fipy and my idea. Please see attached PDF. If
> you need more information, please just ask.

That's fine, thanks. I have modified your file and it runs. Note that you still
need to specify some boundary conditions so that the matrix is not singular.

r.