CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Problems of a structural solver that works in an updated configuration

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 30, 2008, 06:58
Default Hello again, I tried the bo
  #1
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 7
fw407 is on a distinguished road
Hello again,

I tried the boundary conditions I mentioned in (3), and I end up with something like:
.
.
.
.
.

Time = 10

Setting traction on solid patch
ICCG: Solving for DUx, Initial residual = 0, Final residual = 0, No Iterations 0
ICCG: Solving for DUy, Initial residual = 0, Final residual = 0, No Iterations 0
ICCG: Solving for DUz, Initial residual = 0, Final residual = 0, No Iterations 0
ExecutionTime = 0.08 s
.
.
.
.
.
.

I think the problem is due to no load on the boundaries. I am sure the pressure increment is calculated correctly, the problem may lie in setting the pressure("setPressure.H"). could anyone shed some light on this problem?

Kind regards

feng
fw407 is offline   Reply With Quote

Old   July 9, 2008, 06:17
Default Hi everyone, I solved the p
  #2
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 7
fw407 is on a distinguished road
Hi everyone,

I solved the problems I post above, and the solver is now working properly. but there are still problems:

(1) for a simple case which streach a flat plate and calculate the resulting displacement,when the strain is larger than 30%, the computation diverges quikly, this can be solved by taking a very small traction increment, but it takes too long to run it; when I run the static deformable beam case decribled in that paper, the computation diverges quikly even if the strain is very small .

(2) for the flat plate case, the results have a big discrepency(up to 7%) with the analytical solution as the strain increases.

I have looked at these problems several days but made no progress. If someone is interested in this solver, I can email you the solver and the test cases, your comments will be really appreciated!

Kind regards

feng
fw407 is offline   Reply With Quote

Old   July 10, 2008, 23:02
Default hi, I can't help you anythin
  #3
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 7
ivanyao is on a distinguished road
hi,
I can't help you anything,but there is something that I want you to help me,I have go through all the messages in "Parabolic inlet velocity profile".but I don't know how to compile the setParabolicInlet.I am doing a simulation about the compute wind engineering,I do the simulation in rhoturbFoam.my inlet velocity profile:U=U_0*(Z/Z_0)^0.25,"U_0","Z_0"are constants I provide."Z" is the height of the building.I am using OF1.4,and I have downloaded the parabolicVelocity_HJ_17Jan2007.tgz,but I don't know what I should do next.is there anyone help me step by step?that is very important to me.I would very very appreciate it. because no one answer me,and I study in civil. i think you can help me.
thanks,Ivan
ivanyao is offline   Reply With Quote

Old   July 13, 2008, 05:29
Default Hi weihong, Sorry for answe
  #4
Member
 
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 7
fw407 is on a distinguished road
Hi weihong,

Sorry for answering your question so late, in fact I am also new to OpenFOAM and now struggling with C++ basics. but I can give you some ideas, hope that might help.

Please find the following link, I think it has explained clearly how to use it
http://openfoamwiki.net/index.php/Tutorials#How_to_implement_a_new_boundary_cond ition_.28OF-1.4.1.29

Based on my understanding, I think you should take the following steps:

(1)put the interface file(**.H) and implementation file(***.C)at the directory where you put the solvers created on your own, create a folder which contains the the interface file(parabolicVelocity/parabolicVelocityFvPatchVectorField.H) and implementation file(parabolicVelocity/parabolicVelocityFvPatchVectorField.C)

(2) create a "Make" folder which contains a "make" and "options" file, the contents can be found in the link I posted above

(3) use "wmake libso" to compile this utility, and the resulting "*****.so" can be found in the "bin" directory located in the user directory(for example "~/OpenFOAM/snail-1.3/applications/bin").

(4) use the method describable in website of how to use this new B.C

you may run into some errors because of your version of OF, the type of your Linux system and other reasons, I am not so clear, and am still a beginner, if I am wrong, please correct me, your comments are appreciated.

Kind regards

feng
fw407 is offline   Reply With Quote

Old   July 13, 2008, 22:09
Default hi, thank you for answer me!e
  #5
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 7
ivanyao is on a distinguished road
hi,
thank you for answer me!everything goes ok,but when i type "wmake libso" ,it show "/home/ivan/OpenFOAM/OpenFOAM-1.4/wmake/wmake: cannot make.Make directory does not exist",i don't know why.I tell you what is my step.
1.decompressing the file which i download on the website.put it in ~/OpenFOAM/OpenFOAM-1.3/src/finiteVolume/
fields/fvPatchFields/derivedFvPatchFields/ and /home/OpenFOAM/OpenFOAM-1.3/src/
finiteVolume/lnInclude/
2.and then modify the Make/files parabolicVelocityFvPatchVectorField.C
LIB = $(FOAM_ivan_LIBBIN)/libmyBCs and modify the Make/options files ~/OpenFOAM/OpenFOAM-1.3/src/finiteVolume/
fields/fvPatchFields/derivedFvPatchFields/parabolicVelocity and /home/OpenFOAM/OpenFOAM-1.3/src/finiteVolume/lnInclude/ parabolicVelocity
3.type "wmake libso"
my c++ is rusty.if you don't mind,would you try it by yourself?and tell me where i make a mistake,it will not cost you much time.
it can download in the forum:http://www.cfd-online.com/OpenFOAM_D...tml?1215566795,
that is very important to me.
ivanyao is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How solver works? John NUMECA 1 February 24, 2006 09:23
How solver works? John FLUENT 2 January 27, 2006 13:35
How solver works? John CFX 1 January 27, 2006 12:48
How solver works? John Phoenics 1 January 27, 2006 12:25
How solver works? John CD-adapco 1 January 27, 2006 04:10


All times are GMT -4. The time now is 09:54.