def test_pm45(): # {{{1
A = lamina(hf, hr, 100, 45, 0.5)
B = lamina(hf, hr, 100, -45, 0.5)
pm45 = laminate('pm45', [A, B, B, A])
assert math.isclose(pm45.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(pm45.ρ, 1.43, rel_tol=0.01)
assert pm45.vf == 0.5
assert math.isclose(pm45.wf, 0.615, rel_tol=0.01)
assert math.isclose(pm45.Ex, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.Ey, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(pm45.Gxy, 31928.0, rel_tol=0.01)
assert math.isclose(pm45.Gxz, 4608.06, rel_tol=0.01)
assert math.isclose(pm45.Gyz, 4608.06, rel_tol=0.01)
assert math.isclose(pm45.νxy, 0.756, rel_tol=0.01)
assert math.isclose(pm45.νyx, 0.756, rel_tol=0.01)
assert math.isclose(pm45.αx, 6.4269e-06, rel_tol=0.01)
assert math.isclose(pm45.αy, 6.4269e-06, rel_tol=0.01)
assert math.isclose(pm45.tEx, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.tEy, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(pm45.tGxy, 31935.1, rel_tol=0.01)
assert math.isclose(pm45.tGyz, 5529.6, rel_tol=0.01)
assert math.isclose(pm45.tGxz, 5529.6, rel_tol=0.01)
assert math.isclose(pm45.tνxy, 0.75, rel_tol=0.01)
assert math.isclose(pm45.tνxz, 0.102, rel_tol=0.01)
assert math.isclose(pm45.tνyz, 0.102, rel_tol=0.01)
def test_plain_weave(): # {{{1
A = lamina(hf, hr, 100, 0, 0.5)
B = lamina(hf, hr, 100, 90, 0.5)
pw = laminate('pw', [A, B, B, A])
assert math.isclose(pw.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(pw.ρ, 1.43, rel_tol=0.01)
assert pw.vf == 0.5
assert math.isclose(pw.wf, 0.615, rel_tol=0.01)
assert math.isclose(pw.Ex, 68327.5, rel_tol=0.01)
assert math.isclose(pw.Ey, 68327.5, rel_tol=0.01)
assert math.isclose(pw.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(pw.Gxy, 5095.588, rel_tol=0.01)
assert math.isclose(pw.Gxz, 4743.715, rel_tol=0.01)
assert math.isclose(pw.Gyz, 4472.41, rel_tol=0.01)
assert math.isclose(pw.νxy, 0.07, rel_tol=0.01)
assert math.isclose(pw.νyx, 0.07, rel_tol=0.01)
assert math.isclose(pw.αx, 6.427e-06, rel_tol=0.01)
assert math.isclose(pw.αy, 6.427e-06, rel_tol=0.01)
assert math.isclose(pw.tEx, 68335.7, rel_tol=0.01)
assert math.isclose(pw.tEy, 68335.7, rel_tol=0.01)
assert math.isclose(pw.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(pw.tGxy, 5095.588, rel_tol=0.01)
assert math.isclose(pw.tGyz, 5529.675, rel_tol=0.01)
assert math.isclose(pw.tGxz, 5529.675, rel_tol=0.01)
assert math.isclose(pw.tνxy, 0.07, rel_tol=0.01)
assert math.isclose(pw.tνxz, 0.39, rel_tol=0.01)
assert math.isclose(pw.tνyz, 0.39, rel_tol=0.01)
def test_qi(): # {{{1
A = lamina(hf, hr, 100, 0, 0.5)
B = lamina(hf, hr, 100, 90, 0.5)
C = lamina(hf, hr, 100, 45, 0.5)
D = lamina(hf, hr, 100, -45, 0.5)
qi = laminate('qi', [A, B, C, D, D, C, B, A])
assert math.isclose(qi.thickness, 0.9090, rel_tol=0.01)
assert math.isclose(qi.ρ, 1.43, rel_tol=0.01)
assert qi.vf == 0.5
assert math.isclose(qi.wf, 0.615, rel_tol=0.01)
assert math.isclose(qi.Ex, 49236.42, rel_tol=0.01)
assert math.isclose(qi.Ey, 49236.42, rel_tol=0.01)
assert math.isclose(qi.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(qi.Gxy, 18511.8, rel_tol=0.01)
assert math.isclose(qi.Gxz, 4658.93, rel_tol=0.01)
assert math.isclose(qi.Gyz, 4557.19, rel_tol=0.01)
assert math.isclose(qi.νxy, 0.3298, rel_tol=0.01)
assert math.isclose(qi.νyx, 0.3298, rel_tol=0.01)
assert math.isclose(qi.αx, 6.42695e-06, rel_tol=0.01)
assert math.isclose(qi.αy, 6.42695e-06, rel_tol=0.01)
assert math.isclose(qi.tEx, 49242.7, rel_tol=0.01)
assert math.isclose(qi.tEy, 49242.7, rel_tol=0.01)
assert math.isclose(qi.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(qi.tGxy, 18515.3, rel_tol=0.01)
assert math.isclose(qi.tGyz, 5529.6, rel_tol=0.01)
assert math.isclose(qi.tGxz, 5529.6, rel_tol=0.01)
assert math.isclose(qi.tνxy, 0.32977, rel_tol=0.01)
assert math.isclose(qi.tνxz, 0.28244, rel_tol=0.01)
assert math.isclose(qi.tνyz, 0.28244, rel_tol=0.01)
def test_extended3(): # {{{1
f = fiber(240000, 0.2, -0.2e-6, 1.76, "carbon")
r = resin(3000, 0.3, 20e-6, 1.2, "resin")
layers = [
lamina(f, r, 300, 0, 0.50),
"pw 200 45",
lamina(f, r, 100, 45, 0.40),
lamina(f, r, 100, -45, 0.40),
"symmetry",
]
extended = _extended(layers)
assert len(extended) == 5
assert isinstance(extended[0], str) and extended[0] == "pw 200 45"
assert isinstance(extended[3], str) and extended[3] == "unknown"
def test_ud(): # {{{1
la = lamina(hf, hr, 100, 0, 0.5)
ud = laminate('ud', [la, la, la, la])
assert math.isclose(ud.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(ud.ρ, 1.43, rel_tol=0.01)
assert ud.vf == 0.5
assert math.isclose(ud.wf, 0.615, rel_tol=0.01)
assert math.isclose(ud.Ex, 118810, rel_tol=0.01)
assert math.isclose(ud.Ey, 16979.8, rel_tol=0.01)
assert math.isclose(ud.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(ud.Gxy, 5095.588, rel_tol=0.01)
assert math.isclose(ud.Gxz, 4246.32, rel_tol=0.01)
assert math.isclose(ud.Gyz, 4969.80, rel_tol=0.01)
assert math.isclose(ud.νxy, 0.28, rel_tol=0.01)
assert math.isclose(ud.νyx, 0.04, rel_tol=0.01)
assert math.isclose(ud.αx, 2.75e-07, rel_tol=0.01)
assert math.isclose(ud.αy, 4.14e-5, rel_tol=0.01)
assert math.isclose(ud.tEx, 118810.0, rel_tol=0.01)
assert math.isclose(ud.tEy, 16979.80, rel_tol=0.01)
assert math.isclose(ud.tEz, 16979.80, rel_tol=0.01)
assert math.isclose(ud.tGxy, 5095.58, rel_tol=0.01)
assert math.isclose(ud.tGyz, 5963.76, rel_tol=0.01)
assert math.isclose(ud.tGxz, 5095.58, rel_tol=0.01)
assert math.isclose(ud.tνxy, 0.2799, rel_tol=0.01)
assert math.isclose(ud.tνxz, 0.2799, rel_tol=0.01)
assert math.isclose(ud.tνyz, 0.4235, rel_tol=0.01)
def test_lamina(): # {{{1
f = fiber(230000, 0.30, -0.41e-6, 1.76, 'T300')
r = resin(2900, 0.36, 41.4e-6, 1.15, 'Epikote04908')
la = lamina(f, r, 100, 0, 0.5)
areclose(
(la.E1, la.E2, la.G12, la.ν12, la.αx, la.αy, la.ρ),
(116450.0, 10980.86, 3198.53, 0.33, 1.106e-07, 4.14e-05, 1.455)
)
areclose(
(la.Q̅11, la.Q̅12, la.Q̅16, la.Q̅22, la.Q̅26, la.Q̅66),
(117658.22, 3661.28, 0.0, 11094.79, 0.0, 3198.53)
)
