def test_bladestiff2d_lb(): print('Testing linear buckling with BladeStiff2D') spb = StiffPanelBay() spb.a = 2. spb.b = 1. spb.stack = [0, +45, -45, 90, -45, +45] spb.plyt = 1e-3 * 0.125 spb.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9) spb.model = 'plate_clt_donnell_bardell' spb.m = 17 spb.n = 16 Nxx = -50. spb.add_panel(y1=0, y2=spb.b / 2., plyt=spb.plyt, Nxx=Nxx) spb.add_panel(y1=spb.b / 2., y2=spb.b, plyt=spb.plyt, Nxx=Nxx) bb = spb.b / 5. bf = bb stiff = spb.add_bladestiff2d(ys=spb.b / 2., bf=bf, bb=bb, fstack=[0, 90, 90, 0] * 8, fplyt=spb.plyt * 1., flaminaprop=spb.laminaprop, bstack=[0, 90, 90, 0] * 4, bplyt=spb.plyt * 1., blaminaprop=spb.laminaprop, mf=17, nf=16) spb.calc_k0() spb.calc_kG0() eigvals, eigvecs = lb(spb.k0, spb.kG0, silent=True) spb.plot_skin(eigvecs[:, 0], filename='tmp_test_bladestiff2d_lb_skin.png', colorbar=True) spb.plot_stiffener(eigvecs[:, 0], si=0, region='flange', filename='tmp_test_bladestiff2d_lb_stiff_flange.png', colorbar=True) calc = eigvals[0] * Nxx spb.plot_skin(eigvecs[:, 0], filename='tmp_test_bladestiff2d_lb_skin.png', colorbar=True, vec='w', clean=False) assert np.isclose(calc, -759.05689868085778, atol=0.0001, rtol=0.001)
def test_lb_Stiffener2D(): print('Testing linear buckling for StiffPanelBay with a 2D Stiffener') spb = StiffPanelBay() spb.a = 1. spb.b = 0.5 spb.stack = [0, 90, 90, 0] spb.plyt = 1e-3*0.125 spb.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9) spb.model = 'plate_clt_donnell_bardell' spb.mu = 1.3e3 spb.m = 15 spb.n = 16 spb.add_panel(y1=0, y2=spb.b/2., plyt=spb.plyt, Nxx=-1.) spb.add_panel(y1=spb.b/2., y2=spb.b, plyt=spb.plyt, Nxx_cte=1000.) spb.add_bladestiff2d(ys=spb.b/2., m1=14, n1=11, bf=0.05, fstack=[0, 90, 90, 0], fplyt=spb.plyt, flaminaprop=spb.laminaprop) k0 = spb.calc_k0(silent=True) kG = spb.calc_kG0(silent=True) eigvals, eigvecs = lb(k0, kG, silent=True) assert np.isclose(eigvals[0].real, 301.0825234, atol=0.1, rtol=0)
def test_lb_Stiffener2D(): print('Testing linear buckling for StiffPanelBay with a 2D Stiffener') spb = StiffPanelBay() spb.a = 1. spb.b = 0.5 spb.stack = [0, 90, 90, 0] spb.plyt = 1e-3 * 0.125 spb.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9) spb.model = 'plate_clt_donnell_bardell' spb.mu = 1.3e3 spb.m = 15 spb.n = 16 spb.add_panel(y1=0, y2=spb.b / 2., plyt=spb.plyt, Nxx=-1.) spb.add_panel(y1=spb.b / 2., y2=spb.b, plyt=spb.plyt, Nxx_cte=1000.) spb.add_bladestiff2d(ys=spb.b / 2., m1=14, n1=11, bf=0.05, fstack=[0, 90, 90, 0], fplyt=spb.plyt, flaminaprop=spb.laminaprop) k0 = spb.calc_k0(silent=True) kG = spb.calc_kG0(silent=True) eigvals, eigvecs = lb(k0, kG, silent=True) assert np.isclose(eigvals[0].real, 301.0825234, atol=0.1, rtol=0)
def test_tstiff2d_lb(): #NOTE reference values taken from # compmech.panel.assembly.tstiff2d_1stiff_compression a_b_list = [0.5, 2., 10.] ref_values = [-152.607275551, -80.15391295315146, -79.39137361313862] for a_b, ref_value in zip(a_b_list, ref_values): print('Testing linear buckling') spb = StiffPanelBay() spb.b = 1. spb.a = spb.b * a_b spb.stack = [0, +45, -45, 90, -45, +45] spb.plyt = 1e-3 * 0.125 spb.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9) spb.model = 'plate_clt_donnell_bardell' spb.m = 12 spb.n = 13 Nxx = -10. spb.add_panel(y1=0, y2=spb.b / 2., plyt=spb.plyt, Nxx=Nxx) spb.add_panel(y1=spb.b / 2., y2=spb.b, plyt=spb.plyt, Nxx=Nxx) bb = spb.b / 5. if False: stiff = spb.add_tstiff2d(ys=spb.b / 2., bf=bb / 2, bb=bb, fstack=[0, 90, 90, 0] * 8, fplyt=spb.plyt * 1., flaminaprop=spb.laminaprop, bstack=[0, 90, 90, 0] * 4, bplyt=spb.plyt * 1., blaminaprop=spb.laminaprop, mb=13, nb=12, mf=13, nf=12, Nxxf=Nxx) spb.calc_k0() spb.calc_kG0() eigvals, eigvecs = lb(spb.k0, spb.kG0, silent=True) calc = eigvals[0] * Nxx assert np.isclose(calc, ref_value, rtol=0.05)
def test_tstiff2d_lb(): #NOTE reference values taken from # compmech.panel.assembly.tstiff2d_1stiff_compression a_b_list = [0.5, 2., 10.] ref_values = [-152.607275551, -80.15391295315146, -79.39137361313862] for a_b, ref_value in zip(a_b_list, ref_values): print('Testing linear buckling') spb = StiffPanelBay() spb.b = 1. spb.a = spb.b * a_b spb.stack = [0, +45, -45, 90, -45, +45] spb.plyt = 1e-3*0.125 spb.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9) spb.model = 'plate_clt_donnell_bardell' spb.m = 12 spb.n = 13 Nxx = -10. spb.add_panel(y1=0, y2=spb.b/2., plyt=spb.plyt, Nxx=Nxx) spb.add_panel(y1=spb.b/2., y2=spb.b, plyt=spb.plyt, Nxx=Nxx) bb = spb.b/5. if False: stiff = spb.add_tstiff2d(ys=spb.b/2., bf=bb/2, bb=bb, fstack=[0, 90, 90, 0]*8, fplyt=spb.plyt*1., flaminaprop=spb.laminaprop, bstack=[0, 90, 90, 0]*4, bplyt=spb.plyt*1., blaminaprop=spb.laminaprop, mb=13, nb=12, mf=13, nf=12, Nxxf=Nxx) spb.calc_k0() spb.calc_kG0() eigvals, eigvecs = lb(spb.k0, spb.kG0, silent=True) calc = eigvals[0]*Nxx assert np.isclose(calc, ref_value, rtol=0.05)