add icfd model:weak fsi

Author: zhangzhanqun

examples/ICFD/icfd_weak_fsi.py

@@ -0,0 +1,98 @@
+"""
+Weak FSI
+========
+
+This example shows a simple FSI coupling problem using weak/loose/explicit FSI.
+"""
+
+import os
+import sys
+
+
+from ansys.dyna.core.pre.dynasolution import DynaSolution
+from ansys.dyna.core.pre.dynaicfd import (
+    DynaICFD,
+    MatICFD,
+    ICFDPart,
+    ICFDDOF,
+    Curve,
+    ICFDVolumePart,
+    MeshedVolume,
+    ICFDAnalysis,
+    ShellPart,
+    ShellFormulation,
+    PartSet,
+    Curve,
+    DOF,
+    Motion
+)
+from ansys.dyna.core.pre.dynamaterial import MatRigid
+from ansys.dyna.core.pre import examples
+
+
+hostname = "localhost"
+if len(sys.argv) > 1:
+    hostname = sys.argv[1]
+
+solution = DynaSolution(hostname)
+# Import the initial mesh data(nodes and elements)
+fns = []
+path = examples.weak_fsi + os.sep
+fns.append(path + "weak_fsi.k")
+solution.open_files(fns)
+solution.set_termination(termination_time=40)
+icfd = DynaICFD()
+solution.add(icfd)
+
+icfd.set_timestep(tssfac=0.9)
+
+icfdanalysis = ICFDAnalysis()
+icfdanalysis.set_timestep(0.05)
+icfd.add(icfdanalysis)
+
+# define model
+mat = MatICFD(flow_density=1.0, dynamic_viscosity=0.005)
+
+part_inflow = ICFDPart(1)
+part_inflow.set_material(mat)
+part_inflow.set_prescribed_velocity(dof=ICFDDOF.X, motion=Curve(x=[0, 5, 6, 10000], y=[0, 0, 1, 1]))
+part_inflow.set_prescribed_velocity(dof=ICFDDOF.Y, motion=Curve(x=[0, 10000], y=[0, 0]))
+icfd.parts.add(part_inflow)
+
+part_outflow = ICFDPart(2)
+part_outflow.set_material(mat)
+part_outflow.set_prescribed_pressure(pressure=Curve(x=[0, 10000], y=[0, 0]))
+icfd.parts.add(part_outflow)
+
+part_symmetric = ICFDPart(3)
+part_symmetric.set_material(mat)
+part_symmetric.set_free_slip()
+icfd.parts.add(part_symmetric)
+
+part_wall = ICFDPart(4)
+part_wall.set_material(mat)
+part_wall.set_non_slip()
+part_wall.set_fsi()
+part_wall.compute_drag_force()
+part_wall.set_boundary_layer(number=3)
+icfd.parts.add(part_wall)
+
+partvol = ICFDVolumePart(surfaces=[1, 2, 3, 4])
+partvol.set_material(mat)
+icfd.parts.add(partvol)
+# define the volume space that will be meshed,The boundaries
+# of the volume are the surfaces "spids"
+meshvol = MeshedVolume(surfaces=[1, 2, 3, 4])
+icfd.add(meshvol)
+
+#define rigid cylinder
+matrigid = MatRigid(mass_density=1000,young_modulus=2e11,poisson_ratio=0.3)
+cylinder = ShellPart(1)
+cylinder.set_material(matrigid)
+cylinder.set_element_formulation(ShellFormulation.PLANE_STRESS)
+icfd.parts.add(cylinder)
+# Define boundary conddition
+icfd.boundaryconditions.create_imposed_motion(PartSet([1]),Curve(func="2*3.14/10*sin(2*3.14/10*TIME+3.14/2)"),dof=DOF.Y_TRANSLATIONAL,motion=Motion.VELOCITY)
+
+solution.create_database_binary(dt=0.2)
+solution.save_file()

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

@@ -339,6 +339,11 @@ def __init__(self, channel):
                 request_serializer=kwprocess__pb2.DefineCurveRequest.SerializeToString,
                 response_deserializer=kwprocess__pb2.DefineCurveReply.FromString,
                 )
