Commit 323faba3 authored by Jeremy BLEYER's avatar Jeremy BLEYER

Tangent operator blocks

parent 3796acbd
......@@ -33,13 +33,19 @@ ds = Measure("ds", subdomain_data=facets)
material = mf.MFrontNonlinearMaterial("../materials/src/libBehaviour.so",
"StationaryHeatTransfer",
"HeatTransfer",
hypothesis="plane_strain")
problem = mf.MFrontNonlinearProblem(T, material, quadrature_degree=0)
problem.bc = bc
flux = mf.ThermalFlux(T)
problem.define_form(flux)
import mgis.behaviour as mgis_bv
for t in problem.material.behaviour.tangent_operator_blocks:
print('d{}_d{} size {}x{}'.format(t[0].name,t[1].name,
mgis_bv.getVariableSize(t[0], problem.material.hypothesis),
mgis_bv.getVariableSize(t[0], problem.material.hypothesis)))
T.interpolate(Constant(Tl))
problem.solve(T.vector())
......
import mgis.behaviour as mgis_bv
mgis_hypothesis = {"plane_strain": mgis_bv.Hypothesis.PlaneStrain,
"plane_stress": mgis_bv.Hypothesis.PlaneStress,
"3d": mgis_bv.Hypothesis.Tridimensional,
"axisymmetric": mgis_bv.Hypothesis.Axisymmetrical}
"plane_stress": mgis_bv.Hypothesis.PlaneStress,
"3d": mgis_bv.Hypothesis.Tridimensional,
"axisymmetric": mgis_bv.Hypothesis.Axisymmetrical}
class MFrontNonlinearMaterial:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment