Erreur de segmentation (core dumped)

Hi Jabrane,
to let us test your script and help you to solve the problem, please
provide also poutre400_notch.spheres file
cheers
Jan

2015-07-27 10:46 GMT+02:00 Yor1 <email address hidden>:

> New question #269691 on Yade:
> https://answers.launchpad.net/yade/+question/269691
>
> Hi everyone,
>
> I'm trying to determine the uniaxial strain of the beam in 3-points
> bending test.
> When i execute my code i obtain this error message "Erreur de segmentation
> (core dumped)".
> This is a part of my code:
>
> from yade import ymport, utils, pack, export
> from yade import plot
> from pylab import *
> import math
>
> #---------------- SIMULATIONS DEFINED HERE (assembly, material, boundary
> conditions)
>
> #### packing (previously constructed)
> PACKING='poutre400_notch'
> OUT=PACKING+'_flexion_3_points'
>
> #### Simulation Control
> DAMP=0.4 # numerical damping
> saveData=100 # data record interval
> iterMax=200000 # maximum number of iteration (to be adjusted)
> saveVTK=10000 # Vtk files record interval
>
> #### Boundary Conditions
> strainRate=0.002
>
> #### Material microproperties -> Lac du Bonnet granite (cf. A DEM model
> for soft and hard rock, Scholtes & Donze, JMPS 2013)
> intR=1.5# allows near neighbour interaction and coordination number K=13
> (determined with coordinationNumber.py -> to be adjusted for every packing)
> DENS=4000 # could be adapted to match material density:
> dens_DEM=dens_rock*(V_rock/V_particles)=dens_rock*1/(1-poro_DEM) -> poro?
> YOUNG=65e9
> FRICT=10
> ALPHA=0.4
> TENS=8e6
> COH=160e6
> #### material definition
> sphereMat = JCFpmMat(type=1,density=DENS,young=YOUNG,poisson =
> ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH)
> wallMat =
> JCFpmMat(type=0,density=DENS,young=YOUNG,frictionAngle=radians(0))
>
> for mat in (sphereMat,wallMat):
> O.materials.append(mat) # then wallMat will be used if material is not
> specified
>
> #### preprocessing to get dimensions
>
> O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0)))
>
> dim=utils.aabbExtrema()
> xinf=dim[0][0]
> xsup=dim[1][0]
> X=xsup-xinf
> yinf=dim[0][1]
> ysup=dim[1][1]
> Y=ysup-yinf
> zinf=dim[0][2]
> zsup=dim[1][2]
> Z=zsup-zinf
>
> r=1
> #### preprocessing to get dimensions
> R=0
> Rmax=0
> numSpheres=0.
> for o in O.bodies:
> if isinstance(o.shape,Sphere):
> numSpheres+=1
> R+=o.shape.radius
> if o.shape.radius>Rmax:
> Rmax=o.shape.radius
> Rmean=R/numSpheres
>
> O.reset() # all previous lines were for getting dimensions of the packing
> to create walls at the right positions (below) because walls have to be
> genrated after spheres for FlowEngine
>
>
> O.bodies.append(geom.facetCylinder(center=(xinf+X/5.,yinf-0.9*r,Z/2.),radius=r,height=Z,orientation=Quaternion((1,
> 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # bas gauche
> O.bodies.append(geom.facetCylinder(center=(0.8*X,yinf-0.9*r,Z/2.),radius=r,height=Z,orientation=Quaternion((1,
> 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # bas droite
> piston=O.bodies.append(geom.facetCylinder(center=(0.5*X,Y+r,Z/2.),radius=r,height=Z,dynamic=False,orientation=Quaternion((1,
> 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # haut
>
> ### packing
>
> beam=O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat))
>
>
> ### set a color to the spheres
> for o in O.bodies:
> if isinstance(o.shape,Sphere):
> o.shape.color=(0.7,0.5,0.3)
>
> crossSectionArea=X*Z
>
> #---------------- ENGINES DEFINED HERE
>
> #### simulation is defined here (DEM loop, interaction law, servo control,
> recording, etc...)
> ##### simulation piston's motion
>
> for i in range(0,len(piston)):
> O.bodies[piston[i]].state.vel[1]=-1
>
> ##### simulation beam's grains motion
> O.engines=[
> ForceResetter(),
>
> InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb(aabbEnlargeFactor=intR,label='Saabb'),Bo1_Facet_Aabb()]),
> InteractionLoop(
>
> [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR,label='SSgeom'),Ig2_Facet_Sphere_ScGeom()],
>
> [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,label='interactionPhys')],
>
> [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(recordCracks=True,Key=OUT,label='interactionLaw')]
> ),
>
> UniaxialStrainer(strainRate=strainRate,axis=1,asymmetry=0,crossSectionArea=crossSectionArea,blockDisplacements=1,blockRotations=1,setSpeeds=0,stopStrain=0.1,label='strainer'),
>
> GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.4,
> defaultDt=0.1*utils.PWaveTimeStep()),
> NewtonIntegrator(damping=DAMP,label="newton"),
>
> PyRunner(iterPeriod=int(saveData),initRun=True,command='recorder()',label='data'),
>
> VTKRecorder(iterPeriod=int(saveVTK),initRun=True,fileName=OUT+'-',recorders=['spheres','jcfpm','cracks'],Key=OUT,label='vtk')
> ]
>
> Best regards
>
> --
> 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
>

Hi Jan,

My real problem is to determine the bending of the beam.
In my code i mix between bending test and the tensile test. And i think
that it is the problem.
You find the file poutre400_notch.spheres

cheers

2015-07-27 12:36 GMT+02:00 Jan Stránský <
<email address hidden>>:

> Your question #269691 on Yade changed:
> https://answers.launchpad.net/yade/+question/269691
>
> Status: Open => Answered
>
> Jan Stránský proposed the following answer:
> Hi Jabrane,
> to let us test your script and help you to solve the problem, please
> provide also poutre400_notch.spheres file
> cheers
> Jan
>
>
> 2015-07-27 10:46 GMT+02:00 Yor1 <email address hidden>:
>
> > New question #269691 on Yade:
> > https://answers.launchpad.net/yade/+question/269691
> >
> > Hi everyone,
> >
> > I'm trying to determine the uniaxial strain of the beam in 3-points
> > bending test.
> > When i execute my code i obtain this error message "Erreur de
> segmentation
> > (core dumped)".
> > This is a part of my code:
> >
> > from yade import ymport, utils, pack, export
> > from yade import plot
> > from pylab import *
> > import math
> >
> > #---------------- SIMULATIONS DEFINED HERE (assembly, material, boundary
> > conditions)
> >
> > #### packing (previously constructed)
> > PACKING='poutre400_notch'
> > OUT=PACKING+'_flexion_3_points'
> >
> > #### Simulation Control
> > DAMP=0.4 # numerical damping
> > saveData=100 # data record interval
> > iterMax=200000 # maximum number of iteration (to be adjusted)
> > saveVTK=10000 # Vtk files record interval
> >
> > #### Boundary Conditions
> > strainRate=0.002
> >
> > #### Material microproperties -> Lac du Bonnet granite (cf. A DEM model
> > for soft and hard rock, Scholtes & Donze, JMPS 2013)
> > intR=1.5# allows near neighbour interaction and coordination number K=13
> > (determined with coordinationNumber.py -> to be adjusted for every
> packing)
> > DENS=4000 # could be adapted to match material density:
> > dens_DEM=dens_rock*(V_rock/V_particles)=dens_rock*1/(1-poro_DEM) -> poro?
> > YOUNG=65e9
> > FRICT=10
> > ALPHA=0.4
> > TENS=8e6
> > COH=160e6
> > #### material definition
> > sphereMat = JCFpmMat(type=1,density=DENS,young=YOUNG,poisson =
> > ALPHA,frictionAngle=radians(FRICT),tensileStrength=TENS,cohesion=COH)
> > wallMat =
> > JCFpmMat(type=0,density=DENS,young=YOUNG,frictionAngle=radians(0))
> >
> > for mat in (sphereMat,wallMat):
> > O.materials.append(mat) # then wallMat will be used if material is not
> > specified
> >
> > #### preprocessing to get dimensions
> >
> >
> O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0)))
> >
> > dim=utils.aabbExtrema()
> > xinf=dim[0][0]
> > xsup=dim[1][0]
> > X=xsup-xinf
> > yinf=dim[0][1]
> > ysup=dim[1][1]
> > Y=ysup-yinf
> > zinf=dim[0][2]
> > zsup=dim[1][2]
> > Z=zsup-zinf
> >
> > r=1
> > #### preprocessing to get dimensions
> > R=0
> > Rmax=0
> > numSpheres=0.
> > for o in O.bodies:
> > if isinstance(o.shape,Sphere):
> > numSpheres+=1
> > R+=o.shape.radius
> > if o.shape.radius>Rmax:
> > Rmax=o.shape.radius
> > Rmean=R/numSpheres
> >
> > O.reset() # all previous lines were for getting dimensions of the packing
> > to create walls at the right positions (below) because walls have to be
> > genrated after spheres for FlowEngine
> >
> >
> >
> O.bodies.append(geom.facetCylinder(center=(xinf+X/5.,yinf-0.9*r,Z/2.),radius=r,height=Z,orientation=Quaternion((1,
> > 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # bas gauche
> >
> O.bodies.append(geom.facetCylinder(center=(0.8*X,yinf-0.9*r,Z/2.),radius=r,height=Z,orientation=Quaternion((1,
> > 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # bas droite
> >
> piston=O.bodies.append(geom.facetCylinder(center=(0.5*X,Y+r,Z/2.),radius=r,height=Z,dynamic=False,orientation=Quaternion((1,
> > 0, 0), 0),segmentsNumber=20,wire=False,material=wallMat)) # haut
> >
> > ### packing
> >
> >
> beam=O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat))
> >
> >
> > ### set a color to the spheres
> > for o in O.bodies:
> > if isinstance(o.shape,Sphere):
> > o.shape.color=(0.7,0.5,0.3)
> >
> > crossSectionArea=X*Z
> >
> > #---------------- ENGINES DEFINED HERE
> >
> > #### simulation is defined here (DEM loop, interaction law, servo
> control,
> > recording, etc...)
> > ##### simulation piston's motion
> >
> > for i in range(0,len(piston)):
> > O.bodies[piston[i]].state.vel[1]=-1
> >
> > ##### simulation beam's grains motion
> > O.engines=[
> > ForceResetter(),
> >
> >
> InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb(aabbEnlargeFactor=intR,label='Saabb'),Bo1_Facet_Aabb()]),
> > InteractionLoop(
> >
> >
> [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR,label='SSgeom'),Ig2_Facet_Sphere_ScGeom()],
> >
> >
> [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1,label='interactionPhys')],
> >
> >
> [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(recordCracks=True,Key=OUT,label='interactionLaw')]
> > ),
> >
> >
> UniaxialStrainer(strainRate=strainRate,axis=1,asymmetry=0,crossSectionArea=crossSectionArea,blockDisplacements=1,blockRotations=1,setSpeeds=0,stopStrain=0.1,label='strainer'),
> >
> >
> GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=10,timestepSafetyCoefficient=0.4,
> > defaultDt=0.1*utils.PWaveTimeStep()),
> > NewtonIntegrator(damping=DAMP,label="newton"),
> >
> >
> PyRunner(iterPeriod=int(saveData),initRun=True,command='recorder()',label='data'),
> >
> >
> VTKRecorder(iterPeriod=int(saveVTK),initRun=True,fileName=OUT+'-',recorders=['spheres','jcfpm','cracks'],Key=OUT,label='vtk')
> > ]
> >
> > Best regards
> >
> > --
> > 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
> >
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/yade/+question/269691/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/yade/+question/269691
>
> You received this question notification because you asked the question.
>

