solver.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	import numpy as np
	from diapason_engine import DiapasonEngine, DiapasonEngineDamped

	verb = 2

	dampingEta = 0 * 1e4 / 2. / np.pi

	k0 = 1.8e3 # [Hz]
	ktar = 1.8e3 # [Hz]

	theta = 20. * np.pi / 180.
	phi = 10. * np.pi / 180.
	c = 3.e2
	rho = 8e3 * (2. * np.pi) ** 2.
	E = 1.93e11
	nu = .3
	T = 1e6

	###
	if np.isclose(dampingEta, 0.):
	solver = DiapasonEngine(kappa = k0, c = c, rho = rho, E = E, nu = nu,
	T = T, theta = theta, phi = phi, meshNo = 1,
	degree_threshold = 8, verbosity = 0)
	else:
	solver = DiapasonEngineDamped(kappa = k0, c = c, rho = rho, E = E, nu = nu,
	T = T, theta = theta, phi = phi,
	dampingEta = dampingEta, meshNo = 1,
	degree_threshold = 8, verbosity = 0)
	uh = solver.solve(ktar)[0]
	solver.outParaviewTimeDomain(uh, omega = 2. * np.pi * ktar,
	filename = "out/outT{}_{}_".format(ktar, dampingEta),
	forceNewFile = False)