Python UPF Problem #1078

kkytg · 2022-04-28T14:58:00Z

kkytg
Apr 28, 2022

Hello all,

I am studying on UPF's in PyAnsys in Linux Suse.

The current ansys version is 2022R2

I followed the instructions and tried to perform usermat example in PyAnsys User Guide documents as given in https://mapdldocs.pyansys.com/user_guide/extended_examples/Python_UPF/python_upf_examples.html

Then I open Ansys Mechanical APDL Launcher and input my dat file. However I get the following error as shown in Solution Output which I provided below. Please also note that I build Pyansys on Anaconda3 with Python 3.9.

It would be really appreciated if you kindly share your opinions.

Thanks in advance!

 Processing "/upf" found in input file "/home/suse/Documents/upf_studies/usermat.dat"

 Python UPF Detected

PYTHON VERSION : 3.6
**Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'

Current thread 0x00007f0e4f203700 (most recent call first):
/ansys_inc/v221/ansys/bin/ansupf221: line 90:  2800 Aborted                 (core dumped) python${PYTHON_VER} -m pip install ${ANSYS_DIR}/syslib/ansGRPC/User/mapdl-pyupf.tar.gz --user > /dev/null
Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'**

Current thread 0x00007f901edd0700 (most recent call first):

 Ansys Mechanical Enterprise                       


 *------------------------------------------------------------------*
 |                                                                  |
 |   W E L C O M E   T O   T H E   A N S Y S (R)  P R O G R A M     |
 |                                                                  |
 *------------------------------------------------------------------*




 ***************************************************************
 *            ANSYS 2022 R1          LEGAL NOTICES             *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2022 ANSYS, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of ANSYS, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the ANSYS, Inc. online documentation or the ANSYS, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * ANSYS, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * ANSYS, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the ANSYS, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the ANSYS, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

 2022 R1 
     
 Point Releases and Patches installed:   
     
 Ansys, Inc. License Manager 2022 R1 
 Ansys, Inc. Products 2022 R1
 Customization Files for User Programmable Features 2022 R1  
 Mechanical Products 2022 R1 
 Icepak (includes CFD-Post) 2022 R1  
 ACIS Geometry Interface 2022 R1 
 Catia, Version 5 Geometry Interface 2022 R1 
 Parasolid Geometry Interface 2022 R1


          *****  ANSYS COMMAND LINE ARGUMENTS  *****
  INITIAL JOBNAME              = usrmatjhddfg                                                                                                                                                                                                                                            
  BATCH MODE REQUESTED (-b)    = LIST
  INPUT FILE COPY MODE (-c)    = COPY
  DISTRIBUTED MEMORY PARALLEL REQUESTED
       2 PARALLEL PROCESSES REQUESTED WITH SINGLE THREAD PER PROCESS
    TOTAL OF     2 CORES REQUESTED
  MPI OPTION                   = INTELMPI
  INPUT FILE NAME              = /home/suse/Documents/upf_studies/usermat.dat
  START-UP FILE MODE           = READ  
  STOP FILE MODE               = READ  
  LANGUAGE                     = en-us               
  INITIAL DIRECTORY = /home/suse/Documents/upf_studies

 User Link path (ANS_USER_PATH): /ansys_inc/v221/ansys/syslib/ansGRPC/User                                                                                                                                                                                                               

 RELEASE= 2022 R1              BUILD= 22.1      UP20211129   VERSION=LINUX x64   
 CURRENT JOBNAME=usrmatjhddfg0  17:39:47  APR 28, 2022 CP=      0.429
 Note - This ANSYS version was linked by Licensee
     1  /batch,list 
     2  /title,upf-py1s, test usermat.py with 3-D elements  
     3      
     4  /prep7  
     5  /upf,usermat.py 
     6  USRCAL,UANBEG   
     7  tb,user,1,,2
     8  tbdata,1,1e5, 0.3    ! E, Poisson   
     9      
    10  et,1,185
    11      
    12  block,0,1,0,1,0,1   
    13  esize,1 
    14  mshape,0,3D 
    15  vmesh,all   
    16  elist   
    17      
    18  nsel,s,loc,x,0  
    19  d,all,ux
    20      
    21  nsel,s,loc,y,0  
    22  d,all,uy
    23      
    24  nsel,s,loc,z,0  
    25  d,all,uz,   
    26      
    27  allsel,all  
    28  finish  
    29      
    30  /solu   
    31      
    32  time,1  
    33  deltime,0.1 
    34  eresx,no
    35  nsel,s,loc,x,1  
    36  !d,all,ux,-0.01 
    37  sf,all,pres,1000        ! pressure on x-axis
    38  allsel,all  
    39      
    40  outres,all,all  
    41      
    42  solve			       			  
    43      
    44  finish  
    45  /POST1  
    46  set,last
    47  esel,s,elem,,1  
    48  /output 
    49  presol,s
    50  presol,epel 
    51  /com, expected results: Sx=-1000, epel_x=-1e-2  
    52  finish  
    53  /exit,nosave

 Connection to the gRPC Server at address [192.168.1.101:50054]  : OK

 
 ***** /UPF routine = usermat.py  *****
     1  import grpc 
     2  import sys  
     3  import math 
     4  import numpy as np  
     5  from mapdl import * 
     6      
     7  class MapdlUserService( mapdl.MapdlUser_pb2_grpc.MapdlUserServiceServicer ):
     8      
     9      
    10  #   #################################################################   
    11      def UserMat(self, request, context):
    12      
    13          ncomp = request.ncomp   
    14          nDirect = request.nDirect   
    15      
    16          response = mapdl.MapdlUser_pb2.UserMatResponse()
    17      
    18          response.stress[:] = request.stress[:]  
    19          response.ustatev[:] = request.ustatev[:]
    20          response.sedEl = request.sedEl  
    21          response.sedPl = request.sedPl  
    22          response.epseq = request.epseq  
    23          response.epsPl[:] = request.epsPl[:]
    24          response.var0 = request.var0
    25          response.var3 = request.var3
    26          response.var4 = request.var4
    27          response.var5 = request.var5
    28          response.var6 = request.var6
    29          response.var7 = request.var7
    30      
    31          if ncomp > 4:                        # ***    3d, plane strain and axisymmetric example 
    32              usermat3d( request, context, response)  
    33          elif nDirect== 2 and ncomp == 3:     # ***    plane stress example  
    34              usermatps( request, context, response)  
    35          elif ncomp == 3:                     # ***    3d beam example   
    36              usermatbm( request, context, response)  
    37          elif ncomp == 1:                     # ***    1d beam example   
    38              usermat1d( request, context, response)  
    39      
    40          return response 
    41      
    42      
    43      
    44  def usermat3d( request, context, response): 
    45      
    46      ZERO       = 0. 
    47      HALF       = 0.5
    48      THIRD      = 1./3.  
    49      ONE        = 1. 
    50      TWO        = 2. 
    51      SMALL      = 1.e-08 
    52      sqTiny     = 1.e-20 
    53      ONEDM02    = 1.e-02 
    54      ONEDM05    = 1.e-05 
    55      ONEHALF    = 1.5
    56      TWOTHIRD   = 2.0/3.0
    57      mcomp      = 6  
    58      
    59      G = [1., 1., 1., 0., 0. ,0.]
    60      
    61      db.start()                          # Connect to the MAPDL DB gRPC Server   
    62      ncomp = request.ncomp   
    63      
    64      # *** get Young's modulus and Poisson's ratio   
    65      young    = request.prop[0]  
    66      posn     = request.prop[1]  
    67      twoG     = young / (ONE+posn)   
    68      elast1   = young*posn/((1.0+posn)*(1.0-TWO*posn))   
    69      elast2   = HALF*twoG
    70      
    71      #   
    72      # *** calculate elastic stiffness matrix (3d)   
    73      #   
    74      dsdeEl = np.zeros( ( 6, 6)) 
    75      
    76      dsdeEl[0,0] = (elast1+TWO*elast2)*G[0]*G[0] 
    77      dsdeEl[0,1] = elast1*G[0]*G[1]+elast2*TWO*G[3]*G[3] 
    78      dsdeEl[0,2] = elast1*G[0]*G[2]+elast2*TWO*G[4]*G[4] 
    79      dsdeEl[0,3] = elast1*G[0]*G[3]+elast2*TWO*G[0]*G[3] 
    80      dsdeEl[0,4] = elast1*G[0]*G[4]+elast2*TWO*G[0]*G[4] 
    81      dsdeEl[0,5] = elast1*G[0]*G[5]+elast2*TWO*G[3]*G[4] 
    82      
    83      dsdeEl[1,1] = (elast1+TWO*elast2)*G[1]*G[1] 
    84      dsdeEl[1,2] = elast1*G[1]*G[2]+elast2*TWO*G[5]*G[5] 
    85      dsdeEl[1,3] = elast1*G[1]*G[3]+elast2*TWO*G[0]*G[3] 
    86      dsdeEl[1,4] = elast1*G[1]*G[4]+elast2*TWO*G[0]*G[4] 
    87      dsdeEl[1,5] = elast1*G[1]*G[5]+elast2*TWO*G[1]*G[5] 
    88      
    89      dsdeEl[2,2] = (elast1+TWO*elast2)*G[2]*G[2] 
    90      dsdeEl[2,3] = elast1*G[2]*G[3]+elast2*TWO*G[4]*G[5] 
    91      dsdeEl[2,4] = elast1*G[2]*G[4]+elast2*TWO*G[4]*G[2] 
    92      dsdeEl[2,5] = elast1*G[2]*G[5]+elast2*TWO*G[5]*G[2] 
    93      
    94      dsdeEl[3,3] = elast1*G[3]*G[3]+elast2*(G[0]*G[1]+G[3]*G[3]) 
    95      dsdeEl[3,4] = elast1*G[3]*G[4]+elast2*(G[0]*G[5]+G[4]*G[3]) 
    96      dsdeEl[3,5] = elast1*G[3]*G[5]+elast2*(G[3]*G[5]+G[4]*G[1]) 
    97      
    98      dsdeEl[4,4] = elast1*G[4]*G[4]+elast2*(G[0]*G[2]+G[4]*G[4]) 
    99      dsdeEl[4,5] = elast1*G[4]*G[5]+elast2*(G[3]*G[2]+G[4]*G[5]) 
   100      
   101      dsdeEl[5,5] = elast1*G[5]*G[5]+elast2*(G[1]*G[2]+G[5]*G[5]) 
   102      
   103      for i in range( 0, 5):  
   104          for j in range( i+1, 6):
   105              dsdeEl[j,i] = dsdeEl[i,j]   
   106      
   107      Strain = np.zeros( ncomp)   
   108      Strain[0:ncomp] = request.Strain[0:ncomp]   
   109      dStrain = np.zeros( ncomp)  
   110      dStrain[0:ncomp] = request.dStrain[0:ncomp] 
   111      
   112      #   
   113      # *** calculate the stress and  
   114      #     copy elastic moduli dsdeEl to material Jacobian matrix
   115      
   116      strainEl = np.copy(Strain)                  # strainEl = Strain 
   117      strainEl = np.add( strainEl, dStrain)       # strainEl += dStrain   
   118      
   119      dsdePl = np.copy(dsdeEl)
   120      sigElp = np.zeros ( ncomp)  
   121      sigElp = dsdeEl.dot( strainEl)  
   122      
   123      response.stress[:] = sigElp 
   124      dsdePl.shape = (6*6)
   125      response.dsdePl[:] = dsdePl 
   126      
   127      return response 
   128      
   129  if __name__ == '__main__':  
   130      upf.launch( sys.argv[0])
 ***** End of /UPF routine = usermat.py  *****

 RUN SETUP PROCEDURE FROM FILE= /ansys_inc/v221/ansys/apdl/start.ans
 Start GRPC DBServer

 ##############################
 ### START GRPC SERVER      ###
 ##############################

 Server Executable   : MapdlGrpc.DB.Server
 Server listening on : 0.0.0.0:50055

 /INPUT FILE= /ansys_inc/v221/ansys/apdl/start.ans  LINE=       0

 TITLE= 
 upf-py1s, test usermat.py with 3-D elements                                   


 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:40:08  APR 28, 2022 CP=      0.851

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).

          ***** ANSYS ANALYSIS DEFINITION (PREP7) *****

 THE ANSYS PROGRAM WILL CALL USER SUBROUTINE(S):
   UANBEG  

 *** NOTE ***                            CP =       0.852   TIME= 17:40:08
 Calling user subroutines ( UANBEG ) is a non-standard use of the ANSYS  
 program.                                                                

 DEFINE DATA TABLE USER FOR MATERIAL 1 
 USER MATERIALS OPTION 
 WITH A MAXIMUM OF 2 DATA ITEMS
 DATA FOR  USER  TABLE FOR MATERIAL   1 
 LOC=  1 1.00000e+05     3.00000e-01    


 ELEMENT TYPE          1 IS SOLID185     3-D 8-NODE STRUCTURAL SOLID 
  KEYOPT( 1- 6)=        0      0      0        0      0      0
  KEYOPT( 7-12)=        0      0      0        0      0      0
  KEYOPT(13-18)=        0      0      0        0      0      0

 CURRENT NODAL DOF SET IS  UX    UY    UZ  
  THREE-DIMENSIONAL MODEL

 CREATE A HEXAHEDRAL VOLUME WITH
 X-DISTANCES FROM      0.000000000     TO      1.000000000    
 Y-DISTANCES FROM      0.000000000     TO      1.000000000    
 Z-DISTANCES FROM      0.000000000     TO      1.000000000    

      OUTPUT VOLUME =     1

 DEFAULT ELEMENT DIVISIONS PER LINE BASED ON ELEMENT SIZE =   1.00    
 PRODUCE ALL HEXAHEDRAL ELEMENTS IN 3D.

 GENERATE NODES AND ELEMENTS   IN  ALL  SELECTED VOLUMES  

 NUMBER OF VOLUMES MESHED   =         1
 MAXIMUM NODE NUMBER        =         8
 MAXIMUM ELEMENT NUMBER     =         1

 LIST ALL SELECTED ELEMENTS.  (LIST NODES)

 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:40:08  APR 28, 2022 CP=      0.871

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).
 
    ELEM MAT TYP REL ESY SEC        NODES
 
       1   1   1   1   0   1      2     4     3     1     5     6     7     8

 SELECT       FOR ITEM=LOC  COMPONENT=X     BETWEEN  0.0000     AND   0.0000    
   KABS=  0.  TOLERANCE= 0.100000E-05

          4  NODES (OF          8  DEFINED) SELECTED BY  NSEL  COMMAND.

 SPECIFIED CONSTRAINT UX   FOR SELECTED NODES            1 TO           8 BY           1
 REAL=  0.00000000       IMAG=  0.00000000    

 SELECT       FOR ITEM=LOC  COMPONENT=Y     BETWEEN  0.0000     AND   0.0000    
   KABS=  0.  TOLERANCE= 0.100000E-05

          4  NODES (OF          8  DEFINED) SELECTED BY  NSEL  COMMAND.

 SPECIFIED CONSTRAINT UY   FOR SELECTED NODES            1 TO           8 BY           1
 REAL=  0.00000000       IMAG=  0.00000000    

 SELECT       FOR ITEM=LOC  COMPONENT=Z     BETWEEN  0.0000     AND   0.0000    
   KABS=  0.  TOLERANCE= 0.100000E-05

          4  NODES (OF          8  DEFINED) SELECTED BY  NSEL  COMMAND.

 SPECIFIED CONSTRAINT UZ   FOR SELECTED NODES            1 TO           8 BY           1
 REAL=  0.00000000       IMAG=  0.00000000    

 SELECT ALL ENTITIES OF TYPE= ALL  AND BELOW


 ***** ROUTINE COMPLETED *****  CP =         0.873



 *****  ANSYS SOLUTION ROUTINE  *****

 TIME=  1.0000    

 USE INITIAL TIME STEP SIZE OF 0.1000000     FOR ALL  DEGREES OF FREEDOM


 ALL INTEGRATION POINT VALUES WILL BE COPIED TO THE NODES

 SELECT       FOR ITEM=LOC  COMPONENT=X     BETWEEN  1.0000     AND   1.0000    
   KABS=  0.  TOLERANCE= 0.500000E-02

          4  NODES (OF          8  DEFINED) SELECTED BY  NSEL  COMMAND.

 GENERATE SURFACE LOAD PRES ON SURFACE DEFINED BY ALL SELECTED NODES
 VALUES=  1000.00000      0.00000000    

 NUMBER OF PRES ELEMENT FACE LOADS STORED =          1

 SELECT ALL ENTITIES OF TYPE= ALL  AND BELOW

 WRITE ALL  ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL 
   FOR ALL APPLICABLE ENTITIES

 *** NOTE ***                            CP =       1.084   TIME= 17:40:08
 The automatic domain decomposition logic has selected the MESH domain   
 decomposition method with 2 processes per solution.                     

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

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

 ELEMENT TYPE         1 IS SOLID185. IT IS ASSOCIATED WITH ELASTOPLASTIC MATERIALS
 ONLY. KEYOPT(2)=2 IS SUGGESTED.



 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:40:08  APR 28, 2022 CP=      1.087

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).

                       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)
   PLASTIC MATERIAL PROPERTIES INCLUDED. . . . . .YES
   NEWTON-RAPHSON OPTION . . . . . . . . . . . . .PROGRAM CHOSEN   
   GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC  

 *** NOTE ***                            CP =       1.169   TIME= 17:40:27
 This nonlinear analysis defaults to using the full Newton-Raphson       
 solution procedure.  This can be modified using the NROPT command.      

 *** NOTE ***                            CP =       1.169   TIME= 17:40:27
 The conditions for direct assembly have been met.  No .emat or .erot    
 files will be produced.                                                 

 *** NOTE ***                            CP =       1.170   TIME= 17:40:27
 The number of CPU processors specified is greater than the number of    
 elements selected in the model.  The number of Distributed ANSYS        
 processes will be reduced to 1 for this solution.                       

                      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    
   AUTOMATIC TIME STEPPING . . . . . . . . . . . .    ON
      STARTING TIME STEP SIZE. . . . . . . . . . . 0.10000    
      MINIMUM TIME STEP SIZE . . . . . . . . . . . 0.10000E-02
      MAXIMUM TIME STEP SIZE . . . . . . . . . . .  1.0000    
   MAXIMUM NUMBER OF EQUILIBRIUM ITERATIONS. . . .    15
   STEP CHANGE BOUNDARY CONDITIONS . . . . . . . .    NO
   TERMINATE ANALYSIS IF NOT CONVERGED . . . . . .YES (EXIT)  
   CONVERGENCE CONTROLS. . . . . . . . . . . . . .USE DEFAULTS
   COPY INTEGRATION POINT VALUES TO NODE . . . . .   YES
   PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
   DATABASE OUTPUT CONTROLS
      ITEM     FREQUENCY   COMPONENT
       ALL        ALL               


 SOLUTION MONITORING INFO IS WRITTEN TO FILE= usrmatjhddfg.mntr                                                                                                                                                                                                                                                   


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

        1         1  SOLID185      0.256   0.256310
 Time at end of element matrix formulation CP = 1.64710021.              

 ALL CURRENT ANSYS DATA WRITTEN TO FILE NAME= usrmatjhddfg.rdb
  FOR POSSIBLE RESUME FROM THIS POINT
     FORCE CONVERGENCE VALUE  =   50.00      CRITERION=  0.2551    

 *** WARNING ***                         CP =       1.688   TIME= 17:44:27
 There are no active degrees of freedom.                                 
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.5102E-13 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED

 *** WARNING ***                         CP =       1.692   TIME= 17:44:27
 A reference force value times the tolerance is used by the              
 Newton-Raphson method for checking convergence.  The calculated         
 reference FORCE CONVERGENCE VALUE = 0 is less than a threshold.  This   
 threshold is internally calculated.  You can overwrite it by            
 specifying MINREF on the CNVTOL command.  Check results carefully.      
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1

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

        1         1  SOLID185      0.001   0.001339

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

        1         1  SOLID185      0.000   0.000176
 *** LOAD STEP     1   SUBSTEP     1  COMPLETED.    CUM ITER =      1
 *** TIME =  0.100000         TIME INC =  0.100000    
 *** AUTO STEP TIME:  NEXT TIME INC = 0.10000      UNCHANGED

     FORCE CONVERGENCE VALUE  =   100.0      CRITERION=  0.5102    
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.1000E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     2  COMPLETED.    CUM ITER =      2
 *** TIME =  0.200000         TIME INC =  0.100000    
 *** AUTO TIME STEP:  NEXT TIME INC = 0.15000      INCREASED (FACTOR = 1.5000)

     FORCE CONVERGENCE VALUE  =   175.0      CRITERION=  0.8929    
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.1000E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     3  COMPLETED.    CUM ITER =      3
 *** TIME =  0.350000         TIME INC =  0.150000    
 *** AUTO TIME STEP:  NEXT TIME INC = 0.22500      INCREASED (FACTOR = 1.5000)

     FORCE CONVERGENCE VALUE  =   287.5      CRITERION=   1.467    
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.1000E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     4  COMPLETED.    CUM ITER =      4
 *** TIME =  0.575000         TIME INC =  0.225000    
 *** AUTO TIME STEP:  NEXT TIME INC = 0.33750      INCREASED (FACTOR = 1.5000)

     FORCE CONVERGENCE VALUE  =   456.2      CRITERION=   2.328    
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.1000E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     5  COMPLETED.    CUM ITER =      5
 *** TIME =  0.912500         TIME INC =  0.337500    
 *** AUTO TIME STEP:  NEXT TIME INC = 0.87500E-01  DECREASED (FACTOR = 0.2593)

     FORCE CONVERGENCE VALUE  =   500.0      CRITERION=   2.551    
     DISP CONVERGENCE VALUE   =   0.000      CRITERION=  0.1000E-04 <<< CONVERGED
    EQUIL ITER   1 COMPLETED.  NEW TRIANG MATRIX.  MAX DOF INC=   0.000    
     FORCE CONVERGENCE VALUE  =   0.000      CRITERION=  0.5206E-04 <<< CONVERGED
    >>> SOLUTION CONVERGED AFTER EQUILIBRIUM ITERATION   1
 *** LOAD STEP     1   SUBSTEP     6  COMPLETED.    CUM ITER =      6
 *** TIME =   1.00000         TIME INC =  0.875000E-01


 *** ANSYS BINARY FILE STATISTICS
  BUFFER SIZE USED= 16384
        0.062 MB WRITTEN ON ELEMENT SAVED DATA FILE: usrmatjhddfg.esav
        0.062 MB WRITTEN ON ASSEMBLED MATRIX FILE: usrmatjhddfg.full
        0.312 MB WRITTEN ON RESULTS FILE: usrmatjhddfg.rst

 FINISH SOLUTION PROCESSING

 *** NOTE ***                            CP =       1.735   TIME= 17:44:27
 Distributed parallel processing has been reactivated.                   


 ***** ROUTINE COMPLETED *****  CP =         1.735



 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:44:27  APR 28, 2022 CP=      1.741

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).

          ***** ANSYS RESULTS INTERPRETATION (POST1) *****

 *** NOTE ***                            CP =       1.742   TIME= 17:44:27
 Reading results into the database (SET command) will update the current 
 displacement and force boundary conditions in the database with the     
 values from the results file for that load set.  Note that any          
 subsequent solutions will use these values unless action is taken to    
 either SAVE the current values or not overwrite them (/EXIT,NOSAVE).    


 USE LAST SUBSTEP ON RESULT FILE  FOR LOAD CASE 0 

 SET COMMAND GOT LOAD STEP=     1  SUBSTEP=     6  CUMULATIVE ITERATION=     6
   TIME/FREQUENCY=  1.0000    
 TITLE= upf-py1s, test usermat.py with 3-D elements                             

 SELECT       FOR ITEM=ELEM COMPONENT=    
  IN RANGE         1 TO          1 STEP          1

          1  ELEMENTS (OF          1  DEFINED) SELECTED BY  ESEL  COMMAND.

 PRINT S    ELEMENT SOLUTION PER ELEMENT

 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:44:27  APR 28, 2022 CP=      1.750

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).
 
 ***** POST1 ELEMENT NODAL STRESS LISTING *****                                
 
  LOAD STEP=     1  SUBSTEP=     6                                             
   TIME=    1.0000      LOAD CASE=   0                                         
 
  THE FOLLOWING X,Y,Z VALUES ARE IN GLOBAL COORDINATES                         
 
  
  ELEMENT=       1        SOLID185
    NODE     SX           SY           SZ           SXY          SYZ          SXZ     
       2   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       4   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       3   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       1   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       5   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       6   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       7   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       8   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     

 PRINT EPEL ELEMENT SOLUTION PER ELEMENT

 *** ANSYS - ENGINEERING ANALYSIS SYSTEM  RELEASE 2022 R1          22.1     ***
 DISTRIBUTED Ansys Mechanical Enterprise                       

 00000000  VERSION=LINUX x64     17:44:27  APR 28, 2022 CP=      1.844

 upf-py1s, test usermat.py with 3-D elements                                   

 Note - This ANSYS version was linked by /UPF with   1 user supplied routine(s).
 
 ***** POST1 ELEMENT NODAL ELASTIC STRAIN LISTING *****                        
 
  LOAD STEP=     1  SUBSTEP=     6                                             
   TIME=    1.0000      LOAD CASE=   0                                         
 
  THE FOLLOWING X,Y,Z VALUES ARE IN GLOBAL COORDINATES                         
 
  
  ELEMENT=       1        SOLID185
    NODE     EPELX        EPELY        EPELZ        EPELXY       EPELYZ       EPELXZ  
       2   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       4   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       3   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       1   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       5   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       6   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       7   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
       8   0.0000       0.0000       0.0000       0.0000       0.0000       0.0000     
  expected results: Sx=-1000, epel_x=-1e-2

 EXIT THE ANSYS POST1 DATABASE PROCESSOR


 ***** ROUTINE COMPLETED *****  CP =         1.845



 EXIT ANSYS WITHOUT SAVING DATABASE


 NUMBER OF WARNING MESSAGES ENCOUNTERED=          2
 NUMBER OF ERROR   MESSAGES ENCOUNTERED=          0

