Journal of Applied and Computational Mechanics
Volume:8 Issue: 2, Spring 2022

Laminated composite materials have the advantage of desired properties and are vastly replacing the existing traditional materials in Civil Engineering construction. In the present study, it is aimed to extend the study on the analysis of laminated composites plates and shells towards structural optimization. Topology optimization is performed using two different objective functions namely strain energy and fundamental frequency. The results of optimization have shown clearly that the distribution of material is dependent on the laminae. The optimal arrangement of material is obtained after using a cut-off relative density. It is confirmed to be a well-connected grid and is examined in detail. The results have shown that the optimal arrangement of material for a simply supported plate carrying a uniformly distributed load is at the centre of the edges and not towards the corners, Hence, the optimal arrangement of beams using strain energy is to align by joining the centre of the edges similar to plus (+) sign.

An analytical model is developed to study the effects of viscous dissipation and chemical reaction in viscoelastic convection from an inclined plate as a simulation of electro-conductive polymer materials processing. The Jeffery’s viscoelastic model is deployed to describe the non-Newtonian characteristics of the fluid and provides a good approximation for polymers, which constitutes a novelty of the present work. The normalized nonlinear boundary value problem is solved computationally with the Keller-Box implicit finite-difference technique. Extensive solutions for velocity, surface temperature and concentration, skin friction, heat, and transfer rates are visualized numerically and graphically for various thermophysical parameters. Validation is conducted with earlier published work for the case of a vertical plate in the absence of viscous dissipation, chemical reaction, and non-Newtonian effects. The boundary layer flow is accelerated with increasing Deborah number whereas temperatures and concentrations are decelerated slightly. Temperatures and concentration are boosted with increasing inclination parameter whereas velocity is lowered. A reverse trend is seen for increasing Richardson number. Increasing chemical reaction reduces velocity and concentration whereas it enhances temperature. Increasing the viscous dissipation parameter is found to enhance velocity and temperature whereas it suppresses concentration.

Heating capacity of an earth to air heat exchanger EAHE equipped with an exterior fan in the arid region like the Southwest of Algeria is investigated experimentally. In-situ measurement of annual undisturbed subsoil vertical temperature profile is shown that it was 28°C at a depth of 1.5 meters. The EAHE made of 66 meters of PVC tube is demonstrated a heating capacity of 13°C and a cooling capacity of 7°C and a big dependence on local climate conditions. Great potentials and thermal comfort with less energy consumption are represented by earth to air or air-ground heat exchanger in the arid regions like the South of Algeria.

Horn is one of the main components of ultrasonic cutting systems. The most important characteristics of the horn design are its resonant frequency and amplification factor. Closed-form equations can be used only for the design of simple horns and do not apply to more complex shapes like surgical tools. In This paper, a designing technique based on the finite element method and experimental tests is presented. The conventional design methods are improved, and designing a high performance surgical ultrasonic horn for bone cutting tools is facilitated. The new and complex bone cutting tool has both the knife-edge and toothed-edge, which could cut the bone easily and accurately. The investigations of cutting forces applied to the tool edges show less force in the toothed edge than the knife edge.

A novel method to determine the health of the industrial parts is using the ultrasound scattering waves. Any changes in the structure of the scattering object or in the boundary conditions will cause a change in the scattering field. The frequency spectrum of the scattering time signal has valuable information, which is studied by resonant ultrasound spectroscopy (RUS). Since any defect, property changes, or changes in boundary conditions can affect the scattering field. Therefore, the possible defects in the piece are detected using the response of the scattering field. One possible defect in the fiber-reinforced composites is the existence of a crack in the matrix or fibers. In the present study, the effect of crack on the far-field backscattering amplitude spectrum is investigated using the finite element method (FEM). To this end, the effect of the crack’s direction in the cylinder and matrix on the form function is scrutinized. The results show that the Rayleigh frequency modes are more sensitive to the cracks existing in the epoxy matrix than the Whispering-gallery frequency modes. Also, the existence of the crack in the aluminum cylinder has the most effect on the Whispering-gallery frequency modes. Besides, the existence of a horizontal crack in the aluminum cylinder leads to a significant reduction in these frequency modes. The validation of the research is determined by comparing the steel cylinder form function obtained from the finite element method’s information and the analytical and experimental form functions in addition to the comparison of the aluminum cylinder form function and reference form function.

