probelm occur when Using Function :create_matrix in 2D situation

[mrzhaojyi]

Hi, Dear.Kristoffer

I try to generate a 2D mesh with finite thickness and I found phon.mesh_tools.create_matrix.create_matrix(mesh, thickness, order) is able to achieve that. But when I run the code:

create_matrix(mesh, 0.0001) create_matrix is EXPERIMENTAL! Traceback (most recent call last): File "", line 1, in File "/usr/local/lib/python2.7/dist-packages/phon/mesh_tools/create_matrix.py", line 55, in create_matrix create_cohesive_elements(mesh) File "/usr/local/lib/python2.7/dist-packages/phon/mesh_tools/create_cohesive_elements.py", line 73, in create_cohesive_elements grain_id_1 = grains_connected_to_face[0]

IndexError: list index out of range

Could you tell me how to solve it?

[KristofferC]

There was a bug where create_matrix didn't call create_cohesive_elements with the correct mesh_dimension.

However, I have only tested create_matrix for 3D and as the code is right now it will not work for 2D. I added an error check for this now. I will try make it work for 2D today.

[mrzhaojyi]

Thanks a lot!

Best regards JIngyi

On Mon, Jan 19, 2015 at 3:54 AM, KristofferC notifications@github.com wrote:

There was a bug where create_matrix didn't call create_cohesive_elements with the correct mesh_dimension.

However, I have only tested create_matrix for 3D and as the code is right now it will not work for 2D. I added an error check for this now. I will try make it work for 2D today.

— Reply to this email directly or view it on GitHub https://github.com/KristofferC/Phon/issues/5#issuecomment-70462060.

[KristofferC]

I've been busy with some other stuff so haven't gotten around looking at it carefully. Will try get it done tomorrow or in the weekend.

[KristofferC]

Okey, I started fixing some stuff: This is how it looks like now.

Note, that the only thing that is happening is that nodes in interfaces are moved in their normal direction. No remeshing is being done which means that if you either have a too big thickness or that you have a very fine mesh you will get elements with bad quality.

I will later experiment with adding a Laplacian smoothing to each grain which will hopefully help with the element quality.

While doing this I also found a quite severe bug with the abaqus exporting for 2d cohesive elements. Abaqus makes a very weird choice in the way their 2d cohesive elements are ordered which lead to this:

However, that should be fixed now. I will keep the issue open until I have tried the Laplacian smoothing.

[mrzhaojyi]

It looks good. Again, Thanks a lot for your help

On Fri, Jan 23, 2015 at 12:21 PM, KristofferC notifications@github.com wrote:

Okey, I started fixing some stuff: This is how it looks like now. [image: 2d_cohes] https://cloud.githubusercontent.com/assets/1282691/5879040/22bf680a-a32c-11e4-8e0c-581d6ef8f08b.png

Note, that the only thing that is happening is that nodes in interfaces are moved in their normal direction. No remeshing is being done which means that if you either have a too big thickness or that you have a very fine mesh you will get elements with bad quality.

I will later experiment with adding a Laplacian smoothing to each grain which will hopefully help with the element quality.

While doing this I also found a quite severe bug with the abaqus exporting for 2d cohesive elements. Abaqus makes a very weird choice in the way their 2d cohesive elements are ordered which lead to this: [image: abaq] https://cloud.githubusercontent.com/assets/1282691/5879095/910f5054-a32c-11e4-9e1f-6092241f37ae.png

However, that should be fixed now. I will keep the issue open until I have tried the Laplacian smoothing.

— Reply to this email directly or view it on GitHub https://github.com/KristofferC/Phon/issues/5#issuecomment-71229046.