def test_Lee_and_Lee_table4(): print('Testing Lee and Lee Table 4') # Lee and Lee. "Vibration analysis of anisotropic plates with eccentric # stiffeners". Computers & Structures, Vol. 57, No. 1, pp. 99-105, # 1995. models = (('model4', 0.00208, 0.0060, 138.99917796302756), ('model5', 0.00260, 0.0075, 175.00597239286196), ('model7', 0.00364, 0.0105, 205.433509024)) for model, hf, bf, value in models: spb = StiffPanelBay() spb.model = 'plate_clt_donnell_bardell' spb.mu = 1.500e3 # plate material density in kg / m^3 spb.laminaprop = (128.e9, 11.e9, 0.25, 4.48e9, 1.53e9, 1.53e9) spb.stack = [0, -45, +45, 90, 90, +45, -45, 0] plyt = 0.00013 spb.plyt = plyt spb.a = 0.5 spb.b = 0.250 spb.m = 14 spb.n = 15 hf = hf bf = bf n = int(hf / plyt) fstack = [0] * (n // 4) + [90] * (n // 4) + [90] * (n // 4) + [0] * (n // 4) # clamping spb.w1rx = 0. spb.w2rx = 0. spb.w1ry = 0. spb.w2ry = 0. spb.add_panel(y1=0, y2=spb.b / 2.) spb.add_panel(y1=spb.b / 2., y2=spb.b) spb.add_bladestiff1d(mu=spb.mu, ys=spb.b / 2., bb=0., bf=bf, fstack=fstack, fplyt=plyt, flaminaprop=spb.laminaprop) k0 = spb.calc_k0(silent=True) M = spb.calc_kM(silent=True) eigvals, eigvecs = freq(k0, M, silent=True) herz = eigvals[0].real / 2 / np.pi assert np.isclose(herz, value, atol=0.001, rtol=0.001)
def test_Lee_and_Lee_table4(): print('Testing Lee and Lee Table 4') # Lee and Lee. "Vibration analysis of anisotropic plates with eccentric # stiffeners". Computers & Structures, Vol. 57, No. 1, pp. 99-105, # 1995. models = ( ('model4', 0.00208, 0.0060, 138.99917796302756), ('model5', 0.00260, 0.0075, 175.00597239286196), ('model7', 0.00364, 0.0105, 205.433509024)) for model, hf, bf, value in models: spb = StiffPanelBay() spb.model = 'plate_clt_donnell_bardell' spb.mu = 1.500e3 # plate material density in kg / m^3 spb.laminaprop = (128.e9, 11.e9, 0.25, 4.48e9, 1.53e9, 1.53e9) spb.stack = [0, -45, +45, 90, 90, +45, -45, 0] plyt = 0.00013 spb.plyt = plyt spb.a = 0.5 spb.b = 0.250 spb.m = 14 spb.n = 15 hf = hf bf = bf n = int(hf/plyt) fstack = [0]*(n//4) + [90]*(n//4) + [90]*(n//4) + [0]*(n//4) # clamping spb.w1rx = 0. spb.w2rx = 0. spb.w1ry = 0. spb.w2ry = 0. spb.add_panel(y1=0, y2=spb.b/2.) spb.add_panel(y1=spb.b/2., y2=spb.b) spb.add_bladestiff1d(mu=spb.mu, ys=spb.b/2., bb=0., bf=bf, fstack=fstack, fplyt=plyt, flaminaprop=spb.laminaprop) k0 = spb.calc_k0(silent=True) M = spb.calc_kM(silent=True) eigvals, eigvecs = freq(k0, M, silent=True) herz = eigvals[0].real/2/np.pi assert np.isclose(herz, value, atol=0.001, rtol=0.001)
def test_freq_models(): print('Testing frequency analysis for StiffPanelBay with 2 plates') # From Table 4 of # Lee and Lee. "Vibration analysis of anisotropic plates with eccentric # stiffeners". Computers & Structures, Vol. 57, No. 1, pp. 99-105, # 1995. for model in [ 'plate_clt_donnell_bardell', 'cpanel_clt_donnell_bardell', 'kpanel_clt_donnell_bardell' ]: spb = StiffPanelBay() spb.a = 0.5 spb.b = 0.250 spb.plyt = 0.00013 spb.laminaprop = (128.e9, 11.e9, 0.25, 4.48e9, 1.53e9, 1.53e9) spb.stack = [0, -45, +45, 90, 90, +45, -45, 0] spb.model = model spb.r = 1.e6 spb.alphadeg = 0. spb.mu = 1.5e3 spb.m = 9 spb.n = 10 # clamping spb.w1rx = 0. spb.w2rx = 0. spb.w1ry = 0. spb.w2ry = 0. spb.add_panel(0, spb.b / 2., plyt=spb.plyt) spb.add_panel(spb.b / 2., spb.b, plyt=spb.plyt) k0 = spb.calc_k0(silent=True) M = spb.calc_kM(silent=True) eigvals, eigvecs = freq(k0, M, silent=True) ref = [ 85.12907802 - 0.j, 134.16422850 - 0.j, 206.77295186 - 0.j, 216.45992453 - 0.j, 252.24546171 - 0.j ] assert np.allclose(eigvals[:5] / 2 / np.pi, ref, atol=0.1, rtol=0)
def test_freq_models(): print('Testing frequency analysis for StiffPanelBay with 2 plates') # From Table 4 of # Lee and Lee. "Vibration analysis of anisotropic plates with eccentric # stiffeners". Computers & Structures, Vol. 57, No. 1, pp. 99-105, # 1995. for model in ['plate_clt_donnell_bardell', 'cpanel_clt_donnell_bardell', 'kpanel_clt_donnell_bardell']: spb = StiffPanelBay() spb.a = 0.5 spb.b = 0.250 spb.plyt = 0.00013 spb.laminaprop = (128.e9, 11.e9, 0.25, 4.48e9, 1.53e9, 1.53e9) spb.stack = [0, -45, +45, 90, 90, +45, -45, 0] spb.model = model spb.r = 1.e6 spb.alphadeg = 0. spb.mu = 1.5e3 spb.m = 9 spb.n = 10 # clamping spb.w1rx = 0. spb.w2rx = 0. spb.w1ry = 0. spb.w2ry = 0. spb.add_panel(0, spb.b/2., plyt=spb.plyt) spb.add_panel(spb.b/2., spb.b, plyt=spb.plyt) k0 = spb.calc_k0(silent=True) M = spb.calc_kM(silent=True) eigvals, eigvecs = freq(k0, M, silent=True) ref = [85.12907802-0.j, 134.16422850-0.j, 206.77295186-0.j, 216.45992453-0.j, 252.24546171-0.j] assert np.allclose(eigvals[:5]/2/np.pi, ref, atol=0.1, rtol=0)
def test_Lee_and_Lee_table4(): print('Testing Lee and Lee Table 4') # Lee and Lee. "Vibration analysis of anisotropic plates with eccentric # stiffeners". Computers & Structures, Vol. 57, No. 1, pp. 99-105, # 1995. models = ( ('model4', 0.00208, 0.0060, 138.801067988), ('model5', 0.00260, 0.0075, 174.624343202), ('model7', 0.00364, 0.0105, 205.433509024)) for model, hf, bf, value in models: spb = StiffPanelBay() spb.model = 'plate_clt_donnell_bardell' spb.mu = 1.500e3 # plate material density in kg / m^3 spb.laminaprop = (128.e9, 11.e9, 0.25, 4.48e9, 1.53e9, 1.53e9) spb.stack = [0, -45, +45, 90, 90, +45, -45, 0] plyt = 0.00013 spb.plyt = plyt spb.a = 0.5 spb.b = 0.250 spb.m = 14 spb.n = 15 hf = hf bf = bf n = int(hf/plyt) fstack = [0]*(n//4) + [90]*(n//4) + [90]*(n//4) + [0]*(n//4) # clamping spb.w1rx = 0. spb.w2rx = 0. spb.w1ry = 0. spb.w2ry = 0. spb.add_panel(y1=0, y2=spb.b/2.) spb.add_panel(y1=spb.b/2., y2=spb.b) spb.add_bladestiff1d(mu=spb.mu, ys=spb.b/2., bb=0., bf=bf, fstack=fstack, fplyt=plyt, flaminaprop=spb.laminaprop) spb.freq(atype=4, silent=True, reduced_dof=False) assert np.isclose(spb.eigvals[0].real/2/np.pi, value)
spb.n = 22 spb.u1tx = 0. spb.u1rx = 1. spb.u2tx = 0. spb.u2rx = 1. spb.u1ty = 0. spb.u1ry = 1. spb.u2ty = 0. spb.u2ry = 1. spb.v1tx = 0. spb.v1rx = 1. spb.v2tx = 0. spb.v2rx = 1. spb.v1ty = 0. spb.v1ry = 1. spb.v2ty = 0. spb.v2ry = 1. spb.w1tx = 0. spb.w1rx = 1. spb.w2tx = 0. spb.w2rx = 1. spb.w1ty = 0. spb.w1ry = 1. spb.w2ty = 0. spb.w2ry = 1. spb.add_panel(y1=0, y2=spb.b, Nxx=-1.) spb.lb(silent=False)
spb.n = 16 spb.u1tx = 0. spb.u1rx = 1. spb.u2tx = 0. spb.u2rx = 1. spb.u1ty = 0. spb.u1ry = 1. spb.u2ty = 0. spb.u2ry = 1. spb.v1tx = 0. spb.v1rx = 1. spb.v2tx = 0. spb.v2rx = 1. spb.v1ty = 0. spb.v1ry = 1. spb.v2ty = 0. spb.v2ry = 1. spb.w1tx = 0. spb.w1rx = 1. spb.w2tx = 0. spb.w2rx = 1. spb.w1ty = 0. spb.w1ry = 1. spb.w2ty = 0. spb.w2ry = 1. spb.add_panel(y1=0, y2=spb.b, Nxx=-1.) spb.lb(silent=False)