Shahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Numerical Simulation of free falling Cylinders using Overset Grid MethodNumerical Simulation of free falling Cylinders using Overset Grid Method11679010.22044/jsfm.2016.790FAH.Amin EsmaeiliM.SaghafianM.A.BadriJournal Article20131128An incompressible and unsteady flow problem in free falling of cylindrical particles was investigated in the intermediate Reynolds Number. Shape of the body that is used in an initial value problem in a vertical channel for the free falling was considered circular, square and rectangular cylinders. Fluid flow was computed by the Navier-Stokes equations for moderate Reynolds numbers. The particles were moved by effect of hydrodynamic, gravity and buoyancy forces. The governing equations for particle included second Newton law and torque were exerted over the particle. The particle-fluid interaction may treat by introducing a fully two dimensional overset grid Scheme. Overset grid Scheme may allow each component of a flow to be treated accurately and efficiently. The amounts of drag coefficients, free falling velocity, angular velocity and deviation for free falling cylinders were obtained for different Reynolds numbers (0.5≤Re50) and then were compared with experimental results. Free falling of two rectangular and square cylinders in different situations were also treated and the results compared with available data in fair agreement.An incompressible and unsteady flow problem in free falling of cylindrical particles was investigated in the intermediate Reynolds Number. Shape of the body that is used in an initial value problem in a vertical channel for the free falling was considered circular, square and rectangular cylinders. Fluid flow was computed by the Navier-Stokes equations for moderate Reynolds numbers. The particles were moved by effect of hydrodynamic, gravity and buoyancy forces. The governing equations for particle included second Newton law and torque were exerted over the particle. The particle-fluid interaction may treat by introducing a fully two dimensional overset grid Scheme. Overset grid Scheme may allow each component of a flow to be treated accurately and efficiently. The amounts of drag coefficients, free falling velocity, angular velocity and deviation for free falling cylinders were obtained for different Reynolds numbers (0.5≤Re50) and then were compared with experimental results. Free falling of two rectangular and square cylinders in different situations were also treated and the results compared with available data in fair agreement.https://jsfm.shahroodut.ac.ir/article_790_66996bfbaea108433807ca6d34d1386a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Simulation, analysis and optimization of airplane wing leading edge structure against bird strikeSimulation, analysis and optimization of airplane wing leading edge structure against bird strike173279510.22044/jsfm.2016.795FAH.PahangeM.H.AbolbashariJournal Article20150419A bird strike incident is one of the most dangerous threats to flight safety. In this study, Smoothed Particles Hydrodynamics (SPH) method has been used for simulating the bird strike to an airplane wing leading edge structure using Ls-Dyna software. In order to verify the model, first, experiment of bird strike to a flat Aluminum plate has been simulated and strains and deformations on target plate have been compared with experimental results. Simulation outputs are in good agreement with experimental results. Then bird impact on an airplane wing leading edge structure has been investigated. At the next stage, considering dimensions of wing internal structural components like ribs, skin and spar as design variables, it has been tried to minimize structural mass and wing skin deformation simultaneously. To do this, relations between design variables and cost functions have been predicted by Response Surface Method (RSM), then Pareto based multi objective genetic algorithm has been used to minimize structural mass and wing skin deformation due to the bird strike. Finally dimensions of wing internal structural components are determined in such a way that wing’s damage after the collision with a bird becomes minimal.A bird strike incident is one of the most dangerous threats to flight safety. In this study, Smoothed Particles Hydrodynamics (SPH) method has been used for simulating the bird strike to an airplane wing leading edge structure using Ls-Dyna software. In order to verify the model, first, experiment of bird strike to a flat Aluminum plate has been simulated and strains and deformations on target plate have been compared with experimental results. Simulation outputs are in good agreement with experimental results. Then bird impact on an airplane wing leading edge structure has been investigated. At the next stage, considering dimensions of wing internal structural components like ribs, skin and spar as design variables, it has been tried to minimize structural mass and wing skin deformation simultaneously. To do this, relations between design variables and cost functions have been predicted by Response Surface Method (RSM), then Pareto based multi objective genetic algorithm has been used to minimize structural mass and wing skin deformation due to the bird strike. Finally dimensions of wing internal structural components are determined in such a way that wing’s damage after the collision with a bird becomes minimal.https://jsfm.shahroodut.ac.ir/article_795_7a187d001cd50b940459ec2855bd830b.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-947563201609223D Reconstruction of Carbon Nanotube Composite Using Statistical Correlation Functions3D Reconstruction of Carbon Nanotube Composite Using Statistical Correlation Functions334279810.22044/jsfm.2016.798FAM.MahdaviM.BaniassadiM.BaghaniJournal Article20150716Microstructure reconstruction from statistical microstructure descriptors is a strong field of interest for researchers worldwide, due to its importance in material design. A new methodology is presented in this paper to reconstruct microstructure with a large number of representative volume elements. The proposed methodology provides a stable input for a deterministic method which can be used to simulate performance and effective properties. The Monte Carlo technique is used as the basis of the reconstruction methodology in this work. In this paper statistical correlation functions are extracted from the images of the experimental samples. The proposed algorithm reconstructs new samples which has similar statistical correlation functions in comparison to the experimental samples. The information of the geometric distribution of the nanotubes of the composites is stored in a database of the node locations of the unit cylinder segments and the corresponding waviness, Instead of using a discrete image matrix. These node locations are attractive results which can be utilized in the simulation software in order to obtain the properties of the studied composite. In this way, robust microstructures with a large number of representative volume elements were reconstructed for the future evaluation.Microstructure reconstruction from statistical microstructure descriptors is a strong field of interest for researchers worldwide, due to its importance in material design. A new methodology is presented in this paper to reconstruct microstructure with a large number of representative volume elements. The proposed methodology provides a stable input for a deterministic method which can be used to simulate performance and effective properties. The Monte Carlo technique is used as the basis of the reconstruction methodology in this work. In this paper statistical correlation functions are extracted from the images of the experimental samples. The proposed algorithm reconstructs new samples which has similar statistical correlation functions in comparison to the experimental samples. The information of the geometric distribution of the nanotubes of the composites is stored in a database of the node locations of the unit cylinder segments and the corresponding waviness, Instead of using a discrete image matrix. These node locations are attractive results which can be utilized in the simulation software in order to obtain the properties of the studied composite. In this way, robust microstructures with a large number of representative volume elements were reconstructed for the future evaluation.https://jsfm.shahroodut.ac.ir/article_798_ed9853796266c9e21bc9db0cdeddb10a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Analytical and Experimental Investigation for Nonlinear Behavior of Flange Joints Under Axial and Lateral LoadingAnalytical and Experimental Investigation for Nonlinear Behavior of Flange Joints Under Axial and Lateral Loading435480110.22044/jsfm.2016.801FAF.MeisamiFerdowsi university of MashhadM.MoavenianFerdowsi university of MashhadA.AfsharfardFerdowsi university of MashhadS.I.PishbinJournal Article20151124The flange joints are widely used in mechanical and aerospace structures. Hence, the static and dynamic behaviors of these joints are very important to be investigated. In the present study, a flange joint is described using linear and torsional springs. The equivalent stiffness of these springs are theoretically obtained. In doing so, the edge of flange is modeled as a cantilever beam and its deformation under various loading is calculated. Then the nonlinear deformed curvature due to the enforced loading is obtained. Having the curvature, the joint stiffness is calculated separately for axial and lateral loadings. It has been observed that the equivalent linear and torsional springs have bilinear stiffness due to different behaviors in positive and negative loadings. An experimental setup consisting a single bolt flange joint specimen has been designed to investigate the load-deflection behavior. The specimen puts into different loading configurations to obtain moment-rotation and force-deflection curves. A finite element model has been developed and result of the presented analytical model are compared with the experimental, FEM and previous numerical results. Finally, using the presented model, effect of changing the practical parameters of flange, like the bolt preload and flange thickness is investigated.The flange joints are widely used in mechanical and aerospace structures. Hence, the static and dynamic behaviors of these joints are very important to be investigated. In the present study, a flange joint is described using linear and torsional springs. The equivalent stiffness of these springs are theoretically obtained. In doing so, the edge of flange is modeled as a cantilever beam and its deformation under various loading is calculated. Then the nonlinear deformed curvature due to the enforced loading is obtained. Having the curvature, the joint stiffness is calculated separately for axial and lateral loadings. It has been observed that the equivalent linear and torsional springs have bilinear stiffness due to different behaviors in positive and negative loadings. An experimental setup consisting a single bolt flange joint specimen has been designed to investigate the load-deflection behavior. The specimen puts into different loading configurations to obtain moment-rotation and force-deflection curves. A finite element model has been developed and result of the presented analytical model are compared with the experimental, FEM and previous numerical results. Finally, using the presented model, effect of changing the practical parameters of flange, like the bolt preload and flange thickness is investigated.https://jsfm.shahroodut.ac.ir/article_801_789017a7bef18c3780793ec7c64fe062.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Parametric Investigation of Buckling Capacity Improvement of Cracked Columns with Piezoelectric PatchesParametric Investigation of Buckling Capacity Improvement of Cracked Columns with Piezoelectric Patches556480210.22044/jsfm.2016.802FAM.JavidH.BiglariJournal Article20151230In this paper, instability behavior of cracked column equipped with piezoelectric patch is investigated, analytically. Polarization orientation is in such a way that axial force is induced in column. Crack is modeled as a torsional spring that connects two intact portion of column. After considering boundary conditions and continuity at the crack position and edges of piezoelectric patch positions, governing equations of buckling behavior of cracked column are obtained. Effects of effective parameters on the first and second buckling capacity of cracked column such as cracked parameters (position and depth), position, length and voltage of piezoelectric patch is investigated using suitable sketches. The results show that crack decreases buckling capacity, prominently. This decrease depends on position and depth of crack. Using appropriate position and length of piezoelectric patch, buckling capacity of cracked column is improved. The appropriate position of piezoelectric patch is not analogues to the crack position, necessarily and may be determined according to the column boundary conditions and crack position. The results show that positive voltage has opposite effects on the first and second buckling capacity of column.In this paper, instability behavior of cracked column equipped with piezoelectric patch is investigated, analytically. Polarization orientation is in such a way that axial force is induced in column. Crack is modeled as a torsional spring that connects two intact portion of column. After considering boundary conditions and continuity at the crack position and edges of piezoelectric patch positions, governing equations of buckling behavior of cracked column are obtained. Effects of effective parameters on the first and second buckling capacity of cracked column such as cracked parameters (position and depth), position, length and voltage of piezoelectric patch is investigated using suitable sketches. The results show that crack decreases buckling capacity, prominently. This decrease depends on position and depth of crack. Using appropriate position and length of piezoelectric patch, buckling capacity of cracked column is improved. The appropriate position of piezoelectric patch is not analogues to the crack position, necessarily and may be determined according to the column boundary conditions and crack position. The results show that positive voltage has opposite effects on the first and second buckling capacity of column.https://jsfm.shahroodut.ac.ir/article_802_0be80dd1ee64cef01d2f8007f272c352.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Analysis and modeling of aerodynamic pressure wave due to high-speed trains passageAnalysis and modeling of aerodynamic pressure wave due to high-speed trains passage657880410.22044/jsfm.2016.804FAV.SarafraziM.TalaeeIran University of Science and TechnologyJournal Article20160224One of the important subjects in aerodynamic design of high speed train is the slipstream or air flow around the train. This issue not only waste a significant amount of energy and fuel but also reduce the overall efficiency of the train. Due to irregularity the turbulence air flow around the high speed trains aerodynamic drag, noise and vibration are the most phenomenon that concern in this area. The slipstream is one of these aerodynamics features that directily related to the safety of passengers and workers near the line. When high speed train passes near to a standing human body, the produced presuer wave forces on the body surface and may lead to fall.In this numerical simulation we used dynamic mesh method and turbulence equation k-e in computational fluid dynamic software (Fluent). The aim of this study is to investigate the effects of air pressure around the train of moving with 300 km/hr velocity. The pressure waves are calculated in both modes of regular passage and passing with operation of aerodynamic brake. The resules are determined as the safty distance of the passenger from the line due to the different speedOne of the important subjects in aerodynamic design of high speed train is the slipstream or air flow around the train. This issue not only waste a significant amount of energy and fuel but also reduce the overall efficiency of the train. Due to irregularity the turbulence air flow around the high speed trains aerodynamic drag, noise and vibration are the most phenomenon that concern in this area. The slipstream is one of these aerodynamics features that directily related to the safety of passengers and workers near the line. When high speed train passes near to a standing human body, the produced presuer wave forces on the body surface and may lead to fall.In this numerical simulation we used dynamic mesh method and turbulence equation k-e in computational fluid dynamic software (Fluent). The aim of this study is to investigate the effects of air pressure around the train of moving with 300 km/hr velocity. The pressure waves are calculated in both modes of regular passage and passing with operation of aerodynamic brake. The resules are determined as the safty distance of the passenger from the line due to the different speedhttps://jsfm.shahroodut.ac.ir/article_804_850af7bef9871770256a3c667530507a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Effects of design and lubrication parameters on thermo-hydrodynamic performance of noncircular three lobe journal bearings with micropolar lubricantEffects of design and lubrication parameters on thermo-hydrodynamic performance of noncircular three lobe journal bearings with micropolar lubricant799680310.22044/jsfm.2016.803FAA.Rasoolizadeh Shoorokiyazd universityA.Dashti RahmatabadiMahdiZare Mehrjardiardakan universityRezaRashidi Meybodipayam noor universityJournal Article20160213Due to wide range of engineering application of noncircular journal bearings in industry, their performance analysis is important. Shear stresses and friction forces in lubricant film cause variation of temperature condition in the bearing. Parameters such as increasing in rotational speed of shaft, decreasing in clearance space of bearing and low coefficient of heat transfer of lubricant, shaft and bearing cause the intensity changes in lubricant film temperature conditions. In the present work, the effects of design parameters such as speed of shaft, thickness of bearing shell, bearing length and micropolar lubricant characteristics on changes of temperature, lubricant film pressure distribution, load carrying capacity and attitude angle of bearing are investigated. By comparing the performance of noncircular three lobe bearings with and without considering temperature effects, it is found that pressure of lubricant film as well as load carrying capacity of bearing system decrease with increasing temperature, especially at high rotational speed of rotor. Also, the results show that with increasing micropolarity properties of lubricant, the effect of temperature on performance of mentioned bearings can be increased.Due to wide range of engineering application of noncircular journal bearings in industry, their performance analysis is important. Shear stresses and friction forces in lubricant film cause variation of temperature condition in the bearing. Parameters such as increasing in rotational speed of shaft, decreasing in clearance space of bearing and low coefficient of heat transfer of lubricant, shaft and bearing cause the intensity changes in lubricant film temperature conditions. In the present work, the effects of design parameters such as speed of shaft, thickness of bearing shell, bearing length and micropolar lubricant characteristics on changes of temperature, lubricant film pressure distribution, load carrying capacity and attitude angle of bearing are investigated. By comparing the performance of noncircular three lobe bearings with and without considering temperature effects, it is found that pressure of lubricant film as well as load carrying capacity of bearing system decrease with increasing temperature, especially at high rotational speed of rotor. Also, the results show that with increasing micropolarity properties of lubricant, the effect of temperature on performance of mentioned bearings can be increased.https://jsfm.shahroodut.ac.ir/article_803_d6f5ebd4dad3518303d8c7732ab4ae8a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Non-linear finite element low-velocity impact response analysis of a viscoelastic composite plate, employing a layerwise theoryNon-linear finite element low-velocity impact response analysis of a viscoelastic composite plate, employing a layerwise theory9710882210.22044/jsfm.2016.822FAM.H.PachenariA.MozaffariM.ShariyatJournal Article20160410In the present paper, dynamic behavior of a multilayer composite plate with viscoelastic structural damping is investigated against a low-velocity impact by a spherical indenter. Hertz contact law is refined to include effect of the lower layers on the stiffness of the contact region and used in a non-linear form. Voltra hierarchical integral is employed for modeling the viscoelastic material and the layerwise theory and non-linear strain-displacement relations are used to model the plate more accurately. To solve the governing integro-differential equations, a combination of the finite element method, trapezoidal integration method for the Volterra integrals, and the Newmark numerical time integration method is used. In the results section, effects of the various viscoelasticity parameters and the indenter velocity on the time histories of the contact force, indentation, and lateral deflection of the plate are investigated. Results show that due to the damping nature of the viscoelastic materials, the plate rigidity and contact force increase whereas the maximum lateral deflection and the indentation decrease. Furthermore, higher contact forces do not necessarily indicates higher indentations.In the present paper, dynamic behavior of a multilayer composite plate with viscoelastic structural damping is investigated against a low-velocity impact by a spherical indenter. Hertz contact law is refined to include effect of the lower layers on the stiffness of the contact region and used in a non-linear form. Voltra hierarchical integral is employed for modeling the viscoelastic material and the layerwise theory and non-linear strain-displacement relations are used to model the plate more accurately. To solve the governing integro-differential equations, a combination of the finite element method, trapezoidal integration method for the Volterra integrals, and the Newmark numerical time integration method is used. In the results section, effects of the various viscoelasticity parameters and the indenter velocity on the time histories of the contact force, indentation, and lateral deflection of the plate are investigated. Results show that due to the damping nature of the viscoelastic materials, the plate rigidity and contact force increase whereas the maximum lateral deflection and the indentation decrease. Furthermore, higher contact forces do not necessarily indicates higher indentations.https://jsfm.shahroodut.ac.ir/article_822_1c1e8dced775176d039e6e401ef36df2.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Exploring the Effect of Foot Position on Performance of a Modular and Hierarchical Movement Planner in Planning the Sit-to-Stand TransferExploring the Effect of Foot Position on Performance of a Modular and Hierarchical Movement Planner in Planning the Sit-to-Stand Transfer10911879610.22044/jsfm.2016.796FAM.AskarishahiSh.Hadian JaziJournal Article20150531Many researchers in the science of human behavior explore the role and function of the central nervous system (CNS) in planning and controlling of human movements. To this end, researchers have presented several different models. In present research a computer simulation of CNS's performance in designing the Sit-to-Stand transfer, which has been recently presented, is developed. This motion simulator is a modular and hierarchical movement planner (MHMP), based on decomposition hypothesis. In this paper the effect of foot position, as an environmental condition, is explored. The performance of the MHMP is evaluated with experimental captured motion. To this end, sit to stand motion is captured for five different foot position. Among these five motions, three motions are used to train the MHMP and the remaining motions are used for evaluating its performance. Results show that the maximum error of ankle, knee and hip angles are 0.122, 0.069 and 0.092 (radian) respectively. The maximum allowed error, suggested in researches, is 0.17(radian). The results show that the MHMP has a good performance in planning the motion phases under this condition.Many researchers in the science of human behavior explore the role and function of the central nervous system (CNS) in planning and controlling of human movements. To this end, researchers have presented several different models. In present research a computer simulation of CNS's performance in designing the Sit-to-Stand transfer, which has been recently presented, is developed. This motion simulator is a modular and hierarchical movement planner (MHMP), based on decomposition hypothesis. In this paper the effect of foot position, as an environmental condition, is explored. The performance of the MHMP is evaluated with experimental captured motion. To this end, sit to stand motion is captured for five different foot position. Among these five motions, three motions are used to train the MHMP and the remaining motions are used for evaluating its performance. Results show that the maximum error of ankle, knee and hip angles are 0.122, 0.069 and 0.092 (radian) respectively. The maximum allowed error, suggested in researches, is 0.17(radian). The results show that the MHMP has a good performance in planning the motion phases under this condition.https://jsfm.shahroodut.ac.ir/article_796_f1eea3ddc4ee9286507ad7dca9a2156e.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Thermal Response Analysis of Lock-in Thermography Test to Identify the Delamination in Composite MaterialsThermal Response Analysis of Lock-in Thermography Test to Identify the Delamination in Composite Materials11913680810.22044/jsfm.2016.808FAJ.Safari ChaykandiM.AbadiJournal Article20160411The delamination defects in the laminated composite structures are detected using the transient response analysis in the lock-in thermography test. The analysis procedure is implemented for the laminated composite plate subjected to harmonic heat flux in the top surface of the plate, in which the delamination defects with different sizes are located in different depths. The thermal wave propagation due to heat loading is simulated in the plate thickness along the healthy and defective regions and the amplitude and phase differences are determine in such regions. The phase difference of reflected waves is more sensitive to the defects than the amplitude of the waves. The results show that the phase difference between reflected waves from healthy and defective regions increases up to specific frequency value and the optimum frequency is determined using the present simulation. The results also show that the logarithmic phase difference depends linearly on the depth of delamination and such dependency can be used to determine the defect depth. The simulation results of lock-in thermography test are used to determine the optimum parameters required to detect the delamination in the laminated composite materials.The delamination defects in the laminated composite structures are detected using the transient response analysis in the lock-in thermography test. The analysis procedure is implemented for the laminated composite plate subjected to harmonic heat flux in the top surface of the plate, in which the delamination defects with different sizes are located in different depths. The thermal wave propagation due to heat loading is simulated in the plate thickness along the healthy and defective regions and the amplitude and phase differences are determine in such regions. The phase difference of reflected waves is more sensitive to the defects than the amplitude of the waves. The results show that the phase difference between reflected waves from healthy and defective regions increases up to specific frequency value and the optimum frequency is determined using the present simulation. The results also show that the logarithmic phase difference depends linearly on the depth of delamination and such dependency can be used to determine the defect depth. The simulation results of lock-in thermography test are used to determine the optimum parameters required to detect the delamination in the laminated composite materials.https://jsfm.shahroodut.ac.ir/article_808_11426b238c43aa7bb60e1617c9659547.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922The study of the effects of thermal environment on Free Vibration analysis of two-dimensional Functionally Graded Rectangular plates on Pasternak elastic foundationThe study of the effects of thermal environment on Free Vibration analysis of two-dimensional Functionally Graded Rectangular plates on Pasternak elastic foundation13714779210.22044/jsfm.2016.792FAK.KhorshidiA.BakhsheshyH.GhadirianJournal Article20140309This paper is the study of the effects of thermal environment on vibration analysis of a two-dimensional functionally graded rectangular plate on Pasternak elastic foundation. The natural frequencies of the plate are calculated by using the Rayleigh–Ritz method based on minimizing the total energy of the plate.analysis procedure is. The transverse displacement of the plate based on the third-order shear deformation plate theory(TSDT) is approximated by a set of admissible trial functions which is required to satisfy the clamped (CL) and simply supported geometric boundary conditions. For verifying the accuracy of this method, results are compared with those reported in the literature. As we seeit is shown a good conformance is derived from the obtained results and the exact solution. In the numerical results, the effects of volume fraction coefficients, thickness ratios ,aspect ratios of the FG plates, foundation stiffness parameters, temperature and boundary conditions on the natural frequencies are examined and discussed in detail.This paper is the study of the effects of thermal environment on vibration analysis of a two-dimensional functionally graded rectangular plate on Pasternak elastic foundation. The natural frequencies of the plate are calculated by using the Rayleigh–Ritz method based on minimizing the total energy of the plate.analysis procedure is. The transverse displacement of the plate based on the third-order shear deformation plate theory(TSDT) is approximated by a set of admissible trial functions which is required to satisfy the clamped (CL) and simply supported geometric boundary conditions. For verifying the accuracy of this method, results are compared with those reported in the literature. As we seeit is shown a good conformance is derived from the obtained results and the exact solution. In the numerical results, the effects of volume fraction coefficients, thickness ratios ,aspect ratios of the FG plates, foundation stiffness parameters, temperature and boundary conditions on the natural frequencies are examined and discussed in detail.https://jsfm.shahroodut.ac.ir/article_792_fff338110143cbbefd762ad5907d13bb.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Predication of roll force and roll torque in hot strip rolling process using artificial neural networks and finite element methodPredication of roll force and roll torque in hot strip rolling process using artificial neural networks and finite element method14916280010.22044/jsfm.2016.800FAH.Rezaei AshtianiJournal Article20150930This paper introduces an artificial neural network (ANN) application to a hot strip mill to improve the model’s prediction ability for rolling force and rolling torque, as a function of various process parameters. To obtain a data basis for training and validation of the neural network, numerous three dimensional finite element simulations were carried out for different sets of process variables. Experimental data were compared with the finite element predictions to verify the model accuracy. Thus the ABAQUS and MATLAB soft wares are used to simulate the finite element method and an artificial neural network, respectively. The back-propagation algorithm and Levenberg–Marquardt Training function were used in the artificial neural network. The input variables are selected to be, initial temperature of the strip, interface heat-transfer coefficient between strip and work roll, percentage of thickness reduction, initial thickness and rolling speed. The resulted ANN model is feasible for on-line control and rolling schedule optimization and can be easily and rapidly predicted the roll force and roll torque.This paper introduces an artificial neural network (ANN) application to a hot strip mill to improve the model’s prediction ability for rolling force and rolling torque, as a function of various process parameters. To obtain a data basis for training and validation of the neural network, numerous three dimensional finite element simulations were carried out for different sets of process variables. Experimental data were compared with the finite element predictions to verify the model accuracy. Thus the ABAQUS and MATLAB soft wares are used to simulate the finite element method and an artificial neural network, respectively. The back-propagation algorithm and Levenberg–Marquardt Training function were used in the artificial neural network. The input variables are selected to be, initial temperature of the strip, interface heat-transfer coefficient between strip and work roll, percentage of thickness reduction, initial thickness and rolling speed. The resulted ANN model is feasible for on-line control and rolling schedule optimization and can be easily and rapidly predicted the roll force and roll torque.https://jsfm.shahroodut.ac.ir/article_800_33961c58ba1c4978b747aaef6f39009b.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Primary creep analysis of nanocomposite thick-walled sphere under thermal, magnetic and mechanical loading using viscoelastic Bergurs modelPrimary creep analysis of nanocomposite thick-walled sphere under thermal, magnetic and mechanical loading using viscoelastic Bergurs model16317480610.22044/jsfm.2016.806FAH.Mohammadi HooyehA.Naddaf OskoueiM.Sayah BadkhorJournal Article20160303In this study, stresses and strains history of nanocomposite thick walled sphere wich is made of Polyimide reinforced by SiO2 nanoparticles are investigated using the viscoelastic Bergurs model. Loading of sphere includes the thermal and magnetic uniform field under hydrostatic inner pressure. Constitutive model of the problem based on elasticity relationship is obtained. Solving this equation with respect to mechanical boundary conditions leads to stresses and strains at zero time so called thermo-elastic solution in which these stresses are used at beginning creep problem procedure. With differentiating of constitutive differential equation with respect of time and using relationship between Bergurs model and Prandtl–Reuss relations drives a new constitutive equation for creep. History of stresses and strains are provided by solving this equation using numerical strain rate method. Results show that maximum effective stress and creep strains occurs at inner surface. Also the primary creep stage occurs up to 1000 second with quick change of stresses and strains. After this stage the value of stresses and strains change with uniform speed.In this study, stresses and strains history of nanocomposite thick walled sphere wich is made of Polyimide reinforced by SiO2 nanoparticles are investigated using the viscoelastic Bergurs model. Loading of sphere includes the thermal and magnetic uniform field under hydrostatic inner pressure. Constitutive model of the problem based on elasticity relationship is obtained. Solving this equation with respect to mechanical boundary conditions leads to stresses and strains at zero time so called thermo-elastic solution in which these stresses are used at beginning creep problem procedure. With differentiating of constitutive differential equation with respect of time and using relationship between Bergurs model and Prandtl–Reuss relations drives a new constitutive equation for creep. History of stresses and strains are provided by solving this equation using numerical strain rate method. Results show that maximum effective stress and creep strains occurs at inner surface. Also the primary creep stage occurs up to 1000 second with quick change of stresses and strains. After this stage the value of stresses and strains change with uniform speed.https://jsfm.shahroodut.ac.ir/article_806_9628ba3821be80baa4216586c0c8b34f.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Fatigue Life Prediction of Adhesive Joints Based on Initial Stiffness and Stiffness DegradationFatigue Life Prediction of Adhesive Joints Based on Initial Stiffness and Stiffness Degradation17518382110.22044/jsfm.2016.821FAP.AkbarzadehDepartment of Mechanical Engineering, Ferdowsi University of MashhadKh.FarhangdoostDepartment of Mechanical engineering, Ferdwsi University of Mashhad, IranJournal Article20150111Use of adhesively bonded joints in structures has been increased in recent years because of their advantages. In structural applications, fatigue is generally considered to be the most important form of loading in respect to long-term service life. In this paper, experimental tests have been performed in order to examine the fatigue damage process of double lap adhesive joints, composed of E-glass laminates and epoxy adhesive with TiO2 particles, under fatigue loading and then a method was proposed for fatigue life prediction based on initial stiffness. The scatter in results shows that the fatigue life is mostly dependent on initial stiffness. According to this, an exponential equation based on the initial stiffness has been proposed to have an initial estimation of fatigue life of the adhesive joints. In addition to the initial stiffness, stiffness degradation of the joints under fatigue loading influences on their final performance and therefore the initial estimation of fatigue life based on the initial stiffness could be modified. Also to have an online control of the fatigue damage process during the fatigue loading, a damage index using stiffness degradation was introduced.Use of adhesively bonded joints in structures has been increased in recent years because of their advantages. In structural applications, fatigue is generally considered to be the most important form of loading in respect to long-term service life. In this paper, experimental tests have been performed in order to examine the fatigue damage process of double lap adhesive joints, composed of E-glass laminates and epoxy adhesive with TiO2 particles, under fatigue loading and then a method was proposed for fatigue life prediction based on initial stiffness. The scatter in results shows that the fatigue life is mostly dependent on initial stiffness. According to this, an exponential equation based on the initial stiffness has been proposed to have an initial estimation of fatigue life of the adhesive joints. In addition to the initial stiffness, stiffness degradation of the joints under fatigue loading influences on their final performance and therefore the initial estimation of fatigue life based on the initial stiffness could be modified. Also to have an online control of the fatigue damage process during the fatigue loading, a damage index using stiffness degradation was introduced.https://jsfm.shahroodut.ac.ir/article_821_92b1150311751f0f1864df6c667501de.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Reliability Analysis and Evaluating of Rotary Drilling Machine –
Case Study: Mine No.1 of Golgohar Iron Ore MineReliability Analysis and Evaluating of Rotary Drilling Machine –
Case Study: Mine No.1 of Golgohar Iron Ore Mine18519479310.22044/jsfm.2016.793FAS.MohammadiM.GhiasiM.AtaeiE.PourzamaniJournal Article20150318Open pit mines have very high economic and operational risk because of high rate of production and using large scale instruments. So, it is necessary to guarantee the availability of the operational equipment. First stage of mining in these mines is drilling for blasting operation and the availability of related machines has direct effect on all subsequent operations and consequently on overall economy of the operation. Today, operation conditions and machines availability can be assessed and analyzed due to development of reliability methods in different engineering fields. Markov method is one of the potent methods in this context which is based on stochastic process theory. In this study, reliability of rotary drilling machines in Golgohar Mine No. 1, Sirjan, is analyzed using this method. In this way, reliability is calculated by investigating failure dad of these machines during 5 years and modeling different conditions of their availability. So, at least 2 machines will be ready and available for drilling operation with the probability of 95.42% corresponding to 347 days out of total 363 working days.Open pit mines have very high economic and operational risk because of high rate of production and using large scale instruments. So, it is necessary to guarantee the availability of the operational equipment. First stage of mining in these mines is drilling for blasting operation and the availability of related machines has direct effect on all subsequent operations and consequently on overall economy of the operation. Today, operation conditions and machines availability can be assessed and analyzed due to development of reliability methods in different engineering fields. Markov method is one of the potent methods in this context which is based on stochastic process theory. In this study, reliability of rotary drilling machines in Golgohar Mine No. 1, Sirjan, is analyzed using this method. In this way, reliability is calculated by investigating failure dad of these machines during 5 years and modeling different conditions of their availability. So, at least 2 machines will be ready and available for drilling operation with the probability of 95.42% corresponding to 347 days out of total 363 working days.https://jsfm.shahroodut.ac.ir/article_793_a69cc03cfb60dc5819b04f2e9a2b4adc.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Performance improvement of hybrid Darrieus-Savonius wind turbinePerformance improvement of hybrid Darrieus-Savonius wind turbine19521280710.22044/jsfm.2016.807FAA.RoshanM. J.Maghrebifaculty of engineering, ferdowsi university on mashhadJournal Article20160316Vertical axis wind turbine, divided into Darrieus and Savonius, have been considered by many scientists owing to simplicity of design and operation independent of wind direction. Although Darrieus turbine has higher efficiency than Savonius, it faces to essential problem of self-starting. The aim of the present study is the combination of Darrieus and Savonius turbine in order to reach to a model with high starting torque and proper performance range. First, a model involving 3 straight-bladed Darrieus with 2 Savonius blades is simulated. For improving the performance, Savonius blade arc angle and airfoil solidity are changed. By optimizing the arc angle and solidity, the initial hybrid model performance is improved. Two-dimensional simulation using computational fluid dynamics is conducted using moving mesh for rotating part. According to results, the combination of Darrieus and Savonius turbine has favorable impact on self starting and makes performance curve smooth. It was also observed that solidity of 0.75 for Darrieus and angle of 150° for savonius blade arc are the optimized mode for maximum power coefficient, starting torque and performance range. In other words, The hybrid model with maximum performance is achieved by using these values at the same time.Vertical axis wind turbine, divided into Darrieus and Savonius, have been considered by many scientists owing to simplicity of design and operation independent of wind direction. Although Darrieus turbine has higher efficiency than Savonius, it faces to essential problem of self-starting. The aim of the present study is the combination of Darrieus and Savonius turbine in order to reach to a model with high starting torque and proper performance range. First, a model involving 3 straight-bladed Darrieus with 2 Savonius blades is simulated. For improving the performance, Savonius blade arc angle and airfoil solidity are changed. By optimizing the arc angle and solidity, the initial hybrid model performance is improved. Two-dimensional simulation using computational fluid dynamics is conducted using moving mesh for rotating part. According to results, the combination of Darrieus and Savonius turbine has favorable impact on self starting and makes performance curve smooth. It was also observed that solidity of 0.75 for Darrieus and angle of 150° for savonius blade arc are the optimized mode for maximum power coefficient, starting torque and performance range. In other words, The hybrid model with maximum performance is achieved by using these values at the same time.https://jsfm.shahroodut.ac.ir/article_807_bbc23a9e237f66dfa06d3251471b13cf.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Using net present value as a single objective function for thermo-economic optimization of a plate-fin recuperator in microturbineUsing net present value as a single objective function for thermo-economic optimization of a plate-fin recuperator in microturbine21322681010.22044/jsfm.2016.810FAP.MaghsoudiP.HanafizadehJournal Article20160424In this paper, genetic and particle swarm algorithms were applied to optimize a recuperator in a 200 kW microturbine, with the net present value (NPV) considered as the objective function. Examined for the optimization purposes is a plate-fin heat exchanger of offset strip fin configuration as well as counter- and cross-flow arrangements. Fin pitch, fin height, fin offset length, recuperator flow length, recuperator width, and recuperator height were considered as the six design parameters investigated. A sensitivity analysis is performed in order to investigate the effect of design parameters on the objective function. Comparing the results of genetic algorithm(GA) to those of particle swarm algorithm indicated that, the particle swarm optimization(PSO) approach had succeeded to achieve better results from computations’ volume and optimum value perspectives. As the used objective function in this research was a function of total cost and efficiency of the recuperator, the results of the present research were compared to thermo-economic optimization results. Being based on NPV and cycle efficiency criteria, the comparisons indicated superiority of the results of the present study, when viewed from technical and economic perspectives. Results show that in the present study NPV and cycle efficiency have increased up to 11.54% and 11.69%.In this paper, genetic and particle swarm algorithms were applied to optimize a recuperator in a 200 kW microturbine, with the net present value (NPV) considered as the objective function. Examined for the optimization purposes is a plate-fin heat exchanger of offset strip fin configuration as well as counter- and cross-flow arrangements. Fin pitch, fin height, fin offset length, recuperator flow length, recuperator width, and recuperator height were considered as the six design parameters investigated. A sensitivity analysis is performed in order to investigate the effect of design parameters on the objective function. Comparing the results of genetic algorithm(GA) to those of particle swarm algorithm indicated that, the particle swarm optimization(PSO) approach had succeeded to achieve better results from computations’ volume and optimum value perspectives. As the used objective function in this research was a function of total cost and efficiency of the recuperator, the results of the present research were compared to thermo-economic optimization results. Being based on NPV and cycle efficiency criteria, the comparisons indicated superiority of the results of the present study, when viewed from technical and economic perspectives. Results show that in the present study NPV and cycle efficiency have increased up to 11.54% and 11.69%.https://jsfm.shahroodut.ac.ir/article_810_779fe49ed028a4d711e415531462c78a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Optimization of Shazand Arak power plant boiler rifled tubes using genetic algorithmOptimization of Shazand Arak power plant boiler rifled tubes using genetic algorithm22723779110.22044/jsfm.2016.791FAH.MokhtariA.MahmoodiR.Haghighi KhoshkhooJournal Article20140126Increase in pressure drop due to the use of rifled tubes and also increase the heat transfer coefficient in two-phase flow and during boiling length (single-phase region) caused that using evolutionary genetic algorithm to optimize the parameters of rifled tubes, with the goal of increasing the heat transfer coefficient over the single-phase region and reduce the pressure drop in natural circulation system of drum should be dicussed in this paper. Choice of objective function was done with respect to the analysis of the pressure drop in each part of the drum circle and its large changes relative to the change of rifled tubes parameters. The optimization was done using two scenarios based on constant number of tubes (scenario I) and constant mass flux (scenario II). The results showed that the design in case of constant number of tubes which is defined base on no changes in length and width of the Shazand Arak power plant boiler, led to reduce the height of boiler from 60m to 27.23m. Finally, the role of rifled tubes parameters and mass flux in this scenario relative to first scenario in reduce of pressure drop was led to set values of the rifled tubes parameters.Increase in pressure drop due to the use of rifled tubes and also increase the heat transfer coefficient in two-phase flow and during boiling length (single-phase region) caused that using evolutionary genetic algorithm to optimize the parameters of rifled tubes, with the goal of increasing the heat transfer coefficient over the single-phase region and reduce the pressure drop in natural circulation system of drum should be dicussed in this paper. Choice of objective function was done with respect to the analysis of the pressure drop in each part of the drum circle and its large changes relative to the change of rifled tubes parameters. The optimization was done using two scenarios based on constant number of tubes (scenario I) and constant mass flux (scenario II). The results showed that the design in case of constant number of tubes which is defined base on no changes in length and width of the Shazand Arak power plant boiler, led to reduce the height of boiler from 60m to 27.23m. Finally, the role of rifled tubes parameters and mass flux in this scenario relative to first scenario in reduce of pressure drop was led to set values of the rifled tubes parameters.https://jsfm.shahroodut.ac.ir/article_791_7472c477939d04275645f239b7e8905c.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Numerical Simulation of the Controlled Fluid Flow Around the Cylinder Using Compound ActuatorsNumerical Simulation of the Controlled Fluid Flow Around the Cylinder Using Compound Actuators23924781510.22044/jsfm.2016.815FANikiRezazadehJournal Article20130729Abstract <br /> In present paper, the fluid flow around the cylinder has been controlled to change the pattern flow. Flow control methods are classified into two types,that were called active and passive methods.In this work, the compound method is applied, the splitter plate in the wake zone as passive method and EHD actuators as active method have been applied simultaneously. For simulation of fluid flow, Navier Stoks equations for fluid flow and conversation of electric charge and Poisson equations for electric field are calculated numerically by finite volume approach. The length of the splitter plate is equal the diameter of the cylinder and it is connected to the ground as well as the cylinder surface. ،The wire electrodes and the cylinder are used as EHD actuators. Calculations are done for the low Re number, Re=40 and the applied voltage V=10,20 kV. The results show that the presence of the splitter plate at low Re number does not stabilize the wake zone and the vortexes behind the cylinder start to move.<br /> Keywords: Ion Wind; Drag Force; ElectroHydroDynamic Actuator; Splitter Plate.Abstract <br /> In present paper, the fluid flow around the cylinder has been controlled to change the pattern flow. Flow control methods are classified into two types,that were called active and passive methods.In this work, the compound method is applied, the splitter plate in the wake zone as passive method and EHD actuators as active method have been applied simultaneously. For simulation of fluid flow, Navier Stoks equations for fluid flow and conversation of electric charge and Poisson equations for electric field are calculated numerically by finite volume approach. The length of the splitter plate is equal the diameter of the cylinder and it is connected to the ground as well as the cylinder surface. ،The wire electrodes and the cylinder are used as EHD actuators. Calculations are done for the low Re number, Re=40 and the applied voltage V=10,20 kV. The results show that the presence of the splitter plate at low Re number does not stabilize the wake zone and the vortexes behind the cylinder start to move.<br /> Keywords: Ion Wind; Drag Force; ElectroHydroDynamic Actuator; Splitter Plate.https://jsfm.shahroodut.ac.ir/article_815_ad8e4cc644f549fac99193f26bd6fa22.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Analysis the influence of air outlet opening’s location in an industrial environment with high temperature radiant heaters on uniformity of thermal conditions and distribution of the contaminant in presence of asymmetric flow fieldAnalysis the influence of air outlet opening’s location in an industrial environment with high temperature radiant heaters on uniformity of thermal conditions and distribution of the contaminant in presence of asymmetric flow field24926080910.22044/jsfm.2016.809FAM.AminzadehA.SafavinejadA.RZolfaghariJournal Article20160414In the present study, the effect of outlet opening location on providing appropriate and uniform thermal conditions and distribution of carbon dioxide species which is produced by high temperature radiant heaters under nonuniform flow field in an industrial environment have been investigated. For this reason, two dimentional environment with one inlet and one outlet opening have been considered. Two locations of air outlet opening near the floor and near the ceiling and also, two different types of radiant heaters’ arrangement have been considered as: single radiant heater, couple radiant heaters. For the mentioned conditions, continuity equation, momentum equations, energy equation, radiative transfer equation and species concentration equation have been solved by OpenFoam numerical solver. Also energy consumption have been evaluated in the present study. The results show that location of air outlet opening near the ceiling effects significantly on decreasing carbon dioxide concentration. As well as, usage of couple radiant heaters makes temperature distribution more uniform. The results indicate that energy consumption increases in presence of couple radiant heaters for outlet opening which is near the floor 3 percent and near the ceiling 15 percent in comparison with single radiant heater.In the present study, the effect of outlet opening location on providing appropriate and uniform thermal conditions and distribution of carbon dioxide species which is produced by high temperature radiant heaters under nonuniform flow field in an industrial environment have been investigated. For this reason, two dimentional environment with one inlet and one outlet opening have been considered. Two locations of air outlet opening near the floor and near the ceiling and also, two different types of radiant heaters’ arrangement have been considered as: single radiant heater, couple radiant heaters. For the mentioned conditions, continuity equation, momentum equations, energy equation, radiative transfer equation and species concentration equation have been solved by OpenFoam numerical solver. Also energy consumption have been evaluated in the present study. The results show that location of air outlet opening near the ceiling effects significantly on decreasing carbon dioxide concentration. As well as, usage of couple radiant heaters makes temperature distribution more uniform. The results indicate that energy consumption increases in presence of couple radiant heaters for outlet opening which is near the floor 3 percent and near the ceiling 15 percent in comparison with single radiant heater.https://jsfm.shahroodut.ac.ir/article_809_2702e5bfb1f4a48fef800735e18b420a.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Influence of air inlet diffuser layout on indoor air quality and occupants’ thermal comfort in a room with baseboard heating systemInfluence of air inlet diffuser layout on indoor air quality and occupants’ thermal comfort in a room with baseboard heating system26127081210.22044/jsfm.2016.812FAZ.Imani NejadA. R.ZolfaghariM.MaerefatH.Pasdar ShahriJournal Article20160610In the recent years, the baseboard heating system has attracted the attention of many HVAC engineers due to its uniform temperature distribution and low feed water temperature. Despite this, lack of ventilation equipments in the buildings with the mentioned heating system, can cause to decrease the indoor air quality. Therefore, the main goal of this study is to evaluate the effects of fresh air supply diffuser layout on the indoor air quality and occupants’ thermal comfort conditions in a baseboard heated building in winter. For this reason, the thermal sensation and indoor air quality indices have been numerically investigated in a standard case room and under the climatic conditions of 0C as winter outdoor design temperature. The results show that the location of air supply diffuser has a significant influence on air temperature, thermal comfort conditions and indoor air quality. So that, changing the location of ceiling air diffuser from near the exterior wall to the center, can cause the average of air temperature to increase about 1.1C (from 21.4 to 22.5C). Moreover, the results indicate that selecting the proper location for fresh air supply diffuser can lead to increase in the indoor air quality improvement index from 0.44 to 0.67.In the recent years, the baseboard heating system has attracted the attention of many HVAC engineers due to its uniform temperature distribution and low feed water temperature. Despite this, lack of ventilation equipments in the buildings with the mentioned heating system, can cause to decrease the indoor air quality. Therefore, the main goal of this study is to evaluate the effects of fresh air supply diffuser layout on the indoor air quality and occupants’ thermal comfort conditions in a baseboard heated building in winter. For this reason, the thermal sensation and indoor air quality indices have been numerically investigated in a standard case room and under the climatic conditions of 0C as winter outdoor design temperature. The results show that the location of air supply diffuser has a significant influence on air temperature, thermal comfort conditions and indoor air quality. So that, changing the location of ceiling air diffuser from near the exterior wall to the center, can cause the average of air temperature to increase about 1.1C (from 21.4 to 22.5C). Moreover, the results indicate that selecting the proper location for fresh air supply diffuser can lead to increase in the indoor air quality improvement index from 0.44 to 0.67.https://jsfm.shahroodut.ac.ir/article_812_270cfd7f4a62e42d23ca734f14c2408d.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Design and Numerical Investigation of a Small Scale Solar Chimney Using Computational Fluid DynamicsDesign and Numerical Investigation of a Small Scale Solar Chimney Using Computational Fluid Dynamics27128080510.22044/jsfm.2016.805FAJ.PirkandiM.HashemabadiJournal Article20160314In the current study, a solar shimney is analyzed and designed using computational fluid dynamics (CFD). The optimal dimensions of the chimney is determined. The chimney has a rectangular cross section and it is modeled in three-dimenstional space. The chimney has been modeled in various dimensions and with various geometric sections and impact of chimney geometries have been investigated on internal air flow and its performance. In addition to the geometry of the chimney effect, impact of change heat flux received from environment, entrance positive pressure and inlet air temperature were studied on outlet air mass flow rate from the chimney. The results indicated that there is optimum width in order to maximize air flow rate in the chimney. Optimum width is increased with Increase in inlet cross section, entrance pressure and chimney height in addition to increased air mass flow rate and performance improvement. The results of this plan can be used for appropriate design of a solar chimney and finally making it.In the current study, a solar shimney is analyzed and designed using computational fluid dynamics (CFD). The optimal dimensions of the chimney is determined. The chimney has a rectangular cross section and it is modeled in three-dimenstional space. The chimney has been modeled in various dimensions and with various geometric sections and impact of chimney geometries have been investigated on internal air flow and its performance. In addition to the geometry of the chimney effect, impact of change heat flux received from environment, entrance positive pressure and inlet air temperature were studied on outlet air mass flow rate from the chimney. The results indicated that there is optimum width in order to maximize air flow rate in the chimney. Optimum width is increased with Increase in inlet cross section, entrance pressure and chimney height in addition to increased air mass flow rate and performance improvement. The results of this plan can be used for appropriate design of a solar chimney and finally making it.https://jsfm.shahroodut.ac.ir/article_805_c6630d6338f933d6605b1c60f7ef9314.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Elasto-plastic time dependent impact analysis of high speed projectile on water surfaceElasto-plastic time dependent impact analysis of high speed projectile on water surface28129881110.22044/jsfm.2016.811FAH.ForouzaniMalek Ashtar University of technologyB.SaranjamR.KamaliA.Abdollahi FarJournal Article20160426In this paper, elasto-plastic time dependent impact of high-speed projectile on water surface is simulated numerically using Arbitrary Lagrangian- Eulerian (ALE) method. The projectile is considered as elasto-plastic solid and its mesh is generated by Lagrangian approach. The water is also assumed as compressible fluid so its mesh is produced by Eulerian method. Three steps simulation are performed in this research; static, dynamic stress analysis and also impact analysis of full degrees of freedom (DOF) projectile on water surface using ALE method. The effects of fluid compressibility and cavitation are considered in last analysis. In order to validate results, the stress wave propagation produced in the projectile due to water impact is compared with exact ones. The results show that the maximum error compare with exact ones is 5%. Also the magnitude of maximum stress and location/path of fracture in the projectile are compared with experimental data. The good agreement between the predicted and analytical/experimental values shows the accuracy of this numerical algorithm. The impact of projectile on water surface is simulated with different angles. The results show that the safe zone of impact angle for present projectile is ±0.5°.In this paper, elasto-plastic time dependent impact of high-speed projectile on water surface is simulated numerically using Arbitrary Lagrangian- Eulerian (ALE) method. The projectile is considered as elasto-plastic solid and its mesh is generated by Lagrangian approach. The water is also assumed as compressible fluid so its mesh is produced by Eulerian method. Three steps simulation are performed in this research; static, dynamic stress analysis and also impact analysis of full degrees of freedom (DOF) projectile on water surface using ALE method. The effects of fluid compressibility and cavitation are considered in last analysis. In order to validate results, the stress wave propagation produced in the projectile due to water impact is compared with exact ones. The results show that the maximum error compare with exact ones is 5%. Also the magnitude of maximum stress and location/path of fracture in the projectile are compared with experimental data. The good agreement between the predicted and analytical/experimental values shows the accuracy of this numerical algorithm. The impact of projectile on water surface is simulated with different angles. The results show that the safe zone of impact angle for present projectile is ±0.5°.https://jsfm.shahroodut.ac.ir/article_811_e4426375e6692026edfa9d1b3d3d771b.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922Numerical analysis of mechanical properties of mesoporous silica and hydroxyapatite-polypropylene nanocompositesNumerical analysis of mechanical properties of mesoporous silica and hydroxyapatite-polypropylene nanocomposites22930379910.22044/jsfm.2016.799FAA.R.AlbooyehJournal Article20150719This paper presents the numerical study on the mechanical properties of polypropylene (PP) reinforced with 1 wt.% (weight percentage) mesoporous silica (MCM-41), Hydroxy Apatite (HA) and the composite of MCM-41 and HA (MCM41-HA or MH) nanoparticles. Numerical simulation according to macromechanical method (continuum method) was performed for different form of Strain Energy Density Function containing: Marlow, Ogden and Van der Waals using the commercial software ABAQUS/CAE. In this research the hyperelastic behavior of porous polymeric nanocomposits containing: mesoporous silica/ polypropylene (PMCM), hydroxyapatite/polypropylene (PHA) nanocomposites and mesoporous silica- hydroxyapatite /polypropylene (PMH) nanohybrid in the uniaxial tensile test, three point bending test and Izod impact test were studied. In order to verify the numerical results, values obtained from the simulation of mechanical properties were compared With similar functional tests. The performed simulation results are in good accordance with the experimental data and the Marlow model has the best agreement with results of experiment.This paper presents the numerical study on the mechanical properties of polypropylene (PP) reinforced with 1 wt.% (weight percentage) mesoporous silica (MCM-41), Hydroxy Apatite (HA) and the composite of MCM-41 and HA (MCM41-HA or MH) nanoparticles. Numerical simulation according to macromechanical method (continuum method) was performed for different form of Strain Energy Density Function containing: Marlow, Ogden and Van der Waals using the commercial software ABAQUS/CAE. In this research the hyperelastic behavior of porous polymeric nanocomposits containing: mesoporous silica/ polypropylene (PMCM), hydroxyapatite/polypropylene (PHA) nanocomposites and mesoporous silica- hydroxyapatite /polypropylene (PMH) nanohybrid in the uniaxial tensile test, three point bending test and Izod impact test were studied. In order to verify the numerical results, values obtained from the simulation of mechanical properties were compared With similar functional tests. The performed simulation results are in good accordance with the experimental data and the Marlow model has the best agreement with results of experiment.https://jsfm.shahroodut.ac.ir/article_799_d13f23c3a9f648cc0ba82511c34706c6.pdfShahrood University of TechnologyJournal of Solid and Fluid Mechanics2251-94756320160922natural convection heat transfer of a nanofluid in a baffle L-Shaped cavitynatural convection heat transfer of a nanofluid in a baffle L-Shaped cavity31132181410.22044/jsfm.2016.814FAN.AlaviT.ArmaghaniE.IzadpanahJournal Article20160713In this article, free convection heat transfer nanofluids water-alumina in enclosure L-shaped BAFFLED has been investigated. Vertical walls left at warm temperatures and constant horizontal intermediate walls Th and Tc is constant cold temperature. Other walls are well insulated. Baffles at 〖 T〗_c and puts them in a container lot of impact on flow behavior. And could disrupt the current order. The governing equations finite volume method and the SIMPLE algorithm simultaneously solved. The results show: AR increased heat transfer is also increased. By increasing the aspect ratio, nanofluids have a greater impact on increasing the heat transfer. During the baffle increases the heat transfer is increased. To increase as the baffle, nanofluids have less impact on the cooling chamber. The effect of Righly Number and nanofluid volume fraction are also investigated and discussed. the results show: for all AR the natural heat transfer increases when nanofluid volume fraction is increased.In this article, free convection heat transfer nanofluids water-alumina in enclosure L-shaped BAFFLED has been investigated. Vertical walls left at warm temperatures and constant horizontal intermediate walls Th and Tc is constant cold temperature. Other walls are well insulated. Baffles at 〖 T〗_c and puts them in a container lot of impact on flow behavior. And could disrupt the current order. The governing equations finite volume method and the SIMPLE algorithm simultaneously solved. The results show: AR increased heat transfer is also increased. By increasing the aspect ratio, nanofluids have a greater impact on increasing the heat transfer. During the baffle increases the heat transfer is increased. To increase as the baffle, nanofluids have less impact on the cooling chamber. The effect of Righly Number and nanofluid volume fraction are also investigated and discussed. the results show: for all AR the natural heat transfer increases when nanofluid volume fraction is increased.https://jsfm.shahroodut.ac.ir/article_814_cd6fce0966c1d9f6343febe4cb54c4c7.pdf