support randles_batmac and some EM materials

Author: zhangzhanqun

src/ansys/dyna/core/pre/dynabase.py

@@ -113,21 +113,22 @@ class Curve:
 
     """
 
-    def __init__(self, sfo=1, x=[], y=[], func=None):
+    def __init__(self, sfo=1, x=[], y=[], func=None, title=""):
         self.sfo = sfo
         self.abscissa = x
         self.ordinate = y
         self.func = func
+        self.title = title
 
     def create(self, stub=None):
         """Create a curve."""
         if stub is None:
             stub = DynaBase.get_stub()
         if self.func != None:
-            ret = stub.CreateDefineCurveFunction(DefineCurveFunctionRequest(function=self.func))
+            ret = stub.CreateDefineCurveFunction(DefineCurveFunctionRequest(function=self.func, title=self.title))
         else:
             ret = stub.CreateDefineCurve(
-                DefineCurveRequest(sfo=self.sfo, abscissa=self.abscissa, ordinate=self.ordinate)
+                DefineCurveRequest(sfo=self.sfo, abscissa=self.abscissa, ordinate=self.ordinate, title=self.title)
             )
         self.id = ret.id
         logging.info(f"Curve {self.id} defined...")
@@ -1529,6 +1530,7 @@ def __init__(self, analysis_type=AnalysisType.IMPLICIT, initial_timestep_size=0)
         self.defined_dynamic = False
         self.defined_eigenvalue = False
         self.defined_solution = False
+        self.defined_mass_matrix = False
         self.imflag = analysis_type.value
         self.dt0 = initial_timestep_size
         self.stub = DynaBase.get_stub()
@@ -1640,6 +1642,10 @@ def set_solution(
         self.maxref = stiffness_reformation_limit
         self.abstol = absolute_convergence_tolerance
 
+    def set_consistent_mass_matrix(self):
+        """Use the consistent mass matrix in implicit dynamics and eigenvalue solutions."""
+        self.defined_mass_matrix = True
+
     def create(self):
         """Create an implicit analysis."""
         if self.defined == False:
@@ -1662,6 +1668,8 @@ def create(self):
                     nsolver=self.nsolver, ilimit=self.ilimit, maxref=self.maxref, abstol=self.abstol
                 )
             )
+        if self.defined_mass_matrix:
+            self.stub.CreateControlImplicitConsistentMass(ControlImplicitConsistentMassRequest(iflag=1))
 
 
 class ThermalAnalysisType(Enum):

src/ansys/dyna/core/pre/dynaem.py

@@ -732,3 +732,68 @@ def create(self):
         self.id = ret.id
         logging.info(f"EM Isopotential Rogo {self.id} Created...")
         return self.id
+
+
+class RandlesCellType(Enum):
+    USER_DEFINED_EQUIVALENT_CIRCUIT_MODEL = -1
+    RANDLES_CELL_0_ORDER = 0
+    RANDLES_CELL_1_ORDER = 1
+    RANDLES_CELL_2_ORDER = 2
+    RANDLES_CELL_3_ORDER = 3
+
+
+class RandlesCell:
+    """Define parameters for a Randles Cell."""
+
+    def __init__(self, set=None):
+        self.stub = DynaBase.get_stub()
+        self.define_batmac = False
+
+    def set_batmac_model(
+        self,
+        cell_type=RandlesCellType.RANDLES_CELL_1_ORDER,
+        cell_parts=None,
+        area=2,
+        cell_capacity=0,
+        soc_conversion_factor=0,
+        charge_init_state=0,
+        circuit_parameter=None,
+        constant_temperature=0,
+        temperature_from_thermal_solver=False,
+        add_heating_to_thermal_solver=False,
+    ):
+        """define the distributed Randles circuit parameters for a Randles cell when using the batmac model."""
+        self.define_batmac = True
+        self.rdltype = cell_type.value
+        self.rdlarea = area
+        self.psid = cell_parts
+        self.q = cell_capacity
+        self.cq = soc_conversion_factor
+        self.socinit = charge_init_state
+        self.prm = circuit_parameter
+        self.temp = constant_temperature
+        self.frther = temperature_from_thermal_solver
+        self.r0toth = add_heating_to_thermal_solver
+
+    def create(self):
+        """Set parameter for Randles Cell."""
+        if self.define_batmac:
+            sid = 0
+            if self.psid is not None:
+                sid = self.psid.create(self.stub)
+            ret = self.stub.CreateEMRandlesBatmac(
+                EMRandlesBatmacRequest(
+                    rdltype=self.rdltype,
+                    rdlarea=self.rdlarea,
+                    psid=sid,
+                    q=self.q,
+                    cq=self.cq,
+                    socinit=self.socinit,
+                    chargedirparam=self.prm,
+                    temp=self.temp,
+                    frther=self.frther,
+                    r0toth=self.r0toth,
+                )
+            )
+            self.id = ret.rdlid
+            logging.info(f"EM Randles Batmac {self.id} Created...")

src/ansys/dyna/core/pre/dynamaterial.py

@@ -383,6 +383,24 @@ def set_em_resistive_heating_2d(self, material_type=EMMATTYPE.CONDUCTOR, initial
         self.em_material_type = material_type.value
         self.em_initial_conductivity = initial_conductivity
 
+    def set_em_randles_batmac(
+        self, positive_current_collector_conductivity=0, negative_current_collector_conductivity=0
+    ):
+        """Define two conductivities per EM node for special applications (Randles Batmac).
+
+        Parameters
+        ----------
+        positive_current_collector_conductivity : float
+            conductivities of the positive current collector materials
+        negative_current_collector_conductivity : float
+            conductivities of the negative current collector materials
+
+        """
+        self.em_mat_type = 6
+        self.em_material_type = 5
+        self.sigp = positive_current_collector_conductivity
+        self.sign = negative_current_collector_conductivity
+
     def set_thermal_isotropic(
         self, density=0, generation_rate=0, generation_rate_multiplier=0, specific_heat=0, conductivity=0
     ):
@@ -436,6 +454,10 @@ def create(self, stub, matid):
                 stub.CreateEMMat004(
                     EMMat004Request(mid=matid, mtype=self.em_material_type, sigma=self.em_initial_conductivity)
                 )
+            elif self.em_mat_type == 6:
+                stub.CreateEMMat006(
+                    EMMat006Request(mid=matid, mtype=self.em_material_type, sigp=self.sigp, sign=self.sign)
+                )
         if self.em:
             if self.em_eos is not None:
                 eosid = self.em_eos.create(stub)
@@ -477,6 +499,29 @@ def create(self, stub):
         logging.info(f"Material {self.name} Created...")
 
 
+class MatPlasticKinematic:
+    """Define material of modelling isotropic and kinematic hardening plasticity."""
+
+    def __init__(
+        self, mass_density=0, young_modulus=0, poisson_ratio=0.3, yield_stress=0, tangent_modulus=0, hardening=0
+    ):
+        self.ro = mass_density
+        self.e = young_modulus
+        self.pr = poisson_ratio
+        self.sigy = yield_stress
+        self.etan = tangent_modulus
+        self.beta = hardening
+
+    def create(self, stub):
+        """Create plastic kinematic material."""
+        ret = stub.CreateMatPlasticKinematic(
+            MatPlasticKinematicRequest(ro=self.ro, e=self.e, pr=self.pr, sigy=self.sigy, etan=self.etan, beta=self.beta)
+        )
+        self.material_id = ret.mid
+        self.name = "Plastic Kinematic"
+        logging.info(f"Material {self.name} Created...")
+
+
 class MatElasticPlasticThermal(MatAdditional):
     """Defines temperature-dependent material coefficients."""
 
@@ -555,6 +600,32 @@ def create(self, stub):
         logging.info(f"Material {self.name} Created...")
 
 
+class MatCrushableFoam:
+    """Define material of modelling crushable foam."""
+
+    def __init__(
+        self, mass_density=0, young_modulus=0, poisson_ratio=0.3, yield_stress_curve=None, tensile_stress_cutoff=0
+    ):
+        self.ro = mass_density
+        self.e = young_modulus
+        self.pr = poisson_ratio
+        self.lcid = yield_stress_curve
+        self.tsl = tensile_stress_cutoff
+
+    def create(self, stub):
+        """Create crushable foam material."""
+        if self.lcid is not None:
+            cid = self.lcid.create(stub)
+        else:
+            cid = 0
+        ret = stub.CreateMatCrushableFoam(
+            MatCrushableFoamRequest(ro=self.ro, e=self.e, pr=self.pr, lcid=cid, tsl=self.tsl)
+        )
+        self.material_id = ret.mid
+        self.name = "Crushable Foam"
+        logging.info(f"Material {self.name} Created...")
+
+
 class MatThermalIsotropic:
     """Defines isotropic thermal properties."""
 
@@ -581,6 +652,30 @@ def create(self, stub):
         logging.info(f"Material {self.name} Created...")
 
 
+class MatThermalOrthotropic:
+    """Defines orthotropic thermal properties."""
+
+    def __init__(self, specific_heat=0, conductivity_x=0, conductivity_y=0, conductivity_z=0):
+        self.hc = specific_heat
+        self.k1 = conductivity_x
+        self.k2 = conductivity_y
+        self.k3 = conductivity_z
+
+    def create(self, stub):
+        """Create orthotropic thermal material."""
+        ret = stub.CreateMatThermalOrthotropic(
+            MatThermalOrthotropicRequest(
+                hc=self.hc,
+                k1=self.k1,
+                k2=self.k2,
+                k3=self.k3,
+            )
+        )
+        self.material_id = ret.mid
+        self.name = "Orthotropic thermal"
+        logging.info(f"Material {self.name} Created...")
+
+
 class MatRigidDiscrete:
     """Defines a rigid material for shells or solids.