From af1d3a937625ad984505de124e90ad7a341e659e Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Maur=C3=ADcio=20Bonatte?= Date: Sat, 26 Aug 2023 17:37:44 +0100 Subject: [PATCH] Adding tests from Flavio Neto --- secan/test/test.py | 82 ++++++++++++++++++++++++++++++++++++---------- 1 file changed, 65 insertions(+), 17 deletions(-) diff --git a/secan/test/test.py b/secan/test/test.py index affb3fb..e4906aa 100644 --- a/secan/test/test.py +++ b/secan/test/test.py @@ -199,7 +199,7 @@ def setUp(self): self.conc = Concrete(0.85*30e6/1.4) self.steel = SteelIdeal(210e9, 500e6/1.15) - self.r = RectSection(0.2, 0.32, self.conc, (0, 0), n_discret=1000) + self.r = RectSection(0.2, 0.32, self.conc, (0, 0), n_discret=300) self.s = Section([self.r]) xy_pairs = [(x,y) for x in [-.07,.07] for y in [-0.11, 0, 0.11]] @@ -223,6 +223,14 @@ def get_rec_normal(self, e0,k): def get_rec_moment(self, e0,k): return np.array([s.get_moment_resistance(e0=e0, k=k, center=self.s.centroid[1]) for s in self.s.section if isinstance(s, RectSection) ]).sum() + + def get_rebar_stiff(self, e0, k): + return np.array([s.get_stiffness(e0=e0, k=k, center=self.s.centroid[1]) + for s in self.s.section if isinstance(s, Rebar) ]).sum(axis=0) + + def get_rec_stiff(self, e0, k): + return np.array([s.get_stiffness(e0=e0, k=k, center=self.s.centroid[1]) + for s in self.s.section if isinstance(s, RectSection) ]).sum(axis=0) def test_example_3_4(self): # 3.4.1 @@ -273,35 +281,75 @@ def test_example_3_4(self): self.assertAlmostEqual(self.s.get_moment_res(e0=-0.5e-3, k=9.3750e-3)/1e6, 0.0489, delta=self.delta) def test_example_4_5(self): - self.s.residual_e0 = 1e-8 - self.s.max_increment_e0 = 0.0005 - self.s.n_ite_e0 = 3000 self.s.tolerance_check_section = 1e-8 - self.assertAlmostEqual(self.s.check_section(target_normal=-0.8160e6, target_moment=0.0373e6, n_ite=1000), (-0.9027403850211792e-3, 5.985830149126867e-3)) - + # 4.5.1 + rebar_stiff = self.get_rebar_stiff(0, 0) + self.assertAlmostEqual(rebar_stiff[0][0]/1e7, 210*4.71239e-2, delta=self.delta) + self.assertEqual(rebar_stiff[0][1], 0,) + self.assertAlmostEqual(rebar_stiff[1][1]/1e5, 7.9828, delta=self.delta) + + rec_stiff = self.get_rec_stiff(0, 0) + self.assertAlmostEqual(rec_stiff[0][0]/1e9, 1.1657, delta=self.delta) + self.assertAlmostEqual(rec_stiff[0][1], 0, delta=self.delta) + self.assertAlmostEqual(rec_stiff[1][1]/1e6, 9.9474, delta=self.delta) + + stiff = self.s.get_stiff(0, 0) + self.assertAlmostEqual(stiff[0][0]/1e9, 1.2647, delta=self.delta) + self.assertAlmostEqual(stiff[0][1], 0, delta=self.delta) + self.assertAlmostEqual(stiff[1][1]/1e7, 1.0746, delta=self.delta) + self.assertAlmostEqual(np.linalg.det(stiff)/1e16, 1.3590, delta=self.delta) + + target_normal=-0.8160e6 + target_moment=0.0373e6 + e0, k = self.s.check_section(target_normal, target_moment, n_ite=1000) + + self.assertAlmostEqual(e0*1e3, -0.9028, delta=self.delta) + self.assertAlmostEqual(k*1e3, 5.9864, delta=0.0006) ###### + + self.assertAlmostEqual(self.s.get_normal_res(e0, k), target_normal, delta=self.delta) + self.assertAlmostEqual(self.s.get_moment_res(e0, k), target_moment, delta=self.delta) + + # 4.5.2 self.assertAlmostEqual(self.s.check_section(target_normal=0, target_moment=0.0373e6, n_ite=1000), (None, None)) + # 4.5.3 self.assertAlmostEqual(self.s.check_section(target_normal=-1.1657e6, target_moment=0.0373e6, n_ite=1000), (None, None)) - self.assertAlmostEqual(self.s.check_section(target_normal=-0.8160e6, target_moment=0, n_ite=1000), (-0.7884934594071175e-3, 9.412299199057889e-20)) + # 4.5.4 + e0, k = self.s.check_section(target_normal=-0.8160e6, target_moment=0, n_ite=1000) + self.assertAlmostEqual(e0*1e3, -0.7884934594071175, delta=self.delta) + self.assertAlmostEqual(k*1e3, 0, delta=self.delta) + # 4.5.5 self.assertAlmostEqual(self.s.check_section(target_normal=-0.8160e6, target_moment=-0.0560e6, n_ite=1000), (None, None)) - self.assertAlmostEqual(self.s.check_section(target_normal=0.1166e6, target_moment=0.0075e6, n_ite=1000), (1.3784131434898693e-3, 9.843707714720228e-3)) #This should be (1.369521e-3, 9.759756e-3) + # 4.5.6 + target_normal=0.1166e6 + target_moment=0.0075e6 + e0, k = self.s.check_section(target_normal, target_moment, n_ite=10) - e0=1.369521e-3 - k=9.759756e-3 + self.assertAlmostEqual(e0*1e3, 1.369521, delta=0.009) ###### + self.assertAlmostEqual(k*1e3, 9.759756, delta=0.09) ###### - rebar_normal = self.get_rebar_normal(e0,k) - rebar_moment = self.get_rebar_moment(e0,k) + self.assertAlmostEqual(self.s.get_normal_res(e0, k), target_normal, delta=self.delta) + self.assertAlmostEqual(self.s.get_moment_res(e0, k), target_moment, delta=self.delta) - - rec_normal = self.get_rec_normal(e0,k) - rec_moment = self.get_rec_moment(e0,k) + self.assertAlmostEqual(self.get_rebar_normal(e0, k)/1e7, 0.0123, delta=self.delta) + self.assertAlmostEqual(self.get_rebar_moment(e0, k)/1e5, 0.0644, delta=self.delta) + + self.assertAlmostEqual(self.get_rec_normal(e0, k)/1e5, 0.0666, delta=0.15) ###### + self.assertAlmostEqual(self.get_rec_moment(e0, k)/1e5, 0.0102, delta=0.0005) ###### + + rebar_stiff = self.get_rebar_stiff(e0, k) + self.assertAlmostEqual(rebar_stiff[0][0]/1e9, 0.0660, delta=self.delta) + self.assertAlmostEqual(rebar_stiff[0][1]/1e8, -0.0363, delta=self.delta) ###### + self.assertAlmostEqual(rebar_stiff[1][1]/1e7, 0.0399, delta=self.delta) - print(rebar_normal, rebar_moment) - print(rec_normal, rec_moment) + rec_stiff = self.get_rec_stiff(e0, k) + self.assertAlmostEqual(rec_stiff[0][0]/1e9, 0.0682, delta=0.003) ###### + self.assertAlmostEqual(rec_stiff[0][1]/1e9, -0.0102, delta=0.0004) ###### + self.assertAlmostEqual(rec_stiff[1][1]/1e8, 0.0154, delta=0.0004) ###### if __name__ == '__main__':