midas-civil 1.0.8__py3-none-any.whl → 1.1.0__py3-none-any.whl

midas_civil/__init__.py

@@ -1,22 +1,21 @@
 import requests
+from colorama import Fore,Style
 from ._mapi import *
-_version_ = "1.0.8"
+_version_ = "1.1.0"
 
 
-print('')
-print('*'*20,'  MIDAS CIVIL-NX PYTHON LIBRARY v',_version_,' 🐍 ','*'*20)
-print('')
+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(
-                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)
+    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 *

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/_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/_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/_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/_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/_section/_TapdbSecSS.py

@@ -1,5 +1,6 @@
 from ._offsetSS import Offset
 from ._offsetSS import _common
+import math
 
 class _SS_TAPERED_DBUSER(_common):
 
@@ -44,4 +45,131 @@ class _SS_TAPERED_DBUSER(_common):
     
     @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)
+        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/_section/_dbSecSS.py

@@ -1,5 +1,6 @@
 from ._offsetSS import Offset
 from ._offsetSS import _common
+import math
 
 class _SS_DBUSER(_common):
 
@@ -42,16 +43,122 @@ class _SS_DBUSER(_common):
     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]
+    def _centerLine(shape,end):
+        if shape.SHAPE == 'SB' :
+            H,B = shape.PARAMS[:2]
 
-        sect_lin_con = [[1,2],[3,1]]
+            sect_lin_con = [[1,2],[3,1]]
 
-        sect_cg = [0,0]
+            sect_cg_LT = [-B/2,H/2]
+            sect_cg_CC = [0,0]
+            sect_cg_RB = [B/2,-H/2]
 
-        sect_shape = [[0,0],[B,0],[0,-H]]
-        sect_thk = [tf,tw]
-        sect_thk_off = [-tf/2,-tw/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/_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/_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/_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/_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.dist-info → midas_civil-1.1.0.dist-info}/METADATA

@@ -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.dist-info → midas_civil-1.1.0.dist-info}/RECORD

@@ -1,36 +1,36 @@
 midas_civil/_BoundaryChangeAssignment.py,sha256=8lQC0DZ8zAXuPfchBdhqUaLZhZMog51EffQcKeayrlc,9511
-midas_civil/__init__.py,sha256=T0NwsEmr2c1166UlxbogGhDfBzvRGbTIU9XX9c6h-_E,1156
+midas_civil/__init__.py,sha256=BTE7nEPN7fTfwk686EuvD7mBxvyxf0MrEITZt7zDqu0,2328
 midas_civil/_analysiscontrol.py,sha256=r2zp7fdC_m-WfO7KSXU4SzVjOgLWg4Wm0djl7X1ioZo,24456
 midas_civil/_boundary.py,sha256=s5mgZIc1b0-P7AdWuaPAQ5cIbHL5CR1YRKJA7KE4W_U,32958
 midas_civil/_construction.py,sha256=q9C3gGdvky5t7jCDzd6lBRv6pm8cMycI91xfkgDV37Y,40673
 midas_civil/_construction_backup.py,sha256=Pj7V-NYCkkT-aMjKXfs1jKa9klsGh48UXDLwn3BLYTY,18225
-midas_civil/_element.py,sha256=1GcI-9PcRdfvHvtXs2tdMNRvvxyQ7A78Lk8rAp5gGOw,34795
+midas_civil/_element.py,sha256=o6gJcOAKsoUA77uNDFNGUrmMV6DSz1SKP5g2MbtMui4,36212
 midas_civil/_group.py,sha256=ul0gPtAhfuc9HPQkVw2dfV0jqCLvnTbxAlx6hEBUmyI,9882
 midas_civil/_load.py,sha256=twsfqofTzbnMHXShfMRSP61ll7CTBA5qos33od7nQ9E,30170
-midas_civil/_mapi.py,sha256=UMGuzuNE7uTjR6YwPsbfI5xHSoFESVY2U0NM5mfgcW8,5486
-midas_civil/_material.py,sha256=uJEIHJM9OhwTRWUI2mtd_0BQSxdlYhATYJu9P7tNNBA,69511
+midas_civil/_mapi.py,sha256=b-leROZA9Ii-sSawBrBPgGS6ow5JXzEYGh5DvkXAOI8,7749
+midas_civil/_material.py,sha256=x96yHu96X3wE9KF7dW2vzbFg9_58CMQ6Gk8maBfLzNM,69573
 midas_civil/_model.py,sha256=huUlmIZVsBbThszrEIhvb9owOSmdggylXaJeK4g7lxY,16735
 midas_civil/_movingload.py,sha256=a1i5fNc7oCmWTjVXbyBEYzDr5BO7NTXPVCKgvjzxuX4,80000
-midas_civil/_node.py,sha256=3iVj5Hs7DJPbfgM_hIazP-Y1szu3j3ofuh_ftUxoLcY,8827
+midas_civil/_node.py,sha256=qcBar89tEvHfziEvEvreJT8DjJ152CryIETHNKdqr2o,8956
 midas_civil/_result copy.py,sha256=siTMENLIwF_6rvydSjP9aQAWaIlt0pReiqNyDhDevGk,24290
 midas_civil/_result.py,sha256=nzXY2lh_342M45TMq2hBYwSp41_wAHyGFSfvt9h2vXc,7638
