コード例 #1
0
ファイル: test_static.py プロジェクト: sovsep/compmech
def test_static():
    wmin_ref = [
            ['clpt_donnell_bc1', -0.0512249327106],
            ['clpt_donnell_bc2', -0.0500855846496],
            ['clpt_donnell_bc3', -0.0509280039584],
            ['clpt_donnell_bc4', -0.0498127720591],
            ['fsdt_donnell_bc1', -0.0516948646563],
            ['fsdt_donnell_bc2', -0.0505316013923],
            ['fsdt_donnell_bc3', -0.0513959737413],
            ['fsdt_donnell_bc4', -0.0502561654612],
            ]

    for model, wmin in wmin_ref:
        cc = ConeCyl()
        cc.model = model
        cc.m1 = 20
        cc.m2 = 10
        cc.n2 = 11
        cc.name = 'Z33'
        cc.laminaprop = (123.55e3 , 8.708e3,  0.319, 5.695e3, 5.695e3, 5.695e3)
        cc.stack = [0, 0, 19, -19, 37, -37, 45, -45, 51, -51]
        cc.plyt = 0.125
        cc.r2 = 250.
        cc.H = 510.
        cc.add_SPL(10)
        for thetadeg in np.linspace(0, 360, 300, endpoint=False):
            cc.add_force(0., thetadeg, -15., 0, 0, increment=True)
        cs = cc.static()
        cc.uvw(cs[0])
        assert np.isclose(cc.w.min(), wmin, rtol=0.01)
コード例 #2
0
ファイル: test_static.py プロジェクト: sovsep/compmech
def test_NL_static():
    wmin_ref = [
            ['clpt_donnell_bc1', -0.012689461685834305],
            ['clpt_donnell_bc2', -0.011741560192200845],
            ['clpt_donnell_bc3', -0.012634776822892537],
            ['clpt_donnell_bc4', -0.011499181513525969],
            ['fsdt_donnell_bc1', -0.012621439862441628],
            ]

    for model, wmin in wmin_ref:
        cc = ConeCyl()
        cc.model = model
        cc.m1 = 20
        cc.m2 = 5
        cc.n2 = 6
        cc.name = 'Z33'
        cc.laminaprop = (123.55e3 , 8.708e3,  0.319, 5.695e3, 5.695e3, 5.695e3)
        cc.stack = [0, 0, 19, -19, 37, -37, 45, -45, 51, -51]
        cc.plyt = 0.125
        cc.r2 = 250.
        cc.H = 510.
        cc.add_SPL(10, increment=False)
        for thetadeg in np.linspace(0, 360, 300, endpoint=False):
            cc.add_force(0., thetadeg, -15., 0, 0, increment=True)
        cc.analysis.initialInc = 0.5
        cs = cc.static(NLgeom=True)
        cc.uvw(cs[0])
        assert np.isclose(cc.w.min(), wmin, rtol=0.01)
コード例 #3
0
ファイル: dev_utils.py プロジェクト: sovsep/compmech
def test1():
    cc = ConeCyl()
    cc.laminaprop = (123.55e3, 8.708e3, 0.319, 5.695e3, 5.695e3, 5.695e3)
    cc.stack = [0, 0, 19, -19, 37, -37, 45, -45, 51, -51]
    cc.plyt = 0.125
    cc.r2 = 250.
    cc.H = 510.
    cc.alphadeg = 0.

    cc.m1 = 120
    cc.m2 = 15
    cc.n2 = 15

    cc.thetaT = None
    cc.Fc = 1.
    cc.pdC = False

    cc.add_SPL(10.)

    cc.model = 'clpt_donnell_bc1'

    cc.lb()
    cc.Fc = cc.eigvals[0]
    cc.plot(cc.eigvecs[:, 0])
    cc.Fc = 100.
    regions_conditions(cc, 'k0L')
コード例 #4
0
ファイル: dev_utils.py プロジェクト: saullocastro/compmech
def test1():
    cc = ConeCyl()
    cc.laminaprop = (123.55e3 , 8.708e3,  0.319, 5.695e3, 5.695e3, 5.695e3)
    cc.stack = [0, 0, 19, -19, 37, -37, 45, -45, 51, -51]
    cc.plyt = 0.125
    cc.r2 = 250.
    cc.H = 510.
    cc.alphadeg = 0.

    cc.m1 = 120
    cc.m2 = 15
    cc.n2 = 15

    cc.thetaT = None
    cc.Fc = 1.
    cc.pdC = False

    cc.add_SPL(10.)

    cc.model = 'clpt_donnell_bc1'

    cc.lb()
    cc.Fc = cc.eigvals[0]
    cc.plot(cc.eigvecs[:,0])
    cc.Fc = 100.
    regions_conditions(cc, 'k0L')
コード例 #5
0
from compmech.conecyl import ConeCyl

cc = ConeCyl()
cc.model = 'clpt_donnell_bc1'
cc.m1 = 15
cc.m2 = 16
cc.n2 = 17
cc.name = 'Z33'
cc.laminaprop = (123.55e3 , 8.708e3,  0.319, 5.695e3, 5.695e3, 5.695e3)
cc.stack = [0, 0, 19, -19, 37, -37, 45, -45, 51, -51]
cc.plyt = 0.125
cc.r2 = 250.
cc.H = 510.
cc.add_SPL(10)
cc.static()
cc.plot(cc.cs[0], colorbar=True, filename='test.png')