Some researchers need to know how my spherical similar equilateral triangle graph comes from. Now I will publish the code about running in MATLAB, hoping to help you with your research.
I have received the help of netizens in solving the relevant equations and modifying the code. I would like to express my deep thanks here.
$% Similarity of equilateral triangles on sphere
%Small equilateral triangle
[x y z]=sphere(50);
mesh(x,y,z);
hold on
axis equal
K=1;
% K is the slope of the spherical circle, tangent of the angle between the plane of the spherical circle and the plane of the equator
R=1;
C=0.3;
% In the direction perpendicular to the equator, the offset of the center of the circle. If C is zero, the circle is big, otherwise it is small.
alpha =pi/6: pi/50: pi*5/6;
A=(1+K^2*(sin(alpha)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'r','linewidth',2);
grid on;
axis square;
hold on;
%------------
alpha =-pi/2: pi/50: pi/6;
A=(1+K^2*(sin(alpha-pi/3)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'g','linewidth',2);
grid on;
axis square;
hold on;
alpha =-pi*7/6: pi/50: -pi/2;
A=(1+K^2*(sin(alpha-2*pi/3)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'b','linewidth',2);
grid on;
axis square;
hold on;
% Similarity of equilateral triangles on sphere
%Large equilateral triangle
[x y z]=sphere(50);
mesh(x,y,z);
hold on
axis equal
K=1;
% K is the slope of the spherical circle, tangent of the angle between the plane of the spherical circle and the plane of the equator
R=1;
C=0.1;
% In the direction perpendicular to the equator, the offset of the center of the circle. If C is zero, the circle is big, otherwise it is small.
alpha =pi/6: pi/50: pi*5/6;
A=(1+K^2*(sin(alpha)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'r','linewidth',2);
grid on;
axis square;
hold on;
%------------
alpha =-pi/2: pi/50: pi/6;
A=(1+K^2*(sin(alpha-pi/3)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'g','linewidth',2);
grid on;
axis square;
hold on;
alpha =-pi*7/6: pi/50: -pi/2;
A=(1+K^2*(sin(alpha-2*pi/3)).^2).^0.5;
B=C*(1+K^2);
beta = 2*pi+acos(-B./A)-asin(1./A);
[X,Y,Z] = sph2cart(alpha,beta,R);
plot3(X,Y,Z, 'b','linewidth',2);
grid on;
axis square;
hold on;$
It can be seen from this code that since $K$ is constant, the corresponding edges are parallel, and the corresponding angles are equal. Since $K$ is an equilateral triangle, the corresponding edges are also proportional.

Three sides of an equilateral triangle:
$K=1$
$C=0.2$
$\alpha =\cfrac{\pi}{6}: \cfrac{\pi}{50}: \cfrac{5\pi}{6};$
$\beta = 2\pi+acos(-\cfrac{C(1+K^2 )}{\sqrt{1+K^2 (\sin\alpha)^2 } })-asin(\cfrac{1}{\sqrt{1+K^2 (\sin\alpha)^2 }} ); $
$\alpha =-\cfrac{\pi}{2}: \cfrac{\pi}{50}: \cfrac{\pi}{6};$
$\beta = 2\pi+acos(-\cfrac{C(1+K^2 )}{\sqrt{1+K^2 (sin(\alpha-\frac{\pi}{3}))^2 }} )-asin(\cfrac{1}{\sqrt{1+K^2 (sin(\alpha-\frac{\pi}{3}))^2 }} ); $
$\alpha =-\cfrac{7\pi}{6}: \cfrac{\pi}{50}: -\cfrac{\pi}{2};$
$\beta = 2\pi+acos(-\cfrac{C(1+K^2 )}{ \sqrt{1+K^2 (sin(\alpha-\frac{2\pi}{3}))^2 }} )-asin(\cfrac{1}{\sqrt{1+K^2 (sin(\alpha-\frac{2\pi}{3}))^2 } });$
If K does not change, C can be changed to ensure the parallelism of the corresponding sides of the triangle.