example.py

import pathlib

import numpy as np

from microstructure_ve import (
    Heading,
    GridNodes,
    PeriodicBoundaryCondition,
    RectangularElements,
    ElementSet,
    ViscoelasticMaterial,
    Material,
    periodic_assign_intph,
    load_viscoelasticity,
    DisplacementBoundaryCondition,
    Dynamic,
    Step,
    NodeSet,
    Model,
    Simulation,
)

scale = 0.0025
displacement = 0.005
layers = 5

ms_img = np.load("ms.npy")
intph_img = periodic_assign_intph(ms_img, [layers])

base_path = pathlib.Path(__file__).parent
youngs_path = base_path / "PMMA_shifted_R10_data.txt"
freq, youngs_cplx = load_viscoelasticity(youngs_path)
# This is one way to assign a long term modulus, but it is not universal!
# Another strategy is to use 0 for true viscoelastic liquids.
# Pick something physically reasonable for your system.
youngs_plat = youngs_cplx[0].real

heading = Heading("Example RVE simulation")
nodes = GridNodes.from_intph_img(intph_img, scale)
drive_nset = NodeSet("DRIVE", [nodes.virtual_node])
elements = RectangularElements(nodes)
filler_elset, intph_elset, mat_elset = ElementSet.from_intph_image(intph_img)

filler_material = Material(filler_elset, density=2.65e-15, youngs=5e5, poisson=0.15)
intph_material = ViscoelasticMaterial(
    intph_elset,
    density=1.18e-15,
    poisson=0.35,
    shift=-4.0,
    youngs=youngs_plat,
    freq=freq,
    youngs_cplx=youngs_cplx,
    left_broadening=1.8,
    right_broadening=1.5,
)
mat_material = ViscoelasticMaterial(
    mat_elset,
    density=1.18e-15,
    poisson=0.35,
    youngs=youngs_plat,
    freq=freq,
    youngs_cplx=youngs_cplx,
    shift=-6.0,
)
model = Model(
    nodes=nodes,
    nsets=[drive_nset],
    elements=elements,
    materials=[filler_material, intph_material, mat_material],
    bcs=[
        PeriodicBoundaryCondition(
            nodes=nodes, nset=drive_nset, first_dof=1, last_dof=1, displacement=0.0
        )
    ],
)

disp_bc = DisplacementBoundaryCondition(
    drive_nset,
    first_dof=1,
    last_dof=1,
    displacement=displacement,
)
dyn = Dynamic(
    f_initial=1e-7,
    f_final=1e5,
    f_count=30,
    bias=1,
)
step = Step(subsections=[dyn, disp_bc], perturbation=True)

with open("example.inp", mode="w", encoding="ascii") as inp_file_obj:
    Simulation(
        heading=heading,
        model=model,
        steps=[step],
    ).to_inp(inp_file_obj)

# from microstructure_ve import run_job, read_odb
#
# run_job("example", 4)
# read_odb("example", drive_nset)

# import csv
# tsv = csv.reader(open("example-reaction-force.tsv", "r"), dialect=csv.excel_tab)