Optimized source code

Author: codingbychanche

M.png

@@ -1,5 +1,3 @@
-Beam Calc
-=========
 Calculates the supporting forces of a simply supported beam.
 
 Beam and acting loads
@@ -10,39 +8,57 @@ Beam and acting loads
 Sample output
 -------------
 	Beam:
-	Length=4.0
-	Supports sorted by distance from left end of beam:
-	Support A (Left)    xn=0.0
-	Support B (Right)    xn=4.0
-
-	Number of Loads:2
-	Loads sorted by distance from left end of beam:
-	Name:q1 Magnitute at start:-5.0 Point Load:-5.0 distance from left end:0.0 Lengtht:4.0  Force at end:-5.0 Resulting force:-20.0 is acting 2.0 m from left end of beam.
-	Name:F1 Magnitute at start:-2.0 Point Load:-2.0 distance from left end:2.0 Lengtht:0.0  Force at end:-2.0 Resulting force:-2.0 is acting 2.0 m from left end of beam.
-
-	RESULT:
-	Left support:10.707106781186548 N. Right bearing:10.707106781186548 N.
-	Horizontal force at right support:1.414213562373095
-	Term:N=2.0 N x sin(45.0) = 1.41N
-	Term:B (Right)=5.0N x 4.0m + 2.0N x cos(45.0) + (-10.71N) = 10.71N
-	Term:A (Left)=(5.0 N/m x 4.0m x 2.0m+2.0N x cos(45.0) x 2.0m)/4.0m = 10.71N
-
-	=== Shearing forces => Q ===
-	Local maxi-/  minima
-	Points of disconuity in Q
-	x=0.0 m   Q=10.706106781186547 N
-	x=2.0 m   Q=-0.7081067811881641 N
-	x=4.0 m   Q=-5.000000032335805E-4 N
-
-	=== Bending moment => M ===
-	Local maxi-/  minima
-	x=1.9998999999997962 m   M=11.410643101686322 Nm
-	Points of disconuity in M
-	x=0.0 m   M=0.0 Nm
-	x=1.9999999999997962 m   M=11.410289362442864 Nm
-
-	=== Normal forces => N ===
-	x=2.0 m   N[N]=-1.414213562373095 N
+Length=4.0
+Supports sorted by distance from left end of beam:
+Support A (Left)    xn=1.0
+Support B (Right)    xn=3.0
+
+Number of Loads:2
+Loads sorted by distance from left end of beam:
+Name:q1 Magnitute at start:-5.0 Point Load:-5.0 distance from left end:0.0 Lengtht:4.0  Force at end:-5.0 Resulting force:-20.0 is acting 2.0 m from left end of beam.
+Name:F1 Magnitute at start:-2.0 Point Load:-2.0 distance from left end:2.0 Lengtht:0.0  Force at end:-2.0 Resulting force:-2.0 is acting 2.0 m from left end of beam.
+
+RESULT:
+Left support:11.0 N. Right bearing:11.0 N.
+Horizontal force at right support:0.0
+Term:N=0.00N
+Term:B (Right)=5.0N x 4.0m + 2.0N + (-11.00N) = 11.00N
+Term:A (Left)=(5.0 N/m x 4.0m x 1.0m+2.0N x 1.0m)/2.0m = 11.00N
+
+=== Shearing forces => Q ===
+Local maxi-/  minima
+Points of disconuity in Q
+x=0.0 m   Q=-0.01 N
+x=1.0 m   Q=5.995000000000084 N
+x=2.0 m   Q=-1.0099999999998097 N
+x=3.0 m   Q=4.990000000000296 N
+x=3.9999999999996705 m   Q=-0.009999999999619792 N
+Zero points in Q
+x=0.9980000000000008 m
+x=1.9989999999998906 m
+x=2.9989999999997807 m
+x=3.9979999999996707 m
+
+=== Bending moment => M ===
+Local maxi-/  minima
+x=0.0 m   M=0.0 Nm
+x=0.9980000000000008 m   M=-2.5049725149999653 Nm
+x=1.9989999999998906 m   M=0.9970029950001726 Nm
+x=2.9989999999997807 m   M=-2.5269875049995836 Nm
+x=3.9979999999996707 m   M=-0.018966014999237325 Nm
+Points of disconuity in M
+x=0.9990000000000008 m   M=-2.488004504999965 Nm
+x=1.9999999999998905 m   M=0.991985000000173 Nm
+x=2.9999999999997806 m   M=-2.4890234999995826 Nm
+Zero points in M
+x=0.0 m
+x=1.5339999999999419 m
+x=2.4569999999998404 m
+x=3.9999999999996705 m
+
+=== Normal forces => N ===
+x=2.0 m   N[N]=0.0 NBeam Calc
+
 
 ![](LeadingSigns.png)
 

