def elem_dep_data_qt10ms(na, nb, nt):
fens, fes = t4_block(a, b, t, na, nb, nt, orientation='a')
fens, fes = t4_to_t10(fens, fes)
femm = FEMMDeforLinearQT10MS(material=m, material_csys=mcsys, fes=fes)
bfes = mesh_boundary(femm.fes)
fesel = fe_select(fens, bfes, box=[a, a, 0, b, 0, t], inflate=htol)
tsfes = bfes.subset(fesel)
sfemm = FEMMDeforLinear(fes=tsfes, integration_rule=TriRule(npts=3))
return fens, femm, sfemm
def elem_dep_data_qt10ms(na, nb, nts):
print(na, nb, nts)
# fens, fes = t4_block(a, b, t, na, nb, nt, orientation='a')
fens, fes = t4_composite_plate(a, b, ts, na, nb, nts)
fens, fes = t4_to_t10(fens, fes)
bfes = mesh_boundary(fes)
t = sum(ts)
fesel = fe_select(fens, bfes, box=[0, a, 0, b, t, t], inflate=htol)
tsfes = bfes.subset(fesel)
femms = []
for layer in range(len(nts)):
aangle = angles[layer] / 180. * math.pi
print(aangle)
mcsys = CSys(matrix=rotmat(numpy.array([0.0, 0.0, aangle])))
el = fe_select(fens, fes, bylabel=True, label=layer)
femms.append(
FEMMDeforLinearQT10MS(material=m,
material_csys=mcsys,
fes=fes.subset(el)))
sfemm = FEMMDeforLinear(fes=tsfes, integration_rule=TriRule(npts=3))
return fens, femms, sfemm
nu = 0.4999
W = 2.5
H = 5
L = 50
# nW, nL, nH = 20, 20, 20
nW, nL, nH = 4, 20, 4
htol = min(L, H, W) / 1000
magn = -0.2 * 12.2334 / 4
Force = magn * W * H * 2
Force * L**3 / (3 * E * W * H**3 * 2 / 12)
uzex = -12.0935378981478
m = MatDeforTriaxLinearIso(e=E, nu=nu)
start = time.time()
fens, fes = t4_block(W, L, H, nW, nL, nH, orientation='ca')
fens, fes = t4_to_t10(fens, fes)
print('Mesh generation', time.time() - start)
model_data = {}
model_data['fens'] = fens
model_data['regions'] = [{'femm': FEMMDeforLinearQT10MS(material=m, fes=fes)}]
algo_common.plot_mesh(model_data)
model_data['boundary_conditions'] = {}
# Clamped face
cn = fenode_select(fens, box=numpy.array([0, W, 0, 0, 0, H]), inflate=htol)
essential = [
{
'node_list': cn,
'comp': 0,