There are discrepancy at both Gauss point locations and strain values. The Gauss points are derived using shape functions, and the strain-displacement matrix B, which is used to calculate strain, is obtained by differentiating these shape functions. Therefore, I suspect that the primary source of discrepancy lies in the Gauss point locations. I didn't used PRESOL command when I compare MATLAB result and ANSYS's. Â I have also confirmed that the "Integration Point Results" option is set to "averaged."
                      [simple model]
                    [complex model]
I noticed a substantial discrepancy in the x-coordinate of integration points as the complexity of the model increases. This discrepancy does not seem to rely on the number of elements, as the results align well for simpler models.
Is there an option within ANSYS to disable such automatic adjustments to the Gauss point positions or shape functions? Or, if you have any insights or suggestions on how to address this issue, I would greatly appreciate your advice.