osc2.py 2.0 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273
  1. """As osc.py, but without wheels."""
  2. from pysketcher import *
  3. L = 12.
  4. H = L/6
  5. W = L/6
  6. xmax = L
  7. drawing_tool.set_coordinate_system(xmin=-L, xmax=xmax,
  8. ymin=-1, ymax=L+H,
  9. axis=False,
  10. instruction_file='tmp_mpl.py')
  11. x = 0
  12. drawing_tool.set_linecolor('black')
  13. def make_dashpot(x):
  14. d_start = (-L,2*H)
  15. d = Dashpot(start=d_start, total_length=L+x, width=W,
  16. bar_length=3*H/2, dashpot_length=L/2, piston_pos=H+x)
  17. d.rotate(-90, d_start)
  18. return d
  19. def make_spring(x):
  20. s_start = (-L,4*H)
  21. s = Spring(start=s_start, length=L+x, bar_length=3*H/2, teeth=True)
  22. s.rotate(-90, s_start)
  23. return s
  24. d = make_dashpot(0)
  25. s = make_spring(0)
  26. M = Rectangle((0,H), 4*H, 4*H).set_linewidth(4)
  27. left_wall = Rectangle((-L,H),H/10,L-H).set_filled_curves(pattern='/')
  28. ground = Wall(x=[-L/2,L], y=[H,H], thickness=-H/10)
  29. fontsize = 18
  30. text_m = Text('$m$', (2*H, H+2*H), fontsize=fontsize)
  31. text_kx = Text('$s(u)$', (-L/2, H+4*H), fontsize=fontsize)
  32. text_bv = Text('$f(u)$', (-L/2, H), fontsize=fontsize)
  33. x_axis = Axis((2*H, L), H, '$u(t)$', fontsize=fontsize,
  34. label_spacing=(0.04, -0.01))
  35. x_axis_start = Line((2*H, L-H/4), (2*H, L+H/4)).set_linewidth(4)
  36. fig = Composition({
  37. 'spring': s, 'mass': M, 'left wall': left_wall,
  38. 'ground': ground,
  39. 'text_m': text_m, 'text_kx': text_kx,
  40. 'x_axis': x_axis, 'x_axis_start': x_axis_start})
  41. fig.draw()
  42. drawing_tool.display()
  43. drawing_tool.savefig('tmp_oscillator2_spring')
  44. drawing_tool.erase()
  45. fig['dashpot'] = d
  46. fig['text_bv'] = text_bv
  47. # or fig = Composition(dict(fig=fig, dashpot=d, text_bv=text_bv))
  48. fig.draw()
  49. drawing_tool.display()
  50. drawing_tool.savefig('tmp_oscillator2')
  51. drawing_tool.erase()
  52. F_force = Force((4*H, H+2*H), (4*H+H, H+2*H), '$F(t)$',
  53. text_spacing=(0.035, -0.01), text_pos='end', fontsize=fontsize)
  54. fig['F_force'] = F_force
  55. fig.draw()
  56. drawing_tool.savefig('tmp_oscillator2_force')
  57. raw_input()