In this article, the influence of microstructure in the Casson fluid flow through a porous medium is investigated, by extending the Buckingham-Reiner’s one-dimensional model to plane-Poiseuille flow and Hagen-Poiseuille flow geometries. While analyzing the flow characteristics in single-channel/pipes or multiple channels/pipes of different width/radius, four different probability density functions are used to model the pores widths/radii distributions. It is found that when the pressure gradient increases, the Buckingham-Reiner function raises slowly in the plane-Poiseuille flow, whereas in Hagen-Poiseuille flow, it rises rapidly. In all kinds of distribution of pores, the fluid’s mean velocity and porosity of the flow medium are considerably higher in the Hagen-Poiseuille flow than in the plane-Poiseuille flow, and this behavior is reversed for the permeability of the flow medium. The fluid’s mean velocity, porosity, and permeability of the flow medium increases appreciably with the rise of the channel width and pipe radius. The porosity of the flow medium slumps with the rise of the periodof the channels and pipes distribution from 1 to 2, and it decreases very slowly with the further rise of the period H of the channels and pipes from 2 to 11.

Unmanned Aerial Vehicles (UAVs) are becoming popular due to its versatile maneuvering and high pay load carrying capabilities. Military, navy and coastal guard makes crucial use of the amphibious UAVs which includes the working functionalities of both hover craft and multi-rotor systems. Inculcation of these two systems and make it as amphibious UAV for water quality monitoring, sampling and analysis is essential to serve the human-kind for providing clean water. On this note, an amphibious UAV is designed for carrying a water sampler mechanism with an on-board sensor unit. In order to examine the stability of designed UAV under diverse wind load conditions and to examine the aerodynamic performance characteristics, computational fluid dynamic analysis (CFD) is performed. For various flight conditions such as pitch, roll, yaw and hovering, the flow characteristics around the vehicle body is examined. The aerodynamic phenomenon at the rotor section, vortex, turbulent regions, wake and tip vortex are identified. In addition, CFD analysis are conducted to determine the thrust forces during forward and hovering conditions through varying the wind speed 3 to 10 m/sec and speed of rotor 2000 to 5000 rpm. The effect of non-dimensional parameters such as advance ratio and induced inflow ratio on estimating the thrust characteristics are studied. Simulation results suggested that at 5° angle of attack and 8 m/sec wind speed condition, the aerodynamic performance of the vehicle is superior and stable flight is guaranteed. The amphibious UAV with flying and gliding modes for collecting water samples in remote water bodies and also in-situ water quality measurement can be well utilized for water quality monitoring.

The combination of mechanical and bonded joints creates a new connection type, called hybrid joint which has the benefits of both mechanical and bonded joints. In this research, the mechanical behavior of the tubular bonded lap joint between aluminum tubes subjected to pure torsion has been investigated experimentally and numerically, and the results have been compared. The mechanical behavior of the hybrid (rivet/bonded) joint has been investigated numerically and the outcomes have been compared. The adhesive and rivets have cohesive elements and bushing connector elements, respectively. The results from the hybrid joints and the damage mechanism show that the rivets change the interface shear stress and the stress distribution of the joint, and affect the joint’s torque capacity and strengths. It has been observed that for specimens with overlap lengths close to effective length, the hybrid joint is more effective.

The separation between repair mortars and the concrete substrate is one of the serious problems in repairing concrete structures. One of the main causes of this separation is the lack of proper curing and, consequently, excessive shrinkage of the repair mortar, which reduces the bond strength between the concrete substrate and the repair layer and has an adverse effect on the compressive and tensile strength of the repair mortars. In this paper, the mechanical properties, shrinkage of repair mortars, as well as their shear and tensile bond strength is investigated on the concrete substrate of different ages under the curings of "abandoned in the laboratory space," "water-submerged" and "curing agent." In-situ "friction-transfer" and "pull-off" methods are used to measure adhesion. Furthermore, the relationships between compressive strength, tensile strength, and readings are obtained from "friction-transfer" and "pull-off" methods on repair mortars and the stress distribution method used in the above-mentioned methods are presented using nonlinear finite element analysis (Abaqus/CAE). The results indicate a significant effect of curing method on shrinkage and mechanical properties of repair mortars; as a result, effective curing increases the shear and tensile bond strength at the substrate and repair layer joint boundary. It is also observed that there is a linear relationship between the experimental results obtained from the two methods used in this study with a high correlation coefficient, highly consistent with the results obtained from nonlinear finite element analysis. Thus, they can be used as in-situ methods for determining the compressive and tensile strength of repair mortars.

