Dear All,
I would like to raise a question regarding the projection of user-defined variables in Parfeap.
In the attached files, we consider a notched plate composed of two materials: an elastic solid and a user-defined cohesive element that separates the two halves of the specimen. To visualize specific magnitudes of the cohesive elements (such as damage or normal stresses), we need to project values from the integration points to the nodes. This projection is implemented in the function gaussValueProjection47 within elmt47.f. The user-defined magnitudes for the cohesive elements begin from the 13th variable onward (the first 12 correspond to the 6 stress and 6 strain components of the elastic material, to avoid conflicts).
When this procedure is performed using standard FEAP (see the "Feap" folder), the resulting isocontour for damage in Paraview appears smooth, as illustrated in the attached PDF. However, when using ParFEAP (tested with 8 and more partitions; see the "Parfeap" folder), the isocontour for damage displays noticeable discontinuities across partition boundaries, as shown in the attached PNG images.
Our objective is to obtain the same smooth damage isocontour in ParFEAP as observed with FEAP. We suspect this issue may be related to the change in global element numbering (n_el variable) due to the parallel decomposition, and possibly also due to the use of multiple materials, with only one of them being user-defined. Is there any way to look for a tag of the elements?
We would appreciate any guidance to resolve this issue. Many thanks!
