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 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')
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))
cc.add_SPL(0.6, pt=0.25, theta=deg2rad(30))
