add test - stack

Konstantin8105 · Jan 25, 2025 · 4cf5579 · 4cf5579
1 parent c5d8f32
commit 4cf5579
Show file tree

Hide file tree

Showing 2 changed files with 96 additions and 0 deletions.
diff --git a/wind/shape_test.go b/wind/shape_test.go
@@ -43,4 +43,13 @@ func Test(t *testing.T) {
 		}
 		compare.Test(t, "test.cylinder", buf.Bytes())
 	})
+	t.Run("stack", func(t *testing.T) {
+		var buf bytes.Buffer
+		Wsum := Cylinder(&buf, ZoneA, RegionII, LogDecriment15, 0.200, 4.710, 100.100, 2.800,
+			[]float64{0.5, 1.5})
+		for _, z := range []float64{2.8, 5, 10, 20, 50, 100} {
+			fmt.Fprintf(&buf, "Wsum[z = %6.3f m] = %6.1f Pa\n", z, Wsum(z))
+		}
+		compare.Test(t, "test.stack", buf.Bytes())
+	})
 }
diff --git a/wind/test.stack b/wind/test.stack
@@ -0,0 +1,87 @@
+Sketch:
+
+          |<--- d ----->|
+          |             |
+          ***************-------
+          *             *     |
+  Wind    *             *     |
+ ----->   *             *     |
+          *             *     |
+          *             *     |
+          *             *     |
+          *             *     h
+          ***************---  |
+                          |   |
+                          zo  |
+                          |   |
+   ---------- ground -------------
+
+Wind zone: A
+Wind region:  II with value = 300.0 Pa
+Wind log decrement: δ = 0.15
+Natural frequency : [0.5 1.5]
+
+Dimensions:
+b   4.710 m
+d   4.710 m
+zo  2.800 m
+h  100.100 m
+
+Re = 116.677*10^5 for ze=0.8*h = 80.640
+
+Cx∞ =  1.100
+
+Elongation:
+λ  20.658
+λe 41.316
+ϕ   1.000
+Kλ  0.854
+
+Cx  =  0.939
+
+The spatial correlation coefficient of pressure pulsations:
+ρ  4.710
+χ 100.100
+ν  0.712
+
+| z       ze      Kz     ζ      ξ      | Wm     Wp     Wsum   |
+| m       m                            | Pa     Pa     Pa     |
+|                                      |                      |
+|  2.800   2.800   0.750  0.850  2.576 |  211.4  329.3  540.7 |
+|  5.000   5.000   0.750  0.850  2.576 |  211.4  329.3  540.7 |
+| 10.000  10.000   1.000  0.760  2.576 |  281.8  392.6  674.4 |
+| 15.000  15.000   1.129  0.715  2.576 |  318.3  417.2  735.5 |
+| 20.000  20.000   1.231  0.685  2.576 |  347.0  435.6  782.6 |
+| 25.000  25.000   1.316  0.662  2.576 |  371.0  450.4  821.4 |
+| 30.000  30.000   1.390  0.645  2.576 |  391.9  462.9  854.8 |
+| 35.000  35.000   1.456  0.630  2.576 |  410.4  473.7  884.1 |
+| 40.000  40.000   1.516  0.617  2.576 |  427.2  483.3  910.5 |
+| 45.000  45.000   1.570  0.607  2.576 |  442.5  491.9  934.5 |
+| 50.000  50.000   1.621  0.597  2.576 |  456.7  499.8  956.5 |
+| 55.000  55.000   1.668  0.589  2.576 |  470.0  507.0  977.0 |
+| 60.000  60.000   1.712  0.581  2.576 |  482.4  513.6  996.1 |
+| 65.000  65.000   1.753  0.574  2.576 |  494.2  519.8 1014.0 |
+| 70.000  70.000   1.793  0.568  2.576 |  505.3  525.7 1030.9 |
+| 75.000  75.000   1.830  0.562  2.576 |  515.8  531.1 1047.0 |
+| 80.000  80.000   1.866  0.556  2.576 |  525.9  536.3 1062.2 |
+| 85.000  85.000   1.900  0.551  2.576 |  535.6  541.2 1076.8 |
+| 90.000  90.000   1.933  0.547  2.576 |  544.8  545.9 1090.7 |
+| 95.000  95.000   1.965  0.542  2.576 |  553.7  550.3 1104.0 |
+| 100.000 100.000  1.995  0.538  2.576 |  562.3  554.5 1116.9 |
+| 100.100 100.100  1.996  0.538  2.576 |  562.5  554.6 1117.1 |
+
+     Ws on top    |----------->             |--------->
+                  |          /              |         |
+                  |--------->    Ws average |--------->
+                  |        /                |         |
+     Ws on zero   |------->                 |--------->
+                --------------- ground ------------------
+
+| side   | width  | Center of Ws | Ws average |
+|  front |  4.710 | 55.779       |   +927.1   |
+Wsum[z =  2.800 m] =  540.7 Pa
+Wsum[z =  5.000 m] =  540.7 Pa
+Wsum[z = 10.000 m] =  674.4 Pa
+Wsum[z = 20.000 m] =  782.6 Pa
+Wsum[z = 50.000 m] =  956.5 Pa
+Wsum[z = 100.000 m] = 1116.9 Pa