N.png

@@ -9,7 +9,40 @@ Beam and acting loads
 
 Sample output
 -------------
-![](Shot_1.png)
+	Beam:
+	Length=4.0
+	Supports sorted by distance from left end of beam:
+	Support A (Left)    xn=0.0
+	Support B (Right)    xn=4.0
+
+	Number of Loads:2
+	Loads sorted by distance from left end of beam:
+	Name:q1 Magnitute at start:-5.0 Point Load:-5.0 distance from left end:0.0 Lengtht:4.0  Force at end:-5.0 Resulting force:-20.0 is acting 2.0 m from left end of beam.
+	Name:F1 Magnitute at start:-2.0 Point Load:-2.0 distance from left end:2.0 Lengtht:0.0  Force at end:-2.0 Resulting force:-2.0 is acting 2.0 m from left end of beam.
+
+	RESULT:
+	Left support:10.707106781186548 N. Right bearing:10.707106781186548 N.
+	Horizontal force at right support:1.414213562373095
+	Term:N=2.0 N x sin(45.0) = 1.41N
+	Term:B (Right)=5.0N x 4.0m + 2.0N x cos(45.0) + (-10.71N) = 10.71N
+	Term:A (Left)=(5.0 N/m x 4.0m x 2.0m+2.0N x cos(45.0) x 2.0m)/4.0m = 10.71N
+
+	=== Shearing forces => Q ===
+	Local maxi-/  minima
+	Points of disconuity in Q
+	x=0.0 m   Q=10.706106781186547 N
+	x=2.0 m   Q=-0.7081067811881641 N
+	x=4.0 m   Q=-5.000000032335805E-4 N
+
+	=== Bending moment => M ===
+	Local maxi-/  minima
+	x=1.9998999999997962 m   M=11.410643101686322 Nm
+	Points of disconuity in M
+	x=0.0 m   M=0.0 Nm
+	x=1.9999999999997962 m   M=11.410289362442864 Nm
+
+	=== Normal forces => N ===
+	x=2.0 m   N[N]=-1.414213562373095 N
 
 ![](LeadingSigns.png)
 

Q.png

