add analysistype(explicit or implicit) for dynamech

Author: zhangzhanqun

examples/Implicit/camry_rc.py

@@ -25,7 +25,7 @@
     ContactSurface,
     DOF,
     OffsetType,
-    ImplicitAnalysis,
+    AnalysisType
 )
 from ansys.dyna.core.pre.dynamaterial import (
     MatNull,
@@ -80,7 +80,7 @@
 ###############################################################################
 # We can then use the implicit analysis methods in DynaMech class to define
 # the IMILICIT control cards.
-camry = DynaMech()
+camry = DynaMech(analysis=AnalysisType.EXPLICIT)
 camry_solution.add(camry)
 
 ###############################################################################

examples/Thermal/thermal_stress.py

@@ -18,7 +18,8 @@
     ThermalAnalysisType,
     SolidPart,
     SolidFormulation,
-    NodeSet
+    NodeSet,
+    AnalysisType
 )
 from ansys.dyna.core.pre.dynamaterial import MatElasticPlasticThermal
 from ansys.dyna.core.pre import examples
@@ -53,7 +54,7 @@
 ###############################################################################
 # To invoke the transient thermal solver, the thermal analysis type in CONTROL_SOLUTION is
 # being set to 2 by ThermalAnalysisType.TRANSIENT.
-ts = DynaMech()
+ts = DynaMech(analysis=AnalysisType.EXPLICIT)
 solution.add(ts)
 
 tanalysis = ThermalAnalysis()

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

@@ -218,6 +218,12 @@ def __init__(self):
         self.constraints = Constraint()
         self.contacts = ContactGroup()
         self.entities = []
+        self.have_accuracy = False
+        self.have_energy = False
+        self.have_hourglass = False
+        self.have_bulk_viscosity = False
+        self.have_control_shell = False
+
 
     def get_stub():
         """Get the stub of this DynaBase object."""
