I am currently model a unidirectional fibre reinforced composite in 3D RVE in ANSYS. After solve the model, where can I read the Young modulus E1,E2 E3, poison ratio v12,v13,v23, and shear modulus G12, G13, G23 of the whole composite, as I modeled the 2 materials separately.
I believe they are not available readily in a click of a button and there is a long procedure. If you have access to this book http://barbero.cadec-online.com/feacm-ansys/index.html , see it's chapter 6. It describes in detail how to extract the material properties of a composite from a laminated RVE analysis. Also, I think you can find it in other composite textbooks too.
I fully agree with Singh on the book of Barbero. I personally have it and I learnt a lot from it.
Prof. Berbero posted on line two very interesting APDL macros for numerical homogenisation in ANSYS:
- periodic boundary conditions: http://barbero.cadec-online.com/feacm-ansys/Examples/Ch6ex/FEAcomp_Ex603/ceRVE.mac
- stress recovery (the homogenised properties Cij): http://barbero.cadec-online.com/feacm-ansys/Examples/Ch6ex/FEAcomp_Ex603/sRecover.mac
Once you have the homogenised Cij, you should find the inverse Sij (compliance matrix) and from the definition of the compliance matrix you can obtain what you are looking for. The following APDL snippet might be of your interest (check it carefully, though):
!
! COMPLIANCE MATRIX
det_C = C11*(C22*C33-C23*C32)-C12*(C21*C33-C23*C31)+C13*(C21*C32-C22*C31)
S11 = +(C22*C33-C23*C32)/det_C
S12 = -(C12*C33-C13*C32)/det_C
S13 = +(C12*C23-C22*C13)/det_C
!
S21 = -(C21*C33-C31*C23)/det_C
S22 = +(C11*C33-C31*C13)/det_C
S23 = -(C11*C23-C21*C13)/det_C
!
S31 = +(C21*C32-C22*C31)/det_C
S32 = -(C11*C32-C31*C12)/det_C
S33 = +(C11*C22-C21*C12)/det_C
!
S44 = 1./C44
S55 = 1./C55
S66 = 1./C66
!
! ENGINEERING CONSTANTS
E11 = 1./S11
E22 = 1./S22
E33 = 1./S33
!
G12 = C66
G13 = C55
G23 = C44
!
n21 = -S12*E22
n12 = -S21*E11
n31 = -S13*E33
n13 = -S31*E11
n32 = -S23*E33
n23 = -S32*E22
!
*cfopen, 00-res, out, , append
*vwrite, C11, C12, C13, C14, C15, C16
(5(es13.6, 1x), es13.6)
*vwrite, C21, C22, C23, C24, C25, C26
(5(es13.6, 1x), es13.6)
*vwrite, C31, C32, C33, C34, C35, C36
(5(es13.6, 1x), es13.6)
*vwrite, C41, C42, C43, C44, C45, C46
(5(es13.6, 1x), es13.6)
*vwrite, C51, C52, C53, C54, C55, C56
(5(es13.6, 1x), es13.6)
*vwrite, C61, C62, C63, C64, C65, C66
(5(es13.6, 1x), es13.6)
*vwrite,
(" ")
*vwrite, E11
("E11 = ", es13.6, 1x)
*vwrite, E22
("E22 = ", es13.6, 1x)
*vwrite, E33
("E33 = ", es13.6, 1x)
*vwrite, G12
("G12 = ", es13.6, 1x)
*vwrite, G23
("G23 = ", es13.6, 1x)
*vwrite, G13
("G13 = ", es13.6, 1x)
*vwrite, n12
("n12 = ", es13.6, 1x)
*vwrite, n13
("n13 = ", es13.6, 1x)
*vwrite, n23
("n23 = ", es13.6, 1x)
*vwrite,
(/,/)
*cfclos
Hi, please check out this link. It might be helpful.
https://www.google.ca/url?sa=t&rct=j&q=&esrc=s&source=web&cd=2&cad=rja&uact=8&ved=0CCQQFjABahUKEwigz4jRp_XGAhUHfJIKHe-GDnk&url=http%3A%2F%2Fwww-mdp.eng.cam.ac.uk%2Fweb%2Flibrary%2Fenginfo%2Fcueddatabooks%2Fmaterials.pdf&ei=Y_6yVeARh_jJBO-NusgH&usg=AFQjCNFxI28LxP6wY5AbmdmPQb6tQEzPlA&sig2=HkliYxyHxlswd1SOiuzScw&bvm=bv.98717601,d.aWw
Hi Mr Gaetano Giunta,
I have a problem to open the second macro from your links when i try to do the tutorial.
I have followed the file name and saved it as .mac in my directory file. I Can't find any solution to figure it.
Hope you can help.
Thank you
Hi,
try the script in the attached file.
It works fine for me.
Best,
Gaetano
Dear all,
I would like to model a 3D composite repeating unit cell consisting of piezoelectric phases in ANSYS. Are there any related tutorials or videos which teach detailed and comprehensive model analysis of piezoelectric composite repeating unit cell? (I think the key is how to apply periodic bcs for the RUC in ansys).
Thank you
- Badges
- Science topic