Heterogeneous material such as position-dependent Youngs modulus

joergbuchwald · April 22, 2021, 5:33pm

What I would do is defining such data on the mesh beforehand.
I added some functionality to VTUinterface to do that:

joergbuchwald/VTUinterface/blob/master/examples/vtu_create_field.py

'''
This example demonstrates the usage of vtuIO.VTUIO to create point field (arrays)
on the mesh.
'''


import vtuIO	# to read and process (point interpolation) vtu- and pvd-files 


# read file
data=vtuIO.VTUIO("square_1e2_pcs_0_ts_1_t_1.000000.vtu", dim=2)

def get_pressure(x,y,z):
    if x<=0.5:
        return -5.0e3
    else:
        return 5.0e3

data.func_to_field(get_pressure, "p0","p0_field.vtu")

This file has been truncated. show original