Skip to content
Snippets Groups Projects
Commit 1e01db5e authored by Van Dung NGUYEN's avatar Van Dung NGUYEN
Browse files

update

parent 964cdf72
Branches
No related tags found
No related merge requests found
import matplotlib.pyplot as plt
import pandas as pd
import os
import numpy as np
plt.rcParams.update({'font.size': 14})
plt.close("all")
for step in range(1,7):
fig, axes = plt.subplots(1,2)
all_files = [each for each in os.listdir(".") if "Tracking_GPFilterIPValueBounds" in each if f"step{step}" in each]
for ff in all_files:
data = pd.read_csv(ff,sep=";")
if not(data.empty):
axes[0].plot(data["xRef"].values,data["yRef"].values,"r.")
axes[0].plot(data["xCur"].values,data["yCur"].values,"g.")
axes[1].plot(data["xCur"].values,data["yCur"].values,"r.")
plt.show()
......@@ -99,6 +99,8 @@ mysolver.internalPointBuildViewIncrement("U_NORM",IPField.U_NORM, 1, 1)
mysolver.internalPointBuildViewRate("epl",IPField.PLASTICSTRAIN, 1, 1)
mysolver.internalPointBuildViewRate("U_NORM",IPField.U_NORM, 1, 1)
gpFilter = GPFilterIPValueBounds(IPField.PLASTICSTRAIN, 1e-5, 1)
mysolver.archivingIPValueTracking(IPField.PLASTICSTRAIN, 1, gpFilter)
mysolver.archivingForceOnPhysicalGroup("Face", 84, 1)
mysolver.archivingNodeDisplacement(19,1,1)
......@@ -108,9 +110,10 @@ mysolver.solve()
check = TestCheck()
check.equal(-2.859523e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-4)
"""
#newsolver = mysolver.clone(meshfile,1000)
newsolver = nonLinearMechSolver(mysolver)
newsolver.solve()
check = TestCheck()
check.equal(-2.859523e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-4)
"""
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment