Suddenly reported an error and the program cannot run? #3045

COOLEST-GET · 2024-04-30T06:08:31Z

COOLEST-GET
Apr 30, 2024

May I ask why I used it normally a few days ago, but today I suddenly reported an error and the program cannot run?

C:\Software\Python36\python.exe C:\Software\Qt_Project\demo6\ansys_module\ansys_model.py 
CRITICAL -  -  logging - handle_exception - Uncaught exception
Traceback (most recent call last):
  File "C:\Software\Qt_Project\demo6\ansys_module\ansys_model.py", line 89, in <module>
    mapdl.solve()  # 提交求解器进行求解
  File "C:\Software\Python36\lib\site-packages\ansys\mapdl\core\_commands\solution\analysis_options.py", line 3344, in solve
    return self.run(command, **kwargs)
  File "C:\Software\Python36\lib\site-packages\ansys\mapdl\core\mapdl.py", line 2166, in run
    raise MapdlRuntimeError(f"\n\nError in instance {self._name}\n\n" + self._response)
ansys.mapdl.core.errors.MapdlRuntimeError: 

Error in instance GRPC_127.0.0.1:50052

*** NOTE ***                            CP =       4.016   TIME= 14:04:23
 The automatic domain decomposition logic has selected the MESH domain   
 decomposition method with 2 processes per solution.                     

 *****  ANSYS SOLVE    COMMAND  *****

 *** NOTE ***                            CP =       4.016   TIME= 14:04:23
 There is no title defined for this analysis.                            

 *** WARNING ***                         CP =       4.031   TIME= 14:04:23
 No shape testing has been performed on 2 of the 2175 selected elements. 
  To perform shape testing on these elements and see any warning         
 messages about the others, please issue the CHECK command.              

 *** NOTE ***                            CP =       4.031   TIME= 14:04:23
 The model data was checked and warning messages were found.             
  Please review output or errors file (                                  
 C:\Users\xiong\AppData\Local\Temp\ansys_pineqzzlev\file0.err ) for      
 these warning messages.                                                 

 *** SELECTION OF ELEMENT TECHNOLOGIES FOR APPLICABLE ELEMENTS ***
                ---GIVE SUGGESTIONS ONLY---

 ELEMENT TYPE    1 IS SOLID186. KEYOPT(2) IS ALREADY SET AS SUGGESTED.



 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2020 R2          20.2     ***
 DISTRIBUTED ANSYS Mechanical Enterprise                       

 20120530  VERSION=WINDOWS x64   14:04:23  APR 30, 2024 CP=      4.031

                                                                               



                       S O L U T I O N   O P T I O N S

   PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D                  
   DEGREES OF FREEDOM. . . . . . UX   UY   UZ  
   ANALYSIS TYPE . . . . . . . . . . . . . . . . .STATIC (STEADY-STATE)
   GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC  

 *** NOTE ***                            CP =       4.031   TIME= 14:04:23
 Present time 0 is less than or equal to the previous time.  Time will   
 default to 1.                                                           

 *** WARNING ***                         CP =       4.031   TIME= 14:04:23
 No constraints have been defined using the D command.                   

 *** NOTE ***                            CP =       4.031   TIME= 14:04:23
 The step data was checked and warning messages were found.              
  Please review output or errors file (                                  
 C:\Users\xiong\AppData\Local\Temp\ansys_pineqzzlev\file0.err ) for      
 these warning messages.                                                 

 *** NOTE ***                            CP =       4.031   TIME= 14:04:23
 The conditions for direct assembly have been met.  No .emat or .erot    
 files will be produced.                                                 

  
  
     D I S T R I B U T E D   D O M A I N   D E C O M P O S E R
  
  ...Number of elements: 2175
  ...Number of nodes:    4581
  ...Decompose to 2 CPU domains
  ...Element load balance ratio =     1.001


                      L O A D   S T E P   O P T I O N S

   LOAD STEP NUMBER. . . . . . . . . . . . . . . .     1
   TIME AT END OF THE LOAD STEP. . . . . . . . . .  1.0000    
   NUMBER OF SUBSTEPS. . . . . . . . . . . . . . .     1
   STEP CHANGE BOUNDARY CONDITIONS . . . . . . . .    NO
   PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
   DATABASE OUTPUT CONTROLS. . . . . . . . . . . .ALL DATA WRITTEN
                                                  FOR THE LAST SUBSTEP


 SOLUTION MONITORING INFO IS WRITTEN TO FILE= file.mntr                                                                                                                                                                                                                                                           


 Range of element maximum matrix coefficients in global coordinates
 Maximum = 2.09009476E+13 at element 2020.                               
 Minimum = 3.360014652E+11 at element 1143.                              

   *** ELEMENT MATRIX FORMULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1      2175  SOLID186      0.281   0.000129
 Time at end of element matrix formulation CP = 4.3125.                  

 DISTRIBUTED SPARSE MATRIX DIRECT SOLVER.
  Number of equations =       13743,    Maximum wavefront =    339

  Local memory allocated for solver              =     18.056 MB
  Local memory required for in-core solution     =     17.403 MB
  Local memory required for out-of-core solution =     10.310 MB

  Total memory allocated for solver              =     37.281 MB
  Total memory required for in-core solution     =     35.903 MB
  Total memory required for out-of-core solution =     20.849 MB

 *** NOTE ***                            CP =       4.422   TIME= 14:04:23
 The Distributed Sparse Matrix Solver is currently running in the        
 in-core memory mode.  This memory mode uses the most amount of memory   
 in order to avoid using the hard drive as much as possible, which most  
 often results in the fastest solution time.  This mode is recommended   
 if enough physical memory is present to accommodate all of the solver   
 data.                                                                   
 Distributed sparse solver maximum pivot= 3.093863632E+13 at node 4376   
 UZ.                                                                     
 Distributed sparse solver minimum pivot= -311.430652 at node 4444 UY.   
 Distributed sparse solver minimum pivot in absolute value= 0.202319209  
 at node 950 UZ.                                                         

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 There is at least 1 small equation solver pivot term (e.g., at the UZ   
 degree of freedom of node 950).  Please check for an insufficiently     
 constrained model.                                                      

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 The value of UY at node 111 is 2703360.  It is greater than the current 
 limit of 1000000 (which can be reset on the NCNV command).  This        
 generally indicates rigid body motion as a result of an unconstrained   
 model.  Verify that your model is properly constrained.                 

   *** ELEMENT RESULT CALCULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1      2175  SOLID186      0.125   0.000057

   *** NODAL LOAD CALCULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1      2175  SOLID186      0.000   0.000000
 *** LOAD STEP     1   SUBSTEP     1  COMPLETED.    CUM ITER =      1
 *** TIME =   1.00000         TIME INC =   1.00000      NEW TRIANG MATRIX


 *** ANSYS BINARY FILE STATISTICS
  BUFFER SIZE USED= 16384
        3.250 MB WRITTEN ON ASSEMBLED MATRIX FILE: file0.full
        1.688 MB WRITTEN ON RESULTS FILE: file0.rst

