@@ -20,6 +20,11 @@ class Motion(Enum):
2020 DISPLACEMENT = 2
2121
2222
23+ class RWMotion (Enum ):
24+ VELOCITY = 0
25+ DISPLACEMENT = 1
26+
27+
2328class DOF (Enum ):
2429 X_TRANSLATIONAL = 1
2530 Y_TRANSLATIONAL = 2
@@ -140,12 +145,27 @@ class Function:
140145
141146 def __init__ (self , Function = None ):
142147 self .function = Function
148+ self .tabulated = False
149+
150+ def set_tabulated (self , heading = "" , function = "" , x = [], y = []):
151+ self .tabulated = True
152+ self .heading = heading
153+ self .function_name = function
154+ self .x = x
155+ self .y = y
143156
144157 def create (self , stub ):
145158 """Create function."""
159+ if self .tabulated :
160+ ret = stub .CreateDefineFunctionTabulated (
161+ DefineFunctionTabulatedRequest (
162+ heading = self .heading , function = self .function_name , abscissa = self .x , ordinate = self .y
163+ )
164+ )
146165 ret = stub .CreateDefineFunction (DefineFunctionRequest (function = self .function ))
147166 self .id = ret .id
148167 logging .info (f"Function { self .id } )
168+
149169 return self .id
150170
151171
@@ -763,7 +783,7 @@ def create_general_keyword(self, opcode, keyworddata):
763783 \$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
764784
765785 opcode = "INITIAL_VELOCITY"
766- keyworddata = "0\\ n1.480E+01,0.000E+00,0.000E+00,0.000E+00,0.000E+00,0.000E+00"
786+ keyworddata = "0\n 1.480E+01,0.000E+00,0.000E+00,0.000E+00,0.000E+00,0.000E+00"
767787 create_general_keyword(opcode = opcode,keyworddata=keyworddata)
768788 """
769789 ret = self .stub .CreateGeneralKWD (GeneralKWDRequest (opcode = opcode , keyworddata = keyworddata ))
@@ -2224,6 +2244,7 @@ def __init__(self):
22242244 self .spclist = []
22252245 self .imposedmotionlist = []
22262246 self .templist = []
2247+ self .convectionlist = []
22272248
22282249 def create_spc (
22292250 self ,
@@ -2305,6 +2326,39 @@ def create_temperature(
23052326 param = [nodeset , curve , scalefactor ]
23062327 self .templist .append (param )
23072328
2329+ def create_convection (
2330+ self ,
2331+ segmentset = None ,
2332+ convection_heat_transfer_coefficient = None ,
2333+ convection_heat_transfer_coefficient_multiplier = 0.0 ,
2334+ environment_temperature = None ,
2335+ environment_temperature_multiplier = 0.0 ,
2336+ ):
2337+ """Apply a convection boundary condition on SEGMENT_SET for a thermal analysis.
2338+
2339+ Parameters
2340+ ----------
2341+ segmentset : SegmentSet.
2342+ Segment set.
2343+ convection_heat_transfer_coefficient : Curve
2344+ Convection heat transfer coefficient.
2345+ convection_heat_transfer_coefficient_multiplier : float
2346+ Curve multiplier for convection heat transfer coefficient.
2347+ environment_temperature : Curve
2348+ Environment temperature.
2349+ environment_temperature_multiplier : float
2350+ Curve multiplier for environment temperature.
2351+
2352+ """
2353+ param = [
2354+ segmentset ,
2355+ convection_heat_transfer_coefficient ,
2356+ convection_heat_transfer_coefficient_multiplier ,
2357+ environment_temperature ,
2358+ environment_temperature_multiplier ,
2359+ ]
2360+ self .convectionlist .append (param )
2361+
23082362 def create (self ):
23092363 """Create a boundary condition."""
23102364 for obj in self .imposedmotionlist :
@@ -2399,7 +2453,7 @@ def create(self):
23992453 nid = nodeset .create (self .stub )
24002454 option = "SET"
24012455 if curve is not None :
2402- cid = curve .id
2456+ cid = curve .create ( self . stub )
24032457 else :
24042458 cid = 0
24052459 self .stub .CreateBdyTemp (
@@ -2411,6 +2465,19 @@ def create(self):
24112465 )
24122466 )
24132467 logging .info ("Boundary Temperature Created..." )
2468+ for obj in self .convectionlist :
2469+ ss , hlc , hmult , tlc , tmult = obj [0 ], obj [1 ], obj [2 ], obj [3 ], obj [4 ]
2470+ ssid , hlcid , tlcid = 0 , 0 , 0
2471+ if ss is not None :
2472+ ssid = ss .create (self .stub )
2473+ if hlc is not None :
2474+ hlcid = hlc .create (self .stub )
2475+ if tlc is not None :
2476+ tlcid = tlc .create (self .stub )
2477+ self .stub .CreateBdyConvection (
2478+ BdyConvectionRequest (ssid = ssid , pserod = 0 , hlcid = hlcid , hmult = hmult , tlcid = tlcid , tmult = tmult , loc = 0 )
2479+ )
2480+ logging .info ("Boundary Convection Created..." )
24142481
24152482
24162483class InitialCondition :
@@ -2528,13 +2595,11 @@ def __init__(self, tail=Point(0, 0, 0), head=Point(0, 0, 0), radius=1, length=10
25282595 self .lcid = 0
25292596 self .dir = Direction (1 , 0 , 0 )
25302597
2531- def set_motion (self , curve , motion = Motion .VELOCITY , dir = Direction (1 , 0 , 0 )):
2598+ def set_motion (self , curve , motion = RWMotion .VELOCITY , dir = Direction (1 , 0 , 0 )):
25322599 """Set the prescribed motion."""
25332600 curve .create (self .stub )
25342601 self .lcid = curve .id
25352602 self .motion = motion .value
2536- if self .motion == Motion .DISPLACEMENT :
2537- self .motion = 1
25382603 self .dir = dir
25392604
25402605 def create (self ):
@@ -2589,13 +2654,11 @@ def __init__(self, center=Point(0, 0, 0), orient=Point(0, 0, 0), radius=1):
25892654 self .lcid = 0
25902655 self .dir = Direction (1 , 0 , 0 )
25912656
2592- def set_motion (self , curve , motion = Motion .VELOCITY , dir = Direction (1 , 0 , 0 )):
2657+ def set_motion (self , curve , motion = RWMotion .VELOCITY , dir = Direction (1 , 0 , 0 )):
25932658 """Set the prescribed motion."""
25942659 curve .create (self .stub )
25952660 self .lcid = curve .id
25962661 self .motion = motion .value
2597- if self .motion == Motion .DISPLACEMENT :
2598- self .motion = 1
25992662 self .dir = dir
26002663
26012664 def create (self ):
0 commit comments