from abapy.indentation import ContactData from matplotlib import pyplot as plt # Axi contact cd = ContactData() x = [0, 1, 2, 3] alt = [-1, .1, 0, 0] press = [1, 0, 0, 0] cd.add_data(x, altitude=alt, pressure=press) points_axi, alt_axi, press_axi, conn_axi = cd.get_3D_data(axi_repeat=20) # 3D contact cd = ContactData(repeat=3, is_3D=True) k = np.cos(np.radians(60)) p = np.sin(np.radians(60)) x = [0, 4, 10, k * 4, k * 10] y = [0, 0, 0, p * 4, p * 10] alt = [-1, 0, 0, 0, 0] cd.add_data(x, altitude=alt) points_3D, alt_3D, press_3D, conn_3D = cd.get_3D_data() fig = plt.figure() plt.clf() ax1 = fig.add_subplot(121) ax1.set_aspect('equal') plt.tricontourf(points_axi[:, 0], points_axi[:, 1], conn_axi, alt_axi) plt.triplot(points_axi[:, 0], points_axi[:, 1], conn_axi) plt.title('axi contact') ax1 = fig.add_subplot(122) ax1.set_aspect('equal')
import numpy as np from abapy.indentation import ContactData X = np.linspace(-3., 3., 512) Y = np.linspace(-3., 3., 512) X, Y = np.meshgrid(X, Y) # Axi cd = ContactData() x = [0, 1, 2, 10] alt = [-1,.1, 0, 0] press = [1, 0, 0, 0] cd.add_data(x, altitude = alt, pressure = press) Alt_axi, Press_axi = cd.interpolate(X, Y, method ='linear') area = cd.contact_area() # 3D cd = ContactData(repeat = 3, is_3D = True) k = np.cos(np.radians(60)) p = np.sin(np.radians(60)) x = [0, 4, 10, k*4, k*10] y = [0, 0, 0, p*4, p*10] alt = [-1, 0, 0, 0, 0] cd.add_data(x, altitude = alt) Alt_3D, Press_3D = cd.interpolate(X, Y, method ='linear') from matplotlib import pyplot as plt fig = plt.figure() plt.clf()
from abapy.indentation import ContactData from matplotlib import pyplot as plt # Axi contact cd = ContactData() x = [0, 1, 2, 3] alt = [-1,.1, 0, 0] press = [1, 0, 0, 0] cd.add_data(x, altitude = alt, pressure = press) points_axi, alt_axi, press_axi, conn_axi = cd.get_3D_data(axi_repeat = 20) # 3D contact cd = ContactData(repeat = 3, is_3D = True) k = np.cos(np.radians(60)) p = np.sin(np.radians(60)) x = [0, 4, 10, k*4, k*10] y = [0, 0, 0, p*4, p*10] alt = [-1, 0, 0, 0, 0] cd.add_data(x, altitude = alt) points_3D, alt_3D, press_3D, conn_3D = cd.get_3D_data() fig = plt.figure() plt.clf() ax1 = fig.add_subplot(121) ax1.set_aspect('equal') plt.tricontourf(points_axi[:,0], points_axi[:,1], conn_axi, alt_axi)