In this article‎, ‎a fluid dynamic code is implemented to investigate a non-linear model for electro-osmotic flow through a one-dimensional Nano-channel‎. ‎Certain mathematical techniques are simultaneously utilized to convert the coupled system of equations into a non-linear differential correlation‎. ‎This correlation is based on the mole fraction of anion‎. ‎By using a modified homotopy perturbation method‎, ‎the achieved non-linear differential equation is converted into a few linear differential equations‎. ‎The mole fraction of anion across the channel is found by solving the linear differential equations‎. ‎Finally‎, ‎the distribution of the mole fraction of cation‎, ‎electrical potential, and velocity are accurately derived based on the mole fraction of anion‎. ‎The present study confirms that by application of a modified homotopy perturbation method‎, ‎the results are in acceptable agreement with the previously validated data‎. ‎However, using the proposed method here‎, ‎a closed-form of the solution is achieved‎.

Reflection and refraction phenomenon pattern of elastic plane wave at the interface between anisotropic monoclinic elastic half-space and isotropic elastic half-spaces is studied. Closed-form expression for phase velocity is obtained. Reflection and transmission coefficients are obtained using the method of Cramer's rule in determinant form. Also, the energy ratios are calculated in terms of reflection and transmission coefficients. Numerical examples are considered to exhibit all the findings graphically. The energy conservation law is implemented at each angle of incidence to validate the numerical results, and it is found that energy ratios are in good agreement with the energy conservation law.

In This paper, the developed Haar wavelet method for solving boundary value problems is described. As known, the orthogonal Haar basis functions are applied widely for solving initial value problems, but In this study, the method for solving systems of ODEs associated with multipoint boundary conditions is generalized in separated or non-separated forms. In this technique, a system of high-order boundary value problems of ordinary differential equations is reduced to a system of algebraic equations. The experimental results confirm the computational efficiency and simplicity of the proposed method. Also, the implementation of the method for solving the systems arising in the real world for phenomena in fluid mechanics and construction engineering approves the applicability of the approach for a variety of problems.

Chaos is a phenomenon that occurs in some non-linear systems. Therefore, the output of the system will be heavily dependent on the initial conditions. Since the main characteristic of the chaos is an abnormal behavior of the system output, it should be considered in designing control systems. In this paper, controlling chaos phenomenon in a time-variant non-linear gear transmission system is investigated. To do this, a non-linear model for the system is introduced considering the effective parameters of the system, and then it is shown that chaos appears in the system by plotting phase plane of state-space variables. It should be noted that there is a great difference between random and chaotic behavior. In random cases, the model or input contains uncertainty, and therefore, the system behavior and output are not predictable. However, in chaotic behavior, there is only a brief uncertainty in the system model, input or initial conditions, and designing controller based on output prediction could be achieved. Therefore, model predictive control (MPC) algorithms are used to control the chaos, using the output prediction concept. In many cases, perturbation term also can be considered as uncertainty, and therefore, a robust controller family can be used for eliminating chaos. Both generalized predictive controller (GPC) and sliding mode controller (SMC) are used for chaos control here. The simulation results show the efficiency of the proposed algorithms.

In transport as well as manufacturing industries, the two basic aspects are heating and cooling. The use of metal or metallic oxide nanofluids has an effective cooling technique than that of conventional fluids. Therefore, the work is aimed at describing the three-dimensional MHD flow of metal and metallic oxide nanofluids past a stretching/shrinking sheet embedding with a permeable media. Further, thermal properties are enhanced by incorporating heat generation/absorption and radiative heat energy in the heat equation, enhancing the efficiency of temperature profiles. The convective boundary condition for temperature is used, which affects the temperature profile. Suitable similarity transformation is used to transform the governing equations to ordinary differential equations. The approximate analytical solution is obtained for these transformed differential equations employing the Adomian Decomposition Method (ADM). The influences of characterizing parameters are obtained and displayed via graphs, and the computation results of the heat transfer rate for various values of constraints are shown in a table. It is observed that both the momentum and energy profiles decrease with an enhance in the porosity parameter. Also, the fluid temperature decreases with an increasing thermal radiation parameter, but the opposite effect is encountered for the energy generation/absorption parameter.

