Esempio n. 1
0
def test_dynamic_with_TStiff2D():
    print('Testing dynamic analysis with TStiff2D')
    spb = StiffPanelBay()
    spb.a = 2.
    spb.b = 1.
    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.m = 12
    spb.n = 13
    spb.mu = 1.3e3

    spb.add_panel(y1=0, y2=spb.b / 2., plyt=spb.plyt)
    spb.add_panel(y1=spb.b / 2., y2=spb.b, plyt=spb.plyt)

    bb = spb.b / 5.
    bf = bb
    stiff = spb.add_tstiff2d(ys=spb.b / 2.,
                             bf=bf,
                             bb=bb,
                             fstack=[0, 90, 90, 0] * 2,
                             fplyt=spb.plyt * 1.,
                             flaminaprop=spb.laminaprop,
                             bstack=[0, 90, 90, 0] * 1,
                             bplyt=spb.plyt * 1.,
                             blaminaprop=spb.laminaprop,
                             mb=12,
                             nb=13,
                             mf=12,
                             nf=13)

    k0 = spb.calc_k0()
    M = spb.calc_kM()
    eigvals, eigvecs = freq(k0, M, silent=True)
Esempio n. 2
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)
Esempio n. 3
0
def test_static_with_TStiff2D():
    print('Testing static analysis with TStiff2D')
    spb = StiffPanelBay()
    spb.a = 2.
    spb.b = 1.
    spb.stack = [0, +45, -45, 90, -45, +45, 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.m = 12
    spb.n = 13

    spb.add_panel(y1=0, y2=spb.b/2., plyt=spb.plyt)
    spb.add_panel(y1=spb.b/2., y2=spb.b, plyt=spb.plyt)

    bb = spb.b/5.
    bf = bb
    stiff = spb.add_tstiff2d(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,
                     mb=14, nb=13, mf=13, nf=12)

    stiff.flange.forces.append([spb.a/2., bf, 0., 0., 1000.])

    k0 = spb.calc_k0()
    fext = spb.calc_fext()
    inc, cs = static(spb.k0, fext, silent=True)

    wpanelmin = spb.uvw_skin(cs[0])[2].min()
    wbasemin = spb.uvw_stiffener(cs[0], 0, region='base')[2].min()
    wflangemax = spb.uvw_stiffener(cs[0], 0, region='flange')[2].max()
    spb.plot_skin(cs[0], filename='tmp_test_tstiff2d_static_skin.png',
            colorbar=True, vec='w', clean=False)
    spb.plot_stiffener(cs[0], si=0, region='base',
            filename='tmp_test_tstiff2d_static_stiff_base.png', colorbar=True,
            vec='w', clean=False)
    spb.plot_stiffener(cs[0], si=0, region='flange',
            filename='tmp_test_tstiff2d_static_stiff_flange.png', colorbar=True,
            clean=False)
Esempio n. 4
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)