Modifica di Asse Trasversale Cerniera (sezione)

===='''Script Python: Effetto degli Errori di Localizzazione del Centro di Rotazione con on Apertura a 0 mm e Cuspidi Inclinate ''<math>5^\circ</math>'''''====
</blockquote><syntaxhighlight lang="python">
import numpy as np
import matplotlib.pyplot as plt

# Definition of the original data
original_center = np.array([11.56, 61.21])  # Original center of the circle
r = 26.38  # Radius of the original circle
theta = np.linspace(0, 2 * np.pi, 10)  # 10 points along the circumference

# Generation of original points on the circumference
points = np.array([
    original_center + r * np.array([np.cos(angle), np.sin(angle)])
    for angle in theta
])

# Definition of the radius for the position of dental cusps
cusp_radius = 75.0  # Distance of the cusp from the center of rotation

# Generation of the positions of the original dental cusps
cusp_positions = np.array([
    original_center + cusp_radius * np.array([np.cos(angle), np.sin(angle)])
    for angle in theta
])

# Simulation of the errors Δ
errors = np.linspace(0, 10, 10)  # Errors Δ from 0 to 10 mm, 10 steps
max_differences_inclined_cusp = []

# Definition of the angle of the dental cusps (in degrees)
cusp_angle_deg = 5  
cusp_angle_rad = np.deg2rad(cusp_angle_deg)  # Conversion to radians

# Calculation of cusp error for mandibular opening at 0 mm with inclined cusps
for delta in errors:
    # Calculation of vertical error due to cusp inclination
    vertical_error_due_to_inclination = delta * np.tan(cusp_angle_rad)
    
    # Save the maximum error in this specific case
    max_differences_inclined_cusp.append(vertical_error_due_to_inclination)
    
    print(f"Error Δ = {delta:.1f} mm: Max Difference in Cusp Positions with inclined cusps = {vertical_error_due_to_inclination:.2f} mm")

plt.plot(errors, max_differences_inclined_cusp, 'ro-', label=f'Inclined Cusps ({cusp_angle_deg}°)')
plt.xlabel('Error Δ (mm) in HA localization')
plt.ylabel('Max Difference in Cusp Positions (mm)')
plt.title(f'Effect of Rotation Center Localization Errors with Inclined Cusps ({cusp_angle_deg}°)')
plt.grid(True)
plt.legend()
plt.show()

</syntaxhighlight>

<br />