KratosDemStructuresCouplingApplication 10.4.4__py3-none-any.whl

KratosMultiphysics/DemStructuresCouplingApplication/sp_statistics.py

@@ -0,0 +1,143 @@
+import matplotlib.pyplot as plt
+import numpy as np
+import h5py
+
+# Read HDF5 file
+file_name = 'sp_data.hdf5'
+f = h5py.File(file_name, 'r')
+
+# File Attributes
+test_id = f.attrs['test_id']
+internal_radius = f.attrs['internal_radius']
+external_radius = f.attrs['external_radius']
+interface_radius = f.attrs['interface_radius']
+thickness = f.attrs['thickness']
+volume = f.attrs['volume']
+real_probe_height = f.attrs['real_probe_height']
+target_porosity = f.attrs['target_porosity']
+porosity = f.attrs['porosity']
+density = f.attrs['density']
+
+# General factors
+height_factor = real_probe_height / thickness
+porosity_factor = (1.0-target_porosity)/(1.0-porosity)
+gram_factor = 1000.0
+
+# Total initial mass
+initial_radii = np.array(f['0'].get('radius'))
+initial_radii_3 = np.power(initial_radii,3)
+initial_masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * initial_radii_3
+initial_mass = np.sum(initial_masses)
+# print(initial_mass)
+
+# SP
+all_times=[]
+all_pressures=[]
+all_sps=[]
+
+# Time and external applied pressure
+failure_step = 175 # TODO
+max_step = 2*failure_step
+max_step = failure_step
+times = np.zeros(max_step+1)
+for i in range(max_step):
+    times[i+1] = f[str(i)].attrs['time']
+p_rate = 5e10 # TODO
+psi_factor = 0.000145038
+pressures = [p_rate*t*psi_factor for t in times] # Pressure in psi
+t_f=times[failure_step]
+# print(t_f)
+# Estimated time of initial sanding. TODO
+# t_is=times[90]
+# t_is=times[130]
+# print(t_is)
+
+# Maximum number of intact bonds to consider spheres as SP
+max_num_bonds = 6.0 # TODO: this should be an int
+
+for numbonds in np.arange(0.0,max_num_bonds,1.0):
+    sps = np.zeros(max_step+1)
+    for i in range(max_step):
+
+        # Read datasets
+        all_radii = np.array(f[str(i)].get('radius'))
+        continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds'))
+        xs = np.array(f[str(i)].get('x'))
+        ys = np.array(f[str(i)].get('y'))
+
+        # Separate internal and external spheres (weak region and strong region)
+        xs_2 = np.power(xs,2)
+        ys_2 = np.power(ys,2)
+        distance_2 = xs_2 + ys_2
+        # weak_radius = 0.5*(internal_radius+interface_radius)
+        # if times[i+1] < t_is:
+        #     weak_radius = internal_radius
+        # else:
+        #     weak_radius = internal_radius + (interface_radius-internal_radius)/(t_f-t_is)*(times[i+1]-t_is)
+        # TODO: count all the spheres in the domain
+        weak_radius = interface_radius
+        internal_radii = np.where(distance_2<weak_radius**2,all_radii,0.0) # spheres between inner wall (hole) and a certain radius (weak region)
+        external_radii = np.where(distance_2>=weak_radius**2,all_radii,0.0) # spheres between a certain radius and the interface wall (dem-fem wall) (strong region)
+
+        # Eliminate spheres that are free (SP), taking special care of those falling near the hole (internal spheres)
+        cont_internal_radii = np.where(continuum_bonds>numbonds,internal_radii,0.0) # We eliminate spheres with small number of bonds in the weak part of the probe
+        cont_external_radii = np.where(continuum_bonds>0.0,external_radii,0.0)
+
+        total_radii_3 = np.power(cont_internal_radii+cont_external_radii,3)
+
+        # Compute SP
+        masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * total_radii_3
+        current_mass = np.sum(masses)
+        current_sp = initial_mass - current_mass
+        sps[i+1] = current_sp
+
+    all_times.append(times)
+    all_pressures.append(pressures)
+    all_sps.append(sps)
+
+# Experiment data
+with open('ctw16_experiment.txt') as exp_f:
+    exp_times=[]
+    exp_pressures=[]
+    exp_sps=[]
+    for line in exp_f:
+        fields = line.strip().split()
+        if fields:
+            exp_times.append(float(fields[0]))
+            exp_pressures.append(float(fields[1]))
+            exp_sps.append(float(fields[2]))
+all_times.append(exp_times)
+all_pressures.append(exp_pressures)
+all_sps.append(exp_sps)
+
+# Graphs
+graph_name='sp_bonds_t.pdf'
+graph_name='sp_bonds.pdf'
+
+graph_labels=['SP up to 0 intact bonds',
+            'SP up to 1 intact bonds',
+            'SP up to 2 intact bonds',
+            'SP up to 3 intact bonds',
+            'SP up to 4 intact bonds',
+            'SP up to 5 intact bonds',
+            'ctw16 experiment',
+                ]
+
+f = plt.figure()
+
+for name, pressures, productions in zip(graph_labels, all_pressures, all_sps):
+    plt.plot(pressures, productions,label=name)
+#p_is = p_rate*t_is*psi_factor
+#plt.axvline(x=p_is, c='k', ls='--', label='estimated numerical initial sanding')
+
+# for name, times, productions in zip(graph_labels, all_times, all_sps):
+#     plt.plot(times, productions,label=name)
+# plt.axvline(x=t_is, c='k', ls='--', label='estimated numerical initial sanding')
+# plt.axvline(x=t_f, c='k', ls=':', label='numerical collapse')
+
+plt.legend(loc=2, prop={'size': 6})
+plt.xlabel('p (psi)')
+# plt.xlabel('Time (s)')
+plt.ylabel('Sand Production (g)')
+plt.title('SP depending on number of intact bonds')
+f.savefig(graph_name, bbox_inches='tight')

KratosMultiphysics/DemStructuresCouplingApplication/sp_statistics_2.py

@@ -0,0 +1,139 @@
+import matplotlib.pyplot as plt
+import numpy as np
+import h5py
+
+# Read HDF5 file
+file_name = 'sp_data.hdf5'
+f = h5py.File(file_name, 'r')
+
+# File Attributes
+test_id = f.attrs['test_id']
+internal_radius = f.attrs['internal_radius']
+external_radius = f.attrs['external_radius']
+interface_radius = f.attrs['interface_radius']
+thickness = f.attrs['thickness']
+volume = f.attrs['volume']
+real_probe_height = f.attrs['real_probe_height']
+target_porosity = f.attrs['target_porosity']
+porosity = f.attrs['porosity']
+density = f.attrs['density']
+
+# Time and external applied pressure
+failure_step = 300 # TODO
+max_step = failure_step
+times = np.zeros(max_step)
+for i in range(max_step):
+    times[i] = f[str(i)].attrs['time']
+
+delta_probe_radius = (interface_radius-internal_radius)/5.0
+# internal_radius_2 = internal_radius**2
+probe_radius_1 = (internal_radius+delta_probe_radius*1)**2
+probe_radius_2 = (internal_radius+delta_probe_radius*2)**2
+probe_radius_3 = (internal_radius+delta_probe_radius*3)**2
+probe_radius_4 = (internal_radius+delta_probe_radius*4)**2
+interface_radius_2 = interface_radius**2
+avg_num_intact_bonds_0 = np.zeros(max_step)
+avg_num_intact_bonds_1 = np.zeros(max_step)
+avg_num_intact_bonds_2 = np.zeros(max_step)
+avg_num_intact_bonds_3 = np.zeros(max_step)
+avg_num_intact_bonds_4 = np.zeros(max_step)
+
+avg_num_broken_bonds_0 = np.zeros(max_step)
+avg_num_broken_bonds_1 = np.zeros(max_step)
+avg_num_broken_bonds_2 = np.zeros(max_step)
+avg_num_broken_bonds_3 = np.zeros(max_step)
+avg_num_broken_bonds_4 = np.zeros(max_step)
+
+for i in range(max_step):
+
+    # Read datasets
+    continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds'))
+    initial_continuum_bonds = np.array(f[str(i)].get('initial_continuum_bonds'))
+    xs = np.array(f[str(i)].get('x'))
+    ys = np.array(f[str(i)].get('y'))
+
+    xs_2 = np.power(xs,2)
+    ys_2 = np.power(ys,2)
+    distance_2 = xs_2 + ys_2
+
+    num_intact_bonds_0 = np.where(distance_2<probe_radius_1,continuum_bonds,-100.0)
+    index_to_delete = np.where(num_intact_bonds_0<0.0)[0]
+    num_intact_bonds_0 = np.delete(num_intact_bonds_0,index_to_delete)
+    num_init_bonds_0 = np.where(distance_2<probe_radius_1,initial_continuum_bonds,-100.0)
+    index_to_delete = np.where(num_init_bonds_0<0.0)[0]
+    num_init_bonds_0 = np.delete(num_init_bonds_0,index_to_delete)
+    num_broken_bonds_0 = num_init_bonds_0 - num_intact_bonds_0
+    avg_num_intact_bonds_0[i] = np.sum(num_intact_bonds_0)/len(num_intact_bonds_0)
+    avg_num_broken_bonds_0[i] = np.sum(num_broken_bonds_0)/len(num_broken_bonds_0)
+
+    num_intact_bonds_1 = np.where(distance_2>=probe_radius_1,continuum_bonds,-100)
+    num_intact_bonds_1 = np.where(distance_2<probe_radius_2,num_intact_bonds_1,-100)
+    index_to_delete = np.where(num_intact_bonds_1<0.0)[0]
+    num_intact_bonds_1 = np.delete(num_intact_bonds_1,index_to_delete)
+    num_init_bonds_1 = np.where(distance_2>=probe_radius_1,initial_continuum_bonds,-100)
+    num_init_bonds_1 = np.where(distance_2<probe_radius_2,num_init_bonds_1,-100)
+    index_to_delete = np.where(num_init_bonds_1<0.0)[0]
+    num_init_bonds_1 = np.delete(num_init_bonds_1,index_to_delete)
+    num_broken_bonds_1 = num_init_bonds_1 - num_intact_bonds_1
+    avg_num_intact_bonds_1[i] = np.sum(num_intact_bonds_1)/len(num_intact_bonds_1)
+    avg_num_broken_bonds_1[i] = np.sum(num_broken_bonds_1)/len(num_broken_bonds_1)
+
+    num_intact_bonds_2 = np.where(distance_2>=probe_radius_2,continuum_bonds,-100)
+    num_intact_bonds_2 = np.where(distance_2<probe_radius_3,num_intact_bonds_2,-100)
+    index_to_delete = np.where(num_intact_bonds_2<0.0)[0]
+    num_intact_bonds_2 = np.delete(num_intact_bonds_2,index_to_delete)
+    num_init_bonds_2 = np.where(distance_2>=probe_radius_2,initial_continuum_bonds,-100)
+    num_init_bonds_2 = np.where(distance_2<probe_radius_3,num_init_bonds_2,-100)
+    index_to_delete = np.where(num_init_bonds_2<0.0)[0]
+    num_init_bonds_2 = np.delete(num_init_bonds_2,index_to_delete)
+    num_broken_bonds_2 = num_init_bonds_2 - num_intact_bonds_2
+    avg_num_intact_bonds_2[i] = np.sum(num_intact_bonds_2)/len(num_intact_bonds_2)
+    avg_num_broken_bonds_2[i] = np.sum(num_broken_bonds_2)/len(num_broken_bonds_2)
+
+    num_intact_bonds_3 = np.where(distance_2>=probe_radius_3,continuum_bonds,-100)
+    num_intact_bonds_3 = np.where(distance_2<probe_radius_4,num_intact_bonds_3,-100)
+    index_to_delete = np.where(num_intact_bonds_3<0.0)[0]
+    num_intact_bonds_3 = np.delete(num_intact_bonds_3,index_to_delete)
+    num_init_bonds_3 = np.where(distance_2>=probe_radius_3,initial_continuum_bonds,-100)
+    num_init_bonds_3 = np.where(distance_2<probe_radius_4,num_init_bonds_3,-100)
+    index_to_delete = np.where(num_init_bonds_3<0.0)[0]
+    num_init_bonds_3 = np.delete(num_init_bonds_3,index_to_delete)
+    num_broken_bonds_3 = num_init_bonds_3 - num_intact_bonds_3
+    avg_num_intact_bonds_3[i] = np.sum(num_intact_bonds_3)/len(num_intact_bonds_3)
+    avg_num_broken_bonds_3[i] = np.sum(num_broken_bonds_3)/len(num_broken_bonds_3)
+
+    num_intact_bonds_4 = np.where(distance_2>probe_radius_4,continuum_bonds,-100.0)
+    index_to_delete = np.where(num_intact_bonds_4<0.0)[0]
+    num_intact_bonds_4 = np.delete(num_intact_bonds_4,index_to_delete)
+    num_init_bonds_4 = np.where(distance_2>probe_radius_4,initial_continuum_bonds,-100.0)
+    index_to_delete = np.where(num_init_bonds_4<0.0)[0]
+    num_init_bonds_4 = np.delete(num_init_bonds_4,index_to_delete)
+    num_broken_bonds_4 = num_init_bonds_4 - num_intact_bonds_4
+    avg_num_intact_bonds_4[i] = np.sum(num_intact_bonds_4)/len(num_intact_bonds_4)
+    avg_num_broken_bonds_4[i] = np.sum(num_broken_bonds_4)/len(num_broken_bonds_4)
+
+
+# Graphs
+fig1, axs1 = plt.subplots(3, 2, figsize=(7, 7))
+
+graph_name='average_number_intact_bonds_probe_radius.pdf'
+
+axs1[0,0].plot(times,avg_num_intact_bonds_0)
+axs1[1,0].plot(times,avg_num_intact_bonds_1)
+axs1[2,0].plot(times,avg_num_intact_bonds_2)
+axs1[0,1].plot(times,avg_num_intact_bonds_3)
+axs1[1,1].plot(times,avg_num_intact_bonds_4)
+
+fig1.savefig(graph_name)
+
+fig2, axs2 = plt.subplots(3, 2, figsize=(7, 7))
+
+graph_name='average_number_broken_bonds_probe_radius.pdf'
+
+axs2[0,0].plot(times,avg_num_broken_bonds_0)
+axs2[1,0].plot(times,avg_num_broken_bonds_1)
+axs2[2,0].plot(times,avg_num_broken_bonds_2)
+axs2[0,1].plot(times,avg_num_broken_bonds_3)
+axs2[1,1].plot(times,avg_num_broken_bonds_4)
+
+fig2.savefig(graph_name)

KratosMultiphysics/DemStructuresCouplingApplication/stress_failure_check_utility.py

@@ -0,0 +1,36 @@
+import KratosMultiphysics
+import KratosMultiphysics.DemStructuresCouplingApplication as DemFem
+
+class StressFailureCheckUtility():
+    def __init__(self, model_part, test_number):
+
+        if not test_number:
+            return
+
+        self.model_part = model_part
+
+        if test_number == 1: # CTW16
+            minimum_radius = 0.00381
+            maximum_radius = 0.01381 #0.00501
+        elif test_number == 2: # CTW10
+            minimum_radius = 0.0127
+            maximum_radius = 0.0139
+        else: # Blind test
+            minimum_radius = 0.0381
+            maximum_radius = 0.0393
+
+        self.settings = KratosMultiphysics.Parameters( """
+        {
+            "cylinder_center": [0.0,0.0,0.0]
+        }  """ )
+
+        self.settings.AddEmptyValue("min_radius")
+        self.settings["min_radius"].SetDouble(minimum_radius)
+        self.settings.AddEmptyValue("max_radius")
+        self.settings["max_radius"].SetDouble(maximum_radius)
+
+        self.utility = DemFem.StressFailureCheckUtilities(self.model_part, self.settings)
+
+    def ExecuteFinalizeSolutionStep(self):
+
+        self.utility.ExecuteFinalizeSolutionStep()

kratosdemstructurescouplingapplication-10.4.4.dist-info/METADATA

@@ -0,0 +1,17 @@
+Metadata-Version: 2.5
+Name: KratosDemStructuresCouplingApplication
+Version: 10.4.4
+Summary:         KRATOS Multiphysics ("Kratos") is a framework for building parallel,         multi-disciplinary simulation software, aiming at modularity,         extensibility, and high performance. Kratos is written in C++,         and counts with an extensive Python interface.     
+Author-email: Kratos Team <kratos@listas.cimne.upc.edu>
+License: BSD-4-Clause
+Classifier: Programming Language :: C++
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering
+Requires-Python: >=3.8
+Requires-Dist: kratosdemapplication==10.4.4
+Requires-Dist: kratosmultiphysics==10.4.4
+Requires-Dist: kratosporomechanicsapplication==10.4.4
+Requires-Dist: kratosstructuralmechanicsapplication==10.4.4
+Description-Content-Type: text/markdown
+
+## DEM-Structures Coupling Application

kratosdemstructurescouplingapplication-10.4.4.dist-info/RECORD

@@ -0,0 +1,27 @@
+KratosMultiphysics/DemStructuresCouplingApplication/dem_fem_coupling_algorithm.py,sha256=7IDDCeDRlbtail89L0IZH3NBhxuR6sBdRUCA3Ir06d0,13360
+KratosMultiphysics/DemStructuresCouplingApplication/sand_production_post_process_tool.py,sha256=Bap4FcD41gCo2UTDcDPmVDUMjQprs1aZVnzdOEM_Fjw,7151
+KratosMultiphysics/DemStructuresCouplingApplication/dem_structures_coupling_gid_output.py,sha256=6EFyVd8cWJQmthpKRpjsomqYMTQskW4WJbB3Zp6Z4Fw,8749
+KratosMultiphysics/DemStructuresCouplingApplication/__init__.py,sha256=yBBXD5CRVs_ZpazcqKkGfZagku3aqOJo07uBLLMKa8k,308
+KratosMultiphysics/DemStructuresCouplingApplication/sp_statistics.py,sha256=hzFLbpvmFjEVXLYY0rg0iubBVlQdoUJLvwsY3_xF4UU,5052
+KratosMultiphysics/DemStructuresCouplingApplication/dem_poromechanics_coupling_analysis_stage.py,sha256=E4_8t9dIWRKnHzxUHaWxWo21Hxpb6Fm9dduPp1ftrQI,24071
+KratosMultiphysics/DemStructuresCouplingApplication/stress_failure_check_utility.py,sha256=O648a86VVioFDlGW6p025bEofGSc6lRQ3lPXPw-JkOw,1140
+KratosMultiphysics/DemStructuresCouplingApplication/multiaxial_control_module_fem_dem_generalized_2d_utility.py,sha256=6aXqVXZACNWu8vfCqzZ-ielD4E2OxLWTYfAw3ARUczg,6652
+KratosMultiphysics/DemStructuresCouplingApplication/sp_dem_fem_launcher.py,sha256=qVjl6yo3D6u9h9RFd4EFnnXc59Kyya__39ExQrvJFSc,389
+KratosMultiphysics/DemStructuresCouplingApplication/dem_main_script_ready_for_coupling_with_fem.py,sha256=u_YJ0dmcWYze8YB42GX1Zi2ZJRcC-B8M2KwFF4OljhU,1114
+KratosMultiphysics/DemStructuresCouplingApplication/sp_dem_fem_coupling_algorithm.py,sha256=wVO6mSRUJEHLiHwfPQgoEOpYCriF3SW__H9XVXT1IG4,15991
+KratosMultiphysics/DemStructuresCouplingApplication/sp_statistics_2.py,sha256=ggqH79yBPaHqgKREbwwp44BcgtycG7pG_FTR77ZdQ4A,6398
+KratosMultiphysics/DemStructuresCouplingApplication/control_module_fem_dem_utility.py,sha256=4kFvdo_E230dsQ2eQCTOIuPms-_fblyRlodzLGFxL5c,4825
+KratosMultiphysics/DemStructuresCouplingApplication/dem_fem_launcher.py,sha256=RWW69NOz6fMY-eZGFUVeXA4hZc7anPqT7xh5LJqP6oA,382
+KratosMultiphysics/DemStructuresCouplingApplication/control_module_process.py,sha256=obzmvT_ssEvG5h8Y6u_ilprOKxhjssLpVk-Y735zYdc,1358
+KratosMultiphysics/.libs/libKratosDemStructuresCouplingCore.dylib,sha256=raV8uRCZ13tbyq_4jzVlOf8g29Wd0vAb6XSvGTXsPys,2494592
+KratosMultiphysics/.libs/KratosDemStructuresCouplingApplication.so,sha256=RJo8KysypkVpUwUthpFajwFmXBwOTok-EwskrRRaIMg,2226880
+kratosdemstructurescouplingapplication-10.4.4.dist-info/RECORD,,
+kratosdemstructurescouplingapplication-10.4.4.dist-info/WHEEL,sha256=QX0UaoWXQIAv8vp5R3do793urekDoJdJRoGCk8nXIBY,115
+kratosdemstructurescouplingapplication-10.4.4.dist-info/METADATA,sha256=qFboNEJQJ7wBjnUwXeVJwS7yog4yP61rg-XVBVhcYWE,882
+kratosdemstructurescouplingapplication.dylibs/libstdc++.6.dylib,sha256=Evj6w_OS_5KZkCtzPWv2krwVowjlzK5MLjAxRd0k8Gc,3403264
+kratosdemstructurescouplingapplication.dylibs/libKratosStructuralMechanicsCore.dylib,sha256=Qb8OCjNNRNzZkunYM8VNCKhDvseac40mhJm7ogg68yY,21304288
+kratosdemstructurescouplingapplication.dylibs/libKratosDEMCore.dylib,sha256=GS5_Y5GGPYEINWpXLXPln-6lplVdp6_RNtUxcpJfJW0,9988256
+kratosdemstructurescouplingapplication.dylibs/libKratosCore.dylib,sha256=4ezGluhB6_qb9DnuWJCmX2VtpS9wWe16aOe9HLXVMSo,63629792
+kratosdemstructurescouplingapplication.dylibs/libKratosPoromechanicsCore.dylib,sha256=LQQiPBK9v22upHeTQgul2dXFTX1gygPU7R584JD6zr0,15479760
+kratosdemstructurescouplingapplication.dylibs/libgcc_s.1.1.dylib,sha256=LEUPpUC5mgEGvZcGWJjec9KZLBW-qoyUL1QxmuVKXDw,219232
+kratosdemstructurescouplingapplication.dylibs/libgomp.1.dylib,sha256=uQnW1i3Whh3qhRZgpBI1xCtHg5C5iHGCybc19JwuAGs,438880

kratosdemstructurescouplingapplication-10.4.4.dist-info/WHEEL

@@ -0,0 +1,6 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.32.4
+Root-Is-Purelib: true
+Tag: py3-none-any
+Generator: delocate 0.13.0
+

kratosdemstructurescouplingapplication.dylibs/libKratosCore.dylib

@@ -0,0 +1,4 @@
+[diffend] Oversized file quarantined before diffing.
+name: kratosdemstructurescouplingapplication.dylibs/libKratosCore.dylib
+size: 63629792 bytes
+sha256: e1ecc696e841ebfa9bf439ee5890a65f656da52f7059ed7a68e7bd1cb5d5312a