This　paper　adopts　a　meshless　method　that　combines　the　moving　least-squares　approximation　and　Galerkin　weak　form　to　investigate　the　mechanical　properties　of　unidirectional　fiber-reinforced　composites　under　a　lateral　load.　The　degree　of　nonuniformity　is　adopted　to　quantify　the　spatial　distribution　of　randomness　during　simulation,　and　then　the　numerical　implementation　is　developed　to　generate　a　representative　volume　element　(RVE)　model　with　random　distribution　of　fibers.　A　statistical　analysis　is　carried　out　by　using　three　descriptors,　that　is,　inter-fiber　distance,　second-order　intensity　function　and　radial　distribution　function.　Numerical　examples　are　presented　to　illustrate　the　accuracy　of　the　proposed　multiscale　meshless　method,　and　excellent　agreement　is　achieved　in　comparison　with　experiments　and　the　finite　element　method.　Finally,　the　performance　of　the　proposed　numerical　technique　is　evaluated　by　considering　the　effects　of　RVE　size,　node　influence　domain,　degree　of　nonuniformity　and　fiber　volume　fraction　separately.