+--------- D I S T R I B U T E D   A N S Y S   S T A T I S T I C S ------------+

Release: 2022 R1            Build: 22.1       Update: UP20211129   Platform: LINUX x64   
Date Run: 04/28/2022   Time: 17:44     Process ID: 3016
Operating System: SUSE Linux Enterprise Server 12 SP5

Processor Model: Intel(R) Xeon(R) Silver 4116 CPU @ 2.10GHz

Compiler: Intel(R) Fortran Compiler Version 19.0.0  (Build: 20190206)
          Intel(R) C/C++ Compiler Version 19.0.0  (Build: 20190206)
          Intel(R) Math Kernel Library Version 2020.0.0 Product Build 20191122
          BLAS Library supplied by Intel(R) MKL

Number of machines requested            :    1
Total number of cores available         :   16
Number of physical cores available      :   16
Number of processes requested           :    2
Number of threads per process requested :    1
Total number of cores requested         :    2 (Distributed Memory Parallel)               
MPI Type: INTELMPI
MPI Version: Intel(R) MPI Library 2018 Update 3 for Linux* OS


GPU Acceleration: Not Requested

Job Name: usrmatjhddfg0
Input File: /home/suse/Documents/upf_studies/usermat.dat

  Core                Machine Name   Working Directory
 -----------------------------------------------------
     0                        suse   /home/suse/Documents/upf_studies
     1                        suse   /home/suse/Documents/upf_studies
 
