Texture Fiber HCP Prismatic and Pyramidal

151 views
Skip to first unread message

Camilo Garzón

unread,
Sep 21, 2019, 1:32:59 PM9/21/19
to MTEX
Hello,

I've been trying to make prismatic and pyramidal fiber of an hcp system without success, next, my procedure:

First I try with the basal fiber doing this:

1. Define the symmetry and an orientation of a basal pole
csa = crystalsymmetry('hexagonal');
o_basal = orientation.byEuler(0,0,0,csa);


2. Define the fiber with the orientation and a vector 3d to align every direction parallel to the normal to the orientation defined.
f_bsl = fibre(o_basal,vector3d.Z)

---
 size: 1 x 1
mineral: Titanium (Alpha) (622, X||a, Y||b*, Z||c)
o1: (0°,0°,0°)
h: (0001)

---
plotPDF(f_bsl,h,csa,'lineColor','r','linewidth',2)

Luckily, I built the basal fiber right, as you can see in the plot attached. However, I couldn't get the same results for prismatic and pyramidal ones. 

For that reason, I try to build the prismatic and pyramidal fibers, with two orientations.

For prismatic I take as reference the orientation (1 1 -2 0) (1 -1 0 0),

o_pr = orientation.byMiller([1 1 0],[1 -1 0],csa);

--
 Bunge Euler angles in degree
  phi1  Phi phi2 Inv.
   360   90   30    0
 
The idea is to complete the fiber since phi1 = 0*degree to phi2 = 360*degree

o_pr1 = orientation.byEuler(0*degree,90*degree,30*degree,csa);
f_pr = fibre(o_pr1,o_pr);

---
 o1: (0°,90°,30°)
 h: (2-1-10)
>> 

It supposed that it must did the prismatic fiber, but the result is other as you can see in the attached pole figure.



How can I solve it? And also I would like that the point in the center of pole figure appear when it supossed to be there?

Best Regards

Camilo Garzón
basal fiber.png
failed prismatic fiber.png

Rüdiger Kilian

unread,
Sep 22, 2019, 5:56:53 AM9/22/19
to mtex...@googlegroups.com
Hi,
how about

fp = fibre(Miller(1,0,-1,1,cs),zvector)

?
That should be a fibre with (1,0,-1,1) as axis parallel to 001 in specimen coordiantes.

Hope that helps.
Cheers,
Rüdiger
________________________________________
From: mtex...@googlegroups.com <mtex...@googlegroups.com> on behalf of Camilo Garzón <scga...@gmail.com>
Sent: Saturday, September 21, 2019 7:32:59 PM
To: MTEX
Subject: {MTEX} Texture Fiber HCP Prismatic and Pyramidal
--
If you want to reduce the number of emails you get through this forum login to https://groups.google.com/forum/?fromgroups=#!forum/mtexmail, click "My membership" and select "Don't send me email updates". You can still get emails on selected topics by staring them.
---
You received this message because you are subscribed to the Google Groups "MTEX" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mtexmail+u...@googlegroups.com<mailto:mtexmail+u...@googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/mtexmail/7306b354-9fe0-422d-93dc-d7e3381eecce%40googlegroups.com<https://groups.google.com/d/msgid/mtexmail/7306b354-9fe0-422d-93dc-d7e3381eecce%40googlegroups.com?utm_medium=email&utm_source=footer>.

Camilo Garzón

unread,
Sep 23, 2019, 12:39:53 PM9/23/19
to MTEX
Hi Rüdiger,

Thanks for your answer, I tried how you tell me, however, I didn't get the expected results.

f_prismatic = fibre(Miller(1, 1, -2, 0,csa),zvector))
--

 size
: 1 x 1
 mineral
: Titanium (Alpha) (622, X||a, Y||b*, Z||c)

 o1
: (330°,90°,30°)
 h
: (11-20)


f_pyramidal
= fibre(Miller(1,0,0,0,csa),zvector)


--

size
: 1 x 1
 mineral
: Titanium (Alpha) (622, X||a, Y||b*, Z||c)

 o1
: (300°,90°,60°)
 h
: (10-10)

When I plotted, I expected different results in the PDF of what I get, In the attached image the I put the pole figures for prismatic and pyramidal fibre according to with the next article.
Texture Analysis in hexagonal materials /DOI: 
  • 10.1016/S0254-0584(03)00168-8

  • What do you think?

To unsubscribe from this group and stop receiving emails from it, send an email to mtex...@googlegroups.com<mailto:mtexmail+unsubscribe@googlegroups.com>.
Expected vs Obtanied fibre.png

ruediger Kilian

unread,
Sep 28, 2019, 12:12:43 PM9/28/19
to mtex...@googlegroups.com
Hi,
if you look at the pole figures for the {11-20} fibre, you see that (10-10) makes the small circle ring (such as (01-10) would do) while you get the equatorial great circle for (1-100) and (-1100). So the fibre is fine so far but I guess you need to superpose those pole figures to see the same as in that publication for now.

Hope that helps.

Cheers,
Rüdiger

cs = crystalSymmetry('622','X||a')
f = fibre(Miller(1,1,-2,0,cs),zvector)
h = Miller(1,0,-1,0,cs)
plotPDF(f,h.symmetrise,'linewidth',2,'linecolor','r')



João Rodrigues

unread,
Sep 30, 2019, 11:17:41 AM9/30/19
to MTEX
Hi Camilo,

You can calculate the odf and then plot pole figure.

f1 = fibre(Miller(1,1,-2,0,cshex),zvector);
odf = fibreODF(f1,'halfwidth',5*degree);

plotPDF(odf,Miller(-2,1,1,0,cshex));

Best regards,

Camilo Garzón

unread,
Sep 30, 2019, 11:55:21 AM9/30/19
to MTEX
Hi

Thanks, Rüediger and João for your help, the solution of João, confirms clearly what you said Rüediger.
 
Best Regards
Camilo
Prismatic Fibre.PNG
Reply all
Reply to author
Forward
0 new messages