This numerical study is conducted to scrutinize the dual solutions and stability analysis of the flow of Casson nanofluid past a permeable extending/contracting wedge and stagnation point. Momentum, heat and mass transfer behaviors of the Casson nanofluid have been modeled with the use of the Buongiorno nanofluid model. Suitable self-similarity variables are employed to convert the fluid transport equations into ordinary differential equations and the bvp4c MATLAB solver is used to solve the equations. The impacts of active parameters on fluid transport properties are illustrated graphically. The outcomes of the present analysis reveal that the influence of Casson fluid parameter on velocity and temperature distributions obtained from the first and second solutions exhibit the opposite natures. From the stability analysis, it is found that the thermophoresis and Brownian motion effects acquire the same critical point value on Nusselt number. The temperature distribution of the Casson nanofluid is higher over the wedge than stagnation point. The two solutions are found for the limited range of extending/contracting parameter. The detailed stability test is carried out to determine which of the two solutions is physically realizable and stable.

In this article, the nonlinear vibrational behavior of a nano-disk was analyzed using the multiple scales method (MSM). The modified couple stress theory was used to consider the small-scale effect via the application of nonlocal parameter. Employing Hamilton's principle, two coupled nonlinear differential equations were derived based on the nonlinear von-Kármán strain-displacement relation and the classical plate theory. The Galerkin-based procedure was utilized to obtain a Duffing-type nonlinear ordinary differential equation with a cubic nonlinear term and solved by the application of MSM. The effects of nonlocal parameter, aspect ratio, different boundary conditions, and the nonlinear shift frequencies, were obtained on the overall behavior of the nano-disk. Results indicate that increasing the central dimensionless amplitude of the nano-disk, the nonlinear frequency, and the shift index exhibit an increasing behavior, while the increase in the non-dimensional nonlocal parameter, causes a decrease in the nonlinear frequency ratios and the shift index. Additionally, the increase in h/r increases the effect of dimensionless central amplitude on the nonlinear frequencies ratios. Additionally, comparison of the current results with those previously published in the literature shows good agreements. This indicates that the MSM can ease up the solution, and hence, can be applied to the solution of nonlinear nano-disks with high accuracy.

This work is devoted to the direct impact method for determining the deformation diagrams of viscoplastic materials at high strain rates. As the conventional Split Hopkinson Pressure Bar method, the direct impact method is based on the measuring bar technique. The description of the experimental scheme and the traditional experimental data proceeding method are given. The description and the results of numerical analysis of the direct impact scheme are presented. A modified procedure for processing experimental information is proposed which allows to expand the area of correct calculation of strains in the specimen according to the experimental data obtained by the direct impact method. As an illustration the deformation diagrams of copper S101 and aluminum alloy D16T in the strain rate range from 1000 to 10000 s-1 have been obtained using the Split Hopkinson Pressure Bar method and the direct impact method. The use of the direct impact method made it possible to obtain deformation curves at strain rates an order of magnitude higher than the conventional SHPB method. The range of studied plastic deformations is increased by 4 times for the case of copper and 3 times for aluminum alloy.

The direct transcription method that employs global collocation at Legendre-Gauss-Radau points is addressed and applied to infinite-dimensional dynamic optimization problems in engineering. The formulation of these latter is considered referring to a Bolza-type performance index. A reduced unconstrained form of it is particularly studied in the pseudospectral domain and the continuous-to-discrete conversion is thoroughly discussed. An equivalent finite-dimension nonlinear programming problem is therefore obtained and hints on its numerical implementation are given. Eventually, a few benchmark historical problems in engineering are revisited, stated, numerically solved and compared to literature.

Peridynamics is a recent nonlocal theory of continuum mechanics that is suitable to describe fracture problems in solid mechanics. In this paper, a new failure criterion based on the stress field is developed by adopting the damage correspondence model in the ordinary state-based peridynamic theory. The proposed stress tensor-based failure criterion is capable of predicting more accurately crack propagation in the mixed mode I-II fracture problems different from other failure criteria in peridynamics. The effectiveness of the proposed model is demonstrated by performing several examples of mixed-mode dynamic fracture in brittle materials.