Hi Jan,

My real problem is to determine the bending of the beam.
In my code i mix between bending test and the tensile test. And i think that it is the problem.
How can i join the file poutre400_notch.spheres in the mailing list?

cheers.

>
>
> How can i join the file poutre400_notch.spheres in the mailing list?
>

you can not :-) please put the content directly to the mail message if the
file is not big, or put it somewhere on the internet and send the link

In my code i mix between bending test and the tensile test. And i think
> that it is the problem.

yes, thanks for the hint, in UniaxialStrainer, you have to specify posIds
and negIds parameter.. but now I don't understand your question:

I'm trying to determine the uniaxial strain of the beam in 3-points bending
> test.

UniaxialStrainer is for imposing boundary conditions of a sample such that
the result is uniaxial strain (tension or compression). negIds and posIds
are bodies on opposite sides, on whose the displacement (velocity in fact)
is applied. I don't see any reason why to use in in 3 point bending..
Shouldn't your question be "I'm trying to determine normal/axial strain of
the beam in 3-points bending test."? :-)

cheers
Jan

Yes that's it: "I'm trying to determine normal/axial strain of the beam in 3-points bending test."
This is the file poutre400_notch.spheres

#format x_y_z_r
7.38853 1.82535 0.681871 0.394846
13.3911 0.384712 1.60798 0.394391
6.70128 1.39815 1.13805 0.393937
4.21398 2.03294 0.980896 0.393482
3.51109 0.967838 0.962592 0.393027
3.12617 1.68253 0.644577 0.392572
6.1753 1.87984 0.810796 0.392118
10.2448 2.62043 0.38459 0.391663
5.10253 0.387385 1.61296 0.391208
14.6194 2.61659 0.556424 0.390754
11.8063 1.66858 0.385115 0.390299
8.00908 1.47264 0.385189 0.389844
11.4241 2.3472 0.382275 0.389389
5.96522 0.377406 1.03497 0.388935
13.7874 0.384502 0.912233 0.38848
5.27655 0.373497 0.67646 0.388025
13.295 2.62372 1.01981 0.387571
1.93232 1.81002 0.851715 0.387116
10.0417 1.87885 0.373648 0.386661
1.34527 1.41541 0.38249 0.386206
9.19451 2.61569 1.61875 0.385752
5.35657 1.28568 1.62158 0.385297
8.72591 2.05614 1.13578 0.384842
2.00504 1.04568 0.466947 0.384388
3.6127 0.380103 0.380232 0.383933
10.3159 0.787797 1.12922 0.383478
8.91993 0.366926 0.37842 0.383023
1.41064 0.866531 0.911929 0.382569
14.1525 1.68776 1.62061 0.382114
4.37599 0.380681 0.375893 0.381659
8.09145 0.982723 0.972467 0.381204
13.5204 2.09942 1.53031 0.38075
12.1593 2.16604 1.17478 0.380295
11.2382 1.14134 0.377277 0.37984
4.61175 1.10362 0.373962 0.379386
8.58242 1.97678 0.378207 0.378931
11.4121 0.368978 1.62607 0.378476
3.62789 2.62839 0.370954 0.378021
11.5986 0.907031 1.14184 0.377567
12.8811 1.99755 1.03826 0.377112
8.18971 2.62583 1.62962 0.376657
9.11962 0.370993 1.10319 0.376203
3.87484 0.886297 1.62964 0.375748
10.1918 2.16147 1.05822 0.375293
2.04221 0.510807 1.63223 0.374838
6.77811 2.63218 1.63574 0.374384
10.8692 2.62967 0.814603 0.373929
4.65212 1.01356 1.63033 0.373474
2.79203 2.51304 1.62726 0.37302
2.63144 0.969922 1.63464 0.372565
1.78873 2.04877 1.63599 0.37211
6.83048 2.14542 1.08639 0.371655
13.213 0.367339 0.368709 0.371201
14.5824 0.37139 0.787223 0.370746
5.64087 2.16887 0.36137 0.370291
7.12869 1.0213 1.6355 0.369837
14.0406 2.6335 1.04398 0.369382
0.714409 0.363826 0.50124 0.368927
9.20937 1.48469 1.19176 0.368472
12.4245 2.1045 0.419074 0.368018
5.58008 1.36106 0.909748 0.367563
13.3667 1.01011 0.927192 0.367108
10.5069 1.287 0.361293 0.366654
4.00331 2.59383 1.45601 0.366199
5.48906 2.63662 0.906177 0.365744
6.20162 2.63955 0.36943 0.365289
0.790989 1.18798 1.6365 0.364835
2.64419 0.35583 1.23654 0.36438
3.58558 1.67164 1.23625 0.363925
0.957268 2.04523 0.361927 0.363471
9.99456 1.43894 0.970477 0.363016
9.74545 2.11684 1.64035 0.362561
3.28906 2.64795 1.08372 0.362106
1.20546 1.79417 1.00977 0.361652
10.6521 1.31384 1.61748 0.361197
3.81279 1.75986 0.356125 0.360742
14.2333 1.00396 0.896231 0.360288
6.71662 2.13457 0.357108 0.359833
13.9372 1.50976 0.478639 0.359378
11.8072 2.64605 1.64057 0.358923
4.35845 0.347838 1.65181 0.358469
2.43262 1.29518 1.01741 0.358014
6.75198 1.00369 0.351857 0.357559
1.3162 2.6443 1.64482 0.357104
7.03494 0.352121 0.871095 0.35665
6.48933 0.349325 0.409623 0.356195
9.7182 1.21899 0.351018 0.35574
5.30543 2.02087 1.64531 0.355286
7.67859 0.811952 0.353809 0.354831
10.6797 0.353211 1.64982 0.354376
6.47074 1.88061 1.65284 0.353921
7.48337 2.46228 1.65196 0.353467
2.20884 0.350594 0.345639 0.352557
6.17986 2.65055 1.08631 0.352103
5.33205 1.01719 0.337163 0.351648
0.347491 0.977319 0.575552 0.351193
12.7419 0.346244 0.910205 0.350738
13.74 0.851376 0.343585 0.350284
3.88875 1.0549 0.338712 0.349829
12.3563 1.54236 0.855447 0.349374
11.181 1.49151 1.00949 0.34892
8.525 0.963934 0.343474 0.348465
4.91226 2.24803 1.05235 0.34801
9.2872 2.14087 0.339746 0.347555
11.1031 2.65752 1.66756 0.347101
4.90575 1.58755 0.827749 0.346646
8.65368 0.345717 1.65422 0.346191
14.6552 1.29338 0.344222 0.345737
9.96344 1.44852 1.6666 0.345282
11.0859 1.8674 1.65667 0.344827
1.98675 1.28071 1.65687 0.344372
0.647441 2.66083 0.875073 0.343918
8.42044 0.340322 0.982853 0.343463
1.42828 0.881684 1.6361 0.343008
2.80367 1.89813 1.26538 0.342554
12.3154 0.849697 1.12158 0.342099
8.19034 0.341042 0.339692 0.341644
12.2979 1.65381 1.66164 0.341189
8.19094 2.6603 0.819432 0.340735
2.8975 0.328941 0.335234 0.34028
2.6427 2.12067 0.329779 0.339825
3.54995 0.334036 1.32275 0.339371
10.4187 2.66273 1.66316 0.338916
9.07569 1.46894 0.332291 0.338461
0.8713 0.491163 1.66572 0.338006
4.65254 0.33509 1.0319 0.337552
8.69036 2.67301 0.336212 0.337097
13.9336 2.6648 0.335043 0.336642
4.72342 1.70027 1.4729 0.336188
0.329472 1.71014 1.66881 0.335733
7.39052 1.82946 1.40696 0.335278
2.10716 2.67226 1.67366 0.334823
6.06974 0.891341 0.330989 0.334369
0.389165 0.343503 1.223 0.333914
5.51859 2.67161 1.66688 0.333459
8.30609 0.928985 1.66849 0.333004
0.32886 2.13423 1.14806 0.33255
12.6582 2.21801 1.67285 0.332095
9.80468 0.319795 1.23705 0.33164
13.2708 2.6693 0.316649 0.331186
8.22128 1.59056 1.31321 0.330731
6.46458 0.812166 1.67456 0.330276
7.11352 0.326307 1.67086 0.329821
9.75103 0.812421 1.6714 0.329367
8.62285 1.42865 0.819318 0.328912
7.74327 1.33984 1.67745 0.328457
9.19105 0.960823 0.737736 0.328003
12.0524 2.67802 0.318205 0.327548
7.54879 2.68578 0.322645 0.327093
12.663 0.322553 1.67627 0.326638
11.8063 1.56241 1.22506 0.326184
1.83948 2.67896 1.07566 0.325729
14.686 2.67605 1.26315 0.325274
10.7492 2.12541 0.322575 0.32482
13.2684 1.32613 0.32079 0.324365
4.70095 2.32372 1.68033 0.32391
11.7816 2.68126 0.899863 0.323455
10.4064 1.94776 1.68425 0.323001
13.108 2.68355 1.69026 0.322546
14.201 0.995357 1.68007 0.322091
13.5944 1.26418 1.68123 0.321637
8.73098 1.42463 1.68157 0.321182
4.85268 2.09859 0.404876 0.320727
1.49744 2.32769 0.660805 0.320272
4.89517 0.942658 0.98739 0.319818
7.91815 2.16692 0.318002 0.319363
12.9988 1.6685 1.69269 0.318908
2.32678 2.69199 0.312255 0.318454
4.1689 1.97444 1.68541 0.317999
1.11506 2.69548 0.314455 0.317544
14.269 2.06367 0.311693 0.317089
0.313761 2.68379 0.311768 0.316635
14.2771 1.84636 0.95384 0.31618
13.0742 1.93386 0.315021 0.315725
11.7321 1.9839 1.69189 0.315271
2.01953 0.311442 0.9768 0.314816
10.1443 0.722061 0.459971 0.314361
1.36449 1.52832 1.68906 0.313906
8.76167 0.946544 1.21379 0.313452
0.639322 2.275 1.69428 0.312997
11.9486 0.302463 0.301857 0.312542
0.29742 2.0555 0.310955 0.312087
4.25078 2.26417 0.312124 0.311633
4.69746 2.69601 0.30727 0.311178
5.83445 1.51879 0.306405 0.310723
1.05853 0.300028 1.07962 0.310269
12.4104 0.877825 0.308238 0.309814
5.86935 1.69866 1.41006 0.309359
10.7554 0.692136 0.562511 0.308904
9.09782 2.54183 0.816758 0.30845
11.4086 0.30704 0.870314 0.307995
2.95039 2.69321 0.30694 0.30754
4.18838 1.3938 1.30473 0.307086
6.52443 0.763067 0.924632 0.306631
14.6994 0.304312 1.44701 0.306176
9.1085 1.92437 1.69478 0.305721
9.38838 0.298656 1.7116 0.305267
12.6538 1.44027 0.29296 0.304812
12.515 1.05195 1.70113 0.304357
11.9465 0.321634 1.15372 0.303903
0.293469 0.760633 1.69863 0.303448
7.48947 1.27856 1.10281 0.302993
6.57849 0.296307 1.33772 0.302538
12.8611 1.34445 1.19468 0.302084
5.84507 0.299185 0.299097 0.301629
7.20105 2.70757 1.12101 0.301174
10.8213 0.297506 1.01803 0.30072
6.09974 0.2994 1.70331 0.300265
2.71169 2.40828 0.888954 0.29981
0.29012 1.10954 1.20974 0.299355
14.708 1.31103 1.52707 0.298901
7.71549 0.716603 1.7025 0.298446
11.353 0.479152 0.294812 0.297991
9.59046 0.295024 0.295473 0.297537
12.9083 0.78297 1.34578 0.297082
11.9542 0.843539 1.70656 0.296627
13.6598 2.09701 0.294861 0.296172
12.5557 0.292847 0.294513 0.295718
3.34912 1.29032 1.71336 0.295263
2.63009 0.890659 0.283777 0.294808
5.24439 1.64699 0.289597 0.294354
9.08565 0.814069 1.70771 0.293899
8.68296 2.72222 1.1996 0.293444
5.44578 0.782418 1.18457 0.292989
5.91826 2.19769 1.72223 0.292535
13.7192 2.71408 1.70856 0.29208
3.41634 2.71115 1.71663 0.291625
2.28127 2.24466 1.25376 0.291171
3.23756 0.717091 1.71079 0.290716
9.63785 2.71521 0.656104 0.290261
2.0197 2.16957 0.285975 0.289806
5.29333 2.71704 0.28729 0.289352
6.13348 2.72448 1.71696 0.288897
0.756441 0.80873 1.03673 0.288442
8.04585 0.286107 1.48624 0.287987
12.6601 2.71343 0.286125 0.287533
6.72691 2.71668 0.780344 0.287078
9.82249 2.69215 1.19424 0.286623
14.7141 1.43179 0.956626 0.286169
0.609331 1.58737 1.13433 0.285714
13.4045 1.68271 0.79077 0.285259
4.31348 2.66478 0.75802 0.284804
5.95034 0.960022 1.38167 0.28435
3.57704 2.14865 0.809794 0.283895
13.8145 1.445 1.10462 0.28344
12.0593 0.278013 1.72274 0.282986
7.64522 2.39543 0.925748 0.282531
1.73019 2.72198 0.278117 0.282076
3.02194 1.32055 1.24673 0.281621
11.2929 1.28073 1.59186 0.281167
10.9617 0.822281 1.22228 0.280712
1.49472 0.66675 0.279804 0.280257
6.94792 2.72336 0.268334 0.279803
1.44251 0.271717 1.72542 0.279348
12.4525 2.62566 0.804868 0.278893
11.9877 1.0873 0.654146 0.278438
7.92233 2.11758 1.32315 0.277984
7.29996 1.30706 0.272246 0.277529
5.4572 1.98327 0.947553 0.277074
13.9205 0.265808 0.273461 0.27662
1.24441 2.72753 1.0213 0.276165
2.40081 1.57089 0.449279 0.27571
7.13401 0.507679 0.270783 0.275255
8.11589 1.94999 0.834476 0.274801
0.272211 2.72829 1.72584 0.274346
4.31895 0.869411 1.09813 0.273891
14.656 2.29689 1.7266 0.273437
1.52121 0.270705 0.667924 0.272982
0.987436 0.88354 0.265902 0.272527
10.7572 2.27737 1.33822 0.272072
6.16111 1.32644 1.73546 0.271618
4.58148 2.73293 1.23878 0.271163
4.18313 1.42199 0.734915 0.270708
14.113 2.32039 1.73254 0.270254
13.7596 2.09839 0.850767 0.269799
9.8235 2.73426 1.74052 0.269344
14.3428 0.7726 0.268191 0.268889
11.2794 2.2912 1.21013 0.268435
4.43152 1.71447 0.266829 0.26798
14.7337 1.84655 0.591696 0.267525
2.92959 0.706208 0.805452 0.267071
6.39974 1.60524 0.258511 0.266616
3.18285 1.8253 1.73513 0.266161
3.62519 2.17492 1.73784 0.265706
0.265268 1.6346 0.706202 0.265252
11.5986 2.00434 0.904155 0.264797
2.58796 1.59433 1.74038 0.264342
14.7415 0.790338 1.73783 0.263887
3.19117 0.255315 0.853036 0.263433
14.7391 1.86238 1.4233 0.262978
12.7067 2.7377 1.26932 0.262523
9.50596 1.69368 0.682649 0.262069
0.977715 2.31964 1.24215 0.261614
12.7642 1.00075 0.747861 0.261159
9.73279 0.76849 0.857986 0.260704
9.27343 2.7416 0.253816 0.26025
9.39227 1.27142 1.74708 0.259795
10.6036 1.33185 0.97355 0.25934
12.4071 2.74476 1.74187 0.258886
3.86219 1.49939 1.75444 0.258431
7.09704 0.962783 0.858321 0.257976
10.6899 1.7758 1.21279 0.257521
2.01314 0.864088 1.11336 0.257067
11.7743 0.834867 0.245609 0.256612
11.1851 1.75892 0.252961 0.256157
11.0229 0.84639 1.75115 0.255703
10.0865 0.254949 1.74621 0.255248
14.7468 2.16405 1.0037 0.254793
9.362 2.0872 1.1602 0.254338
8.11575 2.74873 0.2439 0.253884
2.96703 0.250061 1.74941 0.253429
7.16278 2.42743 0.650653 0.252974
10.609 0.251512 0.249264 0.25252
8.28658 1.97278 1.74806 0.252065
11.7913 1.36132 1.75775 0.25161
2.42116 2.75888 1.18965 0.251155
0.70927 1.32458 0.245715 0.250701
1.61621 1.3862 1.20558 0.250246
11.0549 2.03596 0.800258 0.249791
10.5555 1.77503 0.729762 0.249337
10.3205 0.823625 1.7569 0.248882
14.2495 0.242249 1.75468 0.248427
6.82519 1.56879 0.536288 0.247972
1.92401 1.63146 0.2447 0.247518
13.768 0.835434 1.35457 0.247063
2.18824 2.37489 0.747252 0.246608
2.41803 0.696803 0.786405 0.246154
6.99396 1.60903 1.75896 0.245699
8.68439 2.27064 1.7549 0.245244
3.99515 0.578111 0.829651 0.244335
10.3372 2.75988 1.09527 0.24388
2.25167 1.73018 1.38834 0.243425
14.238 0.568283 1.3132 0.24297
13.0595 1.09712 1.75755 0.242516
3.33459 1.23869 0.241924 0.242061
3.28798 2.1324 0.234973 0.241606
7.41444 0.745112 1.17484 0.241152
9.35796 0.772332 0.231484 0.240697
0.238684 0.231308 1.76325 0.240242
1.17368 2.04878 1.74824 0.239787
6.32381 1.29297 0.645087 0.239333
10.1186 0.228557 0.234112 0.238878
9.94065 0.237617 0.692382 0.238423
14.3288 2.76967 1.68642 0.237969
1.53666 1.98298 0.229957 0.237514
7.02763 2.08991 1.7637 0.237059
14.7745 0.526154 0.234125 0.236604
5.07714 2.76496 1.32335 0.23615
5.82303 2.20824 1.21055 0.235695
9.65077 2.20755 0.786206 0.23524
9.36656 0.907437 1.26971 0.234786
4.08412 0.231171 1.14187 0.234331
4.91464 2.76939 0.800117 0.233876
10.409 0.232857 0.689301 0.233421
1.23297 0.230409 0.231538 0.232967
12.2159 2.77107 1.23377 0.232512
13.3443 1.55291 1.2838 0.232057
14.4208 0.227064 0.228042 0.231603
4.92558 0.603897 0.226628 0.231148
12.944 0.886472 0.229946 0.230693
1.46219 2.29212 1.20149 0.230238
3.8399 2.77464 0.916925 0.229784
10.8323 2.77072 0.229273 0.229329
0.234816 1.45209 0.247981 0.228874
0.751681 1.5741 0.646898 0.22842
7.81227 1.90448 1.77225 0.227965
11.3241 2.78204 1.16491 0.22751
2.80187 1.37576 0.221474 0.227055
4.98306 2.77986 1.77574 0.226601
1.56088 0.224058 1.24022 0.226146
14.7789 2.77677 1.78643 0.225691
6.05596 1.00883 0.888893 0.225237
0.887072 1.1443 0.675171 0.224782
7.27665 2.21975 0.221066 0.224327
3.19626 0.794248 0.223841 0.223872
3.28909 2.15317 1.389 0.223418
5.68946 0.816126 0.727648 0.222963
0.991612 1.25886 1.09234 0.222508
2.85254 1.15076 0.655179 0.222054
4.183 2.78273 0.219276 0.221599
0.218726 0.220445 0.219771 0.221144
0.22161 2.77266 1.23101 0.220689
0.773054 2.78194 1.7806 0.21978
4.44518 2.78645 1.78362 0.219325
14.7863 2.14386 0.213942 0.21887
6.64346 1.31592 1.78366 0.218416
3.1466 2.32475 0.624784 0.217961
10.981 0.216346 0.533016 0.217506
4.3343 1.47996 1.79386 0.217052
12.2265 0.212115 0.734346 0.216597
3.81891 0.214406 1.78885 0.216142
5.65113 0.769227 1.7701 0.215687
12.9455 1.50668 0.713601 0.215233
11.7996 0.633288 0.669203 0.214778
0.205911 0.211302 0.726681 0.214323
1.67768 0.212982 0.211898 0.213869
13.9296 0.56112 1.79464 0.213414

Jabrane

>
>
> Yes that's it: "I'm trying to determine normal/axial strain of the beam in
> 3-points bending test."
>

ok, then do NOT use UniaxialStrainer, as it is intended to *apply*
something like uniaxial stress (not even strain) boundary conditions :-) to
compute strain from particle displacement, TessalationWrapper might be a
choice, but I have no experience with it, but probably somebody else will
give some more hints :-)
cheers
Jan

Hi Jan

I used the particle displacement to compute the strain.
I choose some particles and i determine their displacements with O.state.pos
I will see if TessalationWrapper can be useful in my case.

cheers
Jabrane