Calling pyansys in C++

[COOLEST-GET]

If you could provide assistance, I would greatly appreciate it.
In PyANSYS, the function I implement is to create a model, encapsulate it as a function, and then perform stress analysis, which is also encapsulated as a function.
I want to call the model building function first in C++, and then operate by calling the stress analysis function multiple times.
The problem I am facing now is that the process of stress analysis after building the model is only successful on the first call, and subsequent calls to the stress analysis function do not run.

Here is my code：
mapdl = launch_mapdl()
def build():
# mapdl = launch_mapdl()
mapdl.clear()
mapdl.prep7()
pv.global_theme.camera['position'] = [0, 1, 0]
b1 = mapdl.block(daoru2()) # 创建中间长方体1
b1x1, b1x2, b1y1, b1y2, b1z1, b1z2 = 70, 130, -30, 30, -30, -20
b2x1, b2x2, b2y1, b2y2, b2z1, b2z2 = -130, -70, -30, 30, 20, 30
mapdl.block(b1x1, b1x2, b1y1, b1y2, b1z1, b1z2) # 创建一边的长方体2
mapdl.block(b2x1, b2x2, b2y1, b2y2, b2z1, b2z2) # 另一边的长方体3

c_top = mapdl.cylind(daoru())     # 4
c_bottom = mapdl.cylind(daoru1())     # 5
xy ,yz= 0, -90
mapdl.wprota(xy ,yz)        #设置了工作平面的旋转角度(xy ,yz ,zx)
pianyix, pianyiy, pianyiz = 40, 0, 30
mapdl.wpoffs(pianyix, pianyiy, pianyiz)      #设置了工作平面的偏移量
c_lft = mapdl.cyl4(0, 0, 17.5, "", "", "", -60)     # 6

mapdl.wpoffs(-80, 0, 0)
c_rgt = mapdl.cyl4(0, 0, 17.5, "", "", "", -60)     # 7

mapdl.wpoffs(40, 0, -30)
mapdl.wprota(0, 90)
b_sub = mapdl.block(-40, 40, -30, 30, -7.5, 7.5)    # 8

mapdl.finish()

mapdl.prep7()   # 访问前处理器

mapdl.vlist()
mapdl.vadd(6, 7, 8)
mapdl.vsbv(1, 9)
mapdl.vadd(2, 3, 4, 5, 6)

# 定义属性
mapdl.et(1, "SOLID186")     # 定义单元类型
mapdl.mp("EX", "1", "2E11")     # 定义弹性模量
mapdl.mp("NUXY", "1", "0.3")    # 定义泊松比

# 网格划分
mapdl.type("1")     # 激活单元类型号1
mapdl.mat("1")      # 激活材料号1
mapdl.smrtsize("3")     # 智能分网，大小为3。1-10中数字越小划分越精细
mapdl.mshape(1)     # 四面体网格
mapdl.vmesh("ALL")  # 执行分网

mapdl.finish()   # 退出当前处理器

pl = mapdl.vplot(quality=5, return_plotter=True)
pl.export_obj("C:/Software/Qt_Project/demo6/ansys_module/model")
# 访问求解器
mapdl.slashsolu()
# 设置边界条件
mapdl.nsel("S", "LOC", "Z", "-45")    # 选择z坐标值为-45的全部结点
mapdl.d("ALL", "ALL", "0")          # 指定边界条件为固定约束
mapdl.allsel("ALL")                 # 全选

def calculate():
# 施加载荷
mapdl.f("ALL", "FZ", "198") # 施加荷载
mapdl.allsel("ALL") # 全选

# 求解
mapdl.antype("STATIC")              # 设置分析类型为静力分析
mapdl.solve()                       # 提交求解器进行求解
mapdl.finish()                      # 退出当前处理器

# 后处理
mapdl.post1()
mapdl.set("last")
mapdl.prnsol("EPTO", "COMP")
mapdl.post_processing.plot_nodal_total_eqv_strain(cpos = 'xz', off_screen=True, savefig='seqv_00.png')
mapdl.view("all", 0, 100, 0)
# mapdl.post_processing.plot_nodal_elastic_component_strain('Z', 'YZ')

# 应变
mapdl.set("LAST")
cmd = mapdl.prnsol("EPTO", "PRIN")

# 保存结果
with open("Result_EPTO.txt", "w") as fo:
    fo.write(cmd)

mapdl.wpoffs(40, 0, 30)
mapdl.nsel('s', 'loc', 57.5, 0, 20)
cmd1 = mapdl.prnsol("EPTO", "PRIN")

with open("Result_1.txt", "w") as fo:
    fo.write(cmd1)

allow_ignore= True
# mapdl.exit()