Journal of Civil and Surveying Engineering44120100321Semi-Analytic Method for a Tri-Material Transversely Isotropic Full-Space in Frequency DomainSemi-Analytic Method for a Tri-Material Transversely Isotropic Full-Space in Frequency Domain20748FAAzizollahBehrestaghiMortezaEskandariJournal Article19700101A linear elastic transversely isotropic full-space containing three different regions, which are an upper half-space, a lower half-space and a layer in between is considered in such a way, each region contains different material. The axes of symmetry of different regions are assumed to be normal to the interface of the regions and thus parallel. An arbitrary load in frequency domain is applied on an arbitrary patch located at the interface of the upper half-space and its underneath layer. Integral formulations are presented for the determination of the displacements and stresses in all regions. By means of the complete displacement potentials introduced by Eskandari-Ghadi (2005), Fourier series in circumferential direction and Hankel transform in radial direction, the displacements and stresses are determined in Fourier-Hankel space. Then, the displacements and stresses at any point are given in the line integral form. The solution is numerically evaluated for: (i) a half-space under an arbitrary surface load, (ii) a half-space under an arbitrary buried force, (iii) a half-space fixed at the top and under an arbitrary buried force, (iv) a half-space containing a layer bonded to the top of a half-space under an arbitrary force applied at the interface of two regions, (v) a full-space under an arbitrary load in it, (vi) a bi-material full-space under an arbitrary force applied at the interface of two half-space, and (vii) a layer of finite thickness fixed at the bottom and under an arbitrary surface load. To confirm the accuracy of the numerical evaluation of the integrals involved, some of the results are compared with existing solutions.A linear elastic transversely isotropic full-space containing three different regions, which are an upper half-space, a lower half-space and a layer in between is considered in such a way, each region contains different material. The axes of symmetry of different regions are assumed to be normal to the interface of the regions and thus parallel. An arbitrary load in frequency domain is applied on an arbitrary patch located at the interface of the upper half-space and its underneath layer. Integral formulations are presented for the determination of the displacements and stresses in all regions. By means of the complete displacement potentials introduced by Eskandari-Ghadi (2005), Fourier series in circumferential direction and Hankel transform in radial direction, the displacements and stresses are determined in Fourier-Hankel space. Then, the displacements and stresses at any point are given in the line integral form. The solution is numerically evaluated for: (i) a half-space under an arbitrary surface load, (ii) a half-space under an arbitrary buried force, (iii) a half-space fixed at the top and under an arbitrary buried force, (iv) a half-space containing a layer bonded to the top of a half-space under an arbitrary force applied at the interface of two regions, (v) a full-space under an arbitrary load in it, (vi) a bi-material full-space under an arbitrary force applied at the interface of two half-space, and (vii) a layer of finite thickness fixed at the bottom and under an arbitrary surface load. To confirm the accuracy of the numerical evaluation of the integrals involved, some of the results are compared with existing solutions.https://jcse.ut.ac.ir/article_20748_3c7e157d3c7385ee96d0618e743160c8.pdfJournal of Civil and Surveying Engineering44120100321An Investigation on Steel Knee Braced Frames for Lateral LoadsAn Investigation on Steel Knee Braced Frames for Lateral Loads20749FAMeysamJalaliSMehdiZahraeiSRasoolMirghaderiJournal Article19700101Hysteretic dampers as the most practical passive control methods have improved the seismic performance of steel structures through energy dissipation induced by large inelastic deformations. Among them in this paper, KBF (Knee Braced Frames), a special form of braced frames having a diagonal brace connected to a knee element instead of beam-column joint, is investigated. In such systems, the knee element designed to have inelastic behavior for dissipating seismic energy caused by strong earthquake. In this article, analytical investigation on cyclic performance of knee braced single span one-story frame specimens are described. Nonlinear finite element method have been used for analysis of KBFs. Results have shown good behavior of knee member before its failure and also indicated a number of important considerations for seismic design of KBFs. Moreover hysteretic behavior have shown no pinching effects. At last, a mathematical model has been derived to predict the bilinear behavior of the KBFs and compared with monotonic pushover analysis. Acceptable conformity between mathematical and numerical model proved the accuracy of the proposed model.Hysteretic dampers as the most practical passive control methods have improved the seismic performance of steel structures through energy dissipation induced by large inelastic deformations. Among them in this paper, KBF (Knee Braced Frames), a special form of braced frames having a diagonal brace connected to a knee element instead of beam-column joint, is investigated. In such systems, the knee element designed to have inelastic behavior for dissipating seismic energy caused by strong earthquake. In this article, analytical investigation on cyclic performance of knee braced single span one-story frame specimens are described. Nonlinear finite element method have been used for analysis of KBFs. Results have shown good behavior of knee member before its failure and also indicated a number of important considerations for seismic design of KBFs. Moreover hysteretic behavior have shown no pinching effects. At last, a mathematical model has been derived to predict the bilinear behavior of the KBFs and compared with monotonic pushover analysis. Acceptable conformity between mathematical and numerical model proved the accuracy of the proposed model.https://jcse.ut.ac.ir/article_20749_24ba949ef9bf43a8f295adf9d1a76afb.pdfJournal of Civil and Surveying Engineering44120100321Estimation of Creep Compliance of Hot Mix Asphalt by Artificial Neural NetworksEstimation of Creep Compliance of Hot Mix Asphalt by Artificial Neural Networks20750FAAbolfazlHasaniAliHeydaripanahJournal Article19700101Creep compliance is one of the fundamental tests of Mechanistic- Empirical flexible pavement design procedure in the AASHTO 2002 design guide. In this research, a new Artificial Neural Network model for estimating the HMA creep compliance with the generalization ability of R=0.949 has been developed successfully using feed forward multi layer perceptron Artificial Neural Networks (ANNs) with Levenberg-Marquardt backpropagation training algorithm. The ANN model has 14 inputs including selected percent passings from aggregate gradation curve, asphalt content, asphalt penetration degrees at 77 and 125 Fahrenheit , asphalt penetration index and test temperature. The model has 8 outputs consist of creep compliances of Hot Mix Asphalt at 1, 2, 5, 10, 20, 50, 100 seconds and Poisson coefficient. All needed data has been extracted from LTPP database 2007 by developing many programs by SQL in Access 2000. As the result, 975 creep compliance tests data and 975 associated aggregate tests data series and 975 associated asphalt tests data series has been extracted successfully for ANNs training and evaluation.Creep compliance is one of the fundamental tests of Mechanistic- Empirical flexible pavement design procedure in the AASHTO 2002 design guide. In this research, a new Artificial Neural Network model for estimating the HMA creep compliance with the generalization ability of R=0.949 has been developed successfully using feed forward multi layer perceptron Artificial Neural Networks (ANNs) with Levenberg-Marquardt backpropagation training algorithm. The ANN model has 14 inputs including selected percent passings from aggregate gradation curve, asphalt content, asphalt penetration degrees at 77 and 125 Fahrenheit , asphalt penetration index and test temperature. The model has 8 outputs consist of creep compliances of Hot Mix Asphalt at 1, 2, 5, 10, 20, 50, 100 seconds and Poisson coefficient. All needed data has been extracted from LTPP database 2007 by developing many programs by SQL in Access 2000. As the result, 975 creep compliance tests data and 975 associated aggregate tests data series and 975 associated asphalt tests data series has been extracted successfully for ANNs training and evaluation.https://jcse.ut.ac.ir/article_20750_c4b24f502bfc5a4e3e4a9ef21032e242.pdfJournal of Civil and Surveying Engineering44120100321Dynamic Response of Unlined Circular Tunnels Subjected to Internal ExplosionDynamic Response of Unlined Circular Tunnels Subjected to Internal Explosion20751FAMohammadRahimianBabakOmidvarBenyaminKyomarsiAlirezaSanaeihaJournal Article19700101https://jcse.ut.ac.ir/article_20751_554cd7594d2b4fcf17d44d1f945c0710.pdfJournal of Civil and Surveying Engineering44120100321Analysis of Blast Shock Waves on Immersed PipesAnalysis of Blast Shock Waves on Immersed Pipes20752FAHamedShahmohammadiSoheilMohammadiJournal Article19700101https://jcse.ut.ac.ir/article_20752_781765eb8582da1fe0798e0463c59e2e.pdfJournal of Civil and Surveying Engineering44120100321Determination of River Regime Equations Based on Stream Power EquationDetermination of River Regime Equations Based on Stream Power Equation20753FAAhmadTahershamsiFarhadImanshoarJournal Article19700101In this paper, river regime equations through developing of unit stream power equation and based on data of Sistan and Karaj Rivers were derived. At the first step, methods and out-comes of former researchers are classified. On the next step, based on dimensional analysis, unit stream power as a function of river flow parameters including suspended sediment load and specific gravity of particles is developed. Using the unit stream power equation and Sistan and Karaj Rivers' data, regime equations are derived for these rivers. Finally, the results derived from regime equations are compared with the field data and regime equations of the other researchers. A good agreement was observed between field data and values resulted from developed regime equations in this research.In this paper, river regime equations through developing of unit stream power equation and based on data of Sistan and Karaj Rivers were derived. At the first step, methods and out-comes of former researchers are classified. On the next step, based on dimensional analysis, unit stream power as a function of river flow parameters including suspended sediment load and specific gravity of particles is developed. Using the unit stream power equation and Sistan and Karaj Rivers' data, regime equations are derived for these rivers. Finally, the results derived from regime equations are compared with the field data and regime equations of the other researchers. A good agreement was observed between field data and values resulted from developed regime equations in this research.https://jcse.ut.ac.ir/article_20753_a4b4a702386d6e53a5084baf7d6993c2.pdfJournal of Civil and Surveying Engineering44120100321Three-Dimensional Slope Stability Analysis with Rotational Failure MechanismThree-Dimensional Slope Stability Analysis with Rotational Failure Mechanism20754FAFarajollahAskariOrangFarzanehMasoudFazeliniaJournal Article19700101This paper deals with three-dimensional evaluation of the slopes on the basis of limit analysis. An upper-bound technique of limit analysis is used in this paper to determine either the bearing capacity of a shallow foundation near a slope or the safety factor of the slope. The theorems of limit analysis (upper and lower bound) provide a powerful tool for solving problems in which limit loads need to be found. According to the upper-bound theorem, for a kinematically admissible velocity field an upper bound of the collapse load can be obtained by equating the power dissipated internally in an increment of displacement to the power expended by the external loads. Such kinematically admissible velocity fields have to comply with the kinematical boundary conditions and compatibility conditions. A rotational mechanism consisting of three slip surfaces is used to determine the factor of safety for a slope or the ultimate limit load of a foundation with eccentric load near a slope. This mechanism includes two lateral surfaces and a surface at bottom. As the soil is assumed to obey the associated flow rule, the angle between the relative velocity vector and velocity discontinuity has to be equal everywhere with soil's friction angle (?).
The geometries of lateral surfaces are represented by a nonlinear differential equation. A surface with log-spiral section is used for bottom surface of the mechanism. To obtain the minimum upper-bound, an algorithm is proposed to optimize the failure mechanism. Results from this algorithm for typical conditions are comparable to those of other existing methods. Dimensionless diagrams for various conditions are presented which can be used to predict the bearing capacity of a foundation with eccentric load located on a slope.This paper deals with three-dimensional evaluation of the slopes on the basis of limit analysis. An upper-bound technique of limit analysis is used in this paper to determine either the bearing capacity of a shallow foundation near a slope or the safety factor of the slope. The theorems of limit analysis (upper and lower bound) provide a powerful tool for solving problems in which limit loads need to be found. According to the upper-bound theorem, for a kinematically admissible velocity field an upper bound of the collapse load can be obtained by equating the power dissipated internally in an increment of displacement to the power expended by the external loads. Such kinematically admissible velocity fields have to comply with the kinematical boundary conditions and compatibility conditions. A rotational mechanism consisting of three slip surfaces is used to determine the factor of safety for a slope or the ultimate limit load of a foundation with eccentric load near a slope. This mechanism includes two lateral surfaces and a surface at bottom. As the soil is assumed to obey the associated flow rule, the angle between the relative velocity vector and velocity discontinuity has to be equal everywhere with soil's friction angle (?).
The geometries of lateral surfaces are represented by a nonlinear differential equation. A surface with log-spiral section is used for bottom surface of the mechanism. To obtain the minimum upper-bound, an algorithm is proposed to optimize the failure mechanism. Results from this algorithm for typical conditions are comparable to those of other existing methods. Dimensionless diagrams for various conditions are presented which can be used to predict the bearing capacity of a foundation with eccentric load located on a slope.https://jcse.ut.ac.ir/article_20754_9a85642014d10213bf5db283800e2456.pdfJournal of Civil and Surveying Engineering44120100321Mitigation Response of Steel Structure using Shape Memory Alloy DampersMitigation Response of Steel Structure using Shape Memory Alloy Dampers20755FAMehdiGhasemiehSArdavanMotahariJournal Article19700101https://jcse.ut.ac.ir/article_20755_b87905c3b8d88228116d38f06f8359b0.pdfJournal of Civil and Surveying Engineering44120100321Application of Seismic Inversion Method for Calculation of Porosity in the Saraje Gas FieldApplication of Seismic Inversion Method for Calculation of Porosity in the Saraje Gas Field20756FAMehdiHazinehMohammadMokhtariJournal Article19700101The direct way to access the petrophysical property of hydrocarbon reservoir is based on well data but due to its high cost, the seismic data has been used to obtain these parameters. For this purpose, many different methods are available and each has its own advantages and limitations. From these methods the seismic inversion is the best tool to study the reservoir characterization. The outcome of the seismic inversion is the transformation of seismic reflection data to pseudo-impedance logs at each CMP. It is therefore, the seismic inversion results facilitate better estimations of reservoir properties such as porosity. In this study the results of two different post-stack inversion methods on seismic data will compared. These are sparse-spike and model-based methods. These results were used for acoustic impedance construction based on 2D seismic data in the Saraje Gas Field. By using this model and creation of linear regression the estimated porosity model has been achieved.The direct way to access the petrophysical property of hydrocarbon reservoir is based on well data but due to its high cost, the seismic data has been used to obtain these parameters. For this purpose, many different methods are available and each has its own advantages and limitations. From these methods the seismic inversion is the best tool to study the reservoir characterization. The outcome of the seismic inversion is the transformation of seismic reflection data to pseudo-impedance logs at each CMP. It is therefore, the seismic inversion results facilitate better estimations of reservoir properties such as porosity. In this study the results of two different post-stack inversion methods on seismic data will compared. These are sparse-spike and model-based methods. These results were used for acoustic impedance construction based on 2D seismic data in the Saraje Gas Field. By using this model and creation of linear regression the estimated porosity model has been achieved.https://jcse.ut.ac.ir/article_20756_a52cfd7f0a3d9b2386f2981c44636630.pdfJournal of Civil and Surveying Engineering44120100321Modal Parameter Identification using Ambient Vibration Testing with Introducing a new SoftwareModal Parameter Identification using Ambient Vibration Testing with Introducing a new Software20757FAJavadYadegariOmidBaharJournal Article19700101The ambient vibration measurement is a kind of output data-only dynamic testing where the traffics and winds are used as agents responsible for natural or environmental excitation. Therefore an experimental modal analysis procedure for ambient vibration testing needs to base itself on output-only data. The modal analysis involving output-only measurements present a challenge that requires the use of special modal identification technique, which can deal with very small magnitude of ambient vibration contaminated by noise. Two complementary modal analysis methods are implemented. They are rather simple peak picking (PP) method in frequency domain and more advanced stochastic subspace identification (SSI) method in time domain. According to these methods Graphical User Interface (GUI) software developed in MATLAB. This program is called SIP and can be used in modal parameter identification of large civil structures such as bridges and towers. This paper presents the application of ambient vibration testing and experimental modal analysis on large civil engineering structures and also gives brief review of program. Also modal parameters of the Babolsar arch steel bridge are extracted using SIP and ambient vibration data from dynamic test carried out on this bridge.The ambient vibration measurement is a kind of output data-only dynamic testing where the traffics and winds are used as agents responsible for natural or environmental excitation. Therefore an experimental modal analysis procedure for ambient vibration testing needs to base itself on output-only data. The modal analysis involving output-only measurements present a challenge that requires the use of special modal identification technique, which can deal with very small magnitude of ambient vibration contaminated by noise. Two complementary modal analysis methods are implemented. They are rather simple peak picking (PP) method in frequency domain and more advanced stochastic subspace identification (SSI) method in time domain. According to these methods Graphical User Interface (GUI) software developed in MATLAB. This program is called SIP and can be used in modal parameter identification of large civil structures such as bridges and towers. This paper presents the application of ambient vibration testing and experimental modal analysis on large civil engineering structures and also gives brief review of program. Also modal parameters of the Babolsar arch steel bridge are extracted using SIP and ambient vibration data from dynamic test carried out on this bridge.https://jcse.ut.ac.ir/article_20757_3d9aeb0190b0670ae014bff1fb3f135b.pdfJournal of Civil and Surveying Engineering44120100321A Combined Model to Calculate and Manage of Leakage in Water Distribution NetworksA Combined Model to Calculate and Manage of Leakage in Water Distribution Networks20758FAAmirhoseinAsadiyaniMasoudTabeshJournal Article19700101Leakage in water distribution networks causes waste of water and economical losses. In this paper a new method is presented to calculate nodal and pipes leakage in water distribution systems using capabilities of EPANET2.10 software. The method estimate the total water losses in the system based on IWA methodology using Minimum Night Flow (MNF) measurements. Recognising the pressure dependency of leakage the total water consumption is divided into two parts, one pressure dependent and the other independent of local pressure. The total leakage is then distributed in different nodes based on nodal pressures and then leakage at each pipe is calculated. Furthermore, the effects of leakage management are evaluated and leakage indices are determined. Based on the outputs of a case study in Ilam Province, the capabilities of the proposed methodology is evaluated. The results show that by combining the new method and hydraulic models the system can be simulated and nodal and pipe leakage are obtained using hydraulic parameters. It can also be observed that leakage performance indices have a great role in leakage management, prioritisation of loss reduction scenarios and budget consumption in water companies. Besides, these indices show the degree of success for leak reduction schemes. Finally, it was found that the infrastructure leakage index (ILI) is the most effective factor on the leakage and the best one to evaluate the leakage performance.Leakage in water distribution networks causes waste of water and economical losses. In this paper a new method is presented to calculate nodal and pipes leakage in water distribution systems using capabilities of EPANET2.10 software. The method estimate the total water losses in the system based on IWA methodology using Minimum Night Flow (MNF) measurements. Recognising the pressure dependency of leakage the total water consumption is divided into two parts, one pressure dependent and the other independent of local pressure. The total leakage is then distributed in different nodes based on nodal pressures and then leakage at each pipe is calculated. Furthermore, the effects of leakage management are evaluated and leakage indices are determined. Based on the outputs of a case study in Ilam Province, the capabilities of the proposed methodology is evaluated. The results show that by combining the new method and hydraulic models the system can be simulated and nodal and pipe leakage are obtained using hydraulic parameters. It can also be observed that leakage performance indices have a great role in leakage management, prioritisation of loss reduction scenarios and budget consumption in water companies. Besides, these indices show the degree of success for leak reduction schemes. Finally, it was found that the infrastructure leakage index (ILI) is the most effective factor on the leakage and the best one to evaluate the leakage performance.https://jcse.ut.ac.ir/article_20758_cea0be3b1f4525585e882ecfb550246a.pdfJournal of Civil and Surveying Engineering44120100321Correction of Records of Important Earthquakes of IRAN with Wavelet Transform and Filtering Method and Comparing between these MethodsCorrection of Records of Important Earthquakes of IRAN with Wavelet Transform and Filtering Method and Comparing between these Methods20759FAHesamVahidifardAnooshiravanAnsariAsadollahNoorzadJournal Article19700101During an earthquake, many stations record the ground motion, but only a few numbers of them could be corrected using conventional high-pass and low-pass filtering method and the others identified as highly contaminated by noise and as a result useless. There are two major problems associated to these noisy records. First, since signal to noise ratio (S/N) is low, it is not possible to discriminate between original signal and noise neither in frequency domain nor in time domain. The consequence is that it is not possible to cancel out the noise using conventional filtering methods. The second problem is about non-stationary characteristics of the noise. In other words, in many cases the characteristics of the noise are varying by time and in these situations, it is not possible to apply frequency domain correction schemes. For correction of acceleration signals contaminated with high-level non-stationary noise, there is an important question that is it possible to have estimation about the state of the noise in different bands of time and frequency? The wavelet multi-resolution analysis which decomposes a signal into different time-frequency components besides a suitable criterion for identification of the noise among each component provides the required mathematical tool for correction of highly noisy acceleration records. In this paper, the strong motion records of great earthquake events with Mw>5.5 in Iranian Plateau is corrected. For this purpose, first, signal to noise ratio is calculated for each record. If the value of this diagram is high enough, the signal is corrected with both methods and otherwise, the de-noising is applied using wavelet procedure. The effect of baseline correction of the final results is also investigated for each record. A comprehensive database containing the information of events, stations and correction parameters is also provided in this study.During an earthquake, many stations record the ground motion, but only a few numbers of them could be corrected using conventional high-pass and low-pass filtering method and the others identified as highly contaminated by noise and as a result useless. There are two major problems associated to these noisy records. First, since signal to noise ratio (S/N) is low, it is not possible to discriminate between original signal and noise neither in frequency domain nor in time domain. The consequence is that it is not possible to cancel out the noise using conventional filtering methods. The second problem is about non-stationary characteristics of the noise. In other words, in many cases the characteristics of the noise are varying by time and in these situations, it is not possible to apply frequency domain correction schemes. For correction of acceleration signals contaminated with high-level non-stationary noise, there is an important question that is it possible to have estimation about the state of the noise in different bands of time and frequency? The wavelet multi-resolution analysis which decomposes a signal into different time-frequency components besides a suitable criterion for identification of the noise among each component provides the required mathematical tool for correction of highly noisy acceleration records. In this paper, the strong motion records of great earthquake events with Mw>5.5 in Iranian Plateau is corrected. For this purpose, first, signal to noise ratio is calculated for each record. If the value of this diagram is high enough, the signal is corrected with both methods and otherwise, the de-noising is applied using wavelet procedure. The effect of baseline correction of the final results is also investigated for each record. A comprehensive database containing the information of events, stations and correction parameters is also provided in this study.https://jcse.ut.ac.ir/article_20759_31302016532bb363509863e4c21fe808.pdf