Latency time from master to core     1 =    2.118 microseconds
 
Communication speed from master to core     1 =  6369.92 MB/sec

Total CPU time for main thread                    :        1.2 seconds
Total CPU time summed for all threads             :        2.1 seconds

Elapsed time spent obtaining a license            :       14.2 seconds
Elapsed time spent pre-processing model (/PREP7)  :        0.0 seconds
Elapsed time spent solution - preprocessing       :       19.4 seconds
Elapsed time spent computing solution             :      239.2 seconds
Elapsed time spent solution - postprocessing      :        0.0 seconds
Elapsed time spent post-processing model (/POST1) :        0.1 seconds
 
Equation solver used                              :            Sparse (symmetric)

Sum of memory used on all processes               :       14.0 MB
Sum of memory allocated on all processes          :     3136.0 MB
Physical memory available                         :         34 GB

Total amount of I/O written to disk               :        0.0 GB
Total amount of I/O read from disk                :        0.0 GB

+------ E N D   D I S T R I B U T E D   A N S Y S   S T A T I S T I C S -------+


 *---------------------------------------------------------------------------*
 |                                                                           |
 |                       DISTRIBUTED ANSYS RUN COMPLETED                     |
 |                                                                           |
 |---------------------------------------------------------------------------|
 |                                                                           |
 | Ansys 2022 R1               Build 22.1         UP20211129     LINUX x64   |
 |                                                                           |
 |---------------------------------------------------------------------------|
 |                                                                           |
 | Database Requested(-db)  1024 MB    Scratch Memory Requested      1024 MB |
 | Maximum Database Used       1 MB    Maximum Scratch Memory Used      8 MB |
 |                                                                           |
 |---------------------------------------------------------------------------|
 |                                                                           |
 |        CP Time      (sec) =          2.053       Time  =  17:44:28        |
 |        Elapsed Time (sec) =        296.000       Date  =  04/28/2022      |
 |                                                                           |
 *---------------------------------------------------------------------------*

>> We shutdown Python Server(s) 

/ansys_inc/v221/commonfiles/launcherQt/source/../../../ansys/bin/mapdl: line 404: kill: (2802) - No such process

Answered by germa89

May 9, 2022

Hi @kkytg

From this:

PYTHON VERSION : 3.6
**Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'

Current thread 0x00007f0e4f203700 (most recent call first):
/ansys_inc/v221/ansys/bin/ansupf221: line 90:  2800 Aborted                 (core dumped) python${PYTHON_VER} -m pip install ${ANSYS_DIR}/syslib/ansGRPC/User/mapdl-pyupf.tar.gz --user > /dev/null
Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'**

It seems to me that the module encodings is not installed in the default python3.6 environment.

You can install that package using the following:

/usr/bin/py…

View full answer

FredAns · 2022-05-06T07:34:21Z

FredAns
May 6, 2022
Maintainer

Hi,
Issue comes from the fact that the Python server fails to start at the very beginning, so MAPDL gets no user routine to talk with, so has zeros for all request it has all along the solve process.
@germa89, did you ever succeed in running Python UPF in PyAnsys ?

1 reply

kkytg May 7, 2022
Author

Any other replies?? Would be really helpful...

Thanks in advance

germa89 · 2022-05-09T12:20:25Z

germa89
May 9, 2022
Maintainer