@@ -280,16 +286,19 @@ def set_accuracy(
         bool
             "True" when successful, "False" when failed
         """
-        ret = self.stub.CreateControlAccuracy(
-            ControlAccuracyRequest(
-                osu=objective_stress_updates.value,
-                inn=invariant_node_number.value,
-                pidosu=partsetid_for_objective_stress_updates,
-                iacc=implicit_accuracy_flag.value,
-                exacc=explicit_accuracy_flag.value,
+        ret = True
+        if not self.have_accuracy:
+            ret = self.stub.CreateControlAccuracy(
+                ControlAccuracyRequest(
+                    osu=objective_stress_updates.value,
+                    inn=invariant_node_number.value,
+                    pidosu=partsetid_for_objective_stress_updates,
+                    iacc=implicit_accuracy_flag.value,
+                    exacc=explicit_accuracy_flag.value,
+                )
             )
-        )
-        logging.info("Control Accuracy Created...")
+            self.have_accuracy = True
+            logging.info("Control Accuracy Created...")
         return ret
 
     def set_energy(
@@ -324,16 +333,19 @@ def set_energy(
         bool
             "True" when successful, "False" when failed
         """
-        ret = self.stub.CreateControlEnergy(
-            ControlEnergyRequest(
-                hgen=hourglass_energy.value,
-                rwen=rigidwall_energy.value,
-                slnten=sliding_interface_energy.value,
-                rylen=rayleigh_energy.value,
-                irgen=initial_reference_geometry_energy.value,
+        ret = True
+        if not self.have_energy:
+            ret = self.stub.CreateControlEnergy(
+                ControlEnergyRequest(
+                    hgen=hourglass_energy.value,
+                    rwen=rigidwall_energy.value,
+                    slnten=sliding_interface_energy.value,
+                    rylen=rayleigh_energy.value,
+                    irgen=initial_reference_geometry_energy.value,
+                )
             )
-        )
-        logging.info("Control Energy Created...")
+            self.have_energy = True
+            logging.info("Control Energy Created...")
         return ret
 
     def set_output(
@@ -385,8 +397,11 @@ def set_hourglass(self, controltype=HourglassControl.STANDARD_VISCOSITY_FORM, co
         bool
             "True" when successful, "False" when failed
         """
-        ret = self.stub.CreateControlHourgalss(ControlHourglassRequest(ihq=controltype.value, qh=coefficient))
-        logging.info("Control Hourglass Created...")
+        ret = True
+        if not self.have_hourglass:
+            ret = self.stub.CreateControlHourgalss(ControlHourglassRequest(ihq=controltype.value, qh=coefficient))
+            self.have_hourglass = True
+            logging.info("Control Hourglass Created...")
         return ret
 
     def set_bulk_viscosity(
@@ -411,14 +426,17 @@ def set_bulk_viscosity(
         bool
             "True" when successful, "False" when failed
         """
-        ret = self.stub.CreateControlBulkViscosity(
-            ControlBulkViscosityRequest(
-                q1=quadratic_viscosity_coeff,
-                q2=linear_viscosity_coeff,
-                type=bulk_viscosity_type.value,
+        ret = True
+        if not self.have_bulk_viscosity:
+            ret = self.stub.CreateControlBulkViscosity(
+                ControlBulkViscosityRequest(
+                    q1=quadratic_viscosity_coeff,
+                    q2=linear_viscosity_coeff,
+                    type=bulk_viscosity_type.value,
+                )
             )
-        )
-        logging.info("Control Bulk Viscosity Created...")
+            self.have_bulk_viscosity = True
+            logging.info("Control Bulk Viscosity Created...")
         return ret
 
     def set_init_temperature(self, temp=0):
@@ -482,20 +500,23 @@ def create_control_shell(
         bool
             "True" when successful, "False" when failed
         """
-        ret = self.stub.CreateControlShell(
-            ControlShellRequest(
-                wrpang=wrpang,
-                esort=esort,
-                irnxx=irnxx,
-                istupd=istupd,
-                theory=theory,
-                bwc=bwc,
-                miter=miter,
-                proj=proj,
-                irquad=irquad,
+        ret = True
+        if not self.have_control_shell:
+            ret = self.stub.CreateControlShell(
+                ControlShellRequest(
+                    wrpang=wrpang,
+                    esort=esort,
+                    irnxx=irnxx,
+                    istupd=istupd,
+                    theory=theory,
+                    bwc=bwc,
+                    miter=miter,
+                    proj=proj,
+                    irquad=irquad,
+                )
             )
-        )
-        logging.info("Control Shell Created...")
+            self.have_control_shell = True
+            logging.info("Control Shell Created...")
         return ret
 
     def create_control_solid(

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

@@ -7,13 +7,18 @@
 
 from .dynabase import *  # noqa : F403
 
+class AnalysisType(Enum):
+    EXPLICIT = 1
+    IMPLICIT = 2
+    NONE = 3
 
 class DynaMech(DynaBase):
     """Define an Mechanical analysis."""
 
-    def __init__(self):
+    def __init__(self,analysis=AnalysisType.EXPLICIT):
         DynaBase.__init__(self)
         self.casetype = CaseType.STRUCTURE
+        self.analysis = analysis.value
 
     def create_control_output(self, npopt=0, neecho=0):
         """Set miscellaneous output parameters.
@@ -239,40 +244,70 @@ def create_hourglass(self, ghid, ihq, qm=0.1, q1=1.5, q2=0.06, qb=1e-9, qw=1e-9)
         logging.info("Hourglass 1 Created...")
         return ret
 
-    def save_file(self):
+    def save_file(self,defaultsetting=1):
         """Save keyword files.
 
         Returns
         -------
         bool
             "True" when successful, "False" when failed
         """
-        self.set_accuracy(
-            objective_stress_updates=Switch.ON,
-            invariant_node_number=InvariantNode.ON_FOR_SHELL_TSHELL_SOLID,
-            implicit_accuracy_flag=Switch.ON,
-        )
-        self.set_bulk_viscosity(bulk_viscosity_type=BulkViscosity.COMPUTE_INTERNAL_ENERGY_DISSIPATED)
-        self.set_energy(
-            hourglass_energy=EnergyFlag.COMPUTED,
-            sliding_interface_energy=EnergyFlag.COMPUTED,
-        )
-        self.set_hourglass(
-            controltype=HourglassControl.FLANAGAN_BELYTSCHKO_INTEGRATION_SOLID,
-            coefficient=0,
-        )
-        if self.parts.get_num_shellpart() > 0:
-            self.create_control_shell(
-                wrpang=0,
-                esort=1,
-                irnxx=0,
-                istupd=4,
-                theory=0,
-                bwc=1,
-                miter=1,
-                proj=1,
-                irquad=0,
+        if self.analysis == 1:
+            self.set_accuracy(
+                objective_stress_updates=Switch.OFF,
+                invariant_node_number=InvariantNode.ON_FOR_SHELL_TSHELL_SOLID,
+                implicit_accuracy_flag=Switch.OFF,
+            )
+            self.set_bulk_viscosity(bulk_viscosity_type=BulkViscosity.COMPUTE_INTERNAL_ENERGY_DISSIPATED)
+            self.set_energy(
+                hourglass_energy=EnergyFlag.COMPUTED,
+                sliding_interface_energy=EnergyFlag.COMPUTED,
+            )
+            self.set_hourglass(
+                controltype=HourglassControl.FLANAGAN_BELYTSCHKO_INTEGRATION_SOLID,
+                coefficient=0,
             )
+            if self.parts.get_num_shellpart() > 0:
+                self.create_control_shell(
+                    wrpang=0,
+                    esort=1,
+                    irnxx=0,
+                    istupd=0,
+                    theory=0,
+                    bwc=1,
+                    miter=1,
+                    proj=1,
+                    irquad=0,
+                )
+        elif self.analysis == 2:
+            self.set_accuracy(
+                objective_stress_updates=Switch.ON,
+                invariant_node_number=InvariantNode.ON_FOR_SHELL_TSHELL_SOLID,
+                implicit_accuracy_flag=Switch.ON,
+            )
+            self.set_bulk_viscosity(bulk_viscosity_type=BulkViscosity.COMPUTE_INTERNAL_ENERGY_DISSIPATED)
+            self.set_energy(
+                hourglass_energy=EnergyFlag.COMPUTED,
+                sliding_interface_energy=EnergyFlag.COMPUTED,
+            )
+            self.set_hourglass(
+                controltype=HourglassControl.FLANAGAN_BELYTSCHKO_INTEGRATION_SOLID,
+                coefficient=0,
+            )
+            if self.parts.get_num_shellpart() > 0:
+                self.create_control_shell(
+                    wrpang=0,
+                    esort=1,
+                    irnxx=0,
+                    istupd=4,
+                    theory=0,
+                    bwc=1,
+                    miter=1,
+                    proj=1,
+                    irquad=0,
+                )
+        else:
+            pass
         DynaBase.save_file(self)
 
 

tests/test_dynamech.py

@@ -12,6 +12,7 @@
     ContactCategory,
     Curve,
     Point,
+    AnalysisType
 )
 
 
@@ -34,7 +35,7 @@ def test_mech(mech_initialfile, resolve_output_path, resolve_standard_path):
     fns = []
     fns.append(mech_initialfile)
     solution.open_files(fns)
-    mech = DynaMech()
+    mech = DynaMech(analysis=AnalysisType.IMPLICIT)
     solution.add(mech)
     solution.set_termination(0.03)
     airbag = Airbag(

tests/test_dynasolution.py

@@ -4,7 +4,8 @@
 from ansys.dyna.core.pre.dynasolution import DynaSolution
 from ansys.dyna.core.pre.dynamech import (
     DynaMech,
-    Transform
+    Transform,
+    AnalysisType
 )
 
 
@@ -48,7 +49,7 @@ def test_elementary_transform(resolve_solution_path, resolve_server_path, resolv
     fns = []
     fns.append(solution_initialfile)
     solution.open_files(fns)
-    obj = DynaMech()
+    obj = DynaMech(analysis=AnalysisType.IMPLICIT)
     solution.add(obj)
 
     obj.set_transform(

tests/test_dynathermal.py

@@ -8,6 +8,7 @@
     SolidPart,
     SolidFormulation,
     NodeSet,
+    AnalysisType
 )
 from ansys.dyna.core.pre.dynamaterial import MatElasticPlasticThermal
 
@@ -32,7 +33,7 @@ def test_thermal(thermal_initialfile, resolve_output_path, resolve_standard_path
     fns.append(thermal_initialfile)
     solution.open_files(fns)
     solution.set_termination(3.0)
-    ts = DynaMech()
+    ts = DynaMech(analysis=AnalysisType.IMPLICIT)
     solution.add(ts)
 
     tanalysis = ThermalAnalysis()