Atom topic feed | site map | contact | login | Protection des données personnelles | Powered by FluxBB | réalisation artaban
You are not logged in.
Dear all,
I am a new user of code_aster. I am planning to use it to calculate the displacement and deformation of a super flexible structure (fishing net). At the beginning phase, a validation simulation is conducted to test the ability of this FEM solver. The comparative study is liked that:
I have set up the mesh and non-linear static solver, but it breaks out during the calculation. (see the attachment)
Could anyone give advice or help me to fix the problem?
Last edited by chenghui62000 (2019-02-23 23:33:26)
Offline
I am a new user of code_aster.
Welcome!
I have set up the mesh and non-linear static solver, but it breaks out during the calculation. (see the attachment)
Could anyone give advice or help me to fix the problem?
I see that you have modeled the net with bar element. Do you think the net is structurally stable in such a case?
Moreover, the load applied is vertical but all the bar elements are horizontal (initially). How is it going to achieve force equilibrium? I think answering these question will lead you to the solution.
I also feel modelling should be done with cables in this case.
Offline
Hello,
I found dbpatankar's answer very interesting so I tried to run the analysis with DYNA_NON_LINE in order to overcome instabilities and understand the problem better.
The deformations seem correct with POU_D_T but not with BARRE elements:
I cannot figure out why the BARRE deformation is so strange, is it because of the vertical loads on horizontal elements? The deformations of the BARRE elements are huge, on the above picture I reduced the forces by 100 in order to be able to display it.
Konyaro
Last edited by konyaro (2019-02-03 07:56:17)
失敗は成功のもと (L'échec est la base de la réussite)
Offline
Thank you for your answering.
I see that you have modeled the net with bar element.
Sorry, I did not see the reasons why the bar elements cannot be used here. I will try to use the calbe elements to reconduct this simulation again.
The deformations seem correct with POU_D_T but not with BARRE elements
It seems very interesting in your results, I was also planning to use the POU_D_T elements. Maybe it can give the same results with the cable element.
The reference simulation can be found in the flowing article. It used a spring element to do simulation the deformation of the flexible system. "Physical modeling for underwater flexible systems dynamic simulation" (sorry I could not put the link in the reply)
By administrator, the link:
https://www.sciencedirect.com/science/a … 1804001532
Offline
hello
i have done a few studies like that
it HAS TO BE DONE with CABLE element
i had like to run your case but i cannot find any .comm file in the archive
can you post it here
after that i may explain you why it does not work with BARRE of POU_D_* elements
jean pierre aubry
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hello Jean-Pierre,
attached the original comm file, the modified ones and the mesh.
Regards,
Konyaro
失敗は成功のもと (L'échec est la base de la réussite)
Offline
attached is a modified .comm file for CABLE
study it
i have put a gravity preload to make it run
the material is very strange to me, very soft
beware of too short steps in STAT_NON_LINE sometimes it is not good at all
with BARRE we have an assembly with no continuous tangent on the nodes
and the elements remain straight
and probably if the structure is in a plane and loaded perpendicular to it
it does not converge at the very step
with POU we have bending stiffness which is not the truth
Last edited by jeanpierreaubry (2019-02-21 14:34:11)
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
nice! I tried to remove the gravity at the end of the simulation and it does converge. The gravity is really needed at the very beginning of the simulation, that was a good idea Jean-Pierre!
失敗は成功のもと (L'échec est la base de la réussite)
Offline
nice! I tried to remove the gravity at the end of the simulation and it does converge. The gravity is really needed at the very beginning of the simulation, that was a good idea Jean-Pierre!
yes the very first steps may be a headache with cable structure without a preload more than ever if the load is pependicular i guess
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Dear all,
Thank you for your help!
Sorry, I am a little to get back this topic because of my Ph.D courses.
1.
i have done a few studies like that
it HAS TO BE DONE with CABLE element
Thank you very much for this information. It saves my time to try different elements. I think that the element's names can be different in software. I had done a similar simulation in ANSYS using beam or bar elements.
2.
the material is very strange to me, very soft
I have changed the density of the material in the new .comm file to make it more realistic (see in the attachment). This material is wildly used in the fishing industry. It is a twine rope in Nylon material, and it is very flexible compared to the metal cable.
3.
with BARRE we have an assembly with to continuous tangent on the nodes
and the elements remain straight
and probably is the structure is in a plane and loaded perpendicular to it
it does not converge at the very step
Thank you for explaining why BARRE cannot be used here. But I am still a little confused that: are the motion equations of the discrete elements is different between CABLE and BARRE? From my point of view, it is may because of the nonlinear behaviour we apply in the cable (see the below code) makes it can solve this problem.
COMPORTEMENT=_F(DEFORMATION='GROT_GDEP',
GROUP_MA=('twines', ),
RELATION='CABLE'),
4.
The gravity is really needed at the very beginning of the simulation,
Thank you very much for your contributions to this topic. Both yours and jean pierre's comm file helps me a lot to understand the simulation process. I think I got the same results as you.
There are some questions I would like to put out here:
5. Is the thermal dilation coefficient is necessary during the material definition? I have deleted this parameter, and the simulation can still give the same results.
6. The initial load (pre-tension, I think) is important for the simulation because it can avoid the singular matrix. But, it won't bring overestimation in the final results?
elemprop = AFFE_CARA_ELEM(CABLE=_F(GROUP_MA=('twines', ),
N_INIT=10.0,
SECTION=3.14159265359e-06),
MODELE=model)
7. There only two forces on the nodes are actually applied to the simulation, I could find the reason why. because I think I have applied all the three node forces in my comm file.
Thank you very much again and best regards,
Hui Cheng
Last edited by chenghui62000 (2019-02-11 20:06:50)
Offline
5. Is the thermal dilation coefficient is necessary during the material definition? I have deleted this parameter, and the simulation can still give the same results.
Useless because no thermal load.
6. The initial load (pre-tension, I think) is important for the simulation because it can avoid the singular matrix. But, it won't bring overestimation in the final results?
According to U4.42.01 §16.3:
"This keyword makes it possible to define an initial pre-tensioning allowing the convergence of
calculation into nonlinear (it is thus useless into linear). It is applied only to the first step of time. "
Therefore it does'nt affect the final result, you can try to modify this value, it has no effect.
7. There only two forces on the nodes are actually applied to the simulation, I could find the reason why. because I think I have applied all the three node forces in my comm file.
The force F3 is applied to the GROUP_NO F2.
Last edited by konyaro (2019-02-13 21:09:47)
失敗は成功のもと (L'échec est la base de la réussite)
Offline
Konyaro, Thank you so much!
The force F3 is applied to the GROUP_NO F2.
Sometimes, I could not find this kind of mistake by myself.
Offline
Hello, Thank you all for your kindly help.
There is one more question for this case: Can I use the DYNA_NON_LINE solver to do the same simulation?
Assume the same boundary condition, it should get the same result?
I have tried to modify the .comm file, but the outcome seems wrong (at least not same with the STAT_NON_LINE simulation).
Could you please help me to check the .comm file?
Thank you in advance.
Hui Cheng
Last edited by chenghui62000 (2019-02-21 08:52:55)
Offline
i will answer your question with a question
what is the difference between DYNA_NON_LINE and STAT_NON_LINE?
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hi, Jean Pierre Aubry,
From my point of view, the distinction between the dynamic and the static analysis on the basis of whether the applied load has enough acceleration in comparison to the structure's natural frequency. If a load is applied sufficiently slowly, the inertia forces (Newton's first law of motion) can be ignored and the analysis can be simplified as static analysis.
Dynamic analysis can be used to find dynamic displacements, time history, and modal analysis. A dynamic analysis is also related to the inertia forces developed by a structure when it is excited by means of dynamic loads applied suddenly (e.g., wind blasts, explosion, water current forces).
Best regards,
Hui Cheng
Offline
to keep going with the argument
you are generally right
but i cannot see any initial speed or acceleration in your comm file
i have never done a DYNA_NON_LINE analysis whose .comm file was a direct translation of STAT_NON_LINE .comm file
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hi Jean Pierre Aubry,
Thank you for your reply. Please forgive my knowledge gap, did you mean I need to define an initial speed or acceleration in the .comm file? or I need to define time-dependent load using AFFE_CHAR_MECA_F? Or I need to change the element? I have seen a case in the forum (a cable-coque connection case), there seems no initial definition for speed or acceleration in the .comm file.
BTW, the reason why I want to use the DYAN_NON_LINE: I want to see how the net is deformed in the real-time domain.
Thank you very much,
Hui Cheng
Last edited by chenghui62000 (2019-02-21 14:32:33)
Offline
what i do usually
perform a first STAT_NON_LINE that will set an initial state with the only static load, not varying with time
for example gravity due to the own weight
then perform the DYNA_NON_LINE with it as ETAT_INIT
dynanl=DYNA_NON_LINE(
.......
ETAT_INIT=(
_F(
VITE= vitini,
#some initial speed given to some nodes,
#if the load condition call for that,
#for example an object falling on
#may not be necessary
EVOL_NOLI= statnl, #the previous STAT_NON_LINE analysis concept name
),
),
......
AFFE_CHAR_MECA_F is not necessary here
there is an example of what AFFE_CHAR_MECA_F can be used for in my book
time steps where the load is applied has to be representative of the truth
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hello,
You can transform your STAT_NON_LINE to DYNA_NON_LINE keeping exactly the same boundary conditions and adding, for instance:
SCHEMA_TEMPS=_F(SCHEMA='HHT',
ALPHA=-0.3,
FORMULATION='DEPLACEMENT',),
All the initial velocities and accelerations will be automatically set to zero.
Without damping your net will be oscillating forever so it will be hard to compare the results with the static analysis.
You can add damping either in the materials or use a dissipative scheme, for instance HHT with a low alpha-value, for instance -0.3. That should stabilize the net after a few oscillations.
Konyaro
Last edited by konyaro (2019-02-22 17:30:18)
失敗は成功のもと (L'échec est la base de la réussite)
Offline
@konyaro
you are absolutely right
may i had that setting a
load FONC_MULT 0,0, 1.1 and a INCREMENT of the same value
or
load FONC_MULT 0,0, 100.1 and a INCREMENT of the same value
will produce the same results in STAT_NON_LINE
but
may produce radically different ones in DYNA_NON_LINE
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Thank you Konyaro and Jean Pierre Aubry,
It is your kindly help makes me feel confident in this software.
Offline