This paper presents a general 2.5D meshless MLPG methodology for the computation of the elastic response of longitudinally invariant structure subjected to the incident wave field. A numerical frequency domain model is established using the Fourier transform in time and longitudinal coordinate domains. This allows for significant reduction of computational effort required. In the MLPG method the Moving-Least Squares (MLS) scheme is employed for the approximation of the spatial variation of displacement field. No finite elements are required for the approximation or integration of unknowns. A small circular subdomain is introduced around each nodal point in the analyzed domain. Local integral equations derived from the governing equations are specified on these subdomains. Continuously non-homogeneous material properties are varying in the cross-section of the analyzed structure. A simple patch test is introduced to assess the accuracy and the convergence of developed numerical model. At the end of the paper, numerical examples illustrate the applicability of the proposed numerical formulation.