|
[Sponsors] |
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 28, 2014, 10:42 |
Exact Stress
|
#221 |
New Member
Join Date: Mar 2014
Posts: 23
Rep Power: 12 |
How can exact stresses distributions (normal & shear stresses) along the left BC be obtained in a cantilever beam under tip shear load, as shown in attached Fig?
|
|
May 28, 2014, 10:47 |
Dual BC
|
#222 |
New Member
Join Date: Mar 2014
Posts: 23
Rep Power: 12 |
The missing paper in the previous question can be found in below.
http://www.sciencedirect.com/science...68874X08000140 |
|
June 2, 2014, 09:40 |
Shell & Plate
|
#223 |
New Member
Join Date: Mar 2014
Posts: 23
Rep Power: 12 |
Hi,
Is there any solver for shell & plate, please? |
|
June 3, 2014, 03:38 |
|
#224 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
Hi,
Currently there are only continuum elements implemented within the OpenFOAM solid mechanics solvers, so there are no shell/plate/beam/membrane/etc elements yet; but it is possible to implement them using the Finite Volume method, see eg Ivankovic, A. & Demirdzic, I. (1994). Finite Volume Approach to Modelling of Plates. Computational Mechanics (PM me if you would like a copy). As regards your previous question about application of BCs to the cantilever case, it is possible to apply all the BCs within the paper but you will need to code some yourself, in particular for the BCs with a fixed zero normal displacement and an applied parabolic profile shear traction: it would be very similar to fixedDisplacementZeroShear. Philip |
|
June 3, 2014, 12:22 |
|
#225 | |
New Member
Join Date: Feb 2014
Posts: 9
Rep Power: 12 |
Quote:
--> FOAM Warning : From function eigenValues(const tensor&) in file primitives/Tensor/tensor/tensor.C at line 170 but i can not find this line rho = rho/J in writeFields.H to fix !!! the file writefields.H is that on the link http://sourceforge.net/p/openfoam-ex.../writeFields.H which does not contian this line!! then the question what should I fix?! |
||
June 3, 2014, 12:46 |
|
#226 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
The change mentioned should be fixed in this file:
$FOAM_SOLVERS/solidMechanics/elasticNonLinTLSoldFoam/writeFields.H You can see it here. Though, I am not sure if this is related to your warning: Quote:
|
||
June 3, 2014, 14:11 |
|
#227 |
New Member
Join Date: Feb 2014
Posts: 9
Rep Power: 12 |
ok Philip it is not the solver which I use, but i use elasticNonLinULSolidFoam and after some time steps when my beam is starting to give the desired deformation I got this warning and then it crashes.
Code:
--> FOAM Warning : From function eigenValues(const tensor&) in file primitives/Tensor/tensor/tensor.C at line 170 complex eigenvalues detected for tensor: (24.0327 -0.573962 7.72431 -84.7865 0.703912 -23.9374 -2646.74 -70.2494 2.14876) Floating point exception (core dumped) |
|
June 5, 2014, 02:58 |
|
#228 |
New Member
Join Date: Mar 2013
Posts: 17
Rep Power: 13 |
Hi Philip,
I have another question regarding the boundary conditions for the solid mechanics solvers. For an FSI run, is it possible to first apply a traction boundary condition, and at the end of the time step a time varying fixed displacement boundary condition? I'm currently working with something similar to a geometric multigrid setup, and at the end of the time step I need to synchronize all the solid mechanics solvers. Thanks a lot! David |
|
June 5, 2014, 03:54 |
a problem of icoFsiElasticNonLinULSolidFoam
|
#229 |
Member
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 13 |
dear Philip,
there is another way to use timeVaringFixedDisplacement boundary in icoFsiElasticNonLinULSolidFoam with no need for changing $FOAM_SRC/solidModels/fvPatchFields . for fluid part I change U to Ufluid , and for solid part I change Usolid to U . after wmake , the timeVaringFixedDisplacement boundary works ok in icoFsiElasticNonLinULSolidFoam. but I found another prolbem ! there is something wrong with moving mesh for fluid part when using timeVaryingUniformFixedValue boundary in motionU Code:
ground { type timeVaryingUniformFixedValue; fileName "$FOAM_CASE/constant/timeVsVelo.txt"; outOfBounds clamp; value uniform (0 0 0); } I have tried this boundary condition in icoDyMFoam , it works good . I don't know what is wrong ? maybe it has something to do with pointer ? I had heard that for strong FSI It's necessary to move mesh(fluid) back to previous point field in each subiteration followed by another update with the new point field, the pointers were switched back into the right position. do you have any suggestion? thanks xuhe |
|
June 9, 2014, 20:33 |
limit cracking to specified patch
|
#230 |
Member
Eric Bryant
Join Date: Sep 2013
Location: Texas
Posts: 44
Rep Power: 13 |
Hi Phil -
Can I limit crack origination to only from a specified patch? Looking at ACP's crackingBiMatDcbLinear tutorial: system/fvSolution Code:
solidMechanics { ... //crackPropagationFromSpecifiedPatches no; crackPropagationFromSpecifiedPatches yes; ... } Code:
crackerFvMeshCoeffs { crackPatch crack; openPatch tractionFree; } I'm seeing cracking from non-openPatch patches, so wondered about these configuration options... Best, Eric |
|
June 12, 2014, 10:36 |
Plate
|
#231 |
New Member
Join Date: Mar 2014
Posts: 23
Rep Power: 12 |
Hi,
Could you please provide the following journal paper. Ivankovic, A. & Demirdzic, I. (1994). Finite Volume Approach to Modelling of Plates. Computational Mechanics Thanks. |
|
June 12, 2014, 10:37 |
|
#232 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
||
June 17, 2014, 04:04 |
|
#233 |
Member
xuhe-openfoam
Join Date: Aug 2013
Location: DaLian,china
Posts: 82
Rep Power: 13 |
hi,
I know we could use sampleDict to get Usolid at patch when posting results from the two fsiFoam in solidMechanics . but I wonder if it is possible to use libsampling.so in controlDict to get Usolid at patch ? I have tried ,but only succeed when sampling fluid field . thanks! |
|
June 17, 2014, 05:09 |
fsi classes
|
#234 |
New Member
Ireneusz Czajka
Join Date: Nov 2013
Posts: 6
Rep Power: 12 |
There are slides by Robert L. Campbell and Bent A. Craven from 6th openfoam workshop. In their pressentation authors have mentioned about classes: structure and fsiInterface. Do you know, how to get such classes? Are they available somewhere in the internet for download?
|
|
June 17, 2014, 06:20 |
|
#235 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
Quote:
Best regards, Philip |
||
June 17, 2014, 06:21 |
|
#236 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
Quote:
Philip |
||
June 26, 2014, 08:29 |
icoFsiElasticNonLinULSolidFoam
|
#237 |
New Member
Petteri
Join Date: Feb 2014
Posts: 11
Rep Power: 12 |
I managed to run icoFsiElasticNonLinULSolidFoam from the foam-extend3.0 for a simple case of submerged hydrofoil under an ice sheet and the results are very similar to what I got with an in-house developed code. However, now that I'm trying to dig in more to the icoFsiElasticNonLinULSolidFoam code and since I don't have much experience with C++ or FSI some questions came up.
-First of all, why is the solid mesh not updated within the sub-iteration (outerCorr) loop in icoFsiElasticNonLinULSolidFoam.C? -What does createStressPointMesh.H do and what is the StressMesh created for? Petteri |
|
June 26, 2014, 09:56 |
Error using PyFoam on solid mechanics solvers
|
#238 |
Member
|
Hi Philip and everyone else,
I am using the solid mechanics solvers (elasticSolidFoam etc) on OpenFoam 2.3.0. The solver runs well but the problem arises when I try to plot the residuals with pyFoam by using pyFoamPlotWatcher.py. I get the following error multiple times: Code:
Time: 5 DICPCG: Solving for Ux, Initial residual = 7.86012e-08, Final residual = 7.86012e-08, No Iterations 0 PyFoam WARNING on line 138 of file /usr/local/lib/python2.7/dist-packages/PyFoam/LogAnalysis/GeneralLineAnalyzer.py : Problem doing timelines could not convert string to float: PyFoam WARNING on line 138 of file /usr/local/lib/python2.7/dist-packages/PyFoam/LogAnalysis/GeneralLineAnalyzer.py : Problem doing timelines could not convert string to float: DICPCG: Solving for Uy, Initial residual = 9.42587e-08, Final residual = 9.42587e-08, No Iterations 0 PyFoam WARNING on line 138 of file /usr/local/lib/python2.7/dist-packages/PyFoam/LogAnalysis/GeneralLineAnalyzer.py : Problem doing timelines could not convert string to float: PyFoam WARNING on line 138 of file /usr/local/lib/python2.7/dist-packages/PyFoam/LogAnalysis/GeneralLineAnalyzer.py : Problem doing timelines could not convert string to float: Time 5, Corrector 0, Solving for U using DICPCG, residual = 9.42587e-08, relative residual = 0 Time 5, Solving for U, Initial residual = 9.42587e-08, Final residual = 9.42587e-08, Relative residual = 0, No outer iterations 0 ExecutionTime = 137.95 s ClockTime = 138 s PyFoam WARNING on line 138 of file /usr/local/lib/python2.7/dist-packages/PyFoam/LogAnalysis/GeneralLineAnalyzer.py : Problem doing timelines could not convert string to float: Max epsilonEq = 1.34163e-08 Max sigmaEq = 3212.97 ExecutionTime = 138.01 s Please note that pyFoam runs perfectly with the other solvers in the original OpenFoam release. The problem is only encountered with the OpenFoam Extend solvers. Also I am able to execute the command pyFoamClearCase.py successfully in OpenFoam extend. I will really appreciate if anyone can help me in this regard Jibran Last edited by Jibran; July 3, 2014 at 08:31. |
|
June 27, 2014, 06:25 |
|
#239 | ||
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
Hi Petteri,
Quote:
UL is described a umber of times in literature e.g. Bathe Finite Element Methods, I can send you some FV specific literature if you would like about UL. Quote:
mesh is the mesh for the fluid. Best regards, Philip |
|||
June 27, 2014, 06:27 |
|
#240 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,091
Rep Power: 34 |
Hi Jibran,
Yes I suspect you are correct that the format of the log is not correct for pyFoam. It is very straight-forward to change output format, if we find what format it would like. Best regards, Philip |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
GPU Linear Solvers for OpenFOAM | gocarts | OpenFOAM Announcements from Other Sources | 37 | August 17, 2022 14:22 |
[Virtualization] OpenFOAM oriented tutorial on using VMware Player - support thread | wyldckat | OpenFOAM Installation | 2 | July 11, 2012 16:01 |
New OpenFOAM Forum Structure | jola | OpenFOAM | 2 | October 19, 2011 06:55 |
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 | wyldckat | OpenFOAM Announcements from Other Sources | 3 | September 8, 2010 06:25 |
OpenFOAM Debian packaging current status problems and TODOs | oseen | OpenFOAM Installation | 9 | August 26, 2007 13:50 |