The dynamic response of an anisotropic multilayer magneto-electro-elastic (MEE) plate due to an external excitation is investigated in this work, using the stiffness matrix approach. A parametric study is performed by varying the stacking sequence, polarization direction, boundary conditions, and interlayer thickness. The proposed method yields the numerical estimation of the dispersion curve and the free vibration of the Lamb waves for different crystallographic orientations and magneto-electric boundary conditions. It is demonstrated that the anisotropy highly affects the dispersion curve and the non-dimensional frequency, and also decreases the interlayer thickness of the magneto-electro-elastic multilayer and raises the phase velocity of the fundamental symmetric Lamb mode vibration. The key outcomes of this research can serve as a reference in the design and analysis of new smart magneto-electro-elastic structures.

The spectral method (p-FEM) is used to solve the problem of a thin-walled structure deformation, such as a stiffened panel. The problem of the continuous conjugation of the membrane function from H1 and the deflection function from H2 was solved by modifying the “boundary” functions. Basis systems were constructed that satisfy not only the essential but also the natural boundary conditions, which made it possible to increase the rate of convergence of the approximate solution. The veracity of the results is confirmed by comparing the obtained spectral solution with the solution obtained by the h-FEM. It has been shown that the exponential rate of convergence characteristic of spectral methods is preserved if the Gibbs phenomenon is avoided. The constructed basis systems can be effectively used for solving various problems of mechanics.

The paper deals with the investigation of the flow regimes of water and crude oil in a Y-microchannel. Systematic experiments allow distinguishing four different types of water-oil flow regimes in studied microchannel, namely, plug regime, droplet regime, parallel regime, and the first-revealed chaotic regime, as well as defining the ranges of parameters corresponding to these regimes. The dependences of the water plug length and frequency of their formation in crude oil on various parameters are obtained, and universal correlations are derived. It was demonstrated that when determining the dimensionless plug lengths and the Strouhal number corresponding to droplet detachment, it is necessary to take into account not only the ratio of the dispersed and carrier phase flow rates but also the capillary number. The systematic numerical simulation of the considered flows is conducted using the VOF method. The conducted comparison of the computation and experimental results has shown that the numerical simulation qualitatively correctly predicts the characteristics of water and oil flows in the Y-type channel under consideration.

An influence of mechanical impacts between variable capacitor electrodes on the electrostatic vibration energy harvester (e-VEH) operation is studied theoretically. The analysis is carried out for two conditioning circuits with parallel and serial load connection. A relationship between e-VEH parameters and external mechanical force characteristics enabling to assess the possibility of operation in a periodic impact mode is obtained. Dependences of the average power generated by the impact-enhanced e-VEH versus the number of collisions between the electrodes and the load resistor value are calculated. The operation of the harvester for two circuits in impact and non-impact modes is compared and analyzed. It is shown that the average power generated by the e-VEH for the impact mode can exceed the power for the non-impact mode by 1–2 orders of magnitude along with a significant decrease of the harvester optimal load resistance.

The Gyroid is a periodic minimal surface explored in different applications, such as architecture and nanotechnology. The general topology is suitable for the construction of porous structures. This paper presents a non-iterative, novel methodology for the generation of a NURBS-based Gyroid volume. The Gyroid fundamental patch is defined with the Weierstrass parameterization. Furthermore, the geometry is manipulated to generate a structured mesh, allowing better element quality for FEM and IGA simulations. The re-parametrization is carried out by a Least-Squares approximation with a parametric NURBS surface, enabling a better definition of the mid-surface normals for the generation of the complete Gyroid volume. Different cases of variational thickness and porosity are presented to validate the versatility of our method.

A two-carrier planetary gear train (PGT) configuration, named S13WN, was developed for a specific application with a negative transmission ratio in the interval between -20 and -21. 72 valid combinations of component PGTs ideal torque ratios have been developed for this configuration, from which kinematically feasible combinations providing minimal dimensions and maximum efficiency within the required interval have been selected. The minimal dimensions of the PGT were achieved with a cylindrical case shape, i.e., with the ratio of the ring gears reference diameters close to unity. Nine other PGT configurations have been synthesized for the same transmission ratio, and their design parameters optimized. The two solutions that offer the most improvements over S13WN have been developed into design concepts.

