add EM model:resistive heating 2d

Author: zhangzhanqun

doc/source/_static/pre/em/resistive_heating_2d.png

@@ -54,8 +54,8 @@
 
 matelastic1 = MatElastic(mass_density=8000, young_modulus=1e11, poisson_ratio=0.33)
 matelastic2 = MatElastic(mass_density=7000, young_modulus=1e11, poisson_ratio=0.33)
-matelastic1.set_electromagnetic_property(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=6e7)
-matelastic2.set_electromagnetic_property(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=4e6,eos=EMEOSTabulated1(Curve(x=[0, 25, 50, 100], y=[4e6, 4e6, 4e5, 4e5])))
+matelastic1.set_em_permeability_equal(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=6e7)
+matelastic2.set_em_permeability_equal(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=4e6,eos=EMEOSTabulated1(Curve(x=[0, 25, 50, 100], y=[4e6, 4e6, 4e5, 4e5])))
 
 matthermaliso1 = MatThermalIsotropic(density=8000,specific_heat=400, conductivity=400)
 matthermaliso2 = MatThermalIsotropic(density=7000,specific_heat=450, conductivity=40)

examples/EM/em_resistive_heating.py

@@ -0,0 +1,74 @@
+"""
+Resistive heating 2D
+====================
+
+This example features a simple metal bar where potential is imposed on both ends. \n
+LS-DYNA version : ls-dyna_smp_d_R13.1_138-g8429c8a10f_winx64_ifort190.exe
+"""
+
+import os
+import sys
+
+from ansys.dyna.core.pre.dynasolution import DynaSolution
+from ansys.dyna.core.pre.dynaem import (
+    DynaEM,
+    PartSet,
+    NodeSet,
+    Curve,
+    ShellPart,
+    ShellFormulation,
+    FEMSOLVER,
+    ThermalAnalysis,
+    ThermalAnalysisType,
+    EMType,
+    Isopotential_ConnType,
+    Isopotential
+)
+from ansys.dyna.core.pre.dynamaterial import MatRigid,MatThermalIsotropic,EMMATTYPE,EMEOSTabulated1
+from ansys.dyna.core.pre import examples
+# sphinx_gallery_thumbnail_path = '_static/pre/em/resistive_heating_2d.png'
+
+hostname = "localhost"
+if len(sys.argv) > 1:
+    hostname = sys.argv[1]
+
+solution = DynaSolution(hostname)
+fns = []
+path = examples.em_resistive_heating_2d + os.sep
+fns.append(path + "em_resistive_heating_2d.k")
+solution.open_files(fns)
+solution.set_termination(termination_time=0.0101)
+solution.create_database_binary(dt=1e-4)
+
+emobj = DynaEM()
+solution.add(emobj)
+
+emobj.set_timestep(tssfac=1,timestep_size_for_mass_scaled=1e-4)
+
+emobj.analysis.set_timestep(timestep=1e-4)
+emobj.analysis.set_em_solver(type=EMType.RESISTIVE_HEATING)
+emobj.analysis.set_solver_fem(solver=FEMSOLVER.DIRECT_SOLVER, relative_tol=1e-3)
+
+tanalysis = ThermalAnalysis()
+tanalysis.set_timestep(initial_timestep=1e-4)
+tanalysis.set_solver(analysis_type=ThermalAnalysisType.TRANSIENT)
+emobj.add(tanalysis)
+
+matrigid = MatRigid(mass_density=1, young_modulus=2e11)
+matrigid.set_em_resistive_heating_2d(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=1e4)
+
+matthermaliso = MatThermalIsotropic(density=100,specific_heat=10, conductivity=7)
+
+part = ShellPart(1)
+part.set_material(matrigid,matthermaliso)
+part.set_element_formulation(ShellFormulation.PLANE_STRESS)
+emobj.parts.add(part)
+
+emobj.boundaryconditions.create_imposed_motion(PartSet([1]), Curve(x=[0,10],y=[10,10]))
+emobj.set_init_temperature(temp=25)
+emobj.connect_isopotential(contype = Isopotential_ConnType.VOLTAGE_SOURCE,isopotential1 = Isopotential(NodeSet([521,517,513,509,525])),value=500)
+emobj.connect_isopotential(contype = Isopotential_ConnType.VOLTAGE_SOURCE,isopotential1 = Isopotential(NodeSet([585,605,625,564,565])))
+
+emobj.create_em_output(mats=2, matf=2, sols=2, solf=2)
+
+solution.save_file()

examples/EM/em_resistive_heating_2d.py

@@ -0,0 +1 @@
+python -m grpc_tools.protoc --python_out=. --grpc_python_out=. -I. .\kwprocess.proto

src/ansys/dyna/core/pre/Server/genpyprocs.bat

