Monte-Carlo Simulation of Particle Reinforced Composites Based on Hybrid Stress Elements

WANG Wei; GUO Ran

doi:10.21656/1000-0887.420016

Volume 42 Issue 8

Aug. 2021

Turn off MathJax

Article Contents

Article Navigation > Applied Mathematics and Mechanics > 2021 > 42(8): 794-802

WANG Wei, GUO Ran. Monte-Carlo Simulation of Particle Reinforced Composites Based on Hybrid Stress Elements[J]. Applied Mathematics and Mechanics, 2021, 42(8): 794-802. doi: 10.21656/1000-0887.420016

Citation:

WANG Wei, GUO Ran. Monte-Carlo Simulation of Particle Reinforced Composites Based on Hybrid Stress Elements[J]. Applied Mathematics and Mechanics, 2021, 42(8): 794-802. doi: 10.21656/1000-0887.420016

Citation:

WANG Wei, GUO Ran. Monte-Carlo Simulation of Particle Reinforced Composites Based on Hybrid Stress Elements[J]. Applied Mathematics and Mechanics, 2021, 42(8): 794-802. doi: 10.21656/1000-0887.420016

PDF( 3785 KB)

Monte-Carlo Simulation of Particle Reinforced Composites Based on Hybrid Stress Elements

doi: 10.21656/1000-0887.420016

WANG Wei,
GUO Ran^,

Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650000, P.R.China

Funds:

The National Natural Science Foundation of China(11572142)

Received Date: 2021-01-15
Rev Recd Date: 2021-03-04

Available Online: 2021-08-14

Abstract

Abstract

Based on the assumed high-order stress field, the hybrid stress finite element method has higher calculation accuracy with sparse grids. The quadtree meshes were used to discretize heterogeneous computing domains with advantages of the displacement coordination conditions for hanging nodes automatically satisfied. Moreover, all quadtree elements can be divided into a limited number of types, and the stiffness matrices of these elements can be pre-computed and stored in the memory, retrieved and scaled as required during computations, which greatly improves the computational efficiency. In view of the randomness of inclusions, the effects of the volume ratio, the number and the aspect ratio of random inclusions on the homogeneous equivalent modulus of the composite were discussed with the Monte-Carlo method and the homogenization method. The results show that, the equivalent elastic modulus of the composite increases with the volume ratio, the number and the aspect ratio of inclusions, and is most sensitive to the volume ratio.
- hybrid stress element method,
- quadtree mesh,
- Monte-Carlo simulation,
- particle reinforced composite

FullText(HTML)

References(18)

References

[2]KAMINSKI M, LESNIAK M. Homogenization of metallic fiber-reinforced composites under stochastic ageing[J].Composite Structures,2012,94(2): 386-393.

ESHELBY J D. The determination of the elastic field of an ellipsoidal inclusion, and related problems[J].Proceedings of the Royal Society of London(Series A): Mathematical and Physical Sciences,1957,241(1226): 376-396.

[3]SAKATA S I, ASHIDA F, IWAHASHI D. Stochastic homogenization analysis of a particle reinforced composite material using an approximate Monte-Carlo simulation with the weighted least square method[J].Journal of Computational Science and Technology,2013,7(1): 1-11.

[4]LEE S P, JIN J W, KANG K W. Probabilistic analysis for mechanical properties of glass/epoxy composites using homogenization method and Monte-Carlo simulation[J].Renewable Energy,2014,65: 219-226.

[5]PIAN T H H. Derivation of element stiffness matrices by assumed stress distributions[J].AIAA Journal,1964,2(7): 1333-1336.

[6]PIAN T H H, TONG P. Relations between incompatible model and hybrid stress model[J].International Journal for Numerical Methods in Engineering,1986,22(1): 173-181.

[7]杨锋, 郭然. 多边形应力杂交单元的接触算法研究[J]. 应用数学和力学, 2019,40(10): 1059-1070.(YANG Feng, GUO Ran. Study on contact algorithms for the polygonal hybrid stress element method[J].Applied Mathematics and Mechanics,2019,40(10): 1059-1070.(in Chinese))

[8]GHOSH S, MUKHOPADHYAY S N. A material based finite element analysis of heterogeneous media involving Dirichlet tessellations[J].Computer Methods in Applied Mechanics and Engineering,1993,104(2): 211-247.

[9]GHOSH S, MOORTHY S. Elastic-plastic analysis of arbitrary heterogeneous materials with the Voronoi cell finite element method[J].Computer Methods in Applied Mechanics and Engineering,1995,121(1/4): 373-409.

[10]GUO R, SHI H, YAO Z. Numerical simulation of thermo-mechanical fatigue properties for particulate reinforced composites[J].Acta Mechanica Sinica,2005,21(2): 160-168.

[11] LI H,GUO R,CHENG H M. Calculation of stress intensity factors of matrix crack tip in particle reinforced composites using the singular Voronoi cell finite element method[J].Theoretical and Applied Fracture Mechanics,2019,101: 269-278.

[12] ZHANG R, WANG T, GUO R. Modeling of interphases in multiple heterogeneities reinforced composites using Voronoi cell finite elements[J]. Acta Mechanica Sinica,2020,36(4): 887-901.

[13]郭然. 颗粒增强复合材料界面脱层和热机疲劳的数值模拟[D]. 博士学位论文. 北京: 清华大学, 2003.(GUO Ran. Modeling of interfacial debonding and characterization of fatigue in particle reinforced composites[D]. PhD Thesis. Beijing: Tsinghua University, 2003.(in Chinese))

[14]YANG Q S, QIN Q H. Micro-mechanical analysis of composite materials by BEM[J].Engineering Analysis With Boundary Elements,2004,28(8): 919-926.

[15]SUKUMAR N, TABARRAEI A. Conforming polygonal finite elements[J]. International Journal for Numerical Methods in Engineering,2004,61(12): 2045-2066.

[16]TABARRAEI A, SUKUMAR N. Extended finite element method on polygonal and quadtree meshes[J].Computer Methods in Applied Mechanics and Engineering,2008,197(5): 425-438.

[17]OOI E T, NATARAJAN S, SONG C, et al. Crack propagation modelling in concrete using the scaled boundary finite element method with hybrid polygon-quadtree meshes[J].International Journal of Fracture,2017,203(1): 135-157.

[18]WENG G J. Some elastic properties of reinforced solids, with special reference to isotropic ones containing spherical inclusions[J]. International Journal of Engineering Science,1984,22(7): 845-856.

Relative Articles

Supplements(0)

Cited By

Proportional views