Hi everyone!
I have modeled a landslide in FLAC3D. I am trying to reproduce the stress state in a small portion (a box) of the landslide body, inside which a tunnel is located, in order to deepen the analysis on the consequences of landslide movements on the tunnel itself. Is it possible to do this either in FLAC3D or in 3DEC?
I have already tried exporting the stress state data on the outer surface of the said box in FLAC3D and applying it to faces along the boundaries of the same box in 3DEC using the command ‘block face apply’, but it seems that in this way the stress values are not maintained during cycling. So I am wondering if there is a better approach or command. Thanks!
Not sure if I understood correctly, but for these kind of problems the way to go is to impose velocity/displacement boundary conditions at the gridpoints, you’re basically doing something called ‘nested models’ or ‘sub-modelling’.
You could even make a table with a history of velocities/displacements and apply them at the desired gridpoints. This forces the boundaries of your small model to move exactly as dictated by the larger landslide simulation.
I’ve never done this myself, but theoretically is the way to go.
Best,
1 Like
Hi Pietro,
I have recently transferred stresses from a big model to a smaller box using python. You can extend the code to do the same with gridpoint information.
Best,
Josef
import itasca as it
import numpy as np
it.command(‘python-reset-state false’)
new_model = ‘new_mesh.f3sav’
stress_model = ‘existing_model.f3sav’
save_name = ‘stress_transfer.f3sav’
zone_info =
get zone location from new model
it.command(“”"
model restore “{}”
“”".format(new_model))
for zp in it.zone.list():
z_id = zp.id()
x, y, z = zp.pos()
temp = {
‘id’ : z_id,
‘x’ : x,
‘y’ : y,
‘z’ : z,
}
zone_info.append(temp)
get stresses and deformations from real model
it.command(“”"
model restore “{}”
“”".format(stress_model))
for ind, item in enumerate(zone_info):
loc = np.array([item[‘x’],item[‘y’],item[‘z’]])
zone = it.zone.near(loc)
zone_info[ind][‘stress’] = zone.stress()
set the stress state to the new model
it.command(“”"
model restore “{}”
“”".format(new_model))
for ind, item in enumerate(zone_info):
loc = np.array([item[‘x’],item[‘y’],item[‘z’]])
zone = it.zone.near(loc)
zone.set_stress(item[‘stress’])
it.command(“”"
model save “{}”
“”".format(save_name))
1 Like