Skip to content
Snippets Groups Projects
Commit 67f2f334 authored by Boman Romain's avatar Boman Romain
Browse files

export loads to TSV files

parent 3424c362
No related branches found
No related tags found
No related merge requests found
Pipeline #7159 passed
...@@ -47,28 +47,11 @@ def parms(d={}): ...@@ -47,28 +47,11 @@ def parms(d={}):
'axis_pt1': ['x', 'y', 'z'], 'axis_pt1': ['x', 'y', 'z'],
'axis_pt2': ['x', 'y', 'z'] 'axis_pt2': ['x', 'y', 'z']
} }
# rigid cylinder (used if p['food']=='cylinder')
# R = 20
# p['cylinder'] = {
# 'origin': [-R-8, 0, -150],
# 'radius': R,
# 'width': 100,
# 'friction': 0.1,
# 'penalty': 1e3
# }
# material properties # material properties
# p['density'] = 1.850e-9 # [T/mm³] - bone: 1.850 kg/l
p['Young'] = 17000. # [MPa] elastic modulus - bone: 17-20 GPa p['Young'] = 17000. # [MPa] elastic modulus - bone: 17-20 GPa
p['Poisson'] = 0.3 # [-] Poisson's ratio p['Poisson'] = 0.3 # [-] Poisson's ratio
# numerical parameters
# p['tolNR'] = 1e-6 # [-] equilibrium tolerance
# p['dt0'] = 1.0 # [s] time step size
# gmsh toolbox
# p['use_gmshOld'] = False # use old gmsh interface
p.update(d) p.update(d)
return p return p
...@@ -130,6 +113,12 @@ def solve(p={}): ...@@ -130,6 +113,12 @@ def solve(p={}):
# store data in a structure # store data in a structure
mgroups[name] = MuscleGroup(nodes, tris, ntags, loads) mgroups[name] = MuscleGroup(nodes, tris, ntags, loads)
# save loads to file
with open(f'loads_{name}.tsv','w') as f:
f.write('node#\tfx\tfy\tfz\n')
for i,load in enumerate(loads):
f.write(f"{int(ntags[i])}\t{load.x[0]}\t{load.x[1]}\t{load.x[2]}\n")
if p['food'] == 'fixteeth': if p['food'] == 'fixteeth':
# teeth are fixed along chosen directions # teeth are fixed along chosen directions
create_group(p['contact_pts'], nods_no, create_group(p['contact_pts'], nods_no,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment