-
Notifications
You must be signed in to change notification settings - Fork 0
/
matrix_to_gid_output.py
97 lines (77 loc) · 3.58 KB
/
matrix_to_gid_output.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
import numpy as np
import matplotlib.pyplot as plt
import KratosMultiphysics as KMP
import KratosMultiphysics.gid_output_process as GOP
import KratosMultiphysics.StructuralMechanicsApplication as SMA
def create_out_mdpa(model_part, file_name):
model_part.AddNodalSolutionStepVariable(KMP.DISPLACEMENT)
model_part.AddNodalSolutionStepVariable(KMP.REACTION)
import_flags = KMP.ModelPartIO.READ
KMP.ModelPartIO(file_name, import_flags).ReadModelPart(model_part)
def print_results_to_gid(model_part, snapshot_matrix, residuals_matrix):
gid_output = GOP.GiDOutputProcess(
model_part,
"PredictDiff",
KMP.Parameters("""
{
"result_file_configuration" : {
"gidpost_flags" : {
"GiDPostMode" : "GiD_PostAscii",
"WriteDeformedMeshFlag" : "WriteDeformed",
"WriteConditionsFlag" : "WriteConditions",
"MultiFileFlag" : "SingleFile"
},
"file_label" : "time",
"output_control_type" : "step",
"output_interval" : 1,
"body_output" : true,
"node_output" : false,
"skin_output" : false,
"plane_output" : [],
"nodal_results" : ["DISPLACEMENT","REACTION"],
"gauss_point_results" : [],
"nodal_nonhistorical_results" : []
},
"point_data_configuration" : []
}"""
)
)
gid_output.ExecuteInitialize()
print(snapshot_matrix.shape)
for ts in range(0, snapshot_matrix.shape[0]):
# for ts in range(1):
# ts=snapshot_matrix.shape[0]-1
# for k in range(1):
model_part.ProcessInfo[KMP.STEP] = ts+1
model_part.ProcessInfo[KMP.TIME] = ts+1
gid_output.ExecuteBeforeSolutionLoop()
gid_output.ExecuteInitializeSolutionStep()
snapshot = snapshot_matrix[ts]
residuals = residuals_matrix[ts][0]
# var_utils.SetSolutionStepValuesVector(model_part.Nodes, KMP.DISPLACEMENT, snapshot, 2)
# var_utils.SetSolutionStepValuesVector(model_part.Nodes, KMP.REACTION, residuals, 2)
i=0
c=2
for node in model_part.Nodes:
node.SetSolutionStepValue(KMP.DISPLACEMENT_X,0,snapshot[i*c+0])
node.SetSolutionStepValue(KMP.DISPLACEMENT_Y,0,snapshot[i*c+1])
node.SetSolutionStepValue(KMP.REACTION_X,0,residuals[i*c+0])
node.SetSolutionStepValue(KMP.REACTION_Y,0,residuals[i*c+1])
i+=1
# conditions_array=model_part.Conditions
# for i, condition in enumerate(conditions_array):
# condition.SetValue(SMA.LINE_LOAD, forces)
gid_output.PrintOutput()
gid_output.ExecuteFinalizeSolutionStep()
gid_output.ExecuteFinalize()
if __name__ == "__main__":
snapshots_matrix=np.load('Scipy_x_snapshots.npy')
# snapshots_matrix=np.load('FOM_snaps_30steps.npy')
reactions_matrix=np.load('Scipy_reactions_snapshots.npy')
# plt.plot(reactions_matrix[0][0])
# plt.show()
# exit()
current_model = KMP.Model()
model_part = current_model.CreateModelPart("main_model_part")
create_out_mdpa(model_part, "datasets_two_forces_dense_extended/dense_2forces_cantilever_thick_dense")
print_results_to_gid(model_part, snapshots_matrix, reactions_matrix)