@@ -19,12 +19,12 @@ public static void main(String[] args) {
 
 		// Create a new beam
 		Beam myBeam = new Beam(4);
-		myBeam.addBearing(new Support("A (Left)", 0, Support.ROLLER_SUPPORT));
-		myBeam.addBearing(new Support("B (Right)", 4, Support.PIN_SUPPORT));
+		myBeam.addBearing(new Support("A (Left)", 1, Support.ROLLER_SUPPORT));
+		myBeam.addBearing(new Support("B (Right)", 3, Support.PIN_SUPPORT));
 
 		// Add load
 		// NAME/ Force/ Distance/ Angle /Length
-		myBeam.addLoad(new Load("F1", -2.0, 1.5, 0, 0));
+		myBeam.addLoad(new Load("F1", -2.0, 2, 0, 0));
 		myBeam.addLoad(new Load("q1", -5, 0, 0,4));
 
 
@@ -89,6 +89,12 @@ public static void main(String[] args) {
 			for (StressResultant r:dis) 
 					System.out.println("x=" + r.getX_m() + " m   Q=" + r.getShearingForce()+" "+r.getUnit());
 
+			List <StressResultant> zero=new ArrayList<>();
+			zero=qTable.getZeroPoints();
+			System.out.println("Zero points in Q");
+			for (StressResultant r:zero) 
+					System.out.println("x=" + r.getX_m() + " m");
+			
 			StressResultantDraw d=new StressResultantDraw("Q",myBeam,qTable,600,1200,10,10,"%.2f");
 			d.draw();
 
@@ -110,6 +116,12 @@ public static void main(String[] args) {
 			for (StressResultant r:dis) 
 					System.out.println("x=" + r.getX_m() + " m   M=" + r.getShearingForce()+" "+r.getUnit());
 
+		    zero=new ArrayList<>();
+			zero=mTable.getZeroPoints();
+			System.out.println("Zero points in M");
+			for (StressResultant r:zero) 
+					System.out.println("x=" + r.getX_m() + " m");
+			
 			StressResultantDraw m=new StressResultantDraw("M",myBeam,mTable,600,1200,10,10,"%.2f");
 			m.draw();
 

Readme.md

@@ -44,26 +44,21 @@ public static StressResultantTable solve(StressResultantTable qTable, Beam beam,
 		double sectionLength_m = mTable.getSectionLength_m();
 		double x = 0;
 
+		mTable.getShearingForceAtIndex(0).setMaxima(true);
+		
 		for (int n = 0; n <= qTable.getLength() - 2; n++) {
 
 			q_N = qTable.getShearingForceAtIndex(n).getShearingForce();
 
 			q1_N = qTable.getShearingForceAtIndex(n + 1).getShearingForce();
 			double deltaQ_N = q1_N - q_N;
 
-			// Checks if leading sign changes or if there is an discontiunuity
-			// if so, disregard....
-			// if (Math.signum(q_N) == Math.signum(q1_N) && Math.abs(deltaQ_N)
-			// <= DISCONTIUNUITY_THRESHOLD) {
 			m_Nm = q_N * x;
 			m1_Nm = q1_N * (x + sectionLength_m);
 			deltaM_Nm = m1_Nm - m_Nm;
-			// } else {
-
-			// }
 
 			m_Nm = mTable.getShearingForceAtIndex(n).getShearingForce();
-			mTable.getShearingForceAtIndex(n+1).setShearingForce(m_Nm + deltaM_Nm);
+			mTable.getShearingForceAtIndex(n + 1).setShearingForce(m_Nm + deltaM_Nm);
 
 			// Check for zero points in Q(x). Zero points are local
 			// maxima in M(x).
@@ -72,18 +67,20 @@ public static StressResultantTable solve(StressResultantTable qTable, Beam beam,
 
 			// Check for disconuinity in Q(x) because
 			// M(x) must also be a diconuinity or a local maxima/ minima
-			if (qTable.getShearingForceAtIndex(n).isDiscontiunuity()) {
-				if (!mTable.getShearingForceAtIndex(n).isMaxima()){
-					mTable.getShearingForceAtIndex(n).setDiscontiunuity(true);
+			if (qTable.getShearingForceAtIndex(n).isDiscontiunuity() && !mTable.getShearingForceAtIndex(n).isMaxima()) {
+				mTable.getShearingForceAtIndex(n).setDiscontiunuity(true);
 				mTable.getShearingForceAtIndex(n).setShearingForceDeltaBy(m_Nm + deltaM_Nm);
-				}
 			}
+			
+			if (Math.signum(mTable.getShearingForceAtIndex(n).getShearingForce()) != Math.signum(mTable.getShearingForceAtIndex(n+1).getShearingForce()))
+				mTable.getShearingForceAtIndex(n).setZeroPoint(true);
+
 			// Next
 			x = x + sectionLength_m / (1 / sectionLength_m);
 		}
-		
-		mTable.getShearingForceAtIndex(mTable.getLength()-1).setShearingForce(0);
-		mTable.getShearingForceAtIndex(mTable.getLength()-1).setZeroPoint(true);
+
+		mTable.getShearingForceAtIndex(mTable.getLength() - 1).setShearingForce(0);
+		mTable.getShearingForceAtIndex(mTable.getLength() - 1).setZeroPoint(true);
 		return mTable;
 	}
 }

Readme.md~

@@ -42,7 +42,7 @@ public class QSolver {
 	 * 
 	 */
 	public static StressResultantTable solve(Beam beam, String unit) {
-		double sectionLength_m = .0001; // Small values lead to more precise results.
+		double sectionLength_m = .001; // Small values lead to more precise results.
 
 		BeamResult result = BeamSolver.getResults(beam, "2f");
 

Shot_0.png

@@ -168,35 +168,11 @@ public void draw() {
 
 				String shFormated;
 
-				/*
-				if (r.isDiscontiunuity() && r.isMaxima()) {
-						graphics.setColor(Color.GRAY);
-						graphics.drawLine((int) getXT(x), padY_px, (int) getXT(x),
-								height_px - padY_px + PADDING_TOP_PX);
-						graphics.setColor(Color.RED);
-						shFormated = String.format(numberFormat, r.getShearingForce());
-						graphics.drawString(shFormated + " " + r.getUnit(), (int) getXT(x), (int) getYT(y));
-				} else {
-					if (r.isMaxima()) {
-						graphics.setColor(Color.BLUE);
-						graphics.drawLine((int) getXT(x), padY_px, (int) getXT(x), height_px - padY_px + PADDING_TOP_PX);
-						shFormated = String.format(numberFormat, r.getShearingForce());
-						graphics.drawString(shFormated + " " + r.getUnit(), (int) getXT(x), (int) getYT(y));
-					}
-				}
-				*/
 				if (r.isDiscontiunuity()|| r.isMaxima() || r.isZeroPoint()) {
-					/*
-					graphics.setColor(Color.GRAY);
-					graphics.drawLine((int) getXT(x), padY_px, (int) getXT(x),
-							height_px - padY_px + PADDING_TOP_PX);
-							*/
 					graphics.setColor(Color.RED);
 					shFormated = String.format(numberFormat, r.getShearingForce());
 					graphics.drawString(shFormated + " " + r.getUnit(), (int) getXT(x), (int) getYT(y));
 			}
-
-				
 				yLast = y;
 				xLast = x;
 			}