44 |
|
|
45 |
fix=VectorConstrainerOverBox(debug=d) |
fix=VectorConstrainerOverBox(debug=d) |
46 |
fix.domain=Link(dom,"domain") |
fix.domain=Link(dom,"domain") |
47 |
fix.value=0. |
fix.top=[False, False, False] |
48 |
fix.top=False |
fix.bottom= [True, True, True] |
49 |
fix.bottom=True |
fix.front= [False, False, False] |
50 |
fix.front=False |
fix.back=[False, False, False] |
51 |
fix.back=False |
fix.left=[False, False, False] |
52 |
fix.left=False |
fix.right=[False, False, False] |
|
fix.right=False |
|
|
# dens_dot.density_rate=Link(hist,"mass_changes") |
|
53 |
|
|
54 |
el=LinearElasticStressChange(debug=d) |
el=LinearElasticStressChange(debug=d) |
55 |
el.domain=Link(dom,"domain") |
el.domain=Link(dom,"domain") |
56 |
el.tag_map=Link(dom,"tag_map") |
el.tag_map=Link(dom,"tag_map") |
57 |
el.density=1. |
el.density=8e3*0 |
58 |
el.lame_lambda=2. |
el.lame_lambda=1.7e11 |
59 |
el.lame_mu=1. |
el.lame_mu=1.7e11 |
60 |
el.location_of_fixed_displacement=Link(fix,"location_of_constraint") |
el.location_of_fixed_displacement=Link(fix,"location_of_constraint") |
61 |
el.density_rate=Link(dens_dot,"density_rate") |
el.density_rate=Link(dens_dot,"density_rate") |
62 |
|
|
|
# hist.domain=Link(dom,"domain") |
|
|
|
|
|
|
|
63 |
s=Simulation([sq, hist, dens_dot, fix, el], debug=d) |
s=Simulation([sq, hist, dens_dot, fix, el], debug=d) |
64 |
s.run() |
s.run() |