-midas_civil/_result_extract.py,sha256=QVo1DJmlkTJ20-FwqgguH70bVrwHbx7TfPFENYfCUbE,5706
+midas_civil/_result_extract.py,sha256=IJLnlzIKOP8yxaLCXqOuRKiYUYMQfzj37fTCBs0iuTM,7884
 midas_civil/_settlement.py,sha256=0ZkKutTjQ7Wg47ohfrc6SgoZY1xyP9_ylGQlCVBH-2U,5124
 midas_civil/_temperature.py,sha256=Sr5fRgUHzMx6tfr1E0Nw7ZLYKUmzoH1h6Sl-APqFkHo,24208
-midas_civil/_tendon.py,sha256=pUtR-74c1GcJyJt6Te38TAlBBOZJtlkERb9Ijf5CkH0,33417
-midas_civil/_thickness.py,sha256=yyFkMCaxcwRsqYJcRAxWqnt59fdqzX4Mqkjz23waw-A,3247
-midas_civil/_utils.py,sha256=c_NgCmqtoWP7jv4RoK7oGyF9edDo_UlJYR-p54h572U,2711
+midas_civil/_tendon.py,sha256=oNeT0e5VIWIDhhUfxdf2KQKXdiCkbN7fWjVYjPzs7mg,33464
+midas_civil/_thickness.py,sha256=rKbU0QKkL9AeFnE-MNVY9ZwwQ-u2Ea-Rq9TeNafsvWQ,3335
+midas_civil/_utils.py,sha256=nTF6btow2TK4mufyKupbjwagjlYp5WVRJmfchYAUdMk,2881
 midas_civil/_view.py,sha256=o7rkfoQzmgAb3dT0ujPIQLVwVlveo3rMRIbZU1UosZo,849
-midas_civil/_section/_TapdbSecSS.py,sha256=ANHBwzegppGmUlm7spEBNDHfdJfTbWX7Q-pvn-Th1tA,1649
+midas_civil/_section/_TapdbSecSS.py,sha256=biwt38Twlmb6DuueZG_KbuMCIv1Sq3TAl4ay-30LnRU,5933
 midas_civil/_section/__init__.py,sha256=subwWGtIRFfz19i2mtLPuQbViFAR8-xZY7eHLBGUWO4,15250
 midas_civil/_section/_compositeSS.py,sha256=mrHOkjHkbuVOLAd1dGEiLruyqmye_mv9H99UJIdWdiQ,9942
-midas_civil/_section/_dbSecSS.py,sha256=Dupfs4un0mwJbreqrUPw7iEstbdFiMiQYXVQtam-5k0,1763
+midas_civil/_section/_dbSecSS.py,sha256=GmVm8gxsLxnDsmJsl73koHYn1GRHVo_KDOoijz_j4_o,5541
 midas_civil/_section/_offsetSS.py,sha256=kHjwp_PCeBKWqYKoQCy0s4d8O5Gb0Vk6e1MWMOGxql4,1377
-midas_civil/_section/_pscSS.py,sha256=PnXR5lJNKYtp8MyFddJlW45CYCRSG6SpA-zojDl4DnE,14254
+midas_civil/_section/_pscSS.py,sha256=aTgFxN0WzT6fBdtPgUicBByFq6KyLMmzoyBD-eVUop8,15214
 midas_civil/_section/_tapPSC1CellSS.py,sha256=bjpuOgIot2yQd2zirE6vVmVlgfco1f83RlI7DwIicdU,6700
 midas_civil/_section/_unSupp.py,sha256=46A_a5CzDpQM7bqglkFzIvq6VG7Lj9J6qflY3nwZlqI,1674
-midas_civil-1.0.8.dist-info/licenses/LICENSE,sha256=zrL4RwZC4rb-by_ZHKXwKdIwcs6ATy59TPZ9HxPHCrs,1071
-midas_civil-1.0.8.dist-info/METADATA,sha256=BW2J0aNKOVhFN8yglhcnfnMj8smztMrx4xbSwEzN6Po,2085
-midas_civil-1.0.8.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-midas_civil-1.0.8.dist-info/top_level.txt,sha256=_NFmrlN5V9OxJ-PAO4s_om8OA8uupXho3QqZcSsnbuI,12
-midas_civil-1.0.8.dist-info/RECORD,,
+midas_civil-1.1.0.dist-info/licenses/LICENSE,sha256=zrL4RwZC4rb-by_ZHKXwKdIwcs6ATy59TPZ9HxPHCrs,1071
+midas_civil-1.1.0.dist-info/METADATA,sha256=YX6pmfkhA9Og-0Kszbc6S6qnP9H8kIJBWY7kg5XzZck,2109
+midas_civil-1.1.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+midas_civil-1.1.0.dist-info/top_level.txt,sha256=_NFmrlN5V9OxJ-PAO4s_om8OA8uupXho3QqZcSsnbuI,12
+midas_civil-1.1.0.dist-info/RECORD,,