I hope this e-mail finds you well. Your input will be of incredible importance in my thesis. I am using the PDE (Partial Differential Equation) ToolBox of Matlab to obtain some stress distributions. Basically regarding the geometry, I am modeling a flaw (opening with the geometry of interest visualized upon running the below code). The major problem with the PDE ToolBox is the meshing. I am basically interested in obtaining the values of stress close to the flaw tip (tip of the opening), Therefore, I want my mesh size close to that region to be very very very small. So far, on my laptop, I am able to use a maximum mesh size (Hmax) equal to 0.135mm (smaller than that, it does not work) but i definitely need to get much smaller than this to have precise results close to the flaw tip. Because, with the actual mesh size that I have, I am not obtaining precise results. Your input here will be of extreme importance. You can find enclosed my code and the biggest challenge/problem as mentioned before is the mesh size. I am very interested in having an extremely small mesh size close to the flaw tip.
%flaw width b taken to be 0.032";
%Radius of the flaw tip equivalent to half of the flaw width b;
%L: Length of the flaw;
%La: Length of the rectangle;
xsemiL = 1*25.4;
ysemiL = 2*25.4;
% semiL= 50;
alpha = 10;
R = [3 4 -xsemiL xsemiL xsemiL -xsemiL -ysemiL -ysemiL ysemiL ysemiL]';
R1 = [3 4 (La/2)*cosd(alpha)-Radius*cosd(90-alpha) (La/2)*cosd(alpha)+Radius*cosd(90-alpha) -(La/2)*cosd(alpha)+Radius*cosd(90-alpha) (-La/2)*cosd(alpha)-Radius*cosd(90-alpha) ...
(La/2)*sind(alpha)+Radius*sind(90-alpha) (La/2)*sind(alpha)-Radius*sind(90-alpha) -(La/2)*sind(alpha)-Radius*sind(90-alpha) -(La/2)*sind(alpha)+Radius*sind(90-alpha)]';
C1 = [1 (La/2)*cosd(alpha) (La/2)*sind(alpha) Radius 0 0 0 0 0 0]';
C2 = [1 (-La/2)*cosd(alpha) (-La/2)*sind(alpha) Radius 0 0 0 0 0 0]';
gdm = [R C1 C2 R1];
ns = char('R','C1','C2','R1');
g = decsg(gdm,sf,ns');
title 'Geometry with Edge, Face and Vertex Labels';
generateMesh(model, 'Hmax', Radius/5)
Please, keep me informed if you need any additional information. I can't wait to receive from you.