Skip to content

Commit 9270bea

Browse files
rhorvath02rhorvath02
committed
Saving misc studies for Abishek
1 parent 58c2aa9 commit 9270bea

File tree

5 files changed

+149
-3
lines changed

5 files changed

+149
-3
lines changed

‎src/_2_misc_studies/DIL_tires.ipynb‎

Lines changed: 18 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -2,10 +2,25 @@
22
"cells": [
33
{
44
"cell_type": "code",
5-
"execution_count": 59,
5+
"execution_count": 4,
66
"id": "a7b1a216",
77
"metadata": {},
8-
"outputs": [],
8+
"outputs": [
9+
{
10+
"ename": "KeyError",
11+
"evalue": "'STRUCTURAL'",
12+
"output_type": "error",
13+
"traceback": [
14+
"\u001b[31m---------------------------------------------------------------------------\u001b[39m",
15+
"\u001b[31mKeyError\u001b[39m Traceback (most recent call last)",
16+
"\u001b[36mCell\u001b[39m\u001b[36m \u001b[39m\u001b[32mIn[4]\u001b[39m\u001b[32m, line 6\u001b[39m\n\u001b[32m 3\u001b[39m \u001b[38;5;28;01mimport\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[34;01mmatplotlib\u001b[39;00m\u001b[34;01m.\u001b[39;00m\u001b[34;01mpyplot\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[38;5;28;01mas\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[34;01mplt\u001b[39;00m\n\u001b[32m 4\u001b[39m \u001b[38;5;28;01mimport\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[34;01mnumpy\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[38;5;28;01mas\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[34;01mnp\u001b[39;00m\n\u001b[32m----> \u001b[39m\u001b[32m6\u001b[39m tire = \u001b[43mMF52\u001b[49m\u001b[43m(\u001b[49m\u001b[43mtire_name\u001b[49m\u001b[43m=\u001b[49m\u001b[33;43m\"\u001b[39;49m\u001b[33;43m\"\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mfile_path\u001b[49m\u001b[43m=\u001b[49m\u001b[33;43m\"\u001b[39;49m\u001b[33;43m../../src/_1_model_inputs/Modified_Round_8_Hoosier_R25B_16x7p5_10_on_7in_12psi_PAC02_UM2.tir\u001b[39;49m\u001b[33;43m\"\u001b[39;49m\u001b[43m)\u001b[49m\n",
17+
"\u001b[36mFile \u001b[39m\u001b[32m~/Documents/Github/simulation_toolkit/.venv/lib/python3.12/site-packages/LHR_tire_toolkit/MF52.py:19\u001b[39m, in \u001b[36mMF52.__init__\u001b[39m\u001b[34m(self, tire_name, file_path)\u001b[39m\n\u001b[32m 16\u001b[39m \u001b[38;5;28mself\u001b[39m._tire_params = _Processor(name = tire_name, file_path = file_path)\n\u001b[32m 18\u001b[39m \u001b[38;5;28mself\u001b[39m._init_coeffs()\n\u001b[32m---> \u001b[39m\u001b[32m19\u001b[39m \u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43m_init_consts\u001b[49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\n",
18+
"\u001b[36mFile \u001b[39m\u001b[32m~/Documents/Github/simulation_toolkit/.venv/lib/python3.12/site-packages/LHR_tire_toolkit/MF52.py:34\u001b[39m, in \u001b[36mMF52._init_consts\u001b[39m\u001b[34m(self)\u001b[39m\n\u001b[32m 32\u001b[39m \u001b[38;5;28mself\u001b[39m._operating_conds = \u001b[38;5;28mself\u001b[39m._tire_params.get_parameters(\u001b[33m\"\u001b[39m\u001b[33mTYRE_CONDITIONS\u001b[39m\u001b[33m\"\u001b[39m)\n\u001b[32m 33\u001b[39m \u001b[38;5;28mself\u001b[39m._vertical_coeffs = \u001b[38;5;28mself\u001b[39m._tire_params.get_parameters(\u001b[33m\"\u001b[39m\u001b[33mVERTICAL\u001b[39m\u001b[33m\"\u001b[39m)\n\u001b[32m---> \u001b[39m\u001b[32m34\u001b[39m \u001b[38;5;28mself\u001b[39m.structural = \u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43m_tire_params\u001b[49m\u001b[43m.\u001b[49m\u001b[43mget_parameters\u001b[49m\u001b[43m(\u001b[49m\u001b[33;43m\"\u001b[39;49m\u001b[33;43mSTRUCTURAL\u001b[39;49m\u001b[33;43m\"\u001b[39;49m\u001b[43m)\u001b[49m\n",
19+
"\u001b[36mFile \u001b[39m\u001b[32m~/Documents/Github/simulation_toolkit/.venv/lib/python3.12/site-packages/tire_model/file_processing/_process_tir.py:12\u001b[39m, in \u001b[36m_Processor.get_parameters\u001b[39m\u001b[34m(self, parameter)\u001b[39m\n\u001b[32m 11\u001b[39m \u001b[38;5;28;01mdef\u001b[39;00m\u001b[38;5;250m \u001b[39m\u001b[34mget_parameters\u001b[39m(\u001b[38;5;28mself\u001b[39m, parameter: \u001b[38;5;28mstr\u001b[39m) -> \u001b[38;5;28mdict\u001b[39m:\n\u001b[32m---> \u001b[39m\u001b[32m12\u001b[39m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43m_tire\u001b[49m\u001b[43m[\u001b[49m\u001b[32;43m1\u001b[39;49m\u001b[43m]\u001b[49m\u001b[43m[\u001b[49m\u001b[43mparameter\u001b[49m\u001b[43m]\u001b[49m\n",
20+
"\u001b[31mKeyError\u001b[39m: 'STRUCTURAL'"
21+
]
22+
}
23+
],
924
"source": [
1025
"from LHR_tire_toolkit.MF52 import MF52\n",
1126
"\n",
@@ -224,7 +239,7 @@
224239
"name": "python",
225240
"nbconvert_exporter": "python",
226241
"pygments_lexer": "ipython3",
227-
"version": "3.13.7"
242+
"version": "3.13.5"
228243
}
229244
},
230245
"nbformat": 4,

‎src/_2_misc_studies/abishek_gamma_Fx.py‎

Lines changed: 43 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
1+
from LHR_tire_toolkit.MF52 import MF52
2+
import matplotlib.pyplot as plt
3+
import numpy as np
4+
5+
# Load tire model
6+
tire = MF52(
7+
tire_name="",
8+
file_path="./src/_1_model_inputs/Modified_Round_8_Hoosier_R25B_16x7p5_10_on_7in_12psi_PAC02_UM2.tir"
9+
)
10+
11+
FZNOMIN = tire.get_FNOMIN()
12+
13+
# Sweeps
14+
kappa_sweep = np.linspace(-0.10, 0.10, 100)
15+
gamma_sweep = np.linspace(-4, 0, 5) * np.pi / 180
16+
17+
Fx_lst_lst = []
18+
19+
for gamma in gamma_sweep:
20+
Fx_lst = []
21+
for kappa in kappa_sweep:
22+
Fx = tire.tire_eval(FZ=FZNOMIN, alpha=1e-6, kappa=kappa, gamma=gamma)[0]
23+
Fx_lst.append(Fx)
24+
25+
Fx_lst_lst.append(Fx_lst)
26+
27+
# Create 3D plot
28+
fig = plt.figure(figsize=(10, 7))
29+
ax = fig.gca()
30+
31+
for Fx_vals, gamma_val in zip(Fx_lst_lst, gamma_sweep):
32+
print(Fx_vals)
33+
ax.plot(kappa_sweep, Fx_vals, label=rf"$\gamma = {gamma_val * 180 / np.pi}$")
34+
35+
ax.grid()
36+
ax.set_xlim([0, max(kappa_sweep)])
37+
ax.set_ylim([0, max(Fx_lst_lst[0]) * 1.125])
38+
ax.legend()
39+
40+
ax.set_title(r"$ F_{x} \: vs \: \kappa$" + f" at Fz = {FZNOMIN}")
41+
ax.set_xlabel(r"Slip Ratio, $\kappa$")
42+
ax.set_ylabel(r"Longitudinal Force, $F_{x}$ ($N$)")
43+
plt.show()

‎src/_2_misc_studies/abishek_gamma_Fy.py‎

Lines changed: 43 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
1+
from LHR_tire_toolkit.MF52 import MF52
2+
import matplotlib.pyplot as plt
3+
import numpy as np
4+
5+
# Load tire model
6+
tire = MF52(
7+
tire_name="",
8+
file_path="./src/_1_model_inputs/Round_8_Hoosier_R25B_16x7p5_10_on_7in_12psi_PAC02_UM2.tir"
9+
)
10+
11+
FZNOMIN = tire.get_FNOMIN()
12+
13+
# Sweeps
14+
alpha_sweep = np.linspace(-20, 20, 100) * np.pi / 180 # Slip angle (rad)
15+
gamma_sweep = np.linspace(-4, 0, 5) * np.pi / 180 # Camber angle (rad)
16+
17+
Fy_lst_lst = []
18+
19+
for gamma in gamma_sweep:
20+
Fy_lst = []
21+
for alpha in alpha_sweep:
22+
Fy = tire.tire_eval(FZ=FZNOMIN, alpha=alpha, kappa=0, gamma=gamma)[1]
23+
Fy_lst.append(Fy)
24+
25+
Fy_lst_lst.append(Fy_lst)
26+
27+
28+
# Create 3D plot
29+
fig = plt.figure(figsize=(10, 7))
30+
ax = fig.gca()
31+
32+
for Fy_vals, gamma_val in zip(Fy_lst_lst, gamma_sweep):
33+
ax.plot(alpha_sweep * 180 / np.pi, Fy_vals, label=rf"$\gamma = {gamma_val * 180 / np.pi}$")
34+
35+
ax.grid()
36+
ax.set_xlim([0, max(alpha_sweep) * 180 / np.pi])
37+
ax.set_ylim([0, max(Fy_lst_lst[0]) * 1.125])
38+
ax.legend()
39+
40+
ax.set_title(r"$ F_{y} \: vs \: \alpha$" + f" at Fz = {FZNOMIN}")
41+
ax.set_xlabel(r"Slip Angle, $\alpha$ ($deg$)")
42+
ax.set_ylabel(r"Lateral Force, $F_{y}$ ($N$)")
43+
plt.show()

‎src/_2_misc_studies/abishek_pneu_trail.py‎

Lines changed: 35 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,35 @@
1+
from LHR_tire_toolkit.MF52 import MF52
2+
import matplotlib.pyplot as plt
3+
import numpy as np
4+
5+
# Load tire model
6+
tire = MF52(
7+
tire_name="",
8+
file_path="./src/_1_model_inputs/Modified_Round_8_Hoosier_R25B_16x7p5_10_on_7in_12psi_PAC02_UM2.tir"
9+
)
10+
11+
# Sweeps
12+
Fz_sweep = np.linspace(100, 1600, 50) # vertical load [N]
13+
alpha_sweep = np.linspace(-30, 30, 50) * np.pi / 180 # slip angle [deg]
14+
15+
# Create meshgrid
16+
Alpha, Fz = np.meshgrid(alpha_sweep, Fz_sweep)
17+
18+
# Compute pneumatic trail
19+
pneu_trail = np.zeros_like(Alpha)
20+
for i in range(Fz.shape[0]):
21+
for j in range(Fz.shape[1]):
22+
pneu_trail[i, j] = tire.get_pneu_trail(FZ=Fz[i, j], alpha=Alpha[i, j])
23+
24+
# Plot
25+
fig = plt.figure(figsize=(10, 7))
26+
ax = fig.add_subplot(111, projection="3d")
27+
surf = ax.plot_surface(Alpha * 180 / np.pi, Fz, pneu_trail, cmap="viridis")
28+
29+
ax.set_xlabel("Slip Angle α [deg]")
30+
ax.set_ylabel("Vertical Load Fz [N]")
31+
ax.set_zlabel("Pneumatic Trail [m]")
32+
ax.set_title("Pneumatic Trail vs Slip Angle & Vertical Load")
33+
34+
fig.colorbar(surf, shrink=0.5, aspect=10, label="Pneumatic Trail [m]")
35+
plt.show()

‎src/_2_misc_studies/tire_eval.py‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
from LHR_tire_toolkit.MF52 import MF52
2+
3+
import matplotlib.pyplot as plt
4+
import numpy as np
5+
6+
tire = MF52(tire_name="", file_path="./src/_1_model_inputs/Modified_Round_8_Hoosier_R25B_16x7p5_10_on_7in_12psi_PAC02_UM2.tir")
7+
8+
tire_load = tire.tire_eval(FZ=600, alpha = 5 * np.pi / 180, kappa = 0.15, gamma = 0)
9+
10+
print(tire_load)

0 commit comments

Comments
 (0)