Hi, I am trying to solve a BVP to obtain the micropolar strain localization tensor in a 2D chiral geometry with some periodic boundary conditions. I have 2 unknowns h111 and h211.
The boundary conditions are:
h111(-0.5,y)= h111(0.5,y)
h211(-0.5,y)= h211(0.5,y)
h111(x,-0.5)= h111(x,0.5)
h211(x,-0.5)= h211(x,0.5)
Below is the code I’ve been using:

However, sigma11, calculated using h111 and h211, does not converge by refining the mesh. Sometimes it decreases and sometimes it increases. What should I do?
Thank you in advance for your answers.

Can you post the code as an attachment as my browser makes copy paste
screw ups. I don’t know if I’ll look at it but curious. Do you have a link to areference
or the eqn you are trying to solve? Probably it does not matter a lot on the mesh
for diagnosis although if you have specific features of the mesh it may help to post it.
Thanks.

Thank you for your answer.
The code and the geometry are attached to the post.
In my equation, dx(h211) equals dy(h111) because the micropolar strain localization tensor is symmetric. Before I considered that this tensor is symmetric, my equation was:

I don’t have time to sort it out as there are a lot of constants in the signal expression
and you didn’t post the equation. I get a constant signal over a range of mesh values
using a simple square mesh,
See if the attached code can help you sort it out. Also did you intend to have a
cross product in there- could be a sign mistake.

Also I would get in the habit of writing expressions with decimal points
as the consts may be evaluated as int/int but probably that is just a small
issue. You can check with this though,