-
Notifications
You must be signed in to change notification settings - Fork 8
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Enable multi-material #119
base: main
Are you sure you want to change the base?
Conversation
57d4112
to
45d0c32
Compare
3f6efe0
to
315f5cb
Compare
simplifies template parameters with single-species tag
332546b
to
8686675
Compare
8686675
to
c089abf
Compare
|
||
KOKKOS_INLINE_FUNCTION void operator()( const int i ) const | ||
{ | ||
c( i, i ) = micromodulus( i ); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should move this to a separate loop
exec_space(), low_corner, high_corner, num_cells, halo_width ); | ||
|
||
std::array<double, 2> same_d = { delta, delta }; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
We shouldn't modify the Kalthoff Winkler example. If this is just for testing, we need to create a new example.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I agree, I will remove when you push a new example
"density" : {"value": 3980, "unit": "kg/m^3"}, | ||
"elastic_modulus" : {"value": 370e+9, "unit": "Pa"}, | ||
"elastic_modulus" : {"value": [370e+9, 100e+9, 10e+9], "unit": "Pa", "note": "11, 22, 12"}, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Why do we have 3 values here? I believe we only need for a bimaterial the elastic modulus of material 1 and material 2.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I moved this to draft. But until we have a bi-material example, I need places to test things are working
@@ -82,7 +82,8 @@ void thermalDeformationExample( const std::string filename ) | |||
// Force model | |||
// ==================================================== | |||
auto force_model = CabanaPD::createForceModel( | |||
model_type{}, CabanaPD::Elastic{}, *particles, delta, K, alpha, temp0 ); | |||
model_type{}, CabanaPD::Elastic{}, CabanaPD::TemperatureDependent{}, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This does not seem related to bimaterial models.
Add type lists to new force models, with averaged properties for cross terms
Closes #109