def test_panel_fkG_num():
for model in ['plate_clt_donnell_bardell', 'cpanel_clt_donnell_bardell']:
print('Checking fkG_num for model {0}'.format(model))
# ssss
p = Panel()
p.a = 8.
p.b = 4.
p.r = 1.e8
p.stack = [0, 90, 90, 0, -45, +45]
p.plyt = 1e-3 * 0.125
p.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9)
p.model = model
p.Nxx = -1.
p.m = 8
p.n = 9
p.nx = 9
p.ny = 9
p.u1ty = 1
p.u2ty = 1
p.u2tx = 1
p.v1tx = 0
p.v2tx = 1
p.v1ty = 1
p.v2ty = 1
num = 1000.
f = 0.
for i in range(int(num)):
if i == 0 or i == num - 2:
fx = p.Nxx * p.b / (num - 1.) / 2.
else:
fx = p.Nxx * p.b / (num - 1.)
p.add_force(x=p.a,
y=i * p.b / (num - 1.),
fx=fx,
fy=0.,
fz=0.,
cte=True)
f += fx
p.static(silent=True)
p.lb(silent=True)
assert np.isclose(p.eigvals[0], 4.5290911349518277, atol=0.01, rtol=0)
p.lb(silent=True, c=p.analysis.cs[0])
assert np.isclose(p.eigvals[0], 4.5345057669315239, atol=0.01, rtol=0)
def test_panel_fkG_num():
for model in ['plate_clt_donnell_bardell',
'cpanel_clt_donnell_bardell']:
print('Checking fkG_num for model {0}'.format(model))
# ssss
p = Panel()
p.a = 8.
p.b = 4.
p.r = 1.e8
p.stack = [0, 90, 90, 0, -45, +45]
p.plyt = 1e-3*0.125
p.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9)
p.model = model
p.Nxx = -1.
p.m = 8
p.n = 9
p.nx = 9
p.ny = 9
p.u1ty = 1
p.u2ty = 1
p.u2tx = 1
p.v1tx = 0
p.v2tx = 1
p.v1ty = 1
p.v2ty = 1
num = 1000.
f = 0.
for i in range(int(num)):
if i == 0 or i == num - 2:
fx = p.Nxx*p.b/(num-1.)/2.
else:
fx = p.Nxx*p.b/(num-1.)
p.add_force(x=p.a, y=i*p.b/(num-1.), fx=fx, fy=0., fz=0.,
cte=True)
f += fx
p.static(silent=True)
p.lb(silent=True)
assert np.isclose(p.eigvals[0], 4.5290911349518277, atol=0.01, rtol=0)
p.lb(silent=True, c=p.analysis.cs[0])
assert np.isclose(p.eigvals[0], 4.5345057669315239, atol=0.01, rtol=0)
def test_lb():
for model in ['plate_clt_donnell_bardell',
'plate_clt_donnell_bardell_w',
'cpanel_clt_donnell_bardell',
'kpanel_clt_donnell_bardell']:
print('Linear buckling for model {0}'.format(model))
# ssss
p = Panel()
p.bc_ssss()
p.m = 12
p.n = 13
p.stack = [0, 90, -45, +45]
p.plyt = 0.125e-3
p.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9)
p.model = model
p.a = 1.
p.b = 0.5
p.r = 1.e8
p.alphadeg = 0.
p.Nxx = -1
p.lb(silent=True)
if '_w' in model:
assert np.allclose(p.eigvals[0], 88.4769618837)
else:
assert np.allclose(p.eigvals[0], 85.2911727144)
p.Nxx = 0
p.Nyy = -1
p.lb(silent=True)
if '_w' in model:
assert np.allclose(p.eigvals[0], 26.4588171556)
else:
assert np.allclose(p.eigvals[0], 25.1756170679)
# ssfs
p = Panel()
p.bc_ssfs()
p.m = 12
p.n = 13
p.stack = [0, 90, -45, +45]
p.plyt = 0.125e-3
p.laminaprop = (142.5e9, 8.7e9, 0.28, 5.1e9, 5.1e9, 5.1e9)
p.model = model
p.a = 1.
p.b = 0.5
p.r = 1.e8
p.alphadeg = 0.
p.Nxx = -1
p.lb(silent=True)
if '_w' in model:
assert np.allclose(p.eigvals[0], 17.1442703121)
else:
assert np.allclose(p.eigvals[0], 15.8423562314)
p.bc_sfss()
p.Nxx = 0
p.Nyy = -1
p.lb(silent=True)
if '_w' in model:
assert np.allclose(p.eigvals[0], 15.8099861083)
else:
assert np.allclose(p.eigvals[0], 13.9421987614)
