Merge branch 'master' of https://github.com/DSLsofMath/BScProj2018

Author: Oskar Lundström

Book/build.py

@@ -129,8 +129,8 @@ def build_index(sources):
         ("Single particle mechanics", "Physics/src/NewtonianMechanics/SingleParticle.lhs")
     ]),
     ("Examples", [
-        ("Box on an incline", "Physics/src/Examples/Box_incline.lhs"),
         ("Teeter", "Physics/src/Examples/Teeter.lhs"),
+        ("Box on an incline", "Physics/src/Examples/Box_incline.lhs"),
     ])
 ]
 

Physics/src/Examples/Box_incline.lhs

@@ -6,24 +6,31 @@ Box on an incline
 
 All vectors are in newton.
 
-
 ![Incline](incline.png){.float-img-left}
 
+Notation:
+fg = gravitational accelleration
+m = mass of box
+
 > fg = V2 0 (-10)
 > m = 2
-> alpha = 30
-> 
+
 > unit_normal :: Double -> Vector2 Double
 > unit_normal a = V2 (cos a) (sin a)
-> 
+
+Force against the incline from the box:
+
 > f_l_ :: Vector2 Double -> Angle -> Vector2 Double
 > f_l_ fa a = scale ((magnitude fa) * (cos a)) (unit_normal (a-(pi/2)))
-> 
+
+The normal against the incline:
+
 > fn :: Vector2 Double -> Angle -> Vector2 Double
 > fn fa a = negate (f_l_ fa a)
 
 Frictionfree incline:
 
+Force resultant:
 
 > fr :: Vector2 Double -> Angle -> Vector2 Double
 > fr fa a = (fn fa a) + fa
@@ -52,24 +59,14 @@ Odd:    Motsatt gravitation ger något underligt? Vi säger fortfarande att norm
 
 Good:   90* lutning - faller med G.
 
-
 *Main> fr (V2 0 10) (pi/3)
 
 (-4.330127018922194 x, 7.499999999999999 y)     
 
-
 *Main> fr (V2 0 10) (pi/4)
 
 (-5.0 x, 4.999999999999999 y)                   
 
-Good:   45* lutning - 5N både i x och y led. Pyth: 5^2 + 5^2 =/= 10^2
-
-                                                        det borde bli:  100 = a^2 + a^2
-                                                                        50 = a^2
-                                                                        5*sqrt(2) = a
-
-                                                        och det är det ju. xD
-
 *Main> fr (V2 0 10) (pi/6)
 
 (-4.330127018922193 x, 2.499999999999999 y)

Physics/src/Examples/Teeter.lhs

@@ -6,50 +6,72 @@ Exam excercise 3, 2017-01-13
 > import Dimensions.Quantity
 > import Prelude hiding (length)
 
-> balk_M = 1.0 # mass
+
+Two boxes, m1 and m2, rests on a beam in balance.
+
+Known values:
+
+> beam_M = 1.0 # mass
 > m1 = 2.0 # mass
 > m2 = 5.0 # mass
 > d = 0.75 # length
-> balk_L = 5.0 # length
+> beam_L = 5.0 # length
 > two = 2.0 # one
 
 ![Teeter](teeter.png){.float-img-left}
 
 Direct implication:
 
-$$ f(x) = \frac{a}{b}$$
+> beam_left_L = (beam_L /# two) +# d
+> beam_right_L = beam_L -# beam_left_L
 
-> balk_left_L = (balk_L /# two) +# d
-> balk_right_L = balk_L -# balk_left_L
+We want to be able to represenet the torques.
 
-We want to be able to represenet the angular momentums.
-Here are some propositions.
+A torque (sv. vridmoment) is defined as:
 
-> m1_vrid = m1 *# balk_left_L 
+$$ \tau = distance\ from\ turning\ point \cdot force $$
 
-> balk_L_vrid = ((balk_left_L /# balk_L) *# balk_M) *# (balk_left_L /# two)
+Since all force values will be composited of a mass and the gravitation, we can ignore the gravitation.
 
+$$ \tau = distance\ from\ turning\ point \cdot mass $$
 
-> balk_H_vrid = ((balk_right_L /# balk_L) *# balk_M) *# (balk_right_L /# two)
 
+> m1_torq = m1 *# beam_left_L 
 
-m2_vrid = m2 * x
+To get the beams torque on one side, we need to divide by 2 because the beam's torque is spread out linearly (the density of the beam is equal everywhere), which means the left parts mass centrum is \emph{half the distance} of the left parts total length.
 
-VL = HL
+$$ beamL_{\tau} = beamL_{M} \cdot \frac{distance}{2} $$
 
-m1_vrid + balk_L_vrid = m2_vrid + balk_H_vrid
+$$ beamL_{\tau} = \frac{beam\ left\ length}{beam\ length} \cdot beam_M \cdot \frac{beam\ left\ length}{2} $$
 
-m1_vrid + balk_L_vrid - balk_H_vrid  = m2_vrid
+> beamL_torq = ((beam_left_L /# beam_L) *# beam_M) *# (beam_left_L /# two)
 
-(m1_vrid + balk_L_vrid - balk_H_vrid) / m2  = x
+> beamR_torq = ((beam_right_L /# beam_L) *# beam_M) *# (beam_right_L /# two)
 
-> x = (m1_vrid +# balk_L_vrid -# balk_H_vrid) /# m2
+We make an expression for $m2_{\tau}$, which involves our unknown distance x.
 
-Security check:
+$$ m2_{\tau} = m2 \cdot x $$
+
+For the teeter to be in balance, both sides torques should be equal.
+
+$$ Left\ side\ torque = Right\ side\ angular\ torque $$
+
+We try to break out $m2_{\tau}$ and then x.
 
-> m2_vrid = m2 *# x
+$$ m1_{\tau} + beamL_{\tau} = m2_{\tau} + beamR_{\tau} $$
+
+$$ m1_{\tau} + beamL_{\tau} - beamR_{\tau} = m2_{\tau} $$
+
+$$ \frac{m1_{\tau} + beamL_{\tau} - beamR_{\tau}}{m2}  = x $$
+
+Our solution:
+
+> x = (m1_torq +# beamL_torq -# beamR_torq) /# m2
+
+Security check:
 
-> vL = m1_vrid +# balk_L_vrid
-> hL = m2_vrid +# balk_H_vrid
+> m2_torq = m2 *# x
 
+> left_side_torque = m1_torq +# beamL_torq
+> right_side_torque = m2_torq +# beamR_torq
 

README.md

@@ -4,6 +4,8 @@ Medlemmar: Oskar, Erik, ~~Daniel~~, Björn, Johan
 
 [Projektbeskrivning](DSLsofMath_andra_kurser.md)
 
+[Homepage](https://dslsofmath.github.io/BScProj2018/index.html)
+
 
 ## Deadlines
 * 9/2  : Inlämning av Planeringsrapport.

Redovisningar/Slutredovisning/Innehåll.md

@@ -62,8 +62,8 @@ Förmedla ett budskap;
 
 sammanfattning?
 
-* Ett projekt som vi valt att kalla "LYAP" och som utfördes av [NAMN]
 * Projektbeskrivning: [...] Ta fram DSLsofMath-inspirerat kompletterande material för andra närliggande kurser [...]
+* Ett projekt som vi valt att kalla "LYAP" och som utfördes av [NAMN]
 * Att lära ut fysik mha DSL
 * Ett läromaterial anpassat till datastudenter
 * Resultatet blev en hemsida/bok i stil med LYAH