 |
 |
 |
|
|
|
[Sponsors] | ||||
December 19, 2019, 02:34
|
MultiMaterial in FSI problems
|
#1 |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7  |
Dear all
In solids4Foam solver for solid physics we have a biMatPlateHoleTotalLag tutorial that uses two material. So this solver can work with multiMaterial. The Question is: How can combine twoElasticBeam tutorial with biMatPlateHoleTotalLag. I change materialProperties of solid regions in fluidSolidInteraction physics to two SOLID1 and SOLID2 and use two properties. Solver start to calculate and check system folder and fvScheme, fvSolution to code them to SOLID1SubSet and SOLID2SubSet, Then check materialProperties in constant folder and PolyMesh of constant folder. It need to solver goes to -region solid folder and check All folder of it. Do you think copy other file to root folder of tutorial is OK? or modifying the code? Thanks |
|
|
|
|
|
December 26, 2019, 01:44
|
|
#2 | |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
As it seen, in the mechanicalProperties with more than one material, the solver use mechanicalModel.C and check line 132
Quote:
|
||
|
|
|
||
|
January 10, 2020, 03:18
|
|
#3 | |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
May it due to definition of IOobject in mechanicalModels?
for exmaple: Quote:
|
||
|
|
|
||
|
January 15, 2020, 05:55
|
|
#4 |
|
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Hi,
To use two or more sets of mechanical properties in the solid region, you just need to specify them in the mechanicalProperties dictionary, e.g. Code:
planeStress no;
mechanical
(
steel
{
type linearElastic;
rho rho [1 -3 0 0 0 0 0] 7854;
E E [1 -1 -2 0 0 0 0] 200e+9;
nu nu [0 0 0 0 0 0 0] 0.3;
}
titanium
{
type linearElastic;
rho rho [1 -3 0 0 0 0 0] 7854;
E E [1 -1 -2 0 0 0 0] 100e+9;
nu nu [0 0 0 0 0 0 0] 0.3;
}
);
Note: linearElastic is for small strain (linearGeometry) solid models only; for finite/large strain (nonlinearGeometry) solid models, you should use an appropriate constitutive law e.g. neoHookeanElastic. Philip |
|
|
|
|
|
|
January 15, 2020, 10:16
|
|
#5 |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Dear Philip
I think this part of code for multiMaterial is only uses for solid physic. for running fluidSolidInteraction physic, the solver try to copy fvScheme and fvSolution in subSet regions (note that in fluidSolidInteraction there isn't any fvSolution file in system folder). Then the solver check constant folder for mechanicalProperties file... In constant folder for solid physic it is not consider, but for fluidSolidInteraction, we have fluid and solid folder in constant folder. Do you sure, this part work with fluidSolidInteraction? Thanks Hojatollah |
|
|
|
|
|
|
January 15, 2020, 11:48
|
|
#6 |
|
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Hi Hojatollah,
Hmnn, yep you are right, it seems that solids4foam looks for mechanicalProperties and fvSchemes in the wrong place for FSI with multi-material solids. A temporary work-around is the following: Code:
$> rm -f system/*Subset $> cd constant $> ln -s solid/mechanicalProperties . $> cd ../system $> ln -s solid/fvSchemes . Philip |
|
|
|
|
|
|
January 16, 2020, 03:11
|
|
#7 | ||
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Hi Philip
it seen, we should add more command to it as: Code:
$> rm -f system/*Subset $> cd constant $> ln -s solid/mechanicalProperties . $> ln -s solid/polyMesh . $> cd ../system $> ln -s solid/fvSchemes . $> ln -s solid/fvSolution . Quote:
Quote:
|
|||
|
|
|
|||
|
January 19, 2020, 23:04
|
|
#8 | |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
here is log of simulation
Quote:
|
||
|
|
|
||
|
January 22, 2020, 07:14
|
|
#9 |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Dear Philip
Do you have any idea about it? Thanks |
|
|
|
|
|
|
April 16, 2021, 10:57
|
|
#10 |
|
Member
Cyrille Bonamy
Join Date: Mar 2015
Location: Grenoble, France
Posts: 86
Rep Power: 11 |
Dear Hojatollah,
You can find a discussion about this issue here : https://bitbucket.org/philip_cardiff...n-fsi-problems But the quick fix is : avoid using the solidPointDisplacement functionObject. Thanks again Philip 
|
|
|
|
|
|
|
June 24, 2021, 18:29
|
|
#11 |
|
New Member
Join Date: Jun 2021
Posts: 7
Rep Power: 5 |
Dear cyss38,
Are there any permissions or something like this to open the link at bitbucket? I can only open the link. After my registration I am not allowed to look at the linked thread or discussion about solids4foam... I'm actually facing the problem with two different materials in FSI. Thank you! |
|
|
|
|
|
|
July 5, 2021, 23:39
|
|
#12 | |
|
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Hi
you should communication with Dr. Philip Cardiff via bigphil id to access to solids4Foam repo. Quote:
|
||
|
|
|
||
|
July 7, 2021, 20:18
|
|
#13 | |
|
New Member
Anas Muhamad Pauzi
Join Date: Nov 2019
Posts: 16
Rep Power: 7 |
Quote:
Is there any solution besides of deleting the pointDisp fun Thank you |
||
|
|
|
||
 |
| Tags |
| multimaterial, solids4foam |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Two way FSI related Query | n.soumya | ANSYS | 0 | July 24, 2018 03:50 |
| How To model FSI problems which have piezoelectric material | mechy | CFX | 2 | July 21, 2012 22:16 |
| How To model FSI problems which have piezoelectric material | mechy | ANSYS | 0 | June 30, 2012 01:40 |
| Needed Benchmark Problems for FSI | Mechstud | Main CFD Forum | 4 | July 26, 2011 13:13 |
| What is weakstrong coupling in FSI problems | hajo | OpenFOAM Running, Solving & CFD | 5 | May 15, 2008 02:45 |
3Likes
Linear Mode
