abaqustools 2.2__tar.gz → 2.32__tar.gz

{abaqustools-2.2/abaqustools.egg-info → abaqustools-2.32}/PKG-INFO

@@ -1,20 +1,14 @@
 Metadata-Version: 2.1
 Name: abaqustools
-Version: 2.2
+Version: 2.32
 Summary: Tools for FEM software ABAQUS.
 Author-email: "O. W. Petersen" <oyvind.w.petersen@ntnu.no>
-Keywords: abaqus,suspension,bridge,fem,odb
+Keywords: abaqus,export,fem,odb
 Requires-Python: >=3
 License-File: LICENSE
 Requires-Dist: numpy
-Requires-Dist: os
-Requires-Dist: sys
-Requires-Dist: time
 Requires-Dist: scipy
 Requires-Dist: nexusformat
-Requires-Dist: datetime
 Requires-Dist: h5py
-Requires-Dist: timeit
-Requires-Dist: warnings
 Requires-Dist: putools
 Requires-Dist: ypstruct

{abaqustools-2.2 → abaqustools-2.32/abaqustools.egg-info}/PKG-INFO

@@ -1,20 +1,14 @@
 Metadata-Version: 2.1
 Name: abaqustools
-Version: 2.2
+Version: 2.32
 Summary: Tools for FEM software ABAQUS.
 Author-email: "O. W. Petersen" <oyvind.w.petersen@ntnu.no>
-Keywords: abaqus,suspension,bridge,fem,odb
+Keywords: abaqus,export,fem,odb
 Requires-Python: >=3
 License-File: LICENSE
 Requires-Dist: numpy
-Requires-Dist: os
-Requires-Dist: sys
-Requires-Dist: time
 Requires-Dist: scipy
 Requires-Dist: nexusformat
-Requires-Dist: datetime
 Requires-Dist: h5py
-Requires-Dist: timeit
-Requires-Dist: warnings
 Requires-Dist: putools
 Requires-Dist: ypstruct

{abaqustools-2.2 → abaqustools-2.32}/abaqustools.egg-info/SOURCES.txt

@@ -15,6 +15,7 @@ abaqustools.egg-info/dependency_links.txt
 abaqustools.egg-info/requires.txt
 abaqustools.egg-info/top_level.txt
 develop/Export_debug.py
+develop/Export_new.py
 develop/nordsenga/nordsenga_truss.py
 develop/test_build_gradually/Runmodel.py
 develop/test_retract/Langenuen_tower_exportstatic.py

{abaqustools-2.2 → abaqustools-2.32}/abaqustools.egg-info/requires.txt

@@ -1,12 +1,6 @@
 numpy
-os
-sys
-time
 scipy
 nexusformat
-datetime
 h5py
-timeit
-warnings
 putools
 ypstruct

abaqustools-2.32/develop/Export_new.py