+        self.CreateDefineCurveFunction = channel.unary_unary(
+                '/kwgrpc.kwC2S/CreateDefineCurveFunction',
+                request_serializer=kwprocess__pb2.DefineCurveFunctionRequest.SerializeToString,
+                response_deserializer=kwprocess__pb2.DefineCurveFunctionReply.FromString,
+                )
         self.CreateDefineVector = channel.unary_unary(
                 '/kwgrpc.kwC2S/CreateDefineVector',
                 request_serializer=kwprocess__pb2.DefineVectorRequest.SerializeToString,
@@ -1170,6 +1175,12 @@ def CreateDefineCurve(self, request, context):
         context.set_details('Method not implemented!')
         raise NotImplementedError('Method not implemented!')
 
+    def CreateDefineCurveFunction(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 CreateDefineVector(self, request, context):
         """Missing associated documentation comment in .proto file."""
         context.set_code(grpc.StatusCode.UNIMPLEMENTED)
@@ -2017,6 +2028,11 @@ def add_kwC2SServicer_to_server(servicer, server):
                     request_deserializer=kwprocess__pb2.DefineCurveRequest.FromString,
                     response_serializer=kwprocess__pb2.DefineCurveReply.SerializeToString,
             ),
+            'CreateDefineCurveFunction': grpc.unary_unary_rpc_method_handler(
+                    servicer.CreateDefineCurveFunction,
+                    request_deserializer=kwprocess__pb2.DefineCurveFunctionRequest.FromString,
+                    response_serializer=kwprocess__pb2.DefineCurveFunctionReply.SerializeToString,
+            ),
             'CreateDefineVector': grpc.unary_unary_rpc_method_handler(
                     servicer.CreateDefineVector,
                     request_deserializer=kwprocess__pb2.DefineVectorRequest.FromString,
@@ -3557,6 +3573,23 @@ def CreateDefineCurve(request,
             options, channel_credentials,
             insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
 
+    @staticmethod
+    def CreateDefineCurveFunction(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/CreateDefineCurveFunction',
+            kwprocess__pb2.DefineCurveFunctionRequest.SerializeToString,
+            kwprocess__pb2.DefineCurveFunctionReply.FromString,
+            options, channel_credentials,
+            insecure, call_credentials, compression, wait_for_ready, timeout, metadata)
+
     @staticmethod
     def CreateDefineVector(request,
             target,

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

@@ -604,6 +604,18 @@ def CreateDefineCurve(self, request, context):
         print(msg)
         return kwprocess_pb2.DefineCurveReply(id=lcid)
 
+    def CreateDefineCurveFunction(self, request, context):
+        lcid = request.lcid
+        lcid = self.kwdproc.get_data(gdt.KWD_DEFINE_CURVE_LASTID) + 1
+        function = request.function
+        card1 = str(lcid)
+        card2 = function
+        newk = "*DEFINE_CURVE_FUNCTION\n" + card1 + "\n" + card2
+        self.kwdproc.newkeyword(newk)
+        msg = "DefineCurveFunction " + str(lcid) + "Created..."
+        print(msg)
+        return kwprocess_pb2.DefineCurveFunctionReply(id=lcid)
+
     def CreateDefineVector(self, request, context):
         title = request.title
         vid = request.vid

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

@@ -104,14 +104,18 @@ def create(self, stub):
 class Curve:
     """Define a curve [for example, load (ordinate value)] as a function of time."""
 
-    def __init__(self, sfo=1, x=[], y=[]):
+    def __init__(self, sfo=1, x=[], y=[],func=None):
         self.sfo = sfo
         self.abscissa = x
         self.ordinate = y
+        self.func = func
 
     def create(self, stub):
         """Create curve."""
-        ret = stub.CreateDefineCurve(DefineCurveRequest(sfo=self.sfo, abscissa=self.abscissa, ordinate=self.ordinate))
+        if self.func!=None:
+            ret = stub.CreateDefineCurveFunction(DefineCurveFunctionRequest( function=self.func))
+        else:
+            ret = stub.CreateDefineCurve(DefineCurveRequest(sfo=self.sfo, abscissa=self.abscissa, ordinate=self.ordinate))
         self.id = ret.id
         logging.info(f"Curve {self.id} defined...")
         return self.id
@@ -733,6 +737,7 @@ class ShellFormulation(Enum):
     BELYTSCHKO_TSAY = 2
     SR_HUGHES_LIU = 6
     FULLY_INTEGRATED_BELYTSCHKO_TSAY_MEMBRANE = 9
+    PLANE_STRESS = 12
 
 
 class IGAFormulation(Enum):

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

@@ -19,16 +19,6 @@ def __init__(self):
         self.timestep = 0
         self.termination = 1e28
 
-    def set_timestep(self, timestep=0):
-        """Set time step for the fluid problem.
-
-        Parameters
-        ----------
-        dt : float
-            Time step for the fluid problem.
-        """
-        self.timestep = timestep
-
     def set_termination(self, termination_time):
         """Set total time of simulation for the fluid problem.
 
@@ -483,6 +473,12 @@ def set_non_slip(self):
         logging.info("ICFD boundary nonslip Created...")
         return ret
 
+    def set_fsi(self):
+        """Define fluid surface will be considered in contact with the solid surfaces for fluid-structure interaction (FSI) analysis."""
+        ret = self.stub.ICFDCreateBdyFSI(ICFDBdyFSIRequest(pid=self.id))
+        logging.info("ICFD boundary FSI Created...")
+        return ret
+
     def compute_drag_force(self):
         """Enable the computation of drag forces over given surface parts of the model."""
         ret = self.stub.ICFDCreateDBDrag(ICFDDBDragRequest(pid=self.id))

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

@@ -24,6 +24,7 @@
 free_convection_flow = os.path.join(_module_path, "icfd", "free_convection_flow")
 dam_break = os.path.join(_module_path, "icfd", "dam_break")
 driven_cavity = os.path.join(_module_path, "icfd", "driven_cavity")
+weak_fsi = os.path.join(_module_path, "icfd", "weak_fsi")
 
 iga_sample = os.path.join(_module_path, "iga", "iga_sample")
 camry_rc = os.path.join(_module_path, "implicit", "camry_rc")