1 /* Test suite for Mechanics.asy.
3 * Copyright (C) 2008-2009 W. Trevor King <wking@drexel.edu>
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 2 of the License, or
8 * (at your option) any later version.
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
15 * You should have received a copy of the GNU General Public License along
16 * with this program; if not, write to the Free Software Foundation, Inc.,
17 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
24 LabeledCircle lc = LabeledCircle(center=(-2u, 1u));
27 lc.center = (-2u, 0.5u);
29 lc.draw_label(Label("e", align=E));
34 lc.draw_label(Label("c", align=E));
35 lc.draw_label(rotate(90)*Label("e------I", align=S));
36 lc.draw_label(rotate(45)*Label("e------I", align=S));
37 lc.center = (-2u, -2u);
41 lc.center = (-2u, -3u);
45 Mass a = Mass(center=(0,0));
46 Mass b = Mass(center=(2u,1u), Label("$m_b$", align=N));
47 Mass c = Mass(center=(1u,-2u), Label("$m_c$", align=E));
49 Distance dab = Distance(b.center(), a.center(), scale=u,
50 L=Label("$r_{ab}$", embed=Shift));
51 Distance dac = Distance(a.center(), c.center(), scale=u, L="$r_{ac}$");
52 Distance ds[] = {dab, dac};
54 b.center(), a.center(), c.center(), radius=.7u, fill=red, L="$\theta_T'$");
56 b.center(), a.center(), c.center(), radius=-.5u, L="$\theta_T$");
58 b.center(), a.center(), c.center(), radius=1.5u,
59 L=Label("$\theta_T''$", position=Relative(0.3)));
60 Angle as[] = {bac1, bac2, bac3};
66 vs.push(Vector(center=(-1u,-4u), phi=90, "Out"));
67 vs.push(Vector(center=(-1u,-5.5u), phi=-90, "In"));
68 vs.push(Vector(center=(0, -4u), mag=1.5u, dir=-90, phi=60,
69 Label("60dg OOP", position=EndPoint)));
71 for (int i=0; i<as.length; i+=1)
73 for (int i=0; i<ms.length; i+=1)
75 for (int i=0; i<ds.length; i+=1)
77 for (int i=0; i<vs.length; i+=1)
80 Pendulum p = makePendulum(
81 pivot=(3.5u,-3u), mass=b, length=4u, angleDeg=-20,
82 angleL="$\rho$", stringL=Label("r", embed=Shift));
84 real len = abs(p.pivot.x-b.center().x);
85 Spring s = Spring(pFrom=b.center()-(2u,0), pTo=b.center(), L="$k_1$");
87 s = Spring(pFrom=b.center(), pTo=(p.pivot.x, b.center().y), L="$k_2$");
90 p.draw(drawVertical=true);
92 Surface s = Surface(pFrom=(0,-7u), pTo=(3.5u, -7u), L="Table");
95 Vector v1 = Vector((-2u, -7u), Label("$v_1$", align=E, position=EndPoint));
97 Vector v2 = Vector(v1.center, mag=10mm, dir=90, Label("$v_2$", align=W));
100 v3.label = Label("$v_3 = v_1 + v_2$", align=E, position=EndPoint);
103 v4.label = Label("$v_4 = v_1 - v_2$", align=S, position=EndPoint);
106 Vector v = Velocity();
107 pair vfc = (5u,0); // vector field center
113 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
114 *shift((-.5,-.5))*unitsquare, grey);
115 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
116 outline=green+dashed);
118 vfc -= (0, vfheight + u);
120 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
121 *shift((-.5,-.5))*unitsquare, grey);
122 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
123 outline=green+dashed);
125 vfc -= (0, vfheight + u);
127 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
128 *shift((-.5,-.5))*unitsquare, grey);
129 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
130 outline=green+dashed);
132 vfc += (vfwidth + u, 2*(vfheight + u));
134 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
135 *shift((-.5,-.5))*unitsquare, grey);
136 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
137 outline=green+dashed);
139 vfc -= (0, vfheight + u);
141 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
142 *shift((-.5,-.5))*unitsquare, grey);
143 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
144 outline=green+dashed);
146 vfc -= (0, vfheight + u);
148 draw(shift(vfc)*xscale(vfwidth-2*vfbuf)*yscale(vfheight-2*vfbuf)
149 *shift((-.5,-.5))*unitsquare, grey);
150 vector_field(vfc, width=vfwidth, height=vfheight, dv=vfdv, buf=vfbuf, v=v,
151 outline=green+dashed);