Hi @kkytg

From this:

PYTHON VERSION : 3.6
**Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'

Current thread 0x00007f0e4f203700 (most recent call first):
/ansys_inc/v221/ansys/bin/ansupf221: line 90:  2800 Aborted                 (core dumped) python${PYTHON_VER} -m pip install ${ANSYS_DIR}/syslib/ansGRPC/User/mapdl-pyupf.tar.gz --user > /dev/null
Fatal Python error: Py_Initialize: Unable to get the locale encoding
ModuleNotFoundError: No module named 'encodings'**

It seems to me that the module encodings is not installed in the default python3.6 environment.

You can install that package using the following:

/usr/bin/python3 -m pip install encodings

I believe the default python3 environment is (in most of linux distributions) in /usr/bin/. However, you might need to change that.

To know where is the default python3 executable, you can write this in a terminal

python3 -c "import sys;print(sys.exec_prefix)"

(Make sure you don't have any virtual environment activated)

After you have that fixed, you might have more error regarding missing modules (Numpy for example), just repeat the approach.

I hope it helps!

P.D. Reminder, you need to use batch mode (-b), only works on Linux, and you need to have intel compiler installed.

/ansys_inc/v221/ansys/bin/ansys221 -b -i input.inp -o out.out

2 replies

kkytg May 11, 2022
Author

Thanks for the help, really appreciated. Since I've tested positive I am not able to get to the office to give it a try, until next week. However, as I remember I could've installed encodings for the default version (which is python3.9 in Anaconda3) Ansys does not recognize that Python version, instead it tries to consider another python, which is python 3.7 somewhere in its own directory.

Hope to respond you soon, and thanks again!

germa89 May 12, 2022
Maintainer

No rush. I hope you recover well!

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Python UPF Problem #1078

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Replies: 2 comments 3 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

Python UPF Problem #1078

kkytg Apr 28, 2022

Replies: 2 comments · 3 replies

FredAns May 6, 2022 Maintainer

kkytg May 7, 2022 Author

germa89 May 9, 2022 Maintainer

kkytg May 11, 2022 Author

germa89 May 12, 2022 Maintainer

kkytg
Apr 28, 2022

Replies: 2 comments 3 replies

FredAns
May 6, 2022
Maintainer

kkytg May 7, 2022
Author

germa89
May 9, 2022
Maintainer

kkytg May 11, 2022
Author

germa89 May 12, 2022
Maintainer