@@ -0,0 +1,115 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Thu Feb  1 15:18:47 2024
+
+@author: oyvinpet
+"""
+
+# -*- coding: utf-8 -*-
+
+#%% Import packages
+
+import h5py    
+import numpy as np    
+import matplotlib.pyplot as plt
+import sys
+
+sys.path.append('C:/Cloud/OD_OWP/Work/Python/Github')
+
+from abaqustools import odbexport
+
+#%% Export from odb
+
+folder_odb=r'C:\Cloud\OD_OWP\Work\Abaqus\Hardanger\ProjectStudents'
+jobname='HardangerWindInitialNew'
+folder_save=folder_odb
+folder_python='C:/Cloud/OD_OWP/Work/Python/Github/abaqustools/odbexport'
+
+# Export modes to h5 file
+odbexport.export.modal(folder_odb,jobname,folder_save,folder_python)
+
+# Export static state for dead loads to h5 file
+# odbexport.export.static(folder_odb,jobname,folder_save,folder_python,stepnumber=-2,framenumber=-1)
+
+#%% Import modal properties
+
+hf=h5py.File(r'C:\Cloud\OD_OWP\Work\Abaqus\Hardanger\ProjectStudents\HardangerWindInitialNew_exportmodal.h5', 'r')
+
+# Frequencies
+f=np.array(hf.get('f'))
+
+# Generalized mass
+gm=np.array(hf.get('gm'))
+
+# Node coordinates all nodes (size [N_nodes,4])
+nodecoord=np.array(hf.get('nodecoord'))
+
+# Mode shape matrix for all DOFs (size [6*N_nodes,N_modes])
+phi=np.array(hf.get('phi'))
+
+# Labels corresponding to each row of phi (each DOF), as a list of strings
+phi_label_temp=np.array(hf.get('phi_label'))
+phi_label=phi_label_temp[:].astype('U10').ravel().tolist()
+
+# The DOFs are as follows:
+# U1=x=longitudinal translation (along bridge)
+# U2=y=horizontal translation
+# U3=z=vertical translation
+# UR1=rotation about x aka. torsion
+# UR2=rotation about y aka. vertical bending
+# UR3=rotation about z aka. horizontal bending
+
+#%% Get node coordinates and mode shape for bridge deck only
+
+# Nodes in bridge deck 
+node_deck=np.arange(1001,1311+1,1)
+
+# Create list index of nodes in bridge deck
+index_node_deck=[]
+
+for k in np.arange(len(node_deck)):
+    index_node_deck.append(np.argwhere(node_deck[k]==nodecoord[:,0])[0,0])
+   
+nodecoord_deck=nodecoord[index_node_deck,:]
+
+   
+# Create list of index of y-DOFs,z-DOFs, and t-DOFs in bridge deck
+index_y=[]
+index_z=[]
+index_t=[]
+
+for k in np.arange(len(node_deck)):
+    str_y=str(node_deck[k]) + '_U2'
+    index_y.append(phi_label.index(str_y))  
+    
+    str_z=str(node_deck[k]) + '_U3'
+    index_z.append(phi_label.index(str_z))  
+    
+    str_t=str(node_deck[k]) + '_UR1'
+    index_t.append(phi_label.index(str_t))
+    
+    
+phi_y=phi[index_y,:]
+phi_z=phi[index_z,:]
+phi_t=phi[index_t,:]
+
+
+#%% Plot single mode
+
+mode_plot=14 # Mode number to plot
+scale_factor=1 # Scale factor
+
+plt.figure()
+
+x=nodecoord_deck[:,1] #x-coordinate of deck nodes
+
+h1=plt.plot(x,phi_y[:,mode_plot]*scale_factor, label='Horizontal')
+h2=plt.plot(x,phi_z[:,mode_plot]*scale_factor, label='Vertical')
+h3=plt.plot(x,phi_t[:,mode_plot]*scale_factor, label='Torsion')
+
+plt.xlabel('x [m]')
+plt.ylabel('Modal deflection [m or rad]')
+
+plt.legend()
+
+plt.show()

{abaqustools-2.2 → abaqustools-2.32}/odbexport/export.py

@@ -4,6 +4,7 @@
 
 import sys
 import os
+import time
 import numpy as np
 import putools
 import h5py
@@ -24,7 +25,7 @@ def static(folder_odb,jobname,folder_save,folder_python,variables=None,stepnumbe
     # framenumber: framenumber(s) to export
     # deletetxt: delete txt files
     # saveh5: save to h5
-    # prefix: prefix for txt files    # saveh5: save to h5
+    # prefix: prefix for txt files
     # postfixh5: postfix for h5 file
     # exportscript: name of python script that is created
 
@@ -61,7 +62,7 @@ def modal(folder_odb,jobname,folder_save,folder_python,variables=None,stepnumber
     # framenumber: framenumber(s) to export
     # deletetxt: delete txt files
     # saveh5: save to h5
-    # prefix: prefix for txt files    # saveh5: save to h5
+    # prefix: prefix for txt files
     # postfixh5: postfix for h5 file
     # exportscript: name of python script that is created    # stepnumber
        
@@ -74,7 +75,7 @@ def modal(folder_odb,jobname,folder_save,folder_python,variables=None,stepnumber
     
     # Variables to export
     if variables is None:
-       variables=['f' , 'gm' , 'phi' , 'phi_sf' , 'nodecoord' , 'elconn' ]
+       variables=['f' , 'gm' , 'phi' , 'phi_sf' , 'nodecoord' , 'elconn']
 
     # Script py file
     if exportscript is None:
@@ -121,13 +122,23 @@ def exportmain(folder_odb,jobname,folder_save,folder_python,variables,stepnumber
     if not file_exist_logic:
         print('***** ' + file_name)
         raise Exception('***** ODB file not found')
+        
+    # Check if python folder and python files are given correctly
+    file_name=folder_python + '\\' + 'odbfunc.py'
+    file_exist_logic=os.path.isfile(file_name)
+    if not file_exist_logic:
+        print('***** ' + folder_python)
+        raise Exception('***** Supplied python (odbexport) folder and odbfunc module not found')    
  
     # Ensure list
     if isinstance(variables,str):
         variables=[variables]
- 
+        
+    # Lower case
+    variables=[x.lower() for x in variables]
+
     # Check input variables
-    variables_allowed= ['u' , 'sf' , 'f' , 'gm' , 'phi' , 'phi_sf' , 'nodecoord' , 'elconn' , 'elconn' , 'elset' ]
+    variables_allowed= ['u' , 'sf' , 'phi' , 'phi_sf' , 'f' , 'gm' , 'nodecoord' , 'elconn' , 'elset' ]
     for k in np.arange(len(variables)):
         if not variables[k] in variables_allowed:
             print('***** Check variable ' + variables[k])
@@ -140,66 +151,113 @@ def exportmain(folder_odb,jobname,folder_save,folder_python,variables,stepnumber
     system_cmd='abaqus cae noGUI=' + folder_save + '/' + exportscript + '.py'
     sys_out=os.popen(system_cmd).read()
     
-    # Collect data for h5
-    h5_var=[]
-    h5_data=[]
-    h5_isnum=[]
-    
-    for k in np.arange(len(variables)):
-        
-        # Base for txt files
-        filename_base=folder_save + '/' + prefix
-        
-        # Read txt file
-        data_temp=np.genfromtxt(filename_base+variables[k]+'.txt', delimiter=',')
-        h5_var.append(variables[k])
-        h5_data.append(data_temp)
-        h5_isnum.append(True)
-                    
-        # Special cases where labels also need be imported
-        if variables[k]=='phi':
-            fid=open(filename_base+'phi_label.txt','r')
-            phi_label=fid.read().splitlines()
-            fid.close()
-            h5_var.append('phi_label')
-            h5_data.append(phi_label)
-            h5_isnum.append(False)
-        elif variables[k]=='phi_sf':
-            fid=open(filename_base+'phi_sf_label.txt','r')
-            phi_sf_label=fid.read().splitlines()
-            fid.close()
-            h5_var.append('phi_sf_label')
-            h5_data.append(phi_sf_label)
-            h5_isnum.append(False)
-        elif variables[k]=='u':
-            fid=open(filename_base+'u_label.txt','r')
-            u_label=fid.read().splitlines()
-            fid.close()
-            h5_var.append('u_label')
-            h5_data.append(u_label)
-            h5_isnum.append(False)
-        elif variables[k]=='sf':
-            fid=open(filename_base+'sf_label.txt','r')
-            sf_label=fid.read().splitlines()
-            fid.close()
-            h5_var.append('sf_label')
-            h5_data.append(sf_label)
-            h5_isnum.append(False)
-        elif variables[k]=='elset':
-            fid=open(filename_base+'elset_label.txt','r')
-            elset_label=fid.read().splitlines()
-            fid.close()
-            h5_var.append('elset_label')
-            h5_data.append(elset_label)
-            h5_isnum.append(False)
+    idx_error=sys_out.find('error')
 
+    if idx_error>0:
+        print(sys_out)
+        raise Exception('***** Export aborted due to errors, see above')
+            
+                
     # Save h5 file        
     if saveh5==True:
-        
-        #if os.path.exists(hf_name):
+
+        # Base for txt files
         hf_name=folder_save + '/' + jobname + postfixh5 + '.h5'
-        exporth5(hf_name,h5_var,h5_data,h5_isnum)
         
+        # Overwrite file
+        if os.path.exists(hf_name):
+            os.remove(hf_name)
+            time.sleep(1)
+            
+        hf = h5py.File(hf_name,'w')
+        dt = h5py.special_dtype(vlen=str)
+            
+        for k in np.arange(len(variables)):
+            
+            # Base for txt files
+            filename_base=folder_save + '/' + prefix      
+            
+            # Read txt file
+            data_from_txt=np.genfromtxt(filename_base+variables[k]+'.txt', delimiter=',')
+            
+            
+            if variables[k]=='u':
+            
+                hf.create_dataset('u',data=data_from_txt)
+                hf['u'].attrs.modify('Type', 'Displacement')
+                hf['u'].attrs.modify('Unit', 'm, rad')
+                
+                fid=open(filename_base+'u_label.txt','r'); label=fid.read().splitlines(); fid.close()  
+                hf.create_dataset('u_label',data=label,dtype=dt)
+                hf['u_label'].attrs.modify('Type', 'DOF labels')
+                
+            elif variables[k]=='sf':
+            
+                hf.create_dataset('sf',data=data_from_txt)
+                hf['sf'].attrs.modify('Type', 'Section force')
+                hf['sf'].attrs.modify('Unit', 'm, rad')
+                
+                fid=open(filename_base+'sf_label.txt','r'); label=fid.read().splitlines(); fid.close()  
+                hf.create_dataset('sf_label',data=label,dtype=dt)
+                hf['sf_label'].attrs.modify('Type', 'Element labels')
+                
+            elif variables[k]=='phi':
+            
+                hf.create_dataset('phi',data=data_from_txt)
+                hf['phi'].attrs.modify('Type', 'Modal displacement')
+                hf['phi'].attrs.modify('Unit', 'm, rad')
+                
+                fid=open(filename_base+'phi_label.txt','r'); label=fid.read().splitlines(); fid.close()  
+                hf.create_dataset('phi_label',data=label,dtype=dt)
+                hf['phi_label'].attrs.modify('Type', 'DOF labels')
+                
+            elif variables[k]=='phi_sf':
+            
+                hf.create_dataset('phi_sf',data=data_from_txt)
+                hf['phi_sf'].attrs.modify('Type', 'Modal section force')
+                hf['phi_sf'].attrs.modify('Unit', 'N, Nm')
+                
+                fid=open(filename_base+'phi_sf_label.txt','r'); label=fid.read().splitlines(); fid.close()  
+                hf.create_dataset('phi_sf_label',data=label,dtype=dt)
+                hf['phi_sf_label'].attrs.modify('Type', 'Element labels')        
+                
+            elif variables[k]=='f':
+            
+                hf.create_dataset('f',data=data_from_txt)
+                hf['f'].attrs.modify('Type', 'Natural frequency')
+                hf['f'].attrs.modify('Unit', 'Hz')
+                
+            elif variables[k]=='gm':
+            
+                hf.create_dataset('gm',data=data_from_txt)
+                hf['gm'].attrs.modify('Type', 'Generalized mass')
+                hf['gm'].attrs.modify('Unit', 'kg')
+                
+            elif variables[k]=='nodecoord':
+            
+                hf.create_dataset('nodecoord',data=data_from_txt)
+                hf['nodecoord'].attrs.modify('Type', 'Node coordinate')
+                hf['nodecoord'].attrs.modify('Unit', 'Node number, m')     
+                
+            elif variables[k]=='elconn':
+            
+                hf.create_dataset('elconn',data=data_from_txt)
+                hf['elconn'].attrs.modify('Type', 'Element connectivity')
+                hf['elconn'].attrs.modify('Unit', 'Element number, Node number')                     
+                
+            elif variables[k]=='elset':
+           
+                hf.create_dataset('elset',data=data_from_txt)
+                hf['elset'].attrs.modify('Type', 'Element sets')
+                hf['elset'].attrs.modify('Unit', 'Element number')      
+           
+                fid=open(filename_base+'elset_label.txt','r'); label=fid.read().splitlines(); fid.close()  
+                hf.create_dataset('elset_label',data=label,dtype=dt)
+                hf['elset_label'].attrs.modify('Type', 'Element labels')    
+                            
+
+    hf.close()
+    
     # Delete txt files identified by prefix
     if deletetxt==True:
         deletefiles(folder_save,prefix,['.txt'])
@@ -208,6 +266,8 @@ def exportmain(folder_odb,jobname,folder_save,folder_python,variables,stepnumber
 
 def exporth5(hf_name,h5_var,h5_data,h5_isnum):
     
+    # NOT USED
+    
     # Save h5 file with data
     #
     # Inputs:
@@ -239,17 +299,17 @@ def deletefiles(foldername,name_match,extensions):
     # extensions: list of extensions allowed
     
     # Find files that match
-    FileNameListDir=os.listdir(foldername)
-    IndexMatch=putools.num.listindexsub(FileNameListDir,name_match)
+    filenames_dir=os.listdir(foldername)
+    idx_match=putools.num.listindexsub(filenames_dir,name_match)
     
-    for k in np.arange(len(IndexMatch)):
+    for k in np.arange(len(idx_match)):
             
-        FileNameRemove=FileNameListDir[IndexMatch[k]]
+        filename_remove=filenames_dir[idx_match[k]]
         
-        # Only selected files
-        match_logic=[list_sub in FileNameRemove for list_sub in extensions]
+        # Only files with specified extensions
+        match_logic=[list_sub in filename_remove for list_sub in extensions]
         if sum(match_logic)>0:
-            os.remove(foldername + '/' + FileNameRemove)
+            os.remove(foldername + '/' + filename_remove)
 
 #%% Export modal script
 

{abaqustools-2.2 → abaqustools-2.32}/odbexport/odbfunc.py

@@ -10,7 +10,7 @@ def p(pstring):
     prettyPrint(pstring)
 
 # This module is imported in Abaqus by "import odbfunc"
-# The folder must first be added to the system path
+# The folder containing odbfunc.py must first be added to the system path
 
 #%%
 
@@ -157,7 +157,8 @@ def exportnodecoord(odb_id,stepnumber,framenumber):
     if 'COORD' not in SelectedFrame.fieldOutputs.keys():
         NodeCoordMatrix=np.array([0 , 0 , 0])
         NodeCoordLabelVector=np.array([0])
-        return NodeCoordMatrix,NodeCoordLabelVector
+        NodeCoord=np.hstack((NodeCoordLabelVector,NodeCoordMatrix))
+        return NodeCoord
         
     t_start=timer()
     CoordValues=SelectedFrame.fieldOutputs['COORD'].values

{abaqustools-2.2 → abaqustools-2.32}/odbexport_example/Export.py

@@ -89,7 +89,7 @@ phi_t=phi[index_t,:]
 
 #%% Plot single mode
 
-mode_plot=14 # Mode number to plot
+mode_plot=0 # Mode number to plot
 scale_factor=1 # Scale factor
 
 plt.figure()

abaqustools-2.32/pyproject.toml

@@ -0,0 +1,15 @@
+[build-system]
+requires = ["setuptools"]
+build-backend = "setuptools.build_meta"
+
+[tool.setuptools.packages.find]
+exclude = ["develop/","__old__/","__pycache__/*",".idea/*",".git/*",".github/*"]
+
+[project]
+name = "abaqustools"
+version = "2.32" # Follow PEP440
+authors = [{ name="O. W. Petersen", email="oyvind.w.petersen@ntnu.no" },]
+description = "Tools for FEM software ABAQUS."
+requires-python = ">=3"
+keywords = ["abaqus", "export","fem","odb"]
+dependencies = ["numpy","scipy","nexusformat","h5py","putools","ypstruct"]

abaqustools-2.2/pyproject.toml

@@ -1,15 +0,0 @@
-[build-system]
-requires = ["setuptools"]
-build-backend = "setuptools.build_meta"
-
-[tool.setuptools.packages.find]
-exclude = ["develop/*","__old__/*","__pycache__/*",".idea/*",".git/*",".github/*"]
-
-[project]
-name = "abaqustools"
-version = "2.2" # Follow PEP440
-authors = [{ name="O. W. Petersen", email="oyvind.w.petersen@ntnu.no" },]
-description = "Tools for FEM software ABAQUS."
-requires-python = ">=3"
-keywords = ["abaqus", "suspension","bridge","fem","odb"]
-dependencies = ["numpy","os","sys","time","scipy","nexusformat","datetime","h5py","timeit","warnings","putools","ypstruct"]