Transient analysis in MAPDL unresponsive even when the same script runs well in APDL

[anshulc94]

I am trying to solve a thermal transient analysis using mapdl. The script that I am using works correctly in APDL GUI if I run it as an input file but the same converted (mapdl) script does not run give an output and just gets stuck with no sulution. Please check the script that I am using below:

from ansys.mapdl.core import launch_mapdl

mapdl = launch_mapdl(loglevel="DEBUG")

# # Thermal Analysis

mapdl.run("/CLEAR     ")                                                    # CLEAR DATABASE
mapdl.run("/UNITS,BIN        ")                                           # SET UNITS
mapdl.title("Thermal Analysis")
mapdl.prep7()

# # DEFINE ELEMENTS

mapdl.et(1, "SOLID87")                                                     # KEYOP3 SETS AXISYMMETRIC
mapdl.mshape(0, "3D")
mapdl.mshkey(0)                                                               # 0 FOR FREE MESHING AND 1 FOR MAPPED MESHING
mapdl.esize(0.01)

# # READ IN MATERIAL LIBRARY

mapdl.run("/MPLIB,READ,'C:\matlib'")
mapdl.run("/MPLIB,WRITE,'C:\matlib'")

# # MATERIALS

mapdl.mat(1)

with mapdl.non_interactive:
    mapdl.run("MPREAD,'OFHC','SI_MPL',,LIB")

# # GEOMETRY

mapdl.k(1, 0, 0, 0.0)
mapdl.k(2, -0.13612575, 6.5, 0.0)
mapdl.k(3, -0.15706815, 7.5, 0.0)
mapdl.l(2, 3)   !Line 1
mapdl.k(4, 0.15706815, 7.5, 0.0)
mapdl.larc(3, 4, 1, 7.5)
mapdl.k(5, 0.13612575, 6.5, 0.0)
mapdl.l(4, 5)   !Line 3
mapdl.larc(5, 2, 1, 6.5)
mapdl.a(2, 3, 4, 5)

mapdl.k(6, 0, 4.0, 0.0)
mapdl.circle(6, 0.1)
mapdl.a(7, 8, 9, 10)
              
mapdl.asba(1, 2)                                                           
mapdl.voffst(3, -10)                                                        # Extrude length 
mapdl.vsel("ALL")
mapdl.vlscale("ALL", "", "", 0.0254, 0.0254, 0.0254, "", 0, 1)              # scale geometry from inches to SI units
    

# # MATERIAL ASSIGNMENT    
    
mapdl.asel("S", "AREA", "", "ALL")
mapdl.aatt(1, "", 1)

# # MESHING

mapdl.vsel("ALL")
mapdl.esize(3.0e-3)
mapdl.vmesh("ALL")
    
# # BOUNDARY CONDITIONS
mapdl.run("/SOLU")  
mapdl.nsel("ALL")
mapdl.ic("ALL", "TEMP", 22)                                                 # INITIAL TEMPERATURE = 22 C
    

# # APPLY HEAT FLUX LOAD

mapdl.asel("S", "AREA", "", 6)
mapdl.sfa(6, "", "HFLUX", 1000000)                                      #HEAT FLUX VALUE = 1E6 W/m2
    
# # APPLY CONVECTION LOAD
    
mapdl.asel("S", "AREA", "", 7)
mapdl.asel("A", "AREA", "", 8)
mapdl.asel("A", "AREA", "", 9)
mapdl.asel("A", "AREA", "", 10)
mapdl.sfa("ALL", "", "CONV", 2000, 22)                        # FILM COEFFICIENT = 2000W/m2-C AT BULK TEMP = 22 C 
mapdl.scopt("NO")
    

# # TRANSFER LOAD TO ELEMENTS
    
mapdl.asel("ALL")                                                           # TO TRANSFER ALL SURFACE LOADS TO FE MODEL USING SBCTRAN
mapdl.sbctran()
mapdl.sflist()                                                                  # CHECK
mapdl.run("/PSF,HFLUX,,2")
# mapdl.eplot()
    

# # RUN SOLVER

mapdl.allsel("ALL")
mapdl.antype(4)                                                             # transient analysis
mapdl.trnopt("FULL")                                                     # full transient analysis
mapdl.kbc(1)                                                       

# # Time stepping
end_time = 0.1
mapdl.time(end_time)                                                    # end time for load step
mapdl.autots("ON")                                                       # use automatic time stepping

mapdl.deltim(0.01,0.001,0.01)                                    # TIME STEP SIZE
    
# mapdl.outres("ERASE")
mapdl.outres("NSOL", "LAST")


mapdl.solve()
print(mapdl.solve())                                                        # PRINT SOLVER COMMAND
out = mapdl.finish()

# # POST PROCESSING

mapdl.run("/POST1")
mapdl.plnsol("TEMP", "", 0)
mapdl.post_processing.plot_nodal_temperature()         # PLOT TEMPERATURE VALUES

mapdl.run("/GOPR")

mapdl.exit()

[akaszynski]

I can't see any glaring issues with this. Can you enable verbose output for the solve method and see if that helps clarify things? Use verbose=True. See Additional Options When Running Commands.

[mikerife]

Hi @anshulc94 Areas 1 and 2 do not overlap each other, so the 'asba' boolean fails. And so there is no area 3 to extrude etc.
Mike