Semi-analytical studies for the porosity action of a hollow sphere made of functionally graded piezoelectric material are presented. A semi-analytical technique for radial displacement is presented. A combination of internal and external pressures, temperature distribution, uniform hygrothermal distribution, and an electric potential variation has been discussed between the inner and outer surfaces of the sphere. The material physical properties for the present porous hollow sphere are varying through the thickness due to the power functions of the radius. Numerical outcomes are validated for radial displacement, electric potential, and stresses for the perfect and porous functionally graded hollow sphere. The effect of different mechanical, piezoelectric, and hygrothermal are investigated.

The paper presents an approach to the strength analysis in steel cylindrical panels reinforced from the concave side with an orthogonal grid of stiffeners. A mathematical model of the Timoshenko (Mindlin – Reissner) type is used. Transverse shears and geometric nonlinearity are taken into account. The stiffeners are introduced in two ways: using the method of refined discrete introduction (proposed by author) and the method of structural anisotropy. Computational algorithm based on the Ritz method and the best parameter continuation method. For strength analysis von Mises criterion is used. The values of the maximum permissible strength loss loads are shown for several variants of structures made of steel S345. The extension of areas of non-fulfillment of strength conditions according to the Mises criterion for the stiffened and unstiffened structures are shown.

Long-term monitoring of the safety and reliability of large dams operation has been attracting increasing attention of researchers. Moreover, special consideration is given to the study of dam displacements that characterize its global behavior. The article discusses specifics of constructing predictive mathematical models for studying the deformation process associated with displacements of the high-head dam crest. The authors present the most successfully designed predictive mathematical models for various combinations of input effective factors, including the results of field observations and the calculated values of the component displacements. These models allow forecasting the control points of the dam body for various time stages of its operation. The advantages of using a mathematical model with separate introduction of the main effective factors into the model are shown, thereby eliminating the effect of their multicollinearity. Using the example of the Sayano-Shushenskaya hydro power plant for certain time stages of the dam operation and structures with different temperature conditions (average, warm and cold in respect to annual temperatures), the authors present the results of forecasting the dam displacements.

The field of magnetohydrodynamics (MHD) encompasses a wide range of physical objects due to their stabilising effects. Thus, this study concerns the numerical investigation of the radiative non-Newtonian fluid flow past a shrinking sheet in the presence of an aligned magnetic field. By adopting proper similarity transformations, the governing partial derivatives of multivariable differential equations are converted to similarity equations of a particular form. The numerical results are obtained by using the bvp4c technique. According to the findings, increases in the suction parameter resulted in higher values of the skin friction and heat transfer rate. The same pattern emerges as the aligned angle and magnetic parameter are considered. On the other hand, the inclusion of the Bingham number, the Reiner–Philippoff fluid, and the thermal radiation parameters deteriorate the heat transfer performance, evidently. The dual solutions are established, which results in a stability analysis that upholds the validity of the first solution.

Herein, the heat transfer performance of the curved surface of a louvered fin heat exchanger using computational fluid dynamics (CFD) is examined. Four new models are used with curved surfaces in different locations. The air inlet velocity is 1–9 m/s. The air and fin wall temperature remain constant at 300 and 353 K, respectively. The result of the reference flat fin is confirmed with experimental results. The results demonstrate that curved fins changed the airflow path and created vortices. The air tends to flow between louver fins, improving its velocity and enhancing heat transfer. The result from the case that individual fin is close to the middle fin demonstrated that louver fin provides a 15% increase compared to that of the reference. However, when the air inlet velocity is high, the performance evaluation criteria from the case that individual fin is close to the first fin, is the highest, which results in a 1% increase from that of the reference. Therefore, increasing heat transfer can compensate the effect of pressure drop because of vortices in the louver fin domain. This study can be applied to the air conditioning system to increase its efficiency and cut the operation cost down.

In this paper, dynamic results of carbon nanotube (CNT)-reinforced composite beams under a moving load are presented. The constitutive equations in motion are obtained by the Lagrange procedure according to Timoshenko beam theory and then solved by using the Ritz method. In the solution of the moving load problem, the Newmark average acceleration method is used in the time history. In the numerical results, the effects of CNTs’ volume fraction, patterns of CNTs, and the velocity of moving load on the dynamic responses of CNT-reinforced composite beam are investigated in detail. It is observed that the reinforcement patterns and volume fraction of CNTs are very effective on the behavior of the moving load. Also, it is found that X-Beam and O-Beam have the biggest and lowest rigidities in all models, respectively.