Run a Job

Learn how to apply boundary conditions in a ramped manner within a Step and run a Job.

• create a Step with ramped boundary conditions

• run a Job and export a XDMF time-series file

This tutorial once again covers the essential high-level parts of creating and solving problems with FElupe. This time, however, the external displacements are applied in a ramped manner. The prescribed displacements of a cube under non-homogenous uniaxial loading will be controlled within a step. The Ogden-Roxburgh pseudo-elastic Mullins softening model is combined with an isotropic hyperelastic Neo-Hookean material formulation, which is further applied on a (nearly) incompressible solid body for a realistic analysis of rubber-like materials. Note that the bulk modulus is now an argument of the (nearly) incompressible solid body instead of the constitutive Neo-Hookean material definition.

import felupe as fem

mesh = fem.Cube(n=6)
region = fem.RegionHexahedron(mesh=mesh)
field = fem.FieldContainer([fem.Field(region=region, dim=3)])

boundaries, loadcase = fem.dof.uniaxial(field, clamped=True)

umat = fem.OgdenRoxburgh(material=fem.NeoHooke(mu=1), r=3, m=1, beta=0)
body = fem.SolidBodyNearlyIncompressible(umat=umat, field=field, bulk=5000)

The ramped prescribed displacements for 20 substeps are created with linsteps. A Step is created with a list of items to be considered (here, one single solid body) and a dict of ramped boundary conditions along with the prescribed values.

move = fem.math.linsteps([0, 2], num=10)
uniaxial = fem.Step(
    items=[body], ramp={boundaries["move"]: move}, boundaries=boundaries
)

This step is now added to a Job. The results are exported after each completed and successful substep as a time-series XDMF-file.

job = fem.Job(steps=[uniaxial])
job.evaluate(filename="result.xdmf", verbose=True)

field.plot("Principal Values of Logarithmic Strain").show()

Static Scene

Interactive Scene

  0%|          | 0/11 [00:00<?, ?substep/s]
 18%|█▊        | 2/11 [00:00<00:00, 10.48substep/s]
 36%|███▋      | 4/11 [00:00<00:00,  8.29substep/s]
 45%|████▌     | 5/11 [00:00<00:00,  7.98substep/s]
 55%|█████▍    | 6/11 [00:00<00:00,  7.74substep/s]
 64%|██████▎   | 7/11 [00:00<00:00,  7.75substep/s]
 73%|███████▎  | 8/11 [00:00<00:00,  7.80substep/s]
 82%|████████▏ | 9/11 [00:01<00:00,  7.62substep/s]
 91%|█████████ | 10/11 [00:01<00:00,  7.47substep/s]
100%|██████████| 11/11 [00:01<00:00,  7.46substep/s]
100%|██████████| 11/11 [00:01<00:00,  7.80substep/s]