进程已结束,退出代码1

Answered by germa89

Apr 30, 2024

Hi @COOLEST-GET

In the last weeks, no new PyMAPDL release has been published, nor new MAPDL release.

Hence I doubt it is a library error (I'm assuming you are not using PyMAPDL from a cloned repo).

Now, it seems you might have changed something in the model because you get:

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 There is at least 1 small equation solver pivot term (e.g., at the UZ   
 degree of freedom of node 950).  Please check for an insufficiently     
 constrained model.                                                      

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 The value of UY at node 111 is 2703360.  It is g…

View full answer

COOLEST-GET · 2024-04-30T06:09:22Z

COOLEST-GET
Apr 30, 2024
Author

Here is all my code.

from ansys.mapdl.core import launch_mapdl


def daoru():
    f = open("C://Software//Qt_Project//demo6//ansys_module//shangyuanzhu.txt", "r", encoding='utf-8')
    str_0 = f.read()
    f.close()
    return str_0


def daoru1():
    f = open("C://Software//Qt_Project//demo6//ansys_module//xiayuanzhu.txt", "r", encoding='utf-8')
    str_1 = f.read()
    f.close()
    return str_1


def daoru2():
    f = open("C://Software//Qt_Project//demo6//ansys_module//changfangti.txt", "r", encoding='utf-8')
    str_2 = f.read()
    f.close()
    return str_2


# 初始化 ANSYS Mapdl
mapdl = launch_mapdl()
mapdl.clear()
mapdl.prep7()

# 创建几何体
b1 = mapdl.block(daoru2())  # 创建中间长方体1
mapdl.block(70, 130, -30, 30, -30, -20)  # 创建一边的长方体2
mapdl.block(-130, -70, -30, 30, 20, 30)  # 另一边的长方体3

c_top = mapdl.cylind(daoru())  # 创建上圆柱体
c_bottom = mapdl.cylind(daoru1())  # 创建下圆柱体

mapdl.wprota(0, -90)  # 设置工作平面旋转角度(xy ,yz ,zx)
mapdl.wpoffs(40, 0, 30)  # 设置工作平面偏移量
c_lft = mapdl.cyl4(0, 0, 17.5, "", "", "", -60)  # 创建左圆柱体
mapdl.wpoffs(-80, 0, 0)
c_rgt = mapdl.cyl4(0, 0, 17.5, "", "", "", -60)  # 创建右圆柱体
mapdl.wpoffs(40, 0, -30)
mapdl.wprota(0, 90)
b_sub = mapdl.block(-40, 40, -30, 30, -7.5, 7.5)  # 创建子长方体

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()  # 退出当前处理器

# 绘制模型
mapdl.eplot(off_screen=True)
pl = mapdl.vplot(quality=5, return_plotter=True, off_screen=True)

# 访问求解器
mapdl.slashsolu()

# 设置边界条件
mapdl.nsel("S", "LOC", "Z", "-45")  # 选择z坐标值为-45的全部结点
mapdl.d("ALL", "ALL", "0")  # 指定边界条件为固定约束
mapdl.allsel("ALL")  # 全选

# 施加加载荷
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()
mapdl.post_processing.plot_nodal_elastic_component_strain('X', 'XZ')

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

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

# 选择左上角节点并输出应变结果
leftup = mapdl.nsel('s', 'loc', -80, 0, 18.5)
cmd1 = mapdl.prnsol("leftup", "PRIN")
with open("Result_1.txt", "w") as fo:
    fo.write(cmd1)

allow_ignore = True

mapdl.finish()

0 replies

germa89 · 2024-04-30T08:07:43Z

germa89
Apr 30, 2024
Maintainer

Hi @COOLEST-GET

In the last weeks, no new PyMAPDL release has been published, nor new MAPDL release.

Hence I doubt it is a library error (I'm assuming you are not using PyMAPDL from a cloned repo).

Now, it seems you might have changed something in the model because you get:

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 There is at least 1 small equation solver pivot term (e.g., at the UZ   
 degree of freedom of node 950).  Please check for an insufficiently     
 constrained model.                                                      

 *** ERROR ***                           CP =       4.484   TIME= 14:04:23
 The value of UY at node 111 is 2703360.  It is greater than the current 
 limit of 1000000 (which can be reset on the NCNV command).  This        
 generally indicates rigid body motion as a result of an unconstrained   
 model.  Verify that your model is properly constrained.

As the message explains, it seems your model is not constrained enough so it is moving freely in space. Check you have applied forces and/or displacements in all three directions.

2 replies

COOLEST-GET Apr 30, 2024
Author

I found a problem. It was when I used a txt document as an input parameter, which caused the program to fail to run. After canceling the functions daoru1 and daoru2, the program was able to run.This is strange

germa89 May 7, 2024
Maintainer

MMhh.. I guess it will depend on what are written in those txt files.

COOLEST-GET · 2024-04-30T08:13:57Z

COOLEST-GET
Apr 30, 2024
Author

I used this approach as a solution in the previous version of this code, but it was not resolved in the version I released. Why do I have to make changes to the code when my SB teammate wrote something

1 reply

germa89 May 7, 2024
Maintainer

I'm not sure to understand what you mean with that comment.

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Suddenly reported an error and the program cannot run? #3045

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Replies: 3 comments 3 replies

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

Suddenly reported an error and the program cannot run? #3045

COOLEST-GET Apr 30, 2024

Replies: 3 comments · 3 replies

COOLEST-GET Apr 30, 2024 Author

germa89 Apr 30, 2024 Maintainer

COOLEST-GET Apr 30, 2024 Author

germa89 May 7, 2024 Maintainer

COOLEST-GET Apr 30, 2024 Author

germa89 May 7, 2024 Maintainer

COOLEST-GET
Apr 30, 2024

Replies: 3 comments 3 replies

COOLEST-GET
Apr 30, 2024
Author

germa89
Apr 30, 2024
Maintainer

COOLEST-GET Apr 30, 2024
Author

germa89 May 7, 2024
Maintainer

COOLEST-GET
Apr 30, 2024
Author

germa89 May 7, 2024
Maintainer