Yade

VTKRecorder for show breaklines in concrete

Asked by rayolau

Hi all,

First thanks in advance. My goal is show the breaklines in a slab concrete, for differents positions and types of de fissure. For that, in my first approach, I modified the uniax scritp, changing the cylinder geometry by a rectangle with a fissure in the central zone. Also, I eliminated PyRunner corresponding to "stopdamage" and corresponding "stopdamage" def function.

For the postprocessing, I downloaded ParaView and add, in the script, the Engine:

VTKRecorder(fileName='postprocessing',recorders=['all'],iterPeriod=10)

The script works well, but in the postprocessing, I can not see the zones with more stress.

I think that I need put vtkExporter.exportSomething(params) in PyRunner or maybe use Micro-stress with TesselationWrapper, but I don't know how. Any ideas??

Thanks everyone!,
Laura.

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jan Stránský
Solved:
Last query:
Last reply:
Revision history for this message
rayolau (rayolau88) said : 		#1

The link to download the script:
https://drive.google.com/file/d/0B6azPyaq-2oNODRtN29ESWxzQjg/view?usp=sharing

Thanks,

Laura.

Revision history for this message
Jan Stránský (honzik) said : 		#2

Hi Laura,
have a look at VTKRecprder docs [1], apart from 'all', you can also specify
'cpm' (not included in 'all', see the docs). Then in paraview you can
choose 'cpmDamage' or 'cpmStress' for particles color. I tried your example
and got most of damage next to the prescribed boundary conditions planes,
so you should probably decrease the strain rate. Sometimes I also export
contact points with information about their damage, the cracks
visualization might (or might not, it depends) look better. I can give you
some example if you are interested.
cheers
Jan

[1]
https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.VTKRecorder.recorders

2014-11-13 13:16 GMT+01:00 rayolau <email address hidden>:

> Question #257425 on Yade changed:
> https://answers.launchpad.net/yade/+question/257425
>
> rayolau posted a new comment:
> The link to download the script:
>
> https://drive.google.com/file/d/0B6azPyaq-2oNODRtN29ESWxzQjg/view?usp=sharing
>
> Thanks,
>
> Laura.
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

Revision history for this message
rayolau (rayolau88) said : 		#3

Hi Jan,
Thanks for reply so fast! I added 'cmp' in recorders:
VTKRecorder(fileName='postprocessing',recorders=['cpm','all'],iterPeriod=10)

and now, in ParaView, when I import .vtu files, I can see the options "cmpdamage and cpmstress" but, as I can see in the figure attached (link), the values for cpmstress are cero. Maybe is a problem of options in the configuration in ParaView (I'm novice with this software).Do you have any idea of why? From what I understand, despite the strain rate, were you able to see the breaklines in the slab.

About how to export contact points with the information about their damage, I would like that you show me an example for how to do. A lot of thanks for your help.

Laura.

figure link
https://drive.google.com/file/d/0B6azPyaq-2oNeGR6TjYwbW85b00/view?usp=sharing

Revision history for this message
Best Jan Stránský (honzik) said : 		#4

Hi Laura,

paraview is very powerful, you can also use python scripting for pictures
creation :-) have a look at [1] for very short yade introduction.

> and now, in ParaView, when I import .vtu files, I can see the options
> "cmpdamage and cpmstress" but, as I can see in the figure attached
> (link), the values for cpmstress are cero.

you have to set time (upper top in PV) to nonzero value, there are also
some arrows for gui move or play. In zeroth step (as in attachment), the
stress is (and should be) zero. Furthermore, stress is saved as tensor, so
you can plot only specific value (0=xx,1=xy,2=xz,3=yx,4=yy,...,8=zz) or use
'Calculator' filter to play with it (to display mean stress as average of
sxx,syy and szz for instance). Just on the left from 'cpmStress' there is
"rescale to data range" icon, quite nice to use. Or one more left there is
'edit color map', where you can set the range manually.

> Maybe is a problem of options
> in the configuration in ParaView (I'm novice with this software).Do you
> have any idea of why?

see above

> From what I understand, despite the strain rate,
> were you able to see the breaklines in the slab.
>

yes, I can see clear breaklines on the slab sides, and slightly in the
middle. For damage visualization I prefer cpmDamage more than cpmStress. I
am not familiar with 'breaklines' term, so it does not need to be true for
your case.

>
> About how to export contact points with the information about their
> damage, I would like that you show me an example for how to do. A lot of
> thanks for your help.
>

sorry, just now I don't have time to do so, I will send you something in 2
weeks. Write me reminding email if I forget..

good luck in learning paraview :-)
cheers
Jan

[1] https://yade-dem.org/doc/user.html#id3

Revision history for this message
rayolau (rayolau88) said : 		#5

Thanks Jan Stránský, that solved my question.

Revision history for this message
rayolau (rayolau88) said : 		#6

Hi Jan,

thanks for your reply. I use successfuly ParaView and I can see the damage.

Cheers,
Laura.