midas-civil 1.0.8__tar.gz → 1.1.0__tar.gz

{midas_civil-1.0.8 → midas_civil-1.1.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: midas_civil
-Version: 1.0.8
+Version: 1.1.0
 Summary: Python library for MIDAS Civil NX
 Home-page: https://github.com/MIDASIT-Co-Ltd/midas-civil-python
 Author: Sumit Shekhar
@@ -14,6 +14,7 @@ Requires-Dist: polars
 Requires-Dist: xlsxwriter
 Requires-Dist: requests
 Requires-Dist: scipy
+Requires-Dist: colorama
 Dynamic: author
 Dynamic: author-email
 Dynamic: description

midas_civil-1.1.0/midas_civil/__init__.py

@@ -0,0 +1,48 @@
+import requests
+from colorama import Fore,Style
+from ._mapi import *
+_version_ = "1.1.0"
+
+
+print('\n╭────────────────────────────────────────────────────────────────────────────────────╮')
+print(Style.BRIGHT+f'│                      MIDAS CIVIL-NX PYTHON LIBRARY v{_version_}  🐍                      │')
+print('╰────────────────────────────────────────────────────────────────────────────────────╯\n'+Style.RESET_ALL)
+
+if NX.version_check:
+    resp =  requests.get("https://pypi.org/pypi/midas_civil/json").json()
+    latest_ver =  resp["info"]["version"]
+    if _version_ != latest_ver:        
+        print(Fore.YELLOW +'╭─ ⚠️   ──────────────────────────────────────────────────────────────────────────────╮')
+        print(f"│    Warning: You are using v{_version_}, but the latest available version is v{latest_ver}.      │")
+        print(f"│    Run 'pip install midas_civil --upgrade' to update.                              │")
+        print('╰────────────────────────────────────────────────────────────────────────────────────╯\n'+Style.RESET_ALL)
+
+
+from ._model import *
+from ._boundary import *
+from ._utils import *
+from ._node import *
+from ._element import *
+from ._load import *
+from ._group import *
+from ._result import *
+
+#--- TESTING IMPORTS ---
+from ._material import *
+
+# from ._section import *
+from ._section import *
+
+from ._result_extract import *
+from ._construction import *
+from ._thickness import *
+from ._temperature import *
+
+from ._tendon import *
+from ._view import *
+
+from ._movingload import *
+from ._settlement import *
+from ._analysiscontrol import *
+from ._BoundaryChangeAssignment import*
+

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_element.py

@@ -8,7 +8,7 @@ from math import hypot
 import math
 
 
-def _interpolateAlignment(pointsArray,n_seg=10,deg=1,mSize=0,includePoint:bool=True,yEcc=0,zEcc=0) -> list:
+def _interpolateAlignment(pointsArray,n_seg=10,deg=1,mSize=0,includePoint:bool=True) -> list:
     ''' Returns point list and beta angle list'''
     pointsArray = np.array(pointsArray)
     x_p, y_p , z_p  = pointsArray[:,0] , pointsArray[:,1] , pointsArray[:,2]
@@ -110,7 +110,7 @@ def _triangleAREA(a:Node,b:Node,c:Node):
     mag = np.linalg.norm(np.cross(v1, v2))
     return float(0.5 * mag) , np.cross(v1, v2)/mag
 
-def _calcVector(deltaLocation):
+def _calcVector(deltaLocation,angle=0): # Returns normalised local X,Y,Z for line
     Z_new = np.array([0.000001,0,1])
     X_new = np.array(deltaLocation)
     Y_new = np.cross(Z_new, X_new)
@@ -121,22 +121,45 @@ def _calcVector(deltaLocation):
     Y_new = Y_new / (np.linalg.norm(Y_new)+0.000001)
     Z_new = Z_new / (np.linalg.norm(Z_new)+0.000001)
 
-    return [X_new,Y_new,Z_new]
+    from scipy.spatial.transform import Rotation as R
+    p_y = np.array(Y_new)
+    p_z = np.array(Z_new)
+
+    axis = np.array(X_new)
+    theta = np.deg2rad(angle)               # or radians directly
+    rot = R.from_rotvec(axis * theta)           # axis-angle as rotation vector
+    
+    rt_y = rot.apply(p_y)                         # rotated point around axis through origin
+    rt_z = rot.apply(p_z)  
+    
+    return [X_new,rt_y,rt_z]
+
+def _rotatePT(pt,axis,deg):
+    from scipy.spatial.transform import Rotation as R
+    p = np.array(pt)
+    axis = np.array(axis)
+    theta = np.deg2rad(deg)               # or radians directly
+    rot = R.from_rotvec(axis * theta)           # axis-angle as rotation vector
+    return rot.apply(p)                         # rotated point around axis through origin
+
+def _pointOffset(pts,yEcc=0,zEcc=0,angle=0):
+    from ._utils import _matchArray
+
+    angle2 = _matchArray(pts,angle)
 
-def _pointOffset(pts,yEcc,zEcc):
     norm = []
-    norm.append(_calcVector(np.subtract(pts[1],pts[0])))
+    norm.append(_calcVector(np.subtract(pts[1],pts[0]),angle2[0]))    # first X- along vector
 
-    for i in range(len(pts)-2):
+    for i in range(len(pts)-2): # Averaged X- along vector for middle points
         X_new1 = np.array(np.subtract(pts[i+1],pts[i]))
         X_new2 = np.array(np.subtract(pts[i+2],pts[i+1]))
 
         X_new1 = X_new1 / (np.linalg.norm(X_new1)+0.000001)
         X_new2 = X_new2 / (np.linalg.norm(X_new2)+0.000001)
 
-        norm.append(_calcVector(np.add(X_new1,X_new2)))
+        norm.append(_calcVector(np.add(X_new1,X_new2),angle2[i+1]))
 
-    norm.append(_calcVector(np.subtract(pts[-1],pts[-2])))
+    norm.append(_calcVector(np.subtract(pts[-1],pts[-2]),angle2[-1])) # last X- along vector
 
     # print(norm)
 
@@ -329,6 +352,12 @@ class Element:
         Element.elements = []
         Element.ids = []
         Element.__elemDIC__={}
+    
+    @staticmethod
+    def clear():
+        Element.elements = []
+        Element.ids = []
+        Element.__elemDIC__={}
 
     # --- Element Type Subclasses ---
 
@@ -440,6 +469,10 @@ class Element:
                     i_loc = points_loc
                 else:
                     i_loc = _interpolateAlignment(points_loc,n_div,deg,0,includePoint)
+
+                from ._utils import _matchArray
+                angle = _matchArray(i_loc,angle)
+
                 for i in i_loc:
                     Enode=Node(i[0],i[1],i[2])
                     beam_nodes.append(Enode.ID)
@@ -447,29 +480,36 @@ class Element:
                 for i in range(len(i_loc)-1):
                     if id == 0 : id_new = 0
                     else: id_new = id+i
-                    beam_obj.append(Element.Beam(beam_nodes[i],beam_nodes[i+1],mat,sect,angle,group,id_new,bLocalAxis))
+                    beam_obj.append(Element.Beam(beam_nodes[i],beam_nodes[i+1],mat,sect,angle[i],group,id_new,bLocalAxis))
                 
                 return beam_obj
         
         @staticmethod
         def PLine2(points_loc:list,n_div:int=0,deg:int=1,includePoint:bool=True,mat:int=1,sect:int=1,angle:float=0, group = "" , id: int = 0,bLocalAxis=False,yEcc=0,zEcc=0):
-
+                
+                
                 beam_nodes =[]
                 beam_obj = []
                 if n_div == 0 :
                     i_loc = points_loc
                 else:
-                    i_loc = _interpolateAlignment(points_loc,n_div,deg,0,includePoint,yEcc,zEcc)
+                    i_loc = _interpolateAlignment(points_loc,n_div,deg,0,includePoint)
+                
+                from ._utils import _matchArray
+                angle = _matchArray(i_loc,angle)
 
-                i_loc2 = _pointOffset(i_loc,yEcc,zEcc)
+                i_loc2 = _pointOffset(i_loc,yEcc,zEcc,angle)
                 for i in i_loc2:
                     Enode=Node(i[0],i[1],i[2])
                     beam_nodes.append(Enode.ID)
+
                 
+                
+
                 for i in range(len(i_loc2)-1):
                     if id == 0 : id_new = 0
                     else: id_new = id+i
-                    beam_obj.append(Element.Beam(beam_nodes[i],beam_nodes[i+1],mat,sect,angle,group,id_new,bLocalAxis))
+                    beam_obj.append(Element.Beam(beam_nodes[i],beam_nodes[i+1],mat,sect,angle[i],group,id_new,bLocalAxis))
                 
                 return beam_obj
 

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_mapi.py

@@ -1,17 +1,18 @@
 import requests
 import sys
+from colorama import Fore, Style
 try:import winreg
 except: pass
+import time
 
 
 
 
 def Midas_help():
     """MIDAS Documnetation : https://midas-rnd.github.io/midasapi-python """
-    print("---"*22)
-    print("|   HELP MANUAL : https://midas-rnd.github.io/midasapi-python/   |")
-    print("---"*22,"\n")
-
+    print('\n╭────────────────────────────────────────────────────────────────────────────────────╮')
+    print("│         HELP MANUAL   :     https://midas-rnd.github.io/midasapi-python/           │")
+    print('╰────────────────────────────────────────────────────────────────────────────────────╯\n')
 
 
 
@@ -94,7 +95,6 @@ class MAPI_KEY:
                 my_key = value[0]
                 MAPI_KEY(my_key)
                 print(f' 🔑   MAPI-KEY is taken from Registry entry.  >>  {my_key[:35]}...')
-                print("-"*85)
             except:
                 print(f"🔑   MAPI KEY is not defined. Click on Apps > API Settings to copy the MAPI Key.\n Define it using MAPI_KEY('xxxx')")
                 sys.exit(0)
@@ -125,6 +125,9 @@ def MidasAPI(method:str, command:str, body:dict={})->dict:
         "MAPI-Key": mapi_key
     }
 