@@ -1695,6 +1695,7 @@ message EMMat001Request {
     int32 mid = 1;
 	int32 mtype = 2;
 	float sigma = 3;
+	int32 eosid = 4;
 }
 
 message EMMat001Reply {
@@ -1713,6 +1714,48 @@ message EMMat002Reply {
 	int32 answer = 1;
 }
 
+message EMMat004Request {
+    int32 mid = 1;
+	int32 mtype = 2;
+	float sigma = 3;
+}
+
+message EMMat004Reply {
+	int32 answer = 1;
+}
+
+message EMIsopotentialRequest {
+    int32 isoid = 1;
+	int32 settype = 2;
+	int32 setid = 3;
+}
+
+message EMIsopotentialReply {
+	int32 id = 1;
+}
+
+message EMIsopotentialConnectRequest {
+    int32 conid = 1;
+	int32 contype = 2;
+	int32 isoid1 = 3;
+	int32 isoid2 = 4;
+	float val = 5;
+}
+
+message EMIsopotentialConnectReply {
+	int32 id = 1;
+}
+
+message EMIsopotentialRogoRequest {
+    int32 isoid = 1;
+	int32 settype = 2;
+	int32 setid = 3;
+}
+
+message EMIsopotentialRogoReply {
+	int32 id = 1;
+}
+
 message EMSolverBemRequest {
     float reltol = 1;
 	int32 maxite = 2;
@@ -2025,6 +2068,10 @@ service kwC2S {
 	rpc CreateEMCircuit(EMCircuitRequest) returns (EMCircuitReply) {}
 	rpc CreateEMMat001(EMMat001Request) returns (EMMat001Reply) {}
 	rpc CreateEMMat002(EMMat002Request) returns (EMMat002Reply) {}
+	rpc CreateEMMat004(EMMat004Request) returns (EMMat004Reply) {}
+	rpc CreateEMIsopotential(EMIsopotentialRequest) returns (EMIsopotentialReply) {}
+	rpc CreateEMIsopotentialConnect(EMIsopotentialConnectRequest) returns (EMIsopotentialConnectReply) {}
+	rpc CreateEMIsopotentialRogo(EMIsopotentialRogoRequest) returns (EMIsopotentialRogoReply) {}
 	rpc CreateEMSolverBem(EMSolverBemRequest) returns (EMSolverBemReply) {}
 	rpc CreateEMSolverFem(EMSolverFemRequest) returns (EMSolverFemReply) {}
 	rpc CreateEMSolverBemMat(EMSolverBemMatRequest) returns (EMSolverBemMatReply) {}

src/ansys/dyna/core/pre/Server/kwprocess.proto

@@ -739,6 +739,26 @@ def __init__(self, channel):
                 request_serializer=kwprocess__pb2.EMMat002Request.SerializeToString,
                 response_deserializer=kwprocess__pb2.EMMat002Reply.FromString,
                 )
+        self.CreateEMMat004 = channel.unary_unary(
+                '/kwgrpc.kwC2S/CreateEMMat004',
+                request_serializer=kwprocess__pb2.EMMat004Request.SerializeToString,
+                response_deserializer=kwprocess__pb2.EMMat004Reply.FromString,
+                )
+        self.CreateEMIsopotential = channel.unary_unary(
+                '/kwgrpc.kwC2S/CreateEMIsopotential',
+                request_serializer=kwprocess__pb2.EMIsopotentialRequest.SerializeToString,
+                response_deserializer=kwprocess__pb2.EMIsopotentialReply.FromString,
+                )
+        self.CreateEMIsopotentialConnect = channel.unary_unary(
+                '/kwgrpc.kwC2S/CreateEMIsopotentialConnect',
+                request_serializer=kwprocess__pb2.EMIsopotentialConnectRequest.SerializeToString,
+                response_deserializer=kwprocess__pb2.EMIsopotentialConnectReply.FromString,
+                )
+        self.CreateEMIsopotentialRogo = channel.unary_unary(
+                '/kwgrpc.kwC2S/CreateEMIsopotentialRogo',
+                request_serializer=kwprocess__pb2.EMIsopotentialRogoRequest.SerializeToString,
+                response_deserializer=kwprocess__pb2.EMIsopotentialRogoReply.FromString,
+                )
         self.CreateEMSolverBem = channel.unary_unary(
                 '/kwgrpc.kwC2S/CreateEMSolverBem',
                 request_serializer=kwprocess__pb2.EMSolverBemRequest.SerializeToString,
@@ -1688,6 +1708,30 @@ def CreateEMMat002(self, request, context):
         context.set_details('Method not implemented!')
         raise NotImplementedError('Method not implemented!')
 
+    def CreateEMMat004(self, request, context):
+        """Missing associated documentation comment in .proto file."""
+        context.set_code(grpc.StatusCode.UNIMPLEMENTED)
+        context.set_details('Method not implemented!')
+        raise NotImplementedError('Method not implemented!')
+
+    def CreateEMIsopotential(self, request, context):
+        """Missing associated documentation comment in .proto file."""
+        context.set_code(grpc.StatusCode.UNIMPLEMENTED)
+        context.set_details('Method not implemented!')
+        raise NotImplementedError('Method not implemented!')
+
+    def CreateEMIsopotentialConnect(self, request, context):
+        """Missing associated documentation comment in .proto file."""
+        context.set_code(grpc.StatusCode.UNIMPLEMENTED)
+        context.set_details('Method not implemented!')
+        raise NotImplementedError('Method not implemented!')
+
+    def CreateEMIsopotentialRogo(self, request, context):
+        """Missing associated documentation comment in .proto file."""
+        context.set_code(grpc.StatusCode.UNIMPLEMENTED)
+        context.set_details('Method not implemented!')
+        raise NotImplementedError('Method not implemented!')
+
     def CreateEMSolverBem(self, request, context):
         """Missing associated documentation comment in .proto file."""
         context.set_code(grpc.StatusCode.UNIMPLEMENTED)
@@ -2483,6 +2527,26 @@ def add_kwC2SServicer_to_server(servicer, server):
                     request_deserializer=kwprocess__pb2.EMMat002Request.FromString,
                     response_serializer=kwprocess__pb2.EMMat002Reply.SerializeToString,
             ),
+            'CreateEMMat004': grpc.unary_unary_rpc_method_handler(
+                    servicer.CreateEMMat004,
+                    request_deserializer=kwprocess__pb2.EMMat004Request.FromString,
+                    response_serializer=kwprocess__pb2.EMMat004Reply.SerializeToString,
+            ),
+            'CreateEMIsopotential': grpc.unary_unary_rpc_method_handler(
+                    servicer.CreateEMIsopotential,
+                    request_deserializer=kwprocess__pb2.EMIsopotentialRequest.FromString,
+                    response_serializer=kwprocess__pb2.EMIsopotentialReply.SerializeToString,
+            ),
+            'CreateEMIsopotentialConnect': grpc.unary_unary_rpc_method_handler(
+                    servicer.CreateEMIsopotentialConnect,
+                    request_deserializer=kwprocess__pb2.EMIsopotentialConnectRequest.FromString,
+                    response_serializer=kwprocess__pb2.EMIsopotentialConnectReply.SerializeToString,
+            ),
+            'CreateEMIsopotentialRogo': grpc.unary_unary_rpc_method_handler(
+                    servicer.CreateEMIsopotentialRogo,
+                    request_deserializer=kwprocess__pb2.EMIsopotentialRogoRequest.FromString,
+                    response_serializer=kwprocess__pb2.EMIsopotentialRogoReply.SerializeToString,
+            ),
             'CreateEMSolverBem': grpc.unary_unary_rpc_method_handler(
                     servicer.CreateEMSolverBem,
                     request_deserializer=kwprocess__pb2.EMSolverBemRequest.FromString,
@@ -5013,6 +5077,74 @@ def CreateEMMat002(request,
             options, channel_credentials,
             insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
 
+    @staticmethod
+    def CreateEMMat004(request,
+            target,
+            options=(),
+            channel_credentials=None,
+            call_credentials=None,
+            insecure=False,
+            compression=None,
+            wait_for_ready=None,
+            timeout=None,
+            metadata=None):
+        return grpc.experimental.unary_unary(request, target, '/kwgrpc.kwC2S/CreateEMMat004',
+            kwprocess__pb2.EMMat004Request.SerializeToString,
+            kwprocess__pb2.EMMat004Reply.FromString,
+            options, channel_credentials,
+            insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
+
+    @staticmethod
+    def CreateEMIsopotential(request,
+            target,
+            options=(),
+            channel_credentials=None,
+            call_credentials=None,
+            insecure=False,
+            compression=None,
+            wait_for_ready=None,
+            timeout=None,
+            metadata=None):
+        return grpc.experimental.unary_unary(request, target, '/kwgrpc.kwC2S/CreateEMIsopotential',
+            kwprocess__pb2.EMIsopotentialRequest.SerializeToString,
+            kwprocess__pb2.EMIsopotentialReply.FromString,
+            options, channel_credentials,
+            insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
+
+    @staticmethod
+    def CreateEMIsopotentialConnect(request,
+            target,
+            options=(),
+            channel_credentials=None,
+            call_credentials=None,
+            insecure=False,
+            compression=None,
+            wait_for_ready=None,
+            timeout=None,
+            metadata=None):
+        return grpc.experimental.unary_unary(request, target, '/kwgrpc.kwC2S/CreateEMIsopotentialConnect',
+            kwprocess__pb2.EMIsopotentialConnectRequest.SerializeToString,
+            kwprocess__pb2.EMIsopotentialConnectReply.FromString,
+            options, channel_credentials,
+            insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
+
+    @staticmethod
+    def CreateEMIsopotentialRogo(request,
+            target,
+            options=(),
+            channel_credentials=None,
+            call_credentials=None,
+            insecure=False,
+            compression=None,
+            wait_for_ready=None,
+            timeout=None,
+            metadata=None):
+        return grpc.experimental.unary_unary(request, target, '/kwgrpc.kwC2S/CreateEMIsopotentialRogo',
+            kwprocess__pb2.EMIsopotentialRogoRequest.SerializeToString,
+            kwprocess__pb2.EMIsopotentialRogoReply.FromString,
+            options, channel_credentials,
+            insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
+
     @staticmethod
     def CreateEMSolverBem(request,
             target,

src/ansys/dyna/core/pre/Server/kwprocess_pb2.py

@@ -2597,7 +2597,8 @@ def CreateEMMat001(self, request, context):
         mid = request.mid
         mtype = request.mtype
         sigma = request.sigma
-        card1 = str(mid) + "," + str(mtype) + "," + str(sigma) + ",,,,1e28"
+        eosid = request.eosid
+        card1 = str(mid) + "," + str(mtype) + "," + str(sigma)+ ","+ str(eosid) + ",,,1e28"
         newk = "*EM_MAT_001\n" + card1
         self.kwdproc.newkeyword(newk)
         msg = "*EM_MAT_001 Created..."
@@ -2610,24 +2611,59 @@ def CreateEMMat002(self, request, context):
         sigma = request.sigma
         eosid = request.eosid
         murel = request.murel
-        card1 = (
-            str(mid)
-            + ","
-            + str(mtype)
-            + ","
-            + str(sigma)
-            + ","
-            + str(eosid)
-            + ","
-            + str(murel)
-            + ",0,1e28"
-        )
+        card1 = str(mid)+ ","+ str(mtype)+ ","+ str(sigma)+ ","+ str(eosid)+ ","+ str(murel)+ ",0,1e28"
         newk = "*EM_MAT_002\n" + card1
         self.kwdproc.newkeyword(newk)
         msg = "*EM_MAT_002 Created..."
         print(msg)
         return kwprocess_pb2.EMMat002Reply(answer=0)
+    
+    def CreateEMMat004(self, request, context):
+        mid = request.mid
+        mtype = request.mtype
+        sigma = request.sigma
+        card1 = str(mid)+ ","+ str(mtype)+ ","+ str(sigma)
+        newk = "*EM_MAT_004\n" + card1
+        self.kwdproc.newkeyword(newk)
+        msg = "*EM_MAT_004 Created..."
+        print(msg)
+        return kwprocess_pb2.EMMat004Reply(answer=0)
 
+    def CreateEMIsopotential(self, request, context):
+        isoid = self.kwdproc.get_data(gdt.KWD_EM_ISOPOTENTIAL_LASTID) + 1
+        settype = request.settype
+        setid = request.setid
+        card1 = str(isoid)+ ","+ str(settype)+ ","+ str(setid)
+        newk = "*EM_ISOPOTENTIAL\n" + card1
+        self.kwdproc.newkeyword(newk)
+        msg = "*EM_ISOPOTENTIAL Created..."
+        print(msg)
+        return kwprocess_pb2.EMIsopotentialReply(id=isoid)
+    
+    def CreateEMIsopotentialConnect(self, request, context):
+        conid = self.kwdproc.get_data(gdt.KWD_EM_ISOPOTCONNECT_LASTID) + 1
+        contype = request.contype
+        isoid1 = request.isoid1
+        isoid2 = request.isoid2
+        val = request.val
+        card1 = str(conid)+ ","+ str(contype)+ ","+ str(isoid1)+ ","+ str(isoid2)+ ","+ str(val)
+        newk = "*EM_ISOPOTENTIAL_CONNECT\n" + card1
+        self.kwdproc.newkeyword(newk)
+        msg = "*EM_ISOPOTENTIAL_CONNECT Created..."
+        print(msg)
+        return kwprocess_pb2.EMIsopotentialConnectReply(id=conid)
+    
+    def CreateEMIsopotentialRogo(self, request, context):
+        isoid = request.isoid
+        settype = request.settype
+        setid = request.setid
+        card1 = str(isoid)+ ","+ str(settype)+ ","+ str(setid)
+        newk = "*EM_ISOPOTENTIAL_ROGO\n" + card1
+        self.kwdproc.newkeyword(newk)
+        msg = "*EM_ISOPOTENTIAL_ROGO Created..."
+        print(msg)
+        return kwprocess_pb2.EMIsopotentialRogoReply(id=isoid)
+    
     def CreateEMSolverBem(self, request, context):
         reltol = request.reltol
         maxite = request.maxite