Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
1 contributor

Users who have contributed to this file

44 lines (40 sloc) 1.32 KB
#!/usr/bin/env python
# -*- coding: UTF-8 -*-
# ea.py -- 2018 ckunte
import numpy as np
import matplotlib.pyplot as plt
def main():
d = np.linspace(60, 200) # Chain diameter (mm)
E = [
(5.6 - 0.0 * d) * 1E10, # R3/R4/R5 Gr stud (2018)
(5.4 - 0.004 * d) * 1E10, # R3 Gr. studless (2018)
(5.45 - 0.0025 * d) * 1E10, # R4 Gr. studless (2018)
(6.0 - 0.0033 * d) * 1E10 # R5 Gr. studless (2018)
] # Elasticities for studless chains
lbl = [
'Stud chain R3/R4/R5 (2018)',
'Studless chain R3 (2018)',
'Studless chain R4 (2018)',
'Studless chain R5 (2018)'
] # labels
for i, j in zip(E, lbl):
# plot: E vs. chain dia.
ax1 = plt.subplot(211)
plt.plot(d, i, label=j, linewidth=2)
plt.setp(ax1.get_xticklabels(), visible=False)
# plot: EA vs. chain dia.
EA = i * 2 * np.pi * (d / 1E3)**2 / (4 * 1E6) # (MN)
ax2 = plt.subplot(212)
plt.plot(d, EA, label=j, linewidth=2)
plt.setp(ax2.get_xticklabels(), fontsize=12)
ax1.legend(loc=0)
ax2.legend(loc=0)
ax1.set_ylabel('Effective elastic modulus, E (N/m$^2$)')
plt.ylabel('Axial stiffness, EA (MN)')
plt.xlabel('Chain diameter, d (mm)')
ax1.grid(True)
ax2.grid(True)
plt.show()
pass
if __name__ == '__main__':
main()
You can’t perform that action at this time.