Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Biomechanical Evaluation of Bone Quality Effect on Stresses at Bone-Implant Interface: A Finite Element Study126612751548110.22055/jacm.2020.32323.2005ENBelaidTAHAROUFaculty of Science & Technology, University of Mascara, B.P.305 - Mamounia Road, Mascara (29000), Algeria0000-0001-7014-5027AliMERDJIFaculty of Science & Technology, University of Mascara, B.P.305 - Mamounia Road, Mascara (29000), AlgeriaLaboratory of Mechanics Physical of Material (LMPM), Faculty of Technology, University of Sidi Bel-Abbes, Algeria0000-0003-1280-9723RajshreeHILLSTROMDepartment of Bioengineering, Tandon School of Engineering, New York University, USA0000-0002-2365-5486AliBENAISSAFaculty of Science & Technology, University of Mascara, B.P.305 - Mamounia Road, Mascara (29000), Algeria0000-0002-7566-7171SandipanROYDepartment of Mechanical Engineering, SRM Institute of Science and Technology, Chennai-603203, IndiaNoureddineDELLAFaculty of Science & Technology, University of Mascara, B.P.305 - Mamounia Road, Mascara (29000), Algeria0000-0002-6759-6916AbdelkrimAIDFaculty of Science & Technology, University of Mascara, B.P.305 - Mamounia Road, Mascara (29000), Algeria0000-0001-6355-958XOsamaMUKDADIDepartment of Mechanical and Aerospace Engineering, West Virginia University, Morgantown, WV 26506, USA0000-0001-9250-0720Journal Article20200124The aim of this study is to evaluate the effect of the alveolar bone quality on von Mises stress at the bone-implant interface during occlusal loading. Four (3D) finite element models of fully osteointegrated 3-mm diameter × 11.5-mm length dental implant indifferent alveolar bone with different cortical bone thickness are created, using SolidWorks computer aided design software. The alveolar bone cortical-spongy bone ratio modelled includes I) 90%-10%, II) 60%-40%, III) 40%-60%, and IV) 10%-90%. These models are then exported to ABAQUS software and stress analyses are run under an occlusal load of 70 N acting on the platform face of the dental implant. Results of this study show that the implants are subjected to similar stress distributions in all models; maximum stress values are confined in the outer cervical plate of the cortical bone around the neck. This could explain bone loss and implant de-osseointegration. Peak stresses are lowest in the model with 90% cortical bone (14.2 MPa) and almost doubled in the model with 10% cortical bone (26.6 MPa). The stress values gradually reduce towards the apical area, demonstrating masticatory force transfer from implant to bone. Furthermore, both cortical and spongy bone structures exhibit highest stress values in the model with thinnest cortical layer. The high interfacial stress concentration near the implant-cortical bone junction could lead to bone failure or implant instability induced by fatigue or overload risk. Results of our study could be a first step towards the development of a clinical pre-operative planning tool for dental implantolgy.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Structural Strength Analysis and Fabrication of a Straight Blade of an H-Darrieus Wind Turbine127612821535210.22055/jacm.2020.31452.1876ENFatehFerroudjiUnité de Recherche en Energies Renouvelables en Milieu Saharien, URERMS, Centre de Développement des Energies Renouvelables, CDER, 01000, Adrar, AlegriaLaboratoire de Mécanique des Structures et Matériaux, Université Batna 2, Algeria0000-0003-1467-9221CherifKhelifiUnité de Recherche en Energies Renouvelables en Milieu Saharien, URERMS, Centre de Développement des Energies Renouvelables, CDER, 01000, Adrar, AlegriaJournal Article20191019Small H-Darrieus wind turbines have become popular in the wind power market because of their many advantages, which include simplicity of design, low construction costs, and they are thought to represent an adequate solution even in unconventional installation regions. The blade is generally considered as the most important component of the wind turbine system because it controls the efficiency of the turbine. The blade structure must be designed to support the difficult environmental conditions (e.g., wind, and snow) encountered during the operational life of the wind turbine. This current study uses three-dimensional (3D) modeling and structural strength analysis to fabricate two straight blades (aluminum and galvanized steel) for a small H-Darrieus wind turbine. The 3D modeling of the blade structure is performed using SolidWorks, a computer-aided design (CAD) software package, and the structural strength analysis uses the Finite Element Analysis (FEA) technique to identify the stiffness, resistance, and reliability of the blade structure. The simulation results obtained indicate that no structural failures are predicted for either of the two structures tested because the factors of safety are larger than one, and the all maximum deflections are within the allowable deformation limits for the materials. Manufacturing processes for the two structures are described.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Unsteady Separated Stagnation Point Flow of Nanofluid past a Moving Flat Surface in the Presence of Buongiorno's Model128312901555010.22055/jacm.2020.32261.1992ENA.RenukaDepartment of Mathematics, Bharathiar University, Coimbatore 641 046, India0000-0001-6387-9954M.MuthtamilselvanDepartment of Mathematics, Bharathiar University, Coimbatore 641 046, India0000-0003-0075-2520Qasem M.Al-MdallalDepartment of Mathematical Sciences, United Arab Emirates University, P. O. Box 15551, Al Ain, Abu Dhabi0000-0002-2853-9337Deog-HeeDohDivision of Mechanical Engineering, College of Engineering, Korea Maritime and Ocean University, Busan 606 781, South Korea00000-0001-6093-2975BahaaeldinAbdallaDepartment of Mathematics and General Sciences, Prince Sultan University, P.O. Box 66833, Riyadh 11586, Saudi Arabia 0000-0001-8867-0612Journal Article20200113This paper explores energy and mass transport behavior of unstable separated stagnation point flow of nanofluid over a moving flat surface along with Buongiorno’s model. Characteristic of Brownian diffusion and thermophoresis are considered. Additionally, characteristics of chemical reaction is taken into account. A parametric investigation is performed to investigate the outcome of abundant parameters such as temperature, velocity and concentration. An appropriate equation is converting into a set of ODEs through employing appropriate transformation. The governing equations has been solved numerically by using the classical fourth-order Runge-Kutta integration technique combined with the conventional shooting procedure after adapting it into an initial value problem. Our findings depict that the temperature field θ(ζ) improves for augmenting values of theromophoresis parameter (Nt) with dual solutions of attached flow without inflection and flow with inflection. Also, the difference of Brownian motion parameter (Nb) with two different solutions of attached flow exists with energy profile. It can be found that an energy profile θ(ζ) elevates due to augmenting values of (Nb). It has been perceived that thermal boundary layer thickness elevates due to large amount of Brownian motion parameter (Nb).Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701ADM Solution for Cu/CuO –Water Viscoplastic Nanofluid Transient Slip Flow from a Porous Stretching Sheet with Entropy Generation, Convective Wall Temperature and Radiative Effects129113051562210.22055/jacm.2020.33137.2167ENTHIRUPATHITHUMMADepartment of Mathematics, B V Raju Institute of Technology, Narsapur, Medak, 502313, Telangana State, India0000-0002-7993-5647S.R.MishraDepartment of Mathematics, Siksha ‘O’ Anusandhan Deemed to be University, Bhubaneswar, 751030, Odisha, India0000-0002-3018-394XOsman AnwarBégFluid Mechanics, Propulsion and Nanosystems, Aeronautical and Mechanical Engineering, Newton Bldg., University of Salford, Manchester, M54WT, UK0000-0002-7993-5647Journal Article20200403A mathematical model is presented for entropy generation in transient hydromagnetic flow of an electroconductive magnetic Casson (non-Newtonian) nanofluid over a porous stretching sheet in a porous medium. The model employed is Cattaneo-Christov heat flux to simulate non-Fourier (thermal relaxation) effects. A Rosseland flux model is implemented to model radiative heat transfer. The Darcy model is employed for the porous media bulk drag effect. Momentum slip is also included to simulate non-adherence of the nanofluid at the wall. The transformed, dimensionless governing equations and boundary conditions (featuring velocity slip and convective temperature) characterizing the flow are solved with the Adomian Decomposition Method (ADM). Bejan’s entropy minimization generation method is employed. Cu-water and CuO-water nanofluids are considered. Extensive visualization of velocity, temperature, and entropy generation number profiles is presented for variation in pertinent parameters. The calculation of skin friction and local Nusselt number are also studied. The ADM computations are validated with simpler models from the literature. The solutions show that with elevation in the volume fraction of nanoparticle and Brinkman number, the entropy generation magnitudes are increased. An increase in Darcy number also upsurges the friction factor and heat transfer at the wall. Increasing volume fraction, unsteadiness, thermal radiation, velocity slip, Casson parameters, Darcy, and Biot numbers are all observed to boost temperatures. However, temperatures are reduced with increasing non-Fourier (thermal relaxation) parameter. The simulations are relevant to the high temperature manufacturing fluid dynamics of magnetic nano liquids, smart coating systems.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Study on Heat Transfer and Pressure Drop in a Mini-Channel with Corrugated Walls130613141573310.22055/jacm.2020.34487.2418ENAbdelazizBegagApplied Energy and Thermal Laboratory (ETAP), Department of Mechanical Engineering, Faculty of Technology, Abou Bekr Belkaid University, B.P 119, Tlemcen, 13000, AlgeriaCarnot de Bourgogne Interdisciplinary Laboratory - Site UTBM Sévenans - ICB UMR 6303 CNRS / Bourgogne Franche Comté University. (UBFC) 90010 Belfort cedex, France0000-0002-0330-2974RachidSaimApplied Energy and Thermal Laboratory (ETAP), Department of Mechanical Engineering, Faculty of Technology, Abou Bekr Belkaid University, B.P 119, Tlemcen, 13000, Algeria0000-0002-7509-6612Hakan F.ÖztopDepartment of Mechanical Engineering, Faculty of Technology, Firat University, 23119 Elazig, TurkeyDepartment of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan0000-0002-2161-0639SaidAbboudiCarnot de Bourgogne Interdisciplinary Laboratory - Site UTBM Sévenans - ICB UMR 6303 CNRS / Bourgogne Franche Comté University. (UBFC) 90010 Belfort cedex, Francehttps://orcid.org/00Journal Article20200727This study presents the numerical results relative to the development of heat transfer and pressure drop inside a corrugated channel, under constant heat flux conditions applied to the walls; the working fluid is air. The test section is a channel with two plates having trapezoidal-shaped corrugations with V-folds. The corrugated plates were placed inside a 12.5 m high channel and tested for three different inclination angles, i.e. 20°, 40° and 60°. The model was simulated for a heat flux of 0.58 kW /m<sup>2</sup>, while the Reynolds numbers were considered within the interval ranging from 600 to 1400. The standard turbulent model (<em>k-</em><em>ε</em>) was employed to simulate the flow and heat transfer developments within the channel. In addition, the governing equations were solved using the finite volume method in a <em>structured uniform grid</em> arrangement. Moreover, the effects of the geometric parameters on heat transfer and flow evolution were discussed as well. It is also worth noting that the corrugated surface had a significant impact on the enhancement of heat transfer and pressure drop due to breakage and destabilization occurring in the thermal boundary layer.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Comparison of Functionally Graded Hip Stem Implants with Various Second-Generation Titanium Alloys131513231549810.22055/jacm.2020.32964.2115ENTawakol A.EnabProduction Engineering and Mechanical Design Department, Faculty of Engineering, Mansoura University, P.O. 35516 Mansoura, Egypt0000-0003-1201-6547NohaFoudaProduction Engineering and Mechanical Design Department, Faculty of Engineering, Mansoura University, P.O. 35516 Mansoura, Egypt0000-0002-1206-7486IbrahimEldesoukyProduction Engineering and Mechanical Design Department, Faculty of Engineering, Mansoura University, P.O. 35516 Mansoura, Egypt0000-0003-3926-5471Journal Article20200316Total hip Arthroplasty (THA) is performed every year at a very high frequency to improve the quality of life of thousands of patients all over the globe. Nevertheless, the expected service life of such surgery remains unsuitable for patients under 50 years old. This is mainly related to stress shielding and the potential adverse tissue reaction to some of the elements of the market-dominant implant materials. In this research, functionally graded (FG) implant designs of several titanium alloys layered with hydroxyapatite (HA) are proposed to provide lower implant stiffness compared to a solid stem to approach the requirements of human bone. Moreover, TNZT (Ti35Nb7Zr5Ta), and TMZF (Ti12Mo6Zr2Fe) second-generation titanium alloys are studied as a replacement for the famous Ti6Al4V alloy to avoid the adverse tissue reactions related to aluminum and vanadium elements. The different FG models are numerically tested using a 3D finite element simulation after virtual implantation in a femur bone under the dynamic load of a patient descending stairs. In the numerical study, the variation in stress distribution and strain energy in a femur bone is assessed for different FG hip stems as well as the axial stiffness of the hip stems. Results indicated an increase in strain energy and von Mises stress in the cortical and cancellous bones using FG hip stems. Additionally, the axial stiffness is reduced for all FG hip stems relative to the commercial Ti6Al4V hip stem.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Analysis and Optimization of Truss Structures, Constrained Handling using Genetic Algorithm132413331555210.22055/jacm.2020.33076.2143ENPal Ranjan SastiCharanDepartment of Mechanical Engineering, National Institute of Technology Durgapur, West Bengal, India0000-0003-0183-8668Nirmal BaranHuiDepartment of Mechanical Engineering, National Institute of Technology Durgapur, West Bengal, India0000-0002-6225-9371J.Paulo DavimDepartment of Mechanical Engineering, University of Aveiro, Aveiro, Portugal0000-0002-5659-3111Journal Article20200328In this study, an attempt is made to minimize the weight of Howe roof and ten member-6 Node trusses, separately. Two constraints, maximum allowable deflection and maximum allowable member stresses have been considered. For the first truss, permissible deflection is not known from the literature; therefore, it is determined using the exhaustive search method. Once magnitudes of the constraints are identified, member cross-sectional areas are varied to get the optimal weight. Both the exhaustive search method and the genetic algorithm have been implemented for this purpose. During the optimization, members tending to form a string may be eliminated from the structure. Doing this, we could further reduce the weights of the trusses and even less than the minimum available in the literature. The second truss is an indeterminate structure, and Maxwell Betti reciprocal theorem is applied to calculate the member forces. Also, further reduction of members is made for this truss, keeping in mind that the truss becomes determinate with the decrease in the member(s).Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Non-similar Solutions of MHD Mixed Convection over an Exponentially Stretching Surface: Influence of Non-uniform Heat Source or Sink133413471413310.22055/jacm.2019.27729.1432ENP.M.PatilDepartment of Mathematics, Karnatak University, Pavate Nagar, Dharwad – 580003, IndiaD.N.LathaDepartment of Mathematics, Karnatak University, Pavate Nagar, Dharwad – 580003, IndiaAli J.ChamkhaFaculty of Engineering, Kuwait College of Science and Technology, Doha, Kuwait0000-0002-8335-3121Journal Article20181129In this paper, an analysis of magnetohydrodynamic (MHD) mixed convection over an exponentially stretching surface in the presence of a non-uniform heat source/sink and suction/injection is presented. The governing boundary layer equations are transformed into a set of non-dimensional equations by using a group of non-similar transformations. The resulting highly non-linear coupled partial differential equations are solved by using the implicit finite difference method in combination with the quasilinearization technique. Numerical results for the velocity, temperature and concentration profiles, as well as the skin friction coefficient, wall heat transfer and mass transfer rates are computed and presented graphically for various parameters. The results indicate that the velocity profile reduces, while the temperature profile increases in presence of the effects of magnetic field and suction at the wall. The velocity ratio parameter increases the skin-friction coefficient and the Schmidt number decreases the wall mass transfer rate. The temperature profile increases for the positive values of Eckert number and space as well as temperature dependent heat source/sink parameters, while the opposite behavior is observed for negative values of same parameters.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Optimum Design of Nano-Scaled Beam Using the Social Spider Optimization (SSO) Algorithm134813611504410.22055/jacm.2019.31406.1870ENBüşraUzunDepartment of Civil Engineering, Faculty of Engineering, UludağUniversity, Bursa, Turkey0000-0002-7636-7170OmerCivalekResearch Center for Interneural Computing, China Medical University, Taiwan0000-0003-1907-9479IbrahimAydogduDepartment of Civil Engineering, Faculty of Engineering, Akdeniz University, Antalya, 07050, Turkey0000-0002-8281-2365Journal Article20191016In this research study, the optimum cross-sectional dimensions of nano-scale beam elements are investigated under different load conditions. Euler-Bernoulli beam model based on nonlocal elasticity theory is utilized for the analysis of the beam. Two types of nano-scaled beams are modeled; carbon nanotubes (CNTs) and Boron nitride nanotubes (BNNTs). The novel meta-heuristic based optimization algorithm called Social Spider Optimization (SSO) algorithm is employed to find the beam designs with the objective of minimizing the cross-sectional area. Furthermore, the obtained optimum cross-sectional dimensions, critical stress and displacement values of the beams are compared according to the material type, beam length, and load conditions.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Analysis of an Edge Crack Isotropic Plate with Void/Inclusions under Different Loading by Implementing XFEM136213821501010.22055/jacm.2019.31268.1848ENAchchheLalDepartment of Mechanical Engineering, S. V. National Institute of Technology, Surat-395007, India0000-0001-7496-2112Manoj B.VaghelaDepartment of Mechanical Engineering, S. V. National Institute of Technology, Surat-395007, IndiaKundanMishraDepartment of Mechanical Engineering, S. V. National Institute of Technology, Surat-395007, India0000-0001-7802-0982Journal Article20191001In the present work, the effect of various discontinuities like voids, soft inclusions and hard inclusions of the mixed-mode stress intensity factor (MMSIF), crack growth and energy release rate (ERR) of an edge crack isotropic plate under different loading like tensile, shear, combine and exponential by various numerical examples is investigated. The basic formulation is based on the extended finite element method (XFEM) through the M interaction approach using the level set method. The effect of single and multi voids and inclusions with position variation on MMSIF and crack growth are also investigated. The presented results would be applicable to enhancing the better fracture resistance of cracked structures and various loading conditions.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Comparative Study of Mixed Convective MHD Cu-Water Nanofluid Flow over a Cone and Wedge using Modified Buongiorno’s Model in Presence of Thermal Radiation and Chemical Reaction via Cattaneo-Christov Double Diffusion Model138314021539510.22055/jacm.2020.32143.1975ENSawan KumarRawatDepartment of Mathematics, Statistics and Computer Science, G.B. Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, IndiaHimanshuUpretiDepartment of Mathematics, Statistics and Computer Science, G.B. Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, IndiaManojKumarDepartment of Mathematics, Statistics and Computer Science, G.B. Pant University of Agriculture and Technology, Pantnagar, Uttarakhand 263145, IndiaJournal Article20200103The steady Cu-water nanofluid flow in presence of magnetic field is investigated numerically under the effects of mixed convection, thermal radiation and chemical reaction. For investigating the nanofluid flow, the flow over two different geometries, cone and wedge have been considered. The Tiwari and Das nanofluid model is implemented together with Buongiorno nanofluid model. Thermal and concentration diffusion are studied using the Cattaneo-Christov double diffusion model. At the boundary of the surface, no slip and zero mass flux condition are implemented to control the nanoparticle volume fraction at surface. Constitutive laws of flow are obtained in form of ordinary differential equations by the use of similarity transformation. The modeled flow problem is solved numerically by the Runge-Kutta-Fehlberg method and shooting scheme. Variation in flow properties due to parameters involved is presented graphically and through tabular values. The effect of thermal radiation and thermal relaxation parameter is to increase heat transfer. The temperature of nanofluid and drag force at surface increases due to enhanced magnetic field. The nanoparticles are found to be concentrated near the surface of cone and wedge but concentration decreases with chemical reaction parameter and Schmidt number as fluid moves towards far field.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Simulations of Unsteady 3D MHD Micropolar Fluid Flow over a Slendering Sheet140314121526010.22055/jacm.2020.31062.1821ENS.R.R.ReddyDepartment of Mathematics, SAS, Vellore Institute of Technology, Vellore-632014, India0000-0001-5501-570XP.B.Anki ReddyDepartment of Mathematics, SAS, Vellore Institute of Technology, Al-Kharj, Saudi Arabia0000-0001-5824-1796Journal Article20190912The purpose of the present analysis is to explore the numerical investigation on the time-dependent 3D magnetohydrodynamic flow of micropolar fluid over a slendering stretchable sheet. The prevailing PDEs are rehabilitated into coupled non-linear ODEs with the aid of appropriate similarity variables and then numerically calculated by applying the 4<sup>th</sup> RKM incorporate with shooting scheme. The contributions of various interesting variables are shown graphically. Emerging physical parameters on velocity, microrotation, and the surface drag coefficient are portrayed graphically. It is noticed that the microrotation profiles highly influenced by the vortex viscosity parameter and the micro-inertia density parameter. It is also concluded that the microrotation profiles (<em>h<sub>2</sub></em>) are promoted by increasing the spin gradient viscosity parameter. Excellent accuracy of the present results is observed with the formerly published as a result of a special case.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701How Wavelike Bumps Mitigate the Vortex-induced Vibration of a Drilling Riser141314241651510.22055/jacm.2021.34951.2527ENLipingTangSchool of Mechatronic Engineering, Southwest Petroleum University, Chengdu, 610500, P.R. China0000-0003-0981-2120HongYaoSchool of Mechatronic Engineering, Southwest Petroleum University, Chengdu, 610500, P.R. ChinaZeminHuangChina Petroleum Pipeline Engineering Corporation, Langfang, 065000, P.R. ChinaLiqiWangSchool of Mechatronic Engineering, Southwest Petroleum University, Chengdu, 610500, P.R. ChinaXiaohuaZhuSchool of Mechatronic Engineering, Southwest Petroleum University, Chengdu, 610500, P.R. China0000-0002-0507-3773Journal Article20200910In this paper, computational fluid dynamics is used to study how wavelike bumps influence the suppression of drilling-riser vortex-induced vibration (VIV). The numerical model involves two-dimensional unsteady incompressible turbulent flow around a cylinder, with the flow characteristics regarded as being constant. The results show that wavelike bumps are effective in mitigating the VIV, but the degree of mitigation does not increase indefinitely with the number of bumps. The mitigation is greatest with either 5 or 7 wavelike bumps, reducing the vibration amplitudes of the cylinder in the in-line and cross-flow directions to negligible levels. To know how equipping a circular cylinder with wavelike bumps affected its VIV response, cases with wavelike bumps of 1, 3, 5, 7, 9, and 11 are studied. Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Four Legged Guará Robot: From Inspiration to Implementation142514341656010.22055/jacm.2021.35212.2613ENAntônioBento FilhoDepartment of Mechanical Engineering, Universidade Federal do Espírito Santo, Vitória-ES, 29075-910, Brazil0000-0002-2705-9294CristianePescador TonettoDepartment of Mechanical Engineering, Universidade Federal do Espírito Santo, Vitória-ES, 29075-910, Brazil0000-0002-0857-5023RafhaelMilanezi De AndradeDepartment of Mechanical Engineering, Universidade Federal do Espírito Santo, Vitória-ES, 29075-910, Brazil0000-0002-2839-3649Journal Article20201002Design of legged robots is still an open problem with several implementation challenges and there is little material available to guide new designs. Moreover, usually the actuator sets of the robot are heavy and bulky and are placed near the joints, which increases the weight and inertia of the legs and consequently results in instability, slower movement and high energy consumption. This paper presents the design, modeling, and control of the Guará, a sixteen-degree-of-freedom (DOF) legged robot, covering all stages of development of a quadruped robot, including experiments with a prototype. The system is composed by four legs, each one with 4 DOF. To reduce the leg’s weight and inertia, here we present a design concept for the legs where the actuators are positioned on the robot platform and timing belt transmissions are used to drive the joints. Moreover, the design of the legs mimics the human leg by presenting the knees faced forward and its locomotion pattern mimics the quadruped wave gait. With this configuration, the robot is able to walk on straight lines or free curved paths. The control system of the robot consists of a high-level state machine and a lower level PID position controller. The electronics are embedded over a flat platform coupled to the legs. The kinematics of the robot was studied and an integrated environment was provided for walking simulation, adjustment, and diagnosis of operation. The experimental results are focused on the kinematics of the legs and the stability of the robot and show a good agreement between the designed and executed movements. This paper is presented as a framework that can be generalized to other systems and can be a useful reference for the design of other legged robots.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Investigation on Flow Transition through a Curved Square Duct with Negative Rotation143514471570710.22055/jacm.2020.33606.2253ENMohammad SanjeedHasanDepartment of Mathematics, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj-8100, Bangladesh0000-0002-5309-4034Shamsun NaherDolonDepartment of Mathematics, Jagannath University, Dhaka-1100, Bangladesh0000-0003-0497-1915Himadri ShekharChakrabortyDepartment of Mathematics, Shahjalal University of Science and Technology, Sylet-3114, BangladeshRabindra NathMondalDepartment of Mathematics, Jagannath University, Dhaka-1100, Bangladesh0000-0001-6596-575XGiulioLorenziniDepartment of Engineering and Architecture, University of Parma, Parma 43124, Italy0000-0002-5676-8575Journal Article20200514Application of the rotational phenomena in the curved ducts plays an important role in many engineering areas, so researchers are attracted to innovate something new in this area nowadays. In this regard, the current paper has performed the fluid flow through the curved duct for an extensive range of negative rotation (-10 ≤ <em>Tr</em> ≤ -1500). The other useful parameters such as Dean number (<em>Dn</em>), Curvature (d), Grashof number (<em>Gr</em>), and Prandtl number (<em>Pr</em>) are considered fixed. The investigations are divided into four parts. In the first portion, linear stability of the flows through the duct is discussed. Then time evolution calculations of the unsteady solutions for different Taylor numbers are demonstrated in the “<em>time vs. heat flux</em>” plane. This inquiry shows that the flow undergoes various instabilities for increasing the Taylor number. Thirdly, two types of flow velocity, axial flow and secondary flow and the temperature profiles are represented. It is obtained that two up to six vortex secondary flows are found for the regular and irregular oscillation and the flow patterns are different for a fixed period for regular oscillation. To show more clarity of the periodic and chaotic flow, power spectrum density is further examined. However, it is observed that the flows are mixed and enhanced heat transfer because of the acting of centrifugal force, Coriolis force, and heating induced buoyancy force on the duct. Finally, the numerical results are compared with the experimental data which shows that the numerical data fully matches with the experimental outcome.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701MHD Double-Diffusive Natural Convection in a Closed Space Filled with Liquid Metal: Mesoscopic Analysis144814651571510.22055/jacm.2020.33719.2273ENArunSathiyamoorthiSchool of Mechanical Engineering, Vellore Institute of Technology, Vellore, Tamilnadu, 632014, India0000-0001-5690-6299SatheeshAnbalaganSchool of Mechanical Engineering, Vellore Institute of Technology, Vellore, Tamilnadu, 632014, India0000-0001-5045-5263Hakan F.ÖztopDepartment of Mechanical Engineering, Technology Faculty, Fırat University, Elazig/TurkeyDepartment of Mechanical Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah, 21511, Saudi Arabia0000-0002-2161-0639NidalAbu-HamdehDepartment of Mechanical Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah, 21511, Saudi Arabia0000-0002-4852-2217Journal Article20200523< p>In this paper, the lattice Boltzmann approach is carried out to study the double-diffusive natural convection in a space encapsulating liquid metal is presented. The Uniform magnetic field is applied horizontally at the square domain and an insulated rectangular block is kept stationary at the center of the cavity. The linear increment of temperature and concentration is used at the left wall and cold temperature is applied at the right wall. Horizontal walls are adiabatic conditions. Horizontal walls are adiabatic conditions. The numerical analysis is performed at the range of Rayleigh number (10<sup>3</sup> ≤ Ra ≤ 10<sup>5</sup>), Lewis number (2 ≤ Le ≤ 10), buoyancy ratio (-2 ≤ N ≤ 2), Hartmann number (0 ≤ Ha ≤ 50) with Prandtl number (Pr) = 0.054. Results show that the increase in Ra tends to maximize heat and mass transfer rate while increasing Ha, decreases the same. The rise in Le diminishes heat transfer marginally but increasing the mass transfer significantly. The effect of N differs with different operating conditions, in general, the rate of heat and mass transfer is found to decrease with a decrease of N value.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Behavior of Nanofluid with Variable Brownian and Thermal Diffusion Coefficients Adjacent to a Moving Vertical Plate146614791648810.22055/jacm.2021.34852.2483ENAhmed SaadRashedDepartment of Physics and Engineering Mathematics, Faculty of Engineering, Zagazig University, Zagazig, 44515, EgyptFaculty of Engineering, Delta University for Science and Technology, Gamasa, Egypt0000-0001-7800-2768Tarek A.MahmoudDepartment of Physics and Engineering Mathematics, Faculty of Engineering, Zagazig University, Zagazig, 44515, Egypt0000-0001-8846-0877Magda M.KassemDepartment of Physics and Engineering Mathematics, Faculty of Engineering, Zagazig University, Zagazig, 44515, Egypt0000-0003-3068-6923Journal Article20200831This work was motivated by studying the behavior of nanofluid adjacent to a moving vertical plate. A non-homogeneous distribution of nanoparticles inside the boundary layer was considered with variable Brownian and thermal diffusion coefficients throughout the layer. Employing group similarity transformation method transformed the governing mathematical model into a system of ordinary differential equations. The resultant system was numerically solved using shooting method. The numerical investigation was carried out for different parameters namely: Prandtl number, <em>Pr</em>, temperature difference ratio, Υ, and the ratio of nanoparticles volumetric fraction difference, Υ<sub>φ</sub>, and the attained results were illustrated graphically to examine their effect on different fluid characteristics. The results showed that increasing <em>Pr</em> values decreased the nanofluid velocity, shear stress, temperature distribution and nanoparticles volumetric fraction, while it increased the heat flux and nanoparticles gradient inside the boundary layer. On the other hand, increasing Υ values increased the nanofluid velocity, shear stress and heat flux but it decreased the temperature distribution. Also, increasing Υ<sub>φ</sub> values decreased the nanofluid velocity, shear stress and temperature distribution but it increased the heat flux. The characteristics of nanofluids were studied to enhance the thermal conductivity and the efficiency of heat transfer systems. A comparison between the obtained results and the previous published results indicated an excellent agreement.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Solving Duffing-Van der Pol Oscillator Equations of Fractional Order by an Accurate Technique148014871657510.22055/jacm.2021.35369.2642ENNourhaneAttiaDynamic of Engines and Vibroacoustic Laboratory, University M'hamed Bougara of Boumerdes, Boumerdes, 35000, Algeria0000-0002-9937-8665DjamilaSebaDynamic of Engines and Vibroacoustic Laboratory, University M'hamed Bougara of Boumerdes, Boumerdes, 35000, Algeria0000-0002-7910-3074AliAkgülArt and Science Faculty, Department of Mathematics, Siirt University, Siirt, TR-56100, Turkey0000-0001-9832-1424AbdelkaderNourDynamic of Engines and Vibroacoustic Laboratory, University M'hamed Bougara of Boumerdes, Boumerdes, 35000, AlgeriaJournal Article20201010In this paper, an accurate technique is used to find an approximate solution to the fractional-order Duffing-Van der Pol (DVP, for short) oscillators equation which is reproducing kernel Hilbert space (RKHS, for short ) method. The numerical results show that the n-term approximation is a rapidly convergent series representation and they present also the high accuracy and effectiveness of this method. The efficiency of the proposed method has been proved by the theoretical predictions and confirmed by the numerical experiments.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Effect of Temperature and Moisture on the Impact Behaviour of Adhesive Joints for the Automotive Industry148815001657310.22055/jacm.2021.36089.2793ENDanielRosendoDepartamento de Engenharia Mecânica, Faculdade de Engenharia da Universidade do Porto (FEUP), 4200-465 Porto, PortugalGuilhermeVianaDepartamento de Engenharia Mecânica, Faculdade de Engenharia da Universidade do Porto (FEUP), 4200-465 Porto, PortugalRicardoCarbasInstituto de Ciência e Inovação em Engenharia Mecânica e Engenharia Industrial (INEGI), 4200-465 Porto, Portugal0000-0002-1933-0865EduardoMarquesInstituto de Ciência e Inovação em Engenharia Mecânica e Engenharia Industrial (INEGI), 4200-465 Porto, Portugal0000- 0002-2750-8184Lucas FMDa SilvaDepartamento de Engenharia Mecânica, Faculdade de Engenharia da Universidade do Porto (FEUP), 4200-465 Porto, Portugal0000- 0003-3272-4591Journal Article20201217This study focuses on evaluating the impact of aluminum adhesive joints as a function of temperature and moisture, in an effort to understand how these conditions affect their mechanical properties and behavior. After preparing the required specimens (using two different adhesives and adherend thicknesses), several tests have been made in order to determine these properties and compare their values to the predictions made using analytical methods. These tests were repeated with several distinct combinations of temperatures and moisture levels so that the effect of these properties can be properly interpreted. It was observed that higher temperatures strongly increase the ductility of the adhesive but mixed with moisture this can degrade them. Moisture can increase the energy absorbed through increased plastic deformation of the adhesive and improve behaviour at low temperatures.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Heat Transfer Improvement in an Open Cubic Cavity using a Hybrid Nanofluid150115131605710.22055/jacm.2020.35213.2598ENSalihaBelloutLEAP Laboratory, Department of Mechanical Engineering, Faculty of Sciences Technology, University of Brothers Mentouri-Constantine 1, Route de Ain El. Bey, 25000 Constantine, Algeria0000-0002-2181-9141RachidBessaïhLEAP Laboratory, Department of Mechanical Engineering, Faculty of Sciences Technology, University of Brothers Mentouri-Constantine 1, Route de Ain El. Bey, 25000 Constantine, Algeria0000-0002-0764-7731Journal Article20200929Numerical simulation of convection heat transfer and entropy generation in an open cubic cavity filled with a hybrid nanofluid is carried out. This configuration is heated uniformly by a constant volumetric heat source q<sub>v</sub>. All the walls are adiabatic. The hybrid nanofluid flow (Al<sub>2</sub>O<sub>3</sub>-Cu/water) penetrates in the cavity at a uniform velocity U<sub>0</sub> and a temperature T<sub>0</sub>. To solve the mathematical equations, we used Ansys-Fluent 14.5 software. Results are validated with other works found in the literature. We present our results in terms of streamlines, isotherms, velocity, temperature, local and average Nusselt numbers profiles, and entropy generation for the Reynolds number (300 < Re < 700), the solid volume fraction (0 < φ < 0.08), and heat source location (1cm < d < 3cm). Results indicate that by increasing Re, φ and d<sub>h</sub>, the heat transfer is improved. Moreover, nanohybrid gives better heat transfer than nanofluid, and the use of nanoparticles contributes to the minimization of entropy generation. Compared with the vertical location of the heat source, the horizontal location gives an increase in heat transfer. The Nu<sub>av</sub> correlations are determined for the nanofluid and hybrid nanofluid. This study may help to enhance the heat transfer of electronic equipment.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701On Stokes' Second Problem for Burgers' Fluid over a Plane Wall151415261604410.22055/jacm.2020.35227.2603ENSafiaAkramMCS, National University of Sciences and Technology, Islamabad, Pakistan0000-0001-6288-6095AsiaAnjumDepartment of Computer Science, National University of Modern Languages, H-9, Islamabad, Pakistan0000-0002-1651-7851MasoodKhanDepartment of Mathematics, Quaid-i-Azam University, Islamabad 44000, Pakistan0000-0003-4404-1456AnwarHussainDepartment of Mathematics, Quaid-i-Azam University, Islamabad 44000, Pakistan0000-0002-2571-8997Journal Article20200929The Stokes' second problem for a Burgers' fluid over a plane wall is considered in this paper. The motion of the fluid is induced by the oscillation of the plane wall between two side walls perpendicular to the plane wall. The exact analytical solutions for the velocity field and the adequate shear stress are established in simple forms by means of integral transforms. The solutions that have been obtained, presented as a sum of the steady and the transient solutions, satisfy all imposed initial and boundary conditions. In the absence of the side walls they reduce to the similar solutions over an infinite plate. Finally, the results for the velocity, as well as a comparison between models, are displayed graphically for pertinent parameters to show interesting aspects of the solutions. It is observed that the velocity and the boundary layer thickness were observed to be decreased in the presence of the side walls. Moreover, the Maxwell fluid was observed to be the speediest and the Newtonian was the slowest.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701On Integrability up to the Boundary of the Weak Solutions to a Class of non-Newtonian Compressible Fluids with Vacuum152715361658210.22055/jacm.2021.35444.2655ENJanMuhammadSchool of Mathematics, Department of Mathematics, CNS, Northwest University, Xi’an, 710127, P. R. China0000-0001-8872-2523AbdulSamadSchool of Mathematics, Department of Mathematics, Northwest University, Xi’an, 710127, P. R. China0000-0002-0887-9860Journal Article20201015In this paper, we study the integrability up to the boundary of the weak solutions of non-Newtonian compressible fluid with a nonlinear constitutive equation in ℝ<sup>3</sup> bounded domain. Galerkin approximation will be used for existence of weak solutions and by applying the bounded linear operator B, introduced by Bogovskii, we prove the square integrability of the density up to the boundary.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Nonlinear Vibration of an Electrostatically Actuated Functionally Graded Microbeam under Longitudinal Magnetic Field153715491660710.22055/jacm.2021.35504.2670ENDang VanHieuTNU, Thai Nguyen University of Technology (TNUT), Thainguyen, Vietnam0000-0001-8789-7284Nguyen ThiHoaTNU, Thai Nguyen University of Technology (TNUT), Thainguyen, Vietnam0000-0003-3658-5998Le QuangDuyTNU, Thai Nguyen University of Technology (TNUT), Thainguyen, Vietnam0000-0001-5013-7150Nguyen ThiKim ThoaTNU, Thai Nguyen University of Technology (TNUT), Thainguyen, Vietnam0000-0003-4876-8169Journal Article20201022In this work, we develop a model of an electrostatically actuated functionally graded (FG) microbeam under a longitudinal magnetic field based on the Euler-Bernoulli beam and nonlocal strain gradient theories to investigate the nonlinear vibration problem. The FG microbeam is placed between two electrodes, a DC voltage applied between the two fixed electrodes causes an electrostatic force to be exerted on the FG microbeam. The FG microbeam is composed of metal and ceramic in which the properties of these materials are assumed to change in the thickness direction according to the simple power-law distribution. The Galerkin method and the Hamiltonian Approach are employed to find the approximate frequency of the FG microbeam. The accuracy of the present solution is verified by comparing the obtained results with the numerical results and the published results in the literature. Effects of the power-law index, the material length scale parameter, the nonlocal parameter, the applied voltage and the magnetic force on the nonlinear vibration behaviour of the FG microbeam are studied and discussed.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701A Novel Approach to Compute the Numerical Solution of Variable Coefficient Fractional Burgers' Equation with Delay155015641648710.22055/jacm.2021.35574.2689ENAmit KVermaDepartment of Mathematics, Indian Institute of Technology Patna, Patna–801106, Bihar, India0000-0001-8768-094XMukesh KumarRawaniDepartment of Mathematics, Indian Institute of Technology Patna, Patna–801106, Bihar, IndiaRavi P.AgarwalDepartment of Mathematics, Texas A&M University-Kingville, Kingsville, TX, USA0000-0003-0634-2370Journal Article20201029In this article, we come up with a novel numerical scheme based on Haar wavelet (HW) along with nonstandard finite difference (NSFD) scheme to solve time-fractional Burgers’ equation with variable diffusion coefficient and time delay. In the solution process, we discretize the fractional time derivative by NSFD formula and spatial derivative by HWs series expansion. We use the quasilinearisation process to linearize the nonlinear term. Also, the convergence of the scheme is discussed. The efficiency and correctness of the proposed scheme are assessed by L∞-error and L2 -error norms.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Dynamically Consistent NSFD Methods for Predator-prey System156515741657710.22055/jacm.2021.36537.2862ENAliShokriDepartment of Mathematics, Faculty of Science University of Maragheh, Maragheh, Iran0000-0003-2699-1490MohammadMehdizadeh KhalsaraeiDepartment of Mathematics, Faculty of Science University of Maragheh, Maragheh, Iran0000-0002-0130-0449MaryamMolayiDepartment of Mathematics, Faculty of Science University of Maragheh, Maragheh, Iran0000-0001-6395-9161Journal Article20210102In this paper, we introduce two nonstandard finite difference (NSFD) methods for solving the mathematical model of the Rosenzweig-MacArthur predator-prey system. These new proposed numerical methods have important features such as positivity and elementary stability. Numerical comparisons between the proposed methods and the other methods such as second-order and forth order Runge-Kutta methods (we refer them RK2 and RK4, respectively), Euler method, and NSFD method presented in [6] indicate that the new methods have better accuracy and convergence.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701On the Effect of the End-effector Point Trajectory on the Joint Jerk of the Redundant Manipulators157515821666010.22055/jacm.2021.35350.2635ENXuan BienDuongAdvanced Technology Center, Le Quy Don Technical University, Hoang Quoc Viet, Hanoi, 84, Vietnam0000-0001-7380-8551Journal Article20201009<span>This paper is focused on investigating the joints jerk of industrial serial redundant manipulators with 6 degrees of freedom (6-DOF) under the variation of the end-effector point (EEP) trajectory in the workspace. The EEP trajectories are initially built in the basic planes because of their simplicity, verification and experimentation are smooth, and most of the actual welded structures are performed on these basic planes. The jerk is determined by solving the inverse kinematics problem of the redundant system. This problem is solved based on the algorithm which is used for adjusting the increments of the generalized coordinate vector (AGV). The efficiency of this algorithm is shown through the error between a given trajectory and the recalculated trajectory through the forward kinematics problem. The result of this study allows us to evaluate the effect of the change of the trajectory on the kinematics characteristics of the robot in general and the jerk of the joints in particular. On the other hand, these results can be used as the basis for planning the EEP trajectory for redundant robots, developing algorithms to reduce joint jerky, increase the life of robot systems, and improve the accuracy of the redundant robot movement.</span>Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Nonlinear Multiscale Modelling and Design using Gaussian Processes158315921667810.22055/jacm.2021.36106.2795ENSumuduHerathDepartment of Civil Engineering, University of Moratuwa, Sri Lanka0000-0002-4462-4902UdithHaputhanthriDepartment of Electronic and Telecommunication Engineering, University of Moratuwa, Sri Lanka0000-0003-1609-3822Journal Article20201220A method for nonlinear material modeling and design using statistical learning is proposed to assist in the mechanical analysis of structural materials. Conventional computational homogenization schemes are proven to underperform in analyzing the complex nonlinear behavior of such microstructures with finite deformations. Also, the higher computational cost of the existing homogenization schemes inspires the inception of a data-driven multiscale computational homogenization scheme. In this paper, a statistical nonlinear homogenization scheme is discussed to mitigate these issues using the Gaussian Process Regression technique. A data-driven model is trained for different strain states of microscale unit cells. In the macroscale, nonlinear response of the macroscopic structure is analyzed, for which the stresses and material responses are predicted by the trained surrogate model.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Free Vibration Analysis of Functionally Graded Porous Nano-plates with Different Shapes Resting on Elastic Foundation159316051666110.22055/jacm.2021.36181.2807ENTrac LuatDoanFaculty of Mechanical Engineering, Le Quy Don Technical University, Hanoi, Vietnam0000-0003-1531-5163Pham BinhLeFaculty of Mechanical Engineering, Le Quy Don Technical University, Hanoi, Vietnam0000-0003-2023-6275Trung ThanhTranFaculty of Mechanical Engineering, Le Quy Don Technical University, Hanoi, Vietnam0000-0002-1730-0590Vu KhacTraiAutomobile Enterprise, 751 One Member Limited Liability Company, Ho Chi Minh City, Vietnam0000-0003-3693-6022Quoc HoaPhamDepartment of Training, Tran Dai Nghia University, Ho Chi Minh City, Vietnam0000-0001-6316-3120Journal Article20201225This paper proposes a finite element method (FEM) based on a nonlocal theory for analyzing the free vibration of the functionally graded porous (FGP) nano-plate with different shapes lying on the elastic foundation (EF). The FGP materials with two-parameter are the power-law index (<em>k</em>) and the porosity volume fraction (ξ) in two cases of even and uneven porosity. The EF includes Winkler stiffness (<em>k</em><sub>1</sub>) and Pasternak stiffness (<em>k</em><sub>2</sub>). Some numerical results in our work are compared with other published to verify accuracy and reliability. Moreover, the influence of geometric parameters, materials on the free vibration of the FGP nano-plates resting on the EF is comprehensively investigated.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Description of Anomalous Behavior of Aluminum Alloys with Hill48 Yield Criterion by Using Different Experimental Inputs and Weight Coefficients160616191665110.22055/jacm.2021.36297.2821ENBoraSenerFaculty of Mechanical Engineering, Department of Mechanical Engineering, Yildiz Technical University, Besiktas, Istanbul, 34349, Turkey0000-0002-8237-1950Journal Article20210105The anomalous behavior of aluminum alloys is modeled with quadratic Hill48 yield criterion in this study. An identification method based on minimization of the error function is applied and the effect of the number of experimental input and weight coefficients used in the identification are investigated. Two highly anisotropic aluminum alloys (AA2090-T3 and AA5182-O) are selected in the study. Firstly, Hill48 parameters are determined with four different experimental data set, then the effect of the weight coefficients for each set is investigated. In-plane variations of plastic properties and yield surfaces of the materials are predicted with determined Hill48 parameters and the most appropriate pair (experimental data set and weight coefficient) are selected by comparison of the predicted results with experiment. Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Investigation of an Unsteady and Anisotropic Turbulent Flow Downstream a 90° Bend Pipe with and without Ribs162016381667110.22055/jacm.2021.36399.2837ENRachidChiremselDepartment of Hydraulics, University of Batna 2, Research Laboratory in Applied Hydraulics, Constantine road N°53.Fesdis, Batna, 05078, Algeria0000-0002-3048-7206AliFourarDepartment of Hydraulics, University of Batna 2, Research Laboratory in Applied Hydraulics, Constantine road N°53.Fesdis, Batna, 05078, Algeria0000-0003-2612-3589FawazMassouhNational Higher School of Arts and Crafts (ParisTech; ENSAM), Laboratory of Fluid Mechanics, France0000-0002-0718-5429ZakaryaChiremselSafety Department, University of Batna 2, IHSI-LRPI, Constantine road N°53.Fesdis, Batna, 05078, Algeria0000-0002-7236-137xJournal Article20210115In this work, a numerical study of the dynamical behavior of unsteady and anisotropic turbulent flow downstream a 90° bended pipe was presented. For this purpose, comparative computations are carried out employing two flow configurations, bend pipe with ribs and bend pipe without ribs with a curvature radius ratio Rc/D=2.0. In the bend pipe with ribs, the pitch ratios Pt/e=40 and the rib height to pipe diameter e/D is 0.1. This model has been utilized to assess the effect of ribs on flow where the presence of the ribs leads to a complex velocity field with regions of flow separation upstream and downstream of the ribs. The Reynolds-Averaged Navier–Stokes (RANS) approach is employed and the computational model is validated by comparisons with the existing experimental data. The simulations are conducted with the commercials CFD software FLUENT for Dean number varying from 5000 to 40000. The result analysis shows that the higher resistance generated by the ribs produced relatively larger velocity gradient (∂<em>U</em>/∂<em>y</em>) compared to the case of bend pipe without ribs where a more uniform mean velocity profile is observed. The turbulence intensities are higher in the ribbed bend pipe compared to those in the non-ribbed case and depend faintly on the Dean number. The levels of the Reynolds shear stresses are significantly enhanced by the ribs compared to the case without ribs. This increasing is explained by significantly higher levels of turbulence production over those ribs produced by large values of ∂<em>U</em>/∂<em>y</em>.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Variational Principles and Solitary Wave Solutions of Generalized Nonlinear Schrödinger Equation in the Ocean163916481669610.22055/jacm.2021.36690.2890ENMeng-ZhuLiuCollege of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, China0000-0001-6699-3496Xiao-QunCaoCollege of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, ChinaCollege of Computer, National University of Defense Technology, Changsha 410073, China0000-0002-6135-0712Xiao-QianZhuCollege of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, ChinaCollege of Computer, National University of Defense Technology, Changsha 410073, China0000-0002-5769-9292Bai-NianLiuCollege of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, ChinaCollege of Computer, National University of Defense Technology, Changsha 410073, China0000-0002-2742-1265Ke-ChengPengCollege of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, China0000-0002-3800-271XJournal Article20210220Internal solitary waves are very common physical phenomena in the ocean, which play an important role in the transport of marine matter, momentum and energy. Because the generalized nonlinear Schrödinger equation can well explain the effects of nonlinearity and dispersion in the ocean, it is more suitable for describing the deep-sea internal wave propagation and evolution than other mathematical models. At first, by designing skillfully the trial-Lagrange functional, different kinds of variational principles are successfully established for a generalized nonlinear Schrödinger equation by the semi-inverse method. Then, the constructed variational principles are proved correct by minimizing the functionals with the calculus of variations. Furthermore, some kinds of internal solitary wave solutions are obtained and demonstrated by semi-inverse variational principle for the generalized nonlinear Schrödinger equation.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Stability Analysis of Articulated Bus in Straight-ahead Running Manoeuvre164916621670210.22055/jacm.2021.36566.2869ENAlessandroDe FeliceDepartment of Engineering Enzo Ferrari, University of Modena and Reggio Emilia, Via Vivarelli 10, Modena, 41125, Italy0000-0001-9142-1334MatteoMercantiniProtezioni Elaborazioni Industriali (PEI) Srl, Via Torretta 32, Calderara di Reno, 40012, Italy0000-0002-8626-0836SilvioSorrentinoDepartment of Engineering Enzo Ferrari, University of Modena and Reggio Emilia, Via Vivarelli 10, Modena, 41125, Italy0000-0003-2136-5538Journal Article20210205A comprehensive study on the stability of a planar linearized single-track model of a two-section pusher articulated bus is presented with the aid of a complete set of stability maps. The two sections of the vehicle model are connected at the hitch point by a revolute joint; an equivalent visco-elastic characteristic function describes its rotational visco-elastic properties, playing a major role in stability control and therefore in passive safety. The equations of motion are derived in analytical form, allowing easy implementation of the non-linear model (eventually including a non-linear viscoelastic characteristic functions of the joint). Stability of the linearized model is then studied in equilibrium configurations by means of sensitivity analysis with respect to the model’s governing parameters. Stability maps are drawn on the basis of sets of parameter values related to straight-ahead running, steady-state manoeuvres. The most important parameters controlling the onset of unstable motions are identified, paying attention to the role played by the equivalent rotational damping coefficient and the equivalent torsional stiffness characterizing the connection joint, with the aim of finding criteria for its design.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Finite Element Modelling and Simulation of the Hysteretic Behaviour of Single- and Bi-metal Cantilever Beams using a Modified Non-linear Beta-damping Model166316751670910.22055/jacm.2021.35420.2651ENHamza BinTariqDepartment of Mechanical and Aerospace Engineering, Nazarbayev University, Qabanbay Batyr 53, Nur-Sultan, 010000, Kazakhstan0000-0002-2166-7897CharlesRajakumarDepartment of Mechanical and Aerospace Engineering, Nazarbayev University, Qabanbay Batyr 53, Nur-Sultan, 010000, KazakhstanDichuanZhangDepartment of Civil and Environmental Engineering, Nazarbayev University, Qabanbay Batyr 53, Nur-Sultan, 010000, Kazakhstan0000-0002-9253-4178ChristosSpitasDepartment of Mechanical and Aerospace Engineering, Nazarbayev University, Qabanbay Batyr 53, Nur-Sultan, 010000, KazakhstanJournal Article20201014This paper explores a novel non-linear hysteresis model obtained from the modification of the conventional Kelvin-Voigt model, to produce a non-viscous hysteretic behaviour that is closer to metal damping. Two case studies are carried out for a vibrating cantilever beam under tip loading (bending), the first considering a single uniform material and the second considering a bimetallic structure. The damping behaviour is studied in the frequency domain (constant damping ratio model vs. Kelvin-Voigt/ beta damping model) and time-domain (proposed modified hysteresis model vs. Kelvin-Voigt/ beta damping model). In the frequency domain, it was found that the Kelvin-Voigt model essentially damps out the displacement response of the modes more than the constant damping ratio model does. In the transient analysis, the Kelvin-Voigt model likewise produced unnaturally rapid damping of the oscillations for both the single- and bi-metal beam, compared to the modified hysteretic damping model, which produced a damping behaviour closer to actual metal behaviour. This was consistent with results obtained in the frequency domain.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701A Modified Thermoelastic Fractional Heat Conduction Model with a Single-Lag and Two Different Fractional-Orders167616861563010.22055/jacm.2020.33790.2287ENAhmed E.AbouelregalDepartment of Mathematics, College of Science and Arts, Jouf University, Al-Qurayyat, Saudi ArabiaDepartment of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt0000-0003-3363-7924HijazAhmadDepartment of Basic Sciences, University of Engineering and Technology Peshawar, Pakistan0000-0002-5438-5407Journal Article20200531<span>Recently, fractional calculus theory has been successfully employed in generalized thermoelasticity theory and several models with fractional order have been introduced. In this work, a fractional thermoelastic modified Fourier's Law with phase lag and two different fractional-orders has been constructed. The previous fractional models of thermoelasticity introduced by Sherief et al. </span><span>[1]</span><span>, Ezzat </span><span>[2]</span><span> and Lord and Shulman </span><span>[3]</span><span> as well as classical coupled thermoelasticity </span><span>[4]</span><span> follow as limiting cases. This proposed model is applied to an infinitely annular cylinder that is subject to time-dependent surface temperatures, and whose surfaces are free of traction. The method of the Laplace transform is employed to get the solutions of the field variables. A numerical technique is utilized to invert the Laplace transforms. Some results are presented in tables and figures to extract the effects of fractional order parameters on all variables studied. The theory's predictions have been checked and compared to previous models.</span>Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701An Automatic Program of Generation of Equation of Motion and Dynamic Analysis for Multi-body Mechanical System using GNU Octave168716971568910.22055/jacm.2020.33826.2293ENYonghuiParkDepartment of Mechanical Engineering, Yuhan University, 590, Gyeongin-ro, Bucheon-si, Gyeonggi-do, Republic of Korea0000-0001-6716-6935Journal Article20200602Multi-body dynamics is used to calculate the physical quantities required for component design, such as calculating the dynamic response of mechanical components and the time history of dynamic loads. Advances in analysis software, including DADS, ADAMS, RecurDyn, and DAFUL, have made it possible to easily calculate dynamic responses by defining relationships between components and operating environments from 3D modeling on user-created components. However, when the understating of dynamic analysis is lacking, it is difficult to apply multi-body dynamics analysis in the design process, and it is difficult to analyze the acquired response data. In this study, we developed an automatic code to derive equations of motion in the matrix format and calculate dynamic responses of multi-body systems using GNU Octave, a free high level language. In particular, the process of defining matrices and vectors such as inertia matrix, stiffness matrix, and external force vector concerning the degrees of freedom of components by using Euler-Lagrange equations is shown to understand the structure and process of dynamic analysis. The code application by explaining how to use the code in a different mechanical system is also shown to help understand the usage method for who wants to study Multi-body dynamics.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Stress Mode Superposition for a Priori Detection of Highly Stressed Areas: Mode Normalisation and Loading Influence169817091674210.22055/jacm.2021.36637.2878ENCarstenStrzalkaTechnical University Berlin, Institute of Mechanics, TU-Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany0000-0002-8490-5899DraganMarinkovicTechnical University Berlin, Institute of Mechanics, TU-Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany0000-0002-3583-9434Manfred W.ZehnTechnical University Berlin, Institute of Mechanics, TU-Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany0000-0002-2412-3676Journal Article20210212From the economic and technical point of view, the reduction of development periods and required resources represent a considerable benefit. For the reduction of numerical effort and processed data in numerical stress analysis, the present paper is focused onto the investigation of an efficient method for the <em>a priori</em> detection of a structural component’s highly stressed areas. Based on the theory of stress mode superposition and the frequency domain solution of the decoupled equations of motion, an analytically consistent approach for <em>a priori</em> mode superposition is presented. In this context, the influence of multiaxial loading and mode normalisation is investigated. Validation is performed on a simplified industrial model of a twist-beam rear axle.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701The Dynamics of the Working Body of the Tubular Conveyor with the Chain Drive171017181675910.22055/jacm.2021.35725.2719ENOlegLyashukTernopil Ivan Puluj National Technical University, 56 Ruska str., Ternopil, 46001 Ukraine0000-0003-4881-8568BogdanSokilHetman Petro Sahaidachnyi National Army Academy, 32 Heroes of Maidan street, Lviv, 79026 Ukraine0000-0002-2839-5273RomanHevkoTernopil Ivan Puluj National Technical University, 56 Ruska str., Ternopil, 46001 Ukraine,0000-0001-8536-4599ViktorAulinCentral Ukrainian National Technical University, 8, Prospekt Universytetskyi, Kropyvnytskyi, 25006 Ukraine0000-0003-2737-120XLeonidSerilkoNational University of Water and Environmental Engineering, 11 Soborna str., 33028 Rivne, Ukraine0000-0002-2396-561XYuriyVovkTernopil Ivan Puluj National Technical University, 56 Ruska str., Ternopil, 46001 Ukraine0000-0001-8983-2580DmytroSerilkoNational University of Water and Environmental Engineering, 11 Soborna str., 33028 Rivne, Ukraine0000-0001-5881-2413AndriyDovbyshTernopil Ivan Puluj National Technical University, 56 Ruska str., Ternopil, 46001 Ukraine0000-0001-9683-9617Journal Article20201112The theoretical calculations of the dynamics of the conveyor drive chain oscillations for different speeds of bulk material movement as a result of motion internal perturbations have been presented in the article. Resonant oscillations have been studied. It has been established that the amplitude of the transition through the resonance is greater for higher speeds of the conveyor drive chain and the maximum dynamic tension in the chain increases with increasing speed of transportation of bulk material. The dependences of the natural frequency of the system "drive chain of the conveyor line - grain" on the parameters of the system and the amplitude, as well as the amplitude of the resonant oscillations of the system on the speed of grain movement have been obtained. The dependence of the resonance amplitude of the system "conveyor chain drive - grain mass" on the speed of the drive chain at certain parameters has been determined. Taking into account the obtained theoretical data, an improved construction of a tubular chain conveyor with various working bodies and stand equipment using the Altivar 71 frequency converter for complex tasks of the electric drive from 0.75 to 630 kW has been developed. The experimental researches have been carried out and the dependences for definition of productivity and a rotation moment at transportation by the tubular scraper conveyor on curvilinear routes for loose material (wheat and peas) have been received.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Alternative Integration Approaches in the Weight Function Method for Crack Problems171917251677210.22055/jacm.2021.37137.2968ENMartin A.EderDepartment of Wind Energy, Technical University of Denmark, Frederiksborgvej 399, 4000 Roskilde, Denmark0000-0002-5306-365XXiaoChenDepartment of Wind Energy, Technical University of Denmark, Frederiksborgvej 399, 4000 Roskilde, Denmark0000-0001-6726-4068Journal Article20210412This study proposes two alternative approaches to complement existing integration strategies used in the weight function method for linear elastic crack problems. The first approach is based on an interpolation type integration scheme and the second approach is based on Gauss quadrature. The proposed approaches enable a computationally efficient numerical integration for computing stress intensity factors in 2D fracture problems. The efficiency is gained through a comparatively low number of integration points facilitated by higher-order approximation. The integration weights only need to be computed once for a given crack length-to-width ratio and can be applied to arbitrary continuous and smooth stress distributions. The proposed approaches show excellent accuracy. In particular, the Gauss quadrature approach exhibits several orders of magnitude more accuracy compared to the most commonly used trapezoidal integration.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Regularization of the Movement of a Material Point Along a Flat Trajectory: Application to Robotics Problems172617361684710.22055/jacm.2021.37193.2990ENDinaraAzimovaL.N. Gumilyov Eurasian National University, 2, Satpaev Str., Nur-Sultan, 010008, Kazakhstan0000-0002-3721-9001BalgaishaMukanovaInstitute of Information and Computational Technologies,125, Pushkin str., Almaty, 050000, Kazakhstan0000-0002-0823-6451MaxatAkhmetzhanovAl Farabi Kazakh National University, 71/23, al-Farabi Ave., Almaty, 050040, Kazakhstan0000-0001-7890-5422Journal Article20210427A control problem of the robot’s end-effector movement along a predefined trajectory is considered. The aim is to reduce the work against resistance forces and improve the comfortability of the motion. The integral of kinetic energy and weighted inertia forces for the whole period of motion is introduced as a cost functional. By applying variational methods, the problem is reduced to a system of quasilinear ordinary differential equations of the fourth order. Numerical examples of solving the problem for movement along straight, circular and elliptical trajectories are presented. For the sake of clarity, the problem is studied for a specific kind of a 3D printer in the 2DoF approximation. However, in the case of negligible masses of moving elements compared the mass of an end-effector, the solution is universal, i.e., it remains the same for given trajectories.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701An Experimental and Numerical Study on the Aerodynamic Performance of Vibrating Wind Turbine Blade with Frequency-Domain Method173717501684810.22055/jacm.2021.37406.3011ENShine WinNaungDepartment of Mechanical and Construction Engineering, Northumbria University, Newcastle, Upon Tyne, NE1 8ST, United Kingdom0000-0001-8441-0725Mahdi ErfanianNakhchiDepartment of Mechanical and Construction Engineering, Northumbria University, Newcastle, Upon Tyne NE1 8ST, United Kingdom0000-0003-1034-6091MohammadRahmatiDepartment of Mechanical and Construction Engineering, Northumbria University, Newcastle Upon Tyne NE1 8ST, United Kingdom0000-0003-4903-5370Journal Article20210511A highly efficient nonlinear frequency-domain solution method is proposed and employed to investigate the aerodynamic and aeromechanical performances of an oscillating wind turbine blade aerofoil in this study. Extensive validations of a frequency-domain method against an experiment as well as a typical time-domain solution method are provided in this paper. An experiment is also designed and conducted to measure pressure distributions over an aerofoil as well as to validate the numerical model. Unsteady pressure distributions and aeroelasticity parameters of the oscillating NACA0012 aerofoil are computed at various angles of attack and Reynolds numbers. Results indicate that the difference of unsteady pressure distributions between the two surfaces of the aerofoil becomes larger as the angle of attack is increased, whereas the flow separation on the suction surface is reduced by raising the Reynolds number. The turbulent flow develops in the downstream region due to the laminar vortex shedding at lower Reynolds numbers. It is also revealed that the Reynolds number has an impact on the aeroelasticity, and the aerodynamic damping value is larger at higher Reynolds numbers. The comparison between the frequency-domain method and the time-domain method shows that the frequency-domain method is not only accurate but also computationally very efficient as the computation time is reduced by 90%.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Weighted Dual Approach to an Equivalent Stiffness-based Load Transfer Model for Jacked Open-ended Pile175117631687510.22055/jacm.2021.37430.3013ENNguyen NgocLinhDepartment of Mechanical Engineering, Thuyloi University (TLU), Hanoi, Vietnam0000-0002-6724-462XNgoc AnhNguyenDepartment of Mechanical Engineering, University of Transport and Communications (UTC), Hanoi, Vietnam0000-0001-7322-3940Kuu VanNguyenDepartment of Mechanical Engineering, Thuyloi University (TLU), Hanoi, Vietnam0000-0002-2073-9596Diem DangNguyenDepartment of Mechanical Engineering, University of Transport and Technology (UTT), Hanoi, Vietnam0000-0002-1295-8238Journal Article20210514This paper presents a new equivalent stiffness-based load transfer model for an open-ended pipe pile. The main idea of this model is to replace the sum of unit stiffnesses corresponding with external and internal unit skin frictions in the basic differential equation of load transfer by a weighted average of equivalent unit stiffnesses using a dual approach of equivalent replacement. The contribution of external and internal skin frictions to equivalent unit stiffnesses is evaluated by normalized dimensionless weighting coefficients in the form of average value with the penetration depth. Application of new load transfer model to a jacked open-ended pile concerning semi-empirical models of external and internal unit skin frictions leads to corresponding explicit expressions of weighting coefficient. A computational example of a jacked open-ended pile is carried out. It is shown that the proposed equivalent stiffness-based load transfer model is an effective tool for analyzing behaviors of the open-ended pile in considering the soil plugging effect.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Analysis and Recognition of Standards in Intelligent Hybrid Systems using Natural Computing176417731690110.22055/jacm.2021.37798.3089ENRodrigo Francisco BorgesLourençoUNESP - Univ. Estadual Paulista, Faculty of Engineering of Ilha Solteira, Department of Mechanical Engineering, Ilha Solteira, 15385-000, Brasil0000-0002-7772-6438RobertoOutaFaculty of Technology of Araçatuba, Department of Biofuels, Av. Prestes Maia, 1764 - Ipanema, Araçatuba, 16052-045, Brasil0000-0002-8649-1722Fábio RobertoChavaretteUNESP - Instituto de Química, Departamento de Engenharia, Física e Matemática, Rua Prof. Francisco Degni, 55 – Quitandinha, Araraquara, 14800-060, Brasil0000-0001-5376-3392Aparecido CarlosGonçalvesUNESP - Univ. Estadual Paulista, Faculty of Engineering of Ilha Solteira, Department of Mechanical Engineering, Ilha Solteira, 15385-000, Brasil0000-0001-5376-3392Journal Article20210621This work shows the application of one of the techniques of bioengineering, the perceptron network in the detection of system failures, and also allows the use of the perceptron network technique in choosing the location of the best sensor to be used in the dynamic system. The application of the perceptron network was adopted because it is considered the best binary linear classifier. This work is considered multidisciplinary and difficult to develop. The final result demonstrates a severe application of pre-processing and processing, until the classification and grouping of signals in the two phases of the work. Through the results found, this work can be considered successful and can be applied in several areas of engineering for structural analysis.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Parametric Optimization of a Cyclogiro Aircraft Design for Efficient Hover with Aeroelastic Considerations177417871670410.22055/jacm.2021.35872.2752ENLouisGagnonInstitute of Aerodynamics and Gas Dynamics, Universität Stuttgart, 70569, Germany0000-0002-1619-3236MarcoMorandiniDepartment of Aerospace Science and Technology, Politecnico di Milano, 20156, Italy0000-0002-1992-9077StéphaneFournierDepartment of Aerospace Engineering, University of Bristol, BS8 1TR, United KingdomJournal Article20201127A minimization procedure is proposed to orient the design of a vertical take-off and landing drone towards sustainability. The vehicle is a novel cycloidal rotor drone and the principal objective is to yield the best ratio of payload to power consumption. The drone blades, rotor arms, and frame are designed for fused deposition modeling additive manufacturing with polylactic acid. 10 variables for the geometry, operation parameters, and material infill percentages are explored in search of the optimum design. A special derivation procedure allows obtaining the symbolic equations for the weight and power consumption of the drone. This permits optimization with a hybrid genetic and gradient method and exploring a broad range of aircraft sizes. 7 constraint equations ensure that the necessary assumptions made for the derivation remain valid and that the structural strength is adequate. For each new configuration, this method allows to quickly find a new optimum design using a desktop computer. Also, modifying the constraints, variables, or objective function is straightforward. Finally, the resulting design has a power loading of 0.0876 N/W.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701A General Purpose Variational Formulation for Boundary Value Problems of Orders Greater than Two178818021690210.22055/jacm.2021.37244.2987ENXuefengLiDepartment of Mathematics and Computer Science, Loyola University, New Orleans, LA 70118, USA0000-0002-5462-6274Journal Article20210424We develop a new general purpose variational formulation, particularly suitable for solving boundary value problems of orders greater than two. The functional related to this variational formulation requires only <em>Η</em><sup>1</sup> regularity in order to be well-defined. Using the finite element method based on this new formulation thus becomes simple even for domains in dimensions greater than one. We prove that a saddle-point solution to the new variational formulation is a weak solution to the associated boundary value problem. We also prove the convergence of the numerical methods used to find approximate solutions to the new formulation. We provide numerical tests to demonstrate the efficacy of this new paradigm.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Performance measure and tool for benchmarking metaheuristic optimization algorithms180318131690710.22055/jacm.2021.37664.3060ENFrançoisSchottPercipio Robotics, Maison des Microtechniques 18, rue Alain Savary, Besançon, France0000-0003-3744-1988DominiqueChamoretICB UMR 6303, CNRS, UBFC, UTBM, Belfort, France0000-0002-9771-774XThomasBaronFEMTO-ST institute, Univ. Bourgogne Franche-Comté, CNRS, ENSMM Time and frequency dept., Besançon, France0000-0003-1519-6539SébastienSalmonMy-OCCS, Besançon, France0000-0002-1714-3513YannMeyerUniv. Savoie Mont Blanc, SYMME, FR-74000 Annecy, France0000-0001-9412-6752Journal Article20210608In the last decade, many new algorithms have been proposed to solve optimization problems. Most of them are meta-heuristic algorithms. The issue of accurate performance measure of algorithms is still under discussion in the scientific community. Therefore, a new scoring strategy via a new benchmark is proposed. The idea of this new tool is to determine a score, a measure of efficiency taking into account both the end value of the optimization and the convergence speed. This measure is based on an aggregate of statistical results of different optimization problems. These problems are judiciously chosen to cover as broad a spectrum of resolution configurations as possible. They are defined by combinations of several parameters: dimensions, objective functions and evaluation limit on dimension ratios. Aggregation methods are chosen and set in order to make the problem weight relevant according to the computed score. This scoring strategy is compared to the CEC one thanks to the results of different algorithms: PSO, CMAES, Genetic Algorithm, Cuttlefish and simulated annealing.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Two-Solid Deposition in Fluid Column using Immersed Boundary-Lattice Boltzmann Method181418251690910.22055/jacm.2021.37140.2969ENM. RizqieArbiePhysics of Earth and Complex Systems Division, Department of Physics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung, 40132, Indonesia0000-0002-1359-8377UmarFauziPhysics of Earth and Complex Systems Division, Department of Physics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung, 40132, Indonesia0000-0002-2221-7257Fourier D.E.LatiefPhysics of Earth and Complex Systems Division, Department of Physics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung, 40132, Indonesia0000-0002-5448-929XEnjang J.MustopaPhysics of Earth and Complex Systems Division, Department of Physics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung, 40132, Indonesia0000-0002-0769-2164Journal Article20210413Solid deposition in fluid may involve solids with different density and size and may happen in quiescent fluid or rather in counter flow. We perform a numerical investigation on the role of density-ratios, size-ratio, and initial configuration on the settling of two circular solids in a fluid channel with or without counter-flow. Through this study, we show how settling dynamics of two solids can be controlled. Numerical experiment based on a coupled Immersed Boundary-Lattice Boltzmann is employed. It is shown that certain parameter set leads to guided deposition while denser solid leaves the less dense one as time progressing. However, certain parameter set leads to periodic close encounters which is robust in the presence of Poiseuille-like counter-flow. In this case, the separation between two solids is bounded during the deposition.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Numerical Investigation of Tissue-Temperature Controlled System in Thermal Ablation: A Finite Element Approach182618351691110.22055/jacm.2021.36020.2819ENMridulSannyalDepartment of Industrial and Production Engineering, Shahjalal University of Science and Technology, Sylhet, 3114, Bangladesh0000-0001-5040-5549Abul Mukid MohammadMukaddesDepartment of Industrial and Production Engineering, Shahjalal University of Science and Technology, Sylhet, 3114, Bangladesh0000-0001-7965-4602Journal Article20210104In thermal ablation, several techniques of treating infected cell in human tissue are being used by the physicians. Transferring heat to the infected cell is one of them. The purpose of this research is to investigate the tissue-temperature controlled system in thermal ablation and compare with two different point heating processes, namely constant and step heating. For this purpose, the finite element model of Penne’s bio-heat equation has been developed to measure the temperature within the two-dimensional tissue model embedded with a small tumor. The tissue temperature-controlled heating was designed to restrict the healthy tissue temperature below the damage threshold temperature. Using the temperature profile, tissue damage index was measured with the help of Arrhenius rate equation. The results show that the tissue temperature-controlled system reduces the temperature of healthy tissue nearby the infected cell to 40% compare to constant and step point heating. This system keeps the healthy tissue within the threshold value (43<sup>o</sup>C) up to 1000s when it is 100s for other two techniques. After 200s, healthy tissue nearby the infected cell start to damage for constant and step point heating. But temperature-controlled system always keep the healthy tissue safe. The results of this research conclude the temperature-controlled system a better heating technique to remove the infected cell. The information published in this paper will be helpful for the physicians and bio-medical engineers to treat the infected cell or to design medical equipment.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701A New Quantum-computing-based Algorithm for Robotic Arms and Rigid Bodies’ Orientation183618461691210.22055/jacm.2021.37611.3048ENNadjetZiouiDepartment of Mechanical Engineering, Université du Québec à Trois-Rivières & Institut de Recherche sur l’Hydrogène, 3351 Boulevard des Forges, Trois-Rivières, QC G8Z 4M3, Canada0000-0003-1579-4053YousraMahmoudiDepartment of Mechanical Engineering, Université du Québec à Trois-Rivières & Institut de Recherche sur l’Hydrogène, 3351 Boulevard des Forges, Trois-Rivières, QC G8Z 4M3, CanadaLaboratoire RECITS, Université de la Science et de la Technologie Houari Boumediene, BP 32 Bab Ezzouar, 16111, Algiers, AlgeriaDépartement Hydraulique Urbaine, École Nationale Supérieure d’Hydraulique Arbaoui Abdallah, 29 Route de Soumaâ, Blida, Algeria 0000-0002-0048-4194AichaMahmoudiProjet Véo, Université de Sherbrooke, 2500 Boulevard de l'Université, Sherbrooke, QC J1K 2R1, Canada0000-0003-4595-2112MohamedTadjineLaboratoire de Commande des Processus, Ecole Nationale Polytechnique, Algiers, Algeria0000-0001-9715-2673SaidBentoubaFaculty of Engineering and Science, Western Norway University of Applied Sciences Bergen, 5063, Norway0000-0003-3859-8717Journal Article20210601Quantum computing model of robotic arm orientation is presented. Spherical and vector coordinates, a homogenous rotation matrix, Pauli gates and quantum rotation operators are used to formulate the orientation model and establish a new algorithm. The quantum algorithm uses a single qubit to compute orientation and has the advantage of operation reversibility. This was validated for a SCARA robot and a five-joints articulated robotic arm. The obtained results show the effectiveness of the proposed methodology.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Thermal Behavior of Monocrystalline Silicon Solar Cells: A Numerical and Experimental Investigation on the Module Encapsulation Materials184718551691710.22055/jacm.2021.37852.3101ENAnaPavlovicUniversity of Bologna (UNIBO), Department of Industrial Engineering, Via Fontanelle 40, Forlì, Italy0000-0003-2158-1820CristianoFragassaUniversity of Bologna (UNIBO), Department of Industrial Engineering, Via Fontanelle 40, Forlì, Italy0000-0003-0046-8810MarcoBertoldiUniversity of Bologna (UNIBO), Department of Electrical, Energy and Information Engineering "Guglielmo Marconi", Viale del Risorgimento 2, Bologna0000-0002-6188-2751VladyslavMikhnychKarlsruhe Institute of Technology (KIT), Karlsruhe School of Optics & Photonics, Schloßpl. 19, 76131 Karlsruhe, Germany0000-0003-4930-848XJournal Article20210629This research outlines the numerical predictions of the heat distribution in solar cells, accompanied by their empirical validation. Finite element thermal models of five laminated silicon solar photovoltaic cells were firstly established using a simulation software (ANSYS®). The flexible laminated solar cells under study are made of a highly transparent frontsheet, a silicon cell between two encapsulants, and a backsheet. Different combinations of layers (i.e., materials and thicknesses) were taken into account in order to analyze their effect on thermal behavior. Thermal properties of materials were derived in accordance with the literature. Similarly, boundary conditions, loads, and heat losses by reflection and convection were also specified. The solar cells were tested using solar lamps under standard conditions (irradiance: 1000W/m<sup>2</sup>; room-temperature: 25°C) with real-time temperatures measured by a thermal imager. This analysis offers an interpretation of how temperature evolves through the solar cell and, consequently, how the design choice can influence the cells’ efficiency.Shahid Chamran University of AhvazJournal of Applied and Computational Mechanics2383-45367320210701Buckling of Shell Panels Made of Fiberglass and Reinforced with an Orthogonal Grid of Stiffeners185618611691810.22055/jacm.2021.37768.3078ENAlexeySemenovDepartment of Computer Science, Saint Petersburg State University of Architecture and Civil Engineering, 4, 2nd Krasnoarmeyskaya st., Saint-Petersburg, 190005, Russia0000-0001-9490-7364Journal Article20210618The paper presents an approach to the stress-strain and buckling analysis in fiberglass cylindrical and conical 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 and the method of structural anisotropy. We use a computational algorithm based on the Ritz method and the best parameter continuation method. We also provide buckling load values and make a comparison between two types of approaches to account for stiffeners, which shows good convergence.