I have some explanations of fipy and my idea. Please see attached PDF. If you need more information, please just ask.
在 2014年11月14日星期五UTC+1下午2时54分04秒，Robert Cimrman写道： >
I have just installed fipy and tried it. Could you explain more the shape of the prestress array? It is (nx ny 4, 3), where nx * ny is the number of cells (elements of the mesh). Why the factor 4 there? The values seem to be constant in each cell.
On 11/14/2014 01:03 PM, Ronghai Wu wrote:
I still cannot fix the bug in my codes nor find other better ways. I am wondering have you tried my codes found the problem?
在 2014年11月11日星期二UTC+1下午4时45分43秒，Robert Cimrman写道：
On 11/11/2014 03:03 PM, Ronghai Wu wrote:
Thanks for your reply. I tried your suggestions with the assistance of
http://sfepy.org/doc-devel/users_guide.html#defining-material-parameters but seems not working. There are some error massages which I cannot fix. Could you please provide more help. I additionally have some explanations and questions regarding my wield codes.
So eta is then defined in the mesh vertices? Or elsewhere?
Yes, it seems fine. The same mesh (only "transposed" - the numbering of nodes goes y-axis first) can be obtained by:
./script/blockgen.py --2d -s [4,4,4] -o out.mesh
that uses "from sfepy.mesh.mesh_generators import gen_block_mesh".
The materials are evaluated in quadrature points for numerical integration and not in the mesh vertices. If your values are in the mesh vertices, you would need to interpolate them to the quadrature points, for example using a field.
Your field is bi-linear, so 3 should be ok.
Yes, those set up the boundary conditions and materials. A stationary problem is solved as a single step of a quasi-static solver.
My complete codes are shown below. also attached. Thank you.
I will look at it.