+    start_time = time.perf_counter()
+
+
     if method == "POST":
         response = requests.post(url=url, headers=headers, json=body)
     elif method == "PUT":
@@ -134,16 +137,24 @@ def MidasAPI(method:str, command:str, body:dict={})->dict:
     elif method == "DELETE":
         response = requests.delete(url=url, headers=headers)
 
-    if response.status_code == 404: 
-        print(f"⚠️  Civil NX model is not connected.  Click on 'Apps> Connect' in Civil NX. \nMake sure the MAPI Key in python code is matching with the MAPI key in Civil NX.\n\n")
-        sys.exit(0)
+    end_time = time.perf_counter()
+    elapsed_time = end_time - start_time
 
     if NX.debug_print:
-        print(method, command, response.status_code , "✅")
+        print(f"{method} |  URL : {command} | STATUS : {response.status_code} | TIME : {elapsed_time:.4f} sec")
+
+
 
     if MAPI_KEY.count == 1:
+        MAPI_KEY.count = 0
+        if response.status_code == 404:
+            print(Fore.RED +'\n╭─ 💀   ─────────────────────────────────────────────────────────────────────────────╮')
+            print(f"│  Civil NX model is not connected.  Click on 'Apps > Connect' in Civil NX.          │")
+            print(f"│  Make sure the MAPI Key in python code is matching with the MAPI key in Civil NX.  │")
+            print('╰────────────────────────────────────────────────────────────────────────────────────╯\n'+Style.RESET_ALL)
+            sys.exit(0)
+
         if NX.user_print:
-            MAPI_KEY.count = 0
             _checkUSER()
     
 
@@ -168,8 +179,25 @@ def _setUNIT(unitJS):
 def _checkUSER():
     try:
         resp =  MidasAPI('GET','/config/ver',{})['VER']
-        print(f"Connected to {resp['NAME']}   🟢")
-        print(f"USER : {resp['USER']}    |   COMPANY : {resp['COMPANY']}")
-        print("-"*85)
+
+        # print(f"{' '*15}Connected to {resp['NAME']}")
+        # print(f"{' '*15}USER : {resp['USER']}          COMPANY : {resp['COMPANY']}")
+
+        ln1 = f"Connected to {resp['NAME']}"
+        ln2 = f"USER : {resp['USER']}          COMPANY : {resp['COMPANY']}"
+
+        lg_ln1 = 63-len(ln1)
+        lg_ln2 = 67-len(ln2)
+
+        line1 = f"│{' '*15} {ln1} {' '*lg_ln1} 🟢 │"
+        line2 = f"│{' '*15} {ln2} {' '*lg_ln2}│"
+        print(Fore.GREEN+'\n╭─ 🔔  ──────────────────────────────────────────────────────────────────────────────╮')
+        print(line1)
+        print(line2)
+        print('╰────────────────────────────────────────────────────────────────────────────────────╯\n'+Style.RESET_ALL)
+
+
+        # print('─'*86)
+
     except:
-        print("")
+        pass

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_material.py

@@ -1620,6 +1620,10 @@ class TDMatLink:
         MidasAPI("DELETE","/db/TMAT")
         TDMatLink.mats={}
 
+    @staticmethod
+    def clear():
+        TDMatLink.mats={}
+
     @staticmethod
     def sync():
         a = TDMatLink.get()

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_node.py

@@ -152,6 +152,13 @@ class Node:
         Node.Grid={}
         Node.__nodeDic__ = {}
 
+    @staticmethod
+    def clear():
+        Node.nodes=[]
+        Node.ids=[]
+        Node.Grid={}
+        Node.__nodeDic__ = {}
+
 
 
 

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_result_extract.py

@@ -135,7 +135,7 @@ class Result :
                 "TABLE_NAME": tableName,
                 "STYLES": {
                     "FORMAT": "Fixed",
-                    "PLACE": 12
+                    "PLACE": 5
                 }
             }
         }
@@ -162,7 +162,7 @@ class Result :
 
 
 
-    # ---------- Result TABLE ------------------------------
+    # ---------- Result TABLE (For ALL TABLES)------------------------------
     @staticmethod
     def ResultTable(tabletype:str,keys=[],loadcase:list=[],cs_stage=[],force_unit='KN',len_unit='M'):
         '''
@@ -177,7 +177,7 @@ class Result :
                 "TABLE_TYPE": tabletype,
                 "STYLES": {
                     "FORMAT": "Fixed",
-                    "PLACE": 12
+                    "PLACE": 5
                 }
             }
         }
@@ -204,3 +204,57 @@ class Result :
         ss_json = MidasAPI("POST","/post/table",js_dat)
         _setUNIT(currUNIT)
         return _JSToDF_ResTable(ss_json)
+    
+
+    class TABLE :
+        @staticmethod
+        def BeamForce_VBM(keys=[],loadcase:list=[],items=['all'],parts=["PartI", "PartJ"],components=['all'],force_unit='KN',len_unit='M'):
+            '''
+                Keys : List{int} -> Element/ Node IDs  |  str -> Structure Group Name
+                Loadcase : Loadcase name followed by type. eg. DeadLoad(ST)
+                Items to display : [ "Axial" , "Shear-y" , "Shear-z" , "Torsion" , "Moment-y" , "Moment-z"]
+                Parts : ["PartI", "Part1/4", "Part2/4", "Part3/4", "PartJ"]
+                Components (colms of tabulart result): [ "Elem", "Load", "Part", "Component", "Axial", "Shear-y", "Shear-z", "Torsion", "Moment-y", "Moment-z" ]
+                
+            '''
+
+            js_dat = {
+                "Argument": {
+                    "TABLE_NAME": "SS_Table",
+                    "TABLE_TYPE": "BEAMFORCEVBM",
+                    "STYLES": {
+                        "FORMAT": "Fixed",
+                        "PLACE": 5
+                    },
+                    "PARTS" : parts
+                }
+            }
+
+
+            if isinstance(keys,list):
+                if keys!=[]:
+                    js_dat["Argument"]['NODE_ELEMS'] = {"KEYS": keys}
+            elif isinstance(keys,str):
+                js_dat["Argument"]['NODE_ELEMS'] = {"STRUCTURE_GROUP_NAME": keys}
+
+
+            if loadcase!=[]: js_dat["Argument"]['LOAD_CASE_NAMES'] = loadcase
+
+            if components!=['all']:
+                if "Elem" not in components: components.append("Elem")
+                if "Load" not in components: components.append("Load")
+                if "Part" not in components: components.append("Part")
+                if "Component" not in components: components.append("Component")
+                js_dat["Argument"]['COMPONENTS'] = components
+            
+            if items!=['all']:
+                js_dat["Argument"]['ITEM_TO_DISPLAY'] = items
+
+
+
+            currUNIT = _getUNIT()
+            Model.units(force=force_unit,length=len_unit)
+            ss_json = MidasAPI("POST","/post/table",js_dat)
+            _setUNIT(currUNIT)
+            return _JSToDF_ResTable(ss_json)
+

midas_civil-1.1.0/midas_civil/_section/_TapdbSecSS.py

@@ -0,0 +1,175 @@
+from ._offsetSS import Offset
+from ._offsetSS import _common
+import math
+
+class _SS_TAPERED_DBUSER(_common):
+
+    """ Create Standard USER DEFINED sections"""
+
+    def __init__(self,Name='',Shape='',params_I:list=[],params_J:list=[],Offset=Offset(),useShear=True,use7Dof=False,id:int=0):  
+        """ Shape = 'SB' 'SR' for rectangle \n For cylinder"""
+        self.ID = id
+        self.NAME = Name
+        self.TYPE = 'TAPERED'
+        self.SHAPE = Shape
+        self.PARAMS_I = params_I
+        self.PARAMS_J = params_J
+        self.OFFSET = Offset
+        self.USESHEAR = useShear
+        self.USE7DOF = use7Dof
+        self.DATATYPE = 2
+    
+    def __str__(self):
+         return f'  >  ID = {self.ID}   |  USER DEFINED STANDARD SECTION \nJSON = {self.toJSON()}\n'
+
+
+    def toJSON(sect):
+        js =  {
+                "SECTTYPE": sect.TYPE,
+                "SECT_NAME": sect.NAME,
+                "SECT_BEFORE": {
+                    "SHAPE": sect.SHAPE,
+                    "TYPE": sect.DATATYPE,
+                    "SECT_I": {
+                        "vSIZE": sect.PARAMS_I
+                    },
+                    "SECT_J": {
+                        "vSIZE": sect.PARAMS_J
+                    }
+                }
+            }
+        js['SECT_BEFORE'].update(sect.OFFSET.JS)
+        js['SECT_BEFORE']['USE_SHEAR_DEFORM'] = sect.USESHEAR
+        js['SECT_BEFORE']['USE_WARPING_EFFECT'] = sect.USE7DOF
+        return js
+    
+    @staticmethod
+    def _objectify(id,name,type,shape,offset,uShear,u7DOF,js):
+        return _SS_TAPERED_DBUSER(name,shape,js['SECT_BEFORE']['SECT_I']['vSIZE'],js['SECT_BEFORE']['SECT_J']['vSIZE'],offset,uShear,u7DOF,id)
+    
+    def _centerLine(shape,end):
+        if end: 
+            shape.PARAMS = shape.PARAMS_J
+            # print(' J end taken')
+        else: 
+            # print(' I end taken')
+            shape.PARAMS = shape.PARAMS_I
+
+        if shape.SHAPE == 'SB' :
+            H,B = shape.PARAMS[:2]
+
+            sect_lin_con = [[1,2],[3,1]]
+
+            sect_cg_LT = [-B/2,H/2]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [B/2,-H/2]
+
+            if H > B :
+                sect_shape = [[0,0],[0,H/2],[0,-H/2],[0,H/4],[0,-H/4]]
+                sect_thk = [B,B,B,B]
+                sect_thk_off = [0,0,0,0]
+            else : 
+                sect_shape = [[0,0],[B/2,0],[-B/2,0],[B/4,0],[-B/4,0]]
+                sect_thk = [H,H,H,H]
+                sect_thk_off = [0,0,0,0]
+
+        elif shape.SHAPE == 'L' :
+            H,B,tw,tf = shape.PARAMS[:4]
+
+            sect_cg_LT = [0,0]
+            sect_cg_CC = [(H*tw*tw+B*B*tf)/(2*(B*tw+H*tf)),-(H*H*tw+B*tf*tf)/(2*(B*tw+H*tf))]
+            sect_cg_RB = [B,-H]
+
+            # sect_shape = [[0.5*tw,-H],[0.5*tw,-0.5*tf],[B,-0.5*tf]]
+            sect_shape = [[0,-H],[0,0],[B,0]]
+            sect_lin_con = [[3,2],[2,1]]
+            sect_thk = [tw,tf]
+            # sect_thk_off = [0,0]
+            sect_thk_off = [tw/2,tf/2]
+        
+        elif shape.SHAPE == 'C' :
+            H,B1,tw,tf1,B2,tf2 = shape.PARAMS[:6]
+            if B2 == 0 : B2 = B1
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [0,0]
+            sect_cg_CC = [(B1+B2)*0.2,-H*0.5]
+            sect_cg_RB = [max(B1,B2),-H]
+
+            # sect_shape = [[0.5*tw,-0.5*tf1],[B1,-0.5*tf1],[0.5*tw,-H+0.5*tf2],[B2,-H+0.5*tf2]]
+            sect_shape = [[0,0],[B1,0],[0,-H],[B2,-H]]
+            sect_lin_con = [[2,1],[1,3],[3,4]]
+            sect_thk = [tf1,tw,tf2]
+            # sect_thk_off = [0,0,0]
+            sect_thk_off = [tf1/2,tw/2,tf2/2]
+
+        elif shape.SHAPE == 'H' :
+            H,B1,tw,tf1,B2,tf2,r1,r2 = shape.PARAMS[:8]
+            if B2 == 0 : B2 = B1
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [-0.5*max(B1,B2),0.5*H]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [0.5*max(B1,B2),-0.5*H]
+
+            sect_shape = [[-0.5*B1,0.5*(H-tf1)],[0,0.5*(H-tf1)],[0.5*B1,0.5*(H-tf1)],[-0.5*B2,-0.5*(H-tf2)],[0,-0.5*(H-tf2)],[0.5*B2,-0.5*(H-tf2)]]
+            sect_lin_con = [[2,1],[3,2],[2,5],[4,5],[5,6]]
+            sect_thk = [tf1,tf1,tw,tf2,tf2]
+            sect_thk_off = [0,0,0,0,0]
+        
+        elif shape.SHAPE == 'T' :
+            H,B,tw,tf = shape.PARAMS[:4]
+
+            sect_cg_LT = [-B*0.5,0]
+            sect_cg_CC = [0,-H*0.3]
+            sect_cg_RB = [B*0.5,-H]
+
+            sect_shape = [[-0.5*B,-0.5*tf],[0,-0.5*tf],[0.5*B,-0.5*tf],[0,-H]]
+            sect_lin_con = [[2,1],[3,2],[2,4]]
+            sect_thk = [tf,tf,tw]
+            sect_thk_off = [0,0,0]
+
+        elif shape.SHAPE == 'B' :
+            H,B,tw,tf1,C,tf2 = shape.PARAMS[:6]
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [-0.5*B,0.5*H]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [0.5*B,-0.5*H]
+
+            # sect_shape = [[0.5*(B-tw),0.5*(H-tf1)],[-0.5*(B-tw),0.5*(H-tf1)],[-0.5*(B-tw),-0.5*(H-tf2)],[0.5*(B-tw),-0.5*(H-tf2)]]
+            sect_shape = [[0.5*B,0.5*H],[-0.5*B,0.5*H],[-0.5*B,-0.5*H],[0.5*B,-0.5*H]]
+
+            sect_lin_con = [[1,2],[2,3],[3,4],[4,1]]
+            sect_thk = [tf1,tw,tf2,tw]
+            # sect_thk_off = [0,0,0,0]
+            sect_thk_off = [0.5*tf1,0.5*tw,0.5*tf2,0.5*tw]
+        
+        elif shape.SHAPE == 'P' :
+            D,tw = shape.PARAMS[:2]
+
+            # R = 0.5*(D-tw)
+            R = 0.5*D
+
+            sect_cg_LT = [-R,R]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [R,-R]
+
+            sect_shape = []
+            sect_lin_con = []
+            sect_thk = []
+            sect_thk_off = []
+
+            n = 16
+            for i in range(n):
+                sect_shape.append([R*math.sin(i*2*math.pi/n),R*math.cos(i*2*math.pi/n)])
+                sect_lin_con.append([i+1,i+2])
+                sect_thk.append(tw)
+                sect_thk_off.append(-0.5*tw)
+            sect_lin_con[-1] = [i+1,1]
+
+
+
+        sect_cg = (sect_cg_LT,sect_cg_CC,sect_cg_RB)
+
+        return sect_shape, sect_thk ,sect_thk_off, sect_cg , sect_lin_con

midas_civil-1.1.0/midas_civil/_section/_dbSecSS.py

@@ -0,0 +1,164 @@
+from ._offsetSS import Offset
+from ._offsetSS import _common
+import math
+
+class _SS_DBUSER(_common):
+
+    """ Create Standard USER DEFINED sections"""
+
+    def __init__(self,Name='',Shape='',parameters:list=[],Offset=Offset(),useShear=True,use7Dof=False,id:int=0):  
+        """ Shape = 'SB' 'SR' for rectangle \n For cylinder"""
+        self.ID = id
+        self.NAME = Name
+        self.TYPE = 'DBUSER'
+        self.SHAPE = Shape
+        self.PARAMS = parameters
+        self.OFFSET = Offset
+        self.USESHEAR = useShear
+        self.USE7DOF = use7Dof
+        self.DATATYPE = 2
+    
+    def __str__(self):
+         return f'  >  ID = {self.ID}   |  USER DEFINED STANDARD SECTION \nJSON = {self.toJSON()}\n'
+
+
+    def toJSON(sect):
+        js =  {
+                "SECTTYPE": "DBUSER",
+                "SECT_NAME": sect.NAME,
+                "SECT_BEFORE": {
+                    "SHAPE": sect.SHAPE,
+                    "DATATYPE": sect.DATATYPE,
+                    "SECT_I": {
+                        "vSIZE": sect.PARAMS
+                    }
+                }
+            }
+        js['SECT_BEFORE'].update(sect.OFFSET.JS)
+        js['SECT_BEFORE']['USE_SHEAR_DEFORM'] = sect.USESHEAR
+        js['SECT_BEFORE']['USE_WARPING_EFFECT'] = sect.USE7DOF
+        return js
+    
+    @staticmethod
+    def _objectify(id,name,type,shape,offset,uShear,u7DOF,js):
+        return _SS_DBUSER(name,shape,js['SECT_BEFORE']['SECT_I']['vSIZE'],offset,uShear,u7DOF,id)
+    
+    def _centerLine(shape,end):
+        if shape.SHAPE == 'SB' :
+            H,B = shape.PARAMS[:2]
+
+            sect_lin_con = [[1,2],[3,1]]
+
+            sect_cg_LT = [-B/2,H/2]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [B/2,-H/2]
+
+            if H > B :
+                sect_shape = [[0,0],[0,H/2],[0,-H/2],[0,H/4],[0,-H/4]]
+                sect_thk = [B,B,B,B]
+                sect_thk_off = [0,0,0,0]
+            else : 
+                sect_shape = [[0,0],[B/2,0],[-B/2,0],[B/4,0],[-B/4,0]]
+                sect_thk = [H,H,H,H]
+                sect_thk_off = [0,0,0,0]
+
+        elif shape.SHAPE == 'L' :
+            H,B,tw,tf = shape.PARAMS[:4]
+
+            sect_cg_LT = [0,0]
+            sect_cg_CC = [(H*tw*tw+B*B*tf)/(2*(B*tw+H*tf)),-(H*H*tw+B*tf*tf)/(2*(B*tw+H*tf))]
+            sect_cg_RB = [B,-H]
+
+            # sect_shape = [[0.5*tw,-H],[0.5*tw,-0.5*tf],[B,-0.5*tf]]
+            sect_shape = [[0,-H],[0,0],[B,0]]
+            sect_lin_con = [[3,2],[2,1]]
+            sect_thk = [tw,tf]
+            # sect_thk_off = [0,0]
+            sect_thk_off = [tw/2,tf/2]
+        
+        elif shape.SHAPE == 'C' :
+            H,B1,tw,tf1,B2,tf2 = shape.PARAMS[:6]
+            if B2 == 0 : B2 = B1
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [0,0]
+            sect_cg_CC = [(B1+B2)*0.2,-H*0.5]
+            sect_cg_RB = [max(B1,B2),-H]
+
+            # sect_shape = [[0.5*tw,-0.5*tf1],[B1,-0.5*tf1],[0.5*tw,-H+0.5*tf2],[B2,-H+0.5*tf2]]
+            sect_shape = [[0,0],[B1,0],[0,-H],[B2,-H]]
+            sect_lin_con = [[2,1],[1,3],[3,4]]
+            sect_thk = [tf1,tw,tf2]
+            # sect_thk_off = [0,0,0]
+            sect_thk_off = [tf1/2,tw/2,tf2/2]
+
+        elif shape.SHAPE == 'H' :
+            H,B1,tw,tf1,B2,tf2,r1,r2 = shape.PARAMS[:8]
+            if B2 == 0 : B2 = B1
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [-0.5*max(B1,B2),0.5*H]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [0.5*max(B1,B2),-0.5*H]
+
+            sect_shape = [[-0.5*B1,0.5*(H-tf1)],[0,0.5*(H-tf1)],[0.5*B1,0.5*(H-tf1)],[-0.5*B2,-0.5*(H-tf2)],[0,-0.5*(H-tf2)],[0.5*B2,-0.5*(H-tf2)]]
+            sect_lin_con = [[2,1],[3,2],[2,5],[4,5],[5,6]]
+            sect_thk = [tf1,tf1,tw,tf2,tf2]
+            sect_thk_off = [0,0,0,0,0]
+        
+        elif shape.SHAPE == 'T' :
+            H,B,tw,tf = shape.PARAMS[:4]
+
+            sect_cg_LT = [-B*0.5,0]
+            sect_cg_CC = [0,-H*0.3]
+            sect_cg_RB = [B*0.5,-H]
+
+            sect_shape = [[-0.5*B,-0.5*tf],[0,-0.5*tf],[0.5*B,-0.5*tf],[0,-H]]
+            sect_lin_con = [[2,1],[3,2],[2,4]]
+            sect_thk = [tf,tf,tw]
+            sect_thk_off = [0,0,0]
+
+        elif shape.SHAPE == 'B' :
+            H,B,tw,tf1,C,tf2 = shape.PARAMS[:6]
+            if tf2 == 0 : tf2 = tf1
+
+            sect_cg_LT = [-0.5*B,0.5*H]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [0.5*B,-0.5*H]
+
+            # sect_shape = [[0.5*(B-tw),0.5*(H-tf1)],[-0.5*(B-tw),0.5*(H-tf1)],[-0.5*(B-tw),-0.5*(H-tf2)],[0.5*(B-tw),-0.5*(H-tf2)]]
+            sect_shape = [[0.5*B,0.5*H],[-0.5*B,0.5*H],[-0.5*B,-0.5*H],[0.5*B,-0.5*H]]
+
+            sect_lin_con = [[1,2],[2,3],[3,4],[4,1]]
+            sect_thk = [tf1,tw,tf2,tw]
+            # sect_thk_off = [0,0,0,0]
+            sect_thk_off = [0.5*tf1,0.5*tw,0.5*tf2,0.5*tw]
+        
+        elif shape.SHAPE == 'P' :
+            D,tw = shape.PARAMS[:2]
+
+            # R = 0.5*(D-tw)
+            R = 0.5*D
+
+            sect_cg_LT = [-R,R]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [R,-R]
+
+            sect_shape = []
+            sect_lin_con = []
+            sect_thk = []
+            sect_thk_off = []
+
+            n = 16
+            for i in range(n):
+                sect_shape.append([R*math.sin(i*2*math.pi/n),R*math.cos(i*2*math.pi/n)])
+                sect_lin_con.append([i+1,i+2])
+                sect_thk.append(tw)
+                sect_thk_off.append(-0.5*tw)
+            sect_lin_con[-1] = [i+1,1]
+
+
+
+        sect_cg = (sect_cg_LT,sect_cg_CC,sect_cg_RB)
+
+        return sect_shape, sect_thk ,sect_thk_off, sect_cg , sect_lin_con

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_section/_pscSS.py

@@ -1,5 +1,6 @@
 from ._offsetSS import Offset
 from ._offsetSS import _common
+import math
 
 
 
@@ -114,6 +115,35 @@ class _SS_PSC_12CELL(_common):
                             vD[0],vD[1],vD[2],vD[3],vD[4],vD[5],vD[6],vD[7],
                             offset,uShear,u7DOF,id)
     
+    def _centerLine(shape,end):
+        if shape.SHAPE == '1CEL':
+            BO1 = shape.BO1
+            BO2 = shape.BO2
+            BO3 = shape.BO3
+            HO1 = shape.HO1
+            HO2 = shape.HO2
+            HO3 = shape.HO3
+            HI1 = shape.HI1
+            HI5 = shape.HI5
+            BI3 = shape.BI3
+
+            sect_lin_con = [[1,2],[2,3],[3,4],[2,5],[3,6],[5,6]]
+
+            sect_cg_LT = [0,0]
+            sect_cg_RB = [2*(BO1+BO2+BO3),-(HO1+HO2+HO3)]
+            sect_cg_CC = [(BO1+BO2+BO3),-(HO1+HO2+HO3)/2]
+
+            sect_shape = [[0,0],[BO1,0],[BO1+2*(BO2+BO3),0],[2*(BO1+BO2+BO3),0],[BO1+BO2,-HO1-HO3-HO2],[BO1+BO2+2*BO3,-HO1-HO3-HO2]]
+            sect_thk = [HO1,HI1,HO1,BO3-BI3,BO3-BI3,HI5]
+            sect_thk_off = [-HO1/2,-HI1/2,-HO1/2,(BO3-BI3)/2,-(BO3-BI3)/2,HI5/2]       
+
+
+
+        sect_cg = (sect_cg_LT,sect_cg_CC,sect_cg_RB)
+
+        return sect_shape, sect_thk ,sect_thk_off, sect_cg , sect_lin_con
+
+    
 
 
 

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_tendon.py

@@ -761,7 +761,7 @@ class Tendon:
                             array_temp.append({
                                     'PT' : [self.P_XYZ[j].X,self.P_XYZ[j].Y,self.P_XYZ[j].Z],
                                     'bFIX' : self.bFIX[j],
-                                    'R' : self.RADIUS[j]
+                                    'RADIUS' : self.RADIUS[j]
                             })
                     else:
                         for j in range(len(self.P_XYZ)):
@@ -772,6 +772,8 @@ class Tendon:
                             })
                     
                     
+                    
+                    
                     # --- 3D Main ----
 
                     json_prof = {

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_thickness.py

@@ -127,6 +127,11 @@ class Thickness(_common):
         Thickness.thick=[]
         Thickness.ids=[]
 
+    @staticmethod
+    def clear():
+        Thickness.thick=[]
+        Thickness.ids=[]
+
 
     @staticmethod
     def sync():

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil/_utils.py

@@ -86,4 +86,12 @@ def zz_add_to_dict(dictionary, key, value):
 def _convItem2List(item):
     if isinstance(item,list):
         return item
-    return [item]
+    return [item]
+
+def _matchArray(A,B):
+    A = _convItem2List(A)
+    B = _convItem2List(B)
+    n = len(A)
+    if len(B) >= n:
+        return B[:n]
+    return B + [B[-1]] * (n - len(B))

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: midas_civil
-Version: 1.0.8
+Version: 1.1.0
 Summary: Python library for MIDAS Civil NX
 Home-page: https://github.com/MIDASIT-Co-Ltd/midas-civil-python
 Author: Sumit Shekhar
@@ -14,6 +14,7 @@ Requires-Dist: polars
 Requires-Dist: xlsxwriter
 Requires-Dist: requests
 Requires-Dist: scipy
+Requires-Dist: colorama
 Dynamic: author
 Dynamic: author-email
 Dynamic: description

{midas_civil-1.0.8 → midas_civil-1.1.0}/midas_civil.egg-info/requires.txt

@@ -3,3 +3,4 @@ polars
 xlsxwriter
 requests
 scipy
+colorama

{midas_civil-1.0.8 → midas_civil-1.1.0}/setup.py

@@ -5,7 +5,7 @@ with open('README.md','r') as f:
 
 
 setup(name='midas_civil',
-    version='1.0.8',
+    version='1.1.0',
     description='Python library for MIDAS Civil NX',
     author='Sumit Shekhar',
     author_email='sumit.midasit@gmail.com',
@@ -15,7 +15,8 @@ setup(name='midas_civil',
         'polars',
         'xlsxwriter',
         'requests',
-        'scipy'
+        'scipy',
+        'colorama'
     ],          
     long_description= description,
     long_description_content_type='text/markdown',

midas_civil-1.0.8/midas_civil/__init__.py

@@ -1,49 +0,0 @@
-import requests
-from ._mapi import *
-_version_ = "1.0.8"
-
-
-print('')
-print('*'*20,'  MIDAS CIVIL-NX PYTHON LIBRARY v',_version_,' 🐍 ','*'*20)
-print('')
-
-if NX.version_check:
-    resp =  requests.get("https://pypi.org/pypi/midas_civil/json").json()
-    latest_ver =  resp["info"]["version"]
-    if _version_ != latest_ver:
-        print(
-                f"⚠️  Warning: You are using v{_version_}, "
-                f"but the latest available version is v{latest_ver}.\n"
-                f" Run 'pip install midas_civil --upgrade' to update."
-            )
-        print("-"*85)
-
-
-from ._model import *
-from ._boundary import *
-from ._utils import *
-from ._node import *
-from ._element import *
-from ._load import *
-from ._group import *
-from ._result import *
-
-#--- TESTING IMPORTS ---
-from ._material import *
-
-# from ._section import *
-from ._section import *
-
-from ._result_extract import *
-from ._construction import *
-from ._thickness import *
-from ._temperature import *
-
-from ._tendon import *
-from ._view import *
-
-from ._movingload import *
-from ._settlement import *
-from ._analysiscontrol import *
-from ._BoundaryChangeAssignment import*
-

midas_civil-1.0.8/midas_civil/_section/_TapdbSecSS.py

@@ -1,47 +0,0 @@
-from ._offsetSS import Offset
-from ._offsetSS import _common
-
-class _SS_TAPERED_DBUSER(_common):
-
-    """ Create Standard USER DEFINED sections"""
-
-    def __init__(self,Name='',Shape='',params_I:list=[],params_J:list=[],Offset=Offset(),useShear=True,use7Dof=False,id:int=0):  
-        """ Shape = 'SB' 'SR' for rectangle \n For cylinder"""
-        self.ID = id
-        self.NAME = Name
-        self.TYPE = 'TAPERED'
-        self.SHAPE = Shape
-        self.PARAMS_I = params_I
-        self.PARAMS_J = params_J
-        self.OFFSET = Offset
-        self.USESHEAR = useShear
-        self.USE7DOF = use7Dof
-        self.DATATYPE = 2
-    
-    def __str__(self):
-         return f'  >  ID = {self.ID}   |  USER DEFINED STANDARD SECTION \nJSON = {self.toJSON()}\n'
-
-
-    def toJSON(sect):
-        js =  {
-                "SECTTYPE": sect.TYPE,
-                "SECT_NAME": sect.NAME,
-                "SECT_BEFORE": {
-                    "SHAPE": sect.SHAPE,
-                    "TYPE": sect.DATATYPE,
-                    "SECT_I": {
-                        "vSIZE": sect.PARAMS_I
-                    },
-                    "SECT_J": {
-                        "vSIZE": sect.PARAMS_J
-                    }
-                }
-            }
-        js['SECT_BEFORE'].update(sect.OFFSET.JS)
-        js['SECT_BEFORE']['USE_SHEAR_DEFORM'] = sect.USESHEAR
-        js['SECT_BEFORE']['USE_WARPING_EFFECT'] = sect.USE7DOF
-        return js
-    
-    @staticmethod
-    def _objectify(id,name,type,shape,offset,uShear,u7DOF,js):
-        return _SS_TAPERED_DBUSER(name,shape,js['SECT_BEFORE']['SECT_I']['vSIZE'],js['SECT_BEFORE']['SECT_J']['vSIZE'],offset,uShear,u7DOF,id)

midas_civil-1.0.8/midas_civil/_section/_dbSecSS.py

@@ -1,57 +0,0 @@
-from ._offsetSS import Offset
-from ._offsetSS import _common
-
-class _SS_DBUSER(_common):
-
-    """ Create Standard USER DEFINED sections"""
-
-    def __init__(self,Name='',Shape='',parameters:list=[],Offset=Offset(),useShear=True,use7Dof=False,id:int=0):  
-        """ Shape = 'SB' 'SR' for rectangle \n For cylinder"""
-        self.ID = id
-        self.NAME = Name
-        self.TYPE = 'DBUSER'
-        self.SHAPE = Shape
-        self.PARAMS = parameters
-        self.OFFSET = Offset
-        self.USESHEAR = useShear
-        self.USE7DOF = use7Dof
-        self.DATATYPE = 2
-    
-    def __str__(self):
-         return f'  >  ID = {self.ID}   |  USER DEFINED STANDARD SECTION \nJSON = {self.toJSON()}\n'
-
-
-    def toJSON(sect):
-        js =  {
-                "SECTTYPE": "DBUSER",
-                "SECT_NAME": sect.NAME,
-                "SECT_BEFORE": {
-                    "SHAPE": sect.SHAPE,
-                    "DATATYPE": sect.DATATYPE,
-                    "SECT_I": {
-                        "vSIZE": sect.PARAMS
-                    }
-                }
-            }
-        js['SECT_BEFORE'].update(sect.OFFSET.JS)
-        js['SECT_BEFORE']['USE_SHEAR_DEFORM'] = sect.USESHEAR
-        js['SECT_BEFORE']['USE_WARPING_EFFECT'] = sect.USE7DOF
-        return js
-    
-    @staticmethod
-    def _objectify(id,name,type,shape,offset,uShear,u7DOF,js):
-        return _SS_DBUSER(name,shape,js['SECT_BEFORE']['SECT_I']['vSIZE'],offset,uShear,u7DOF,id)
-    
-    def _centerLine(shape):
-        H = shape.PARAMS[0]
-        B = shape.PARAMS[1]
-        tw = shape.PARAMS[2]
-        tf = shape.PARAMS[3]
-
-        sect_lin_con = [[1,2],[3,1]]
-
-        sect_cg = [0,0]
-
-        sect_shape = [[0,0],[B,0],[0,-H]]
-        sect_thk = [tf,tw]
-        sect_thk_off = [-tf/2,-tw/2]