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)
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)
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')
def convergence_nx_nt(): '''It was selected a hypothetical cone with the laminate of cylinder Z33 and a semi-vertex angle of 35. The data generated from this function is used to define the dictionary `nx_nt_table` inside `ConeCyl.rebuild()`. For simplification, it is assumed `m1=m2=n2`. ''' 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 = 35. cc.linear_kinematics='clpt_donnell' cc.NL_kinematics='donnell_numerical' cc.Fc = 200000 cc.initialInc = 0.2 cc.minInc = 1.e-3 cc.modified_NR = False cc.compute_every_n = 10 cc.pd = False PLs = [5, 20, 30, 45, 60, 70, 80] PLs = [90] for n in [20, 25, 30, 35]: cc.m1 = cc.m2 = cc.n2 = n for nt in [60, 80, 100, 120]: cc.nx = cc.nt = nt print cProfile.runctx('curves = cc.SPLA(PLs, NLgeom=True)', globals(), locals(), 'tester.prof')
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')
def convergence_nx_nt(): '''It was selected a hypothetical cone with the laminate of cylinder Z33 and a semi-vertex angle of 35. The data generated from this function is used to define the dictionary `nx_nt_table` inside `ConeCyl.rebuild()`. For simplification, it is assumed `m1=m2=n2`. ''' 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 = 35. cc.linear_kinematics = 'clpt_donnell' cc.NL_kinematics = 'donnell_numerical' cc.Fc = 200000 cc.initialInc = 0.2 cc.minInc = 1.e-3 cc.modified_NR = False cc.compute_every_n = 10 cc.pd = False PLs = [5, 20, 30, 45, 60, 70, 80] PLs = [90] for n in [20, 25, 30, 35]: cc.m1 = cc.m2 = cc.n2 = n for nt in [60, 80, 100, 120]: cc.nx = cc.nt = nt print( cProfile.runctx('curves = cc.SPLA(PLs, NLgeom=True)', globals(), locals(), 'tester.prof'))
from numpy import deg2rad from compmech.conecyl import ConeCyl if __name__ == '__main__': cc = ConeCyl() 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.alphadeg = 0 cc.m1 = cc.m2 = cc.n2 = 40 cc.pdC = False cc.Fc = 2000. cc.add_SPL(0.4, pt=0.25, theta=deg2rad(-30)) cc.add_SPL(0.4, pt=0.25, theta=deg2rad(-60)) cc.add_SPL(0.4, pt=0.25, theta=deg2rad(-90)) cc.add_SPL(0.3, pt=0.25, theta=deg2rad(-120)) cc.add_SPL(0.3, pt=0.5, theta=deg2rad(-120)) cc.add_SPL(0.3, pt=0.75, theta=deg2rad(-120)) cc.add_SPL(0.4, pt=0.75, theta=deg2rad(-90)) cc.add_SPL(0.4, pt=0.75, theta=deg2rad(-60)) cc.add_SPL(0.4, pt=0.75, theta=deg2rad(-30)) cc.add_SPL(0.5, pt=0.75, theta=deg2rad(30)) cc.add_SPL(0.4, pt=0.50, theta=deg2rad(30))