---

 Investigations on vibration behaviors of railway track systems were attempted in this research. This was made by conducting a comprehensive field investigation into the free vibration of track systems and response of tracks to train moving loads. In-situ modal analysis was used in a railway track field as an efficient method of investigating dynamic properties of railway track systems. Natural frequencies and mode shapes of the track system in different insitu track conditions were obtained for the fist time. The sensitivity of the natural frequencies of the track to the types of sleepers, fastening systems, ballast conditions, and rail joints were studied. Efficiency of rail welded joints in CWR tracks and the effects of replacing timber sleepers with concrete sleepers on dynamic behavior of a track were investigated. Advantages of flexible sleeper fastening system from the aspects of serviceability and passenger riding comfort were discussed. The effects of the track accumulative loading as a main indicator of ballast degradation on track dynamic behavior were studied. Rail deflections were calculated by using auto-spectra obtained from vibrations of the track under trainloads, leading to the development of a new mathematical expression for the calculation of the rail dynamic amplification factor.

---

In this paper, the problem of layout optimization for X-bracing of steel frames is studied using the ant system (AS). A new design method is employed to share the gravity and the lateral loads between the main frame and the bracings according to the requirements of the IBC2006 code. An algorithm is developed which is called optimum steel designer (OSD). An optimization method based on an approximate analysis is also developed for layout optimization of braced frames. This method is called the approximate optimum steel designer (AOSD) and uses a simple deterministic optimization algorithm leading to the optimum patterns and it is much faster than the OSD. Several numerical examples are treated by the proposed methods. Efficiency and accuracy of the methods are then discussed. A comparison is also made with Genetic algorithm for one of the frames.

---

In the current study, an effort is made to determine three dimensional bearing capacity of rectangular foundations using Discrete

Element Method. The soil mass is modeled as discrete blocks connected with Winkler springs. Different factors affect the geometry

of failure surface. Six independent angles are used to define the failure surface. By trial and error, the optimum shape of failure

surface beneath the foundation can be found. The paper includes the derivation of the governing equations for this DEM based

formulation in three dimensional state as well as parametric sensitivity analyses and comparison with other methods. Moreover,

using the current method, bearing capacity coefficients are presented for various friction angles and foundation aspect ratios.

---

This paper describes a numerical model and its finite element implementation that used to compute the cavitation effects on

seismic behavior of concrete dam and reservoir systems. The system is composed of two sub-systems, namely, the reservoir and

the dam. The water is considered as bilinear compressible and inviscid and the equation of motion of fluid domain is expressed

in terms of the pressure variable alone. A bilinear state equation is used to model the pressure–density relationship of a cavitated

fluid. A standard displacement finite element formulation is used for the structure. The Structural damping of the dam material

and the radiation damping of the water and damping from foundation soil and banks have been incorporated in the analysis. The

solution of the coupled system is accomplished by solving the two sub-systems separately with the interaction effects at the damreservoir

interface enforced by a developed iterative scheme. The developed method is validated by testing it against problem for

which, there is existing solution and the effects of cavitation on dynamic response of Konya gravity dam and Morrow Point arch

dam subjected to the first 6 s of the May 1940 El-Centro, California earthquake, is considered. Obtained results show that impact

forces caused by cavitation have a small effect on the dynamic response of dam-reservoir system.

---

It has been realized that the raft (mat) foundations are capable of bearing very large loads when they are assisted with a pile
group. The contribution of both raft and piles to carry the surcharge loads is taken into account, considering the stiffness and
strength of involved elements in the system, i.e. piles, raft and surrounding soil. The piles are usually required not to ensure the
overall stability of the foundation but to act as settlement reducers. There is an alternative design in which, the piles are nonconnected
from the raft to reduce the settlement, which are then known to be "settlement reducer non-connected piles" to increase
the system stiffness. In this paper, two and three dimensional finite element analysis of connected and non-connected pile-raft
systems are performed on three case studies including a 12-storey residential building in Iran, a 39-storey twin towers in
Indonesia, and the Messeturm tower, 256m high, in Frankfurt, Germany. The analyses include the investigation of the effect of
different parameters, e.g. piles spacing, embedment length, piling configuration and raft thickness to optimize the design. The role
of each parameter is also investigated. The parametric study results and comparison to a few field measurements indicate that
by concentrating the piles in the central area of the raft foundation the optimum design with the minimum total length of piles is
achieved, which is considered as control parameter for optimum design. This can be considered as a criterion for project cost
efficiency. On the other hand, non-connected piled-raft systems can significantly reduce the settlements and raft internal bending
moments by increasing the subsoil stratum stiffness. Finally, the comparison indicates that simple and faster 2D analysis has
almost similar results to the time consuming and complicated 3D analysis.

---

A series of tests and also numerical analyses were conducted to explore the mechanical behavior of a mixture of coarse gravelsize
particles floating in a matrix of silt, sand or clay. The research is a step forward in an ongoing investigation on behavior of
composite clay, which is used as the core material of some large embankment dams all over the world. After providing the reader
with an overall image about behavior of such materials through the literature, the paper focuses on a predominant feature of the
composite soil behavior: increase of non-deformable solid inclusions in a mixture leads to formation of heterogeneity of stress
field, excess pore water pressure and strain distribution along the specimens. This paper mainly probes formation of such
heterogeneity by the aid of special experiments and also numerical analyses. In addition to loading details, it is clarified through
the paper that position of inclusions relative to loading direction also affects heterogeneity of stress/strain and excess pore water
pressure distribution through the mixture. Despite the former, the latter redistributes with a rate proportional to material
hydraulic conductivity.

---

Presented is a method of three-dimensional stability analysis of convex slopes in plan view based on the Lower-bound theorem of
the limit analysis approach. The method’s aim is to determine the factor of safety of such slopes using numerical linear finite
element and lower bound limit analysis method to produce some stability charts for three dimensional (3D) homogeneous convex
slopes. Although the conventional two and three dimension limit equilibrium method (LEM) is used more often in practice for
evaluating slope stability, the accuracy of the method is often questioned due to the underlying assumptions that it makes. The
rigorous limit analysis results in this paper together with results of other researchers were found to bracket the slope stability
number to within ±10% or better and therefore can be used to benchmark for solutions from other methods. It was found that using
a two dimensional (2D) analysis to analyze a 3D problem will leads to a significant difference in the factors of safety depending
on the slope geometries. Numerical 3D results of proposed algorithm are presented in the form of some dimensionless graphs which
can be a convenient tool to be used by practicing engineers to estimate the initial stability for excavated or man-made slopes

---

Endurance Time (ET) method is a response history based analysis procedure that can be used for estimating the seismic response

of structures at different excitation levels in each response history. This seismic analysis method utilizes specific intensifying

acceleration functions to analyze seismic behaviors. One of the potential applications of the ET method is in the seismic

assessment of shell structures. In this study, a procedure for linear seismic analysis of shell structures is proposed and

applications of this method is investigated for several cases of shell structures. These structures are analyzed under three ET

acceleration functions in one direction and the results are compared to time history analysis considering seven actual earthquake

records. Moreover, the results of the ET method are compared to response spectrum analysis method. The outcomes of the study

reveal that the ET method predicts the linear seismic performance of shell structures with acceptable precision and significant

reduction in analysis time. Furthermore, it is conluded that scattering of results of three ET analysis is very low and one analysis

can be used instead of three. Finally, the comparison between THA and RSM results verify that response spectuarm method is a

conservative method which occasionally encounters problems to evaluate bending stresses of shell structures

---

In order to lighten the prestressed concrete solid members, nowadays, it is possible to make use of the advantage of HPC (fc'>60

MPa) as well as replacing the solid section with a PSC thin-walled section for certain members such as circular and box columns.

Using the strength theory of ACI, a numerical procedure along with a computer program was developed for the analysis of such

sections subjected to axial compression or tension load and bending moments. The program solves for all possible variables such

as, concrete compressive strength (fc'= 60-100 MPa), type of prestressed steel, concrete cover, ratio of wall thickness to the section

dimensions and the PS steel arrangements to satisfy the given loading cases, thus leading to an optimal cost solution. However,

since the cross section is thin-walled circular or box and the PS steel is located at discrete points along the periphery of a circle

or rectangle, the equations of equilibrium are complex for hand computations (especially for circular section) but suitable for

computer program. So, by use of MATLAB software the interaction diagrams were also drawn for the analysis of such sections

for all mentioned variables. The use of prestressed thin-walled column diagrams is a safe and easy tool for the analysis of such

columns. Finally, the accuracy of the proposed method is demonstrated by comparing its results to those of the available

experimental values and is indicate that the proposed method predict very well the capacity of prestressed thin-walled column.

---

The nonlinear static procedures (NSPs) proposed by design codes do not lead to reliable results especially for tall buildings.

They generally provide inconsistent estimates of inelastic seismic demands, especially for the top floors due to their inabilities in

considering the higher modes effects. In this paper, a new enhanced pushover procedure is proposed which is based on the

envelope of the structural responses resulting from two separate pushover analyses as a combination rule. Also, the suggested

pushover analyses are performed using a newly proposed modal load pattern, i.e., the Modal Spectra Combination (MSC), and

the ASCE41-06 required first mode load pattern. The MSC load pattern is consisted of a number of mode shapes combined with

appropriate weighting factors that depend on their modal participation factors, modal frequencies and design spectral values. A

number of 2-D steel moment resisting frame models with different number of stories are used to investigate the efficiency of the

proposed method. The inter-story drifts and the maximum plastic beam moment and curvature responses are used as a measure

to compare the results obtained from the nonlinear time-history analyses (NL-THA) and some other NSPs. The results obtained

through rigorous nonlinear dynamic analyses show that the application of the proposed method leads to acceptable results for

steel MRF systems in comparison to other available enhanced NSPs. The OpenSees program is used for numerical analysis.

---

Double- unsymmetric-plan medium-rise buildings subjected to bi-directional seismic excitation are complex structures where higher-mode effects in plan and elevation are important in estimating the seismic responses using nonlinear static or pushover analysis. Considering two horizontal components of the ground motions makes the problem more intricate. This paper presents a method for nonlinear static analysis of double unsymmetric-plan low- and medium-rise buildings subjected to the two horizontal components of ground motions. To consider bi-directional seismic excitation in pushover analyses, the proposed method utilizes an iterative process until displacements at a control node (centre of mass at the roof level) progressively reach the predefined target displacements in both horizontal directions. In the case of medium-rise buildings, continuous implementation of modal pushover analyses is used to take higher-mode effects into account. To illustrate the applicability and to appraise the accuracy of the proposed method, it is applied to the 4- and 10-storey torsionally-stiff and torsionally-flexible buildings as representative of low- and medium-rise buildings, respectively. For the purpose of comparison, modal pushover analysis (MPA) is also implemented considering the two horizontal components of the ground motions. The results indicate that the proposed method and the MPA procedure can compute the seismic demands of double unsymmetric-plan low- and medium-rise buildings with reasonable accuracy however, seismic responses resulting from the proposed method deteriorate at the flexible edge of the torsionally-flexible buildings

---

Seismic ground motion is profoundly affected by geometrical and mechanical properties of soil deposits overlaying bedrock.

Local seismic ground response of saturated soil deposits was studied in literature by applying the effects of soil stress state

and index properties on the strain-dependent normalized shear modulus reduction, G/G0, and damping ratio, D, curves in an

equivalent linear analysis. However, experimental investigations revealed that, G0, G/G0, and D of unsaturated soils are

influenced by stress state as well as suction. This study presents the results of linear and equivalent linear seismic ground response

analysis of unsaturated soil deposits incorporating suction effects on G/G0 and D curves. Seismic ground response analyses were

done with the computer program EERA for three sets of soil profiles, which are included in saturated, constant and linearly

variable suction unsaturated soil deposits. The results of current study present the magnitude of variation in natural frequency,

amplification ratio and spectral acceleration of unsaturated soil deposits.

---

Pile foundations are frequently used in industrial projects in southwest lowlands of Iran. Although high setup of shaft resistance

is usually reported in the area, no reliable formulation or guidelines are available for considering the increased capacity in design

applications. Therefore, the pile design practices are usually not optimized. The main objective of this paper is presenting a site

specific formulation for setup effects of a utility plant in southwest Iran in which a good database of prestressed concrete driven

piles is available. Fajr-II Petrochemical site in PetZone of Mahshahr accommodating a utility plant is selected as the database of

the current study. The setup factor (A) and the reference time (t0) are evaluated through processing of a relatively large database

of this well-supervised piling project. As the main portion of variations of driven piles capacity with time is related to shaft, only

shaft resistance variations are considered in this research. The shaft capacity variations are derived from signal matching analysis

on PDA tests. Reliability of PDA tests has been confirmed through comparing with the static load test results. Influence of driving

the surrounding piles on setup factor is also investigated. The results show that the average setup factor (A) and the reference time

(t0) of 0.30 and 0.01 day, respectively, are proper values for estimating the long term capacity in this region. Evaluation of the

results indicates that driving 8 piles around the test pile has increased the “A” factor average of 40% resultingin increase of the

shaft capacity about 19% in one month and 22% in one year, in comparison with the tested piles with no surrounding piles driven.

---

In the process of structural analysis we often come to structures that can be analyzed with simpler methods than the standard approaches. For these structures, known as regular structures, the matrices involved are in canonical forms and their eigen-solution can be performed in a simple manner. However, by adding or removing some elements or nodes, such methods cannot be utilized. Here, an efficient method is developed for the analysis of irregular structures in the form a regular structure with additional or missing nodes or with additional or missing supports. In this method, the saving in computational time is considerable. The power of the method becomes more apparent when the analysis should be repeated very many times as it is the case in optimal design or non-linear analysis.

---

High performance and reliability of refurbish able knee braced steel frames has been confirmed in previous researches trying to get an optimal design for its configuration. Buckling of diagonal member which affects the hysteretic behavior of KBF under cyclic loadings has not been foreseen in previous evaluations of this system. This deficiency can be improved by utilization of adjustable rotary friction damper device (FDD) as knee element. Diagonal element buckling can be prevented considering a suitable value for FDD sliding threshold moment Mf. Lower values of Mf Lower energy dissipation rate in FDD and this leads to an optimization problem. Nonlinear time history analyses have been performed in addition to lateral cyclic loading analyses to evaluate the response of single story KBF subjected to seismic excitation. Optimal Mf in FDD has been chosen according to these analyses results. Roof displacement and acceleration, base shear and diagonal element’s buckling status have been compared in optimally designed KBF and FDD utilized KBF (FKBF) with different configurations. Nonlinear dynamic analyses have been performed for one, four, eight and twelve story frames under different seismic records with several PGAs. More than 60% displacement response reduction has been earned for the FKBF without considerable increase in base shear.

---

Dynamic analysis of the seismic performance of power substation equipment is time-consuming, expensive and uses responses that are sensitive to ground motion. This research proposes a method to derive input waves for dynamic analysis in place of original records from seismic events in Iran. In this study, a power transformer, current transformer, circuit breaker and disconnect switch are analyzed using fifty records from the far-field and near-field earthquake ground motions. Statistical analysis is done on the maximum acceleration and displacement responses to obtain their pushover curves. Sinusoidal waves were created using the fundamental frequencies of the equipments and PGA of 0.1g through 0.5 g as the amplitude. The results are compared with the original records and show that the proposed input waves provide a reasonable fit for an extensive range of near-field and far-field ground motion results.

---

Setback in elevation of a structure is a special irregularity with considerable effect on its seismic performance. This paper addresses multistory Reinforced Concrete (RC) frame buildings, regular and irregular in elevation. Several multistory Reinforced Concrete Moment Resisting Frames (RCMRFs) with different types of setbacks, as well as the regular frames in elevation, are designed according to the provisions of the Iranian national building code and Iranian seismic code for the high ductility class. Inelastic dynamic time-history analysis is performed on all frames subjected to ten input motions. The assessment of the seismic performance is done based on both global and local criteria. Results show that when setback occurs in elevation, the requirements of the life safety level are not satisfied. It is also shown that the elements near the setback experience the maximum damage. Therefore it is necessary to strengthen these elements by appropriate method to satisfy the life safety level of the frames.

---

Soil reinforced with fiber shows characteristics of a composite material, in which fiber inclusion has a significant effect on soil permeability. Concerning to the higher void ratio of carpet fibers, at first stages it may be expected that an increase in fiber content of the reinforced soil would result in an increase in permeability of the mixture. However, the present article demonstrates that fiber inclusion will decrease the permeability of sand-fiber composite.A series of constant head permeability tests have been carried out to show the effects and consequently, a new system of phase relationships was introduced to calculate the dry mass for the sand portion of the composite. Monte Carlo simulation technique adopted with finite element theory was employed to back calculate the hydraulic conductivity of individual porous fibers from the laboratory test results. It was observed that the permeability coefficient of the porous fibers are orders of magnitude less than the skeletal sand portion due to the fine sand particle entrapment and also the fiber volume change characteristics.

---

In saturated soils, heating induces thermal expansion of both grains and the pore fluid. Lower thermal expansion coefficient of aggregates results in the increase of pore pressure and reduction of the effective stress besides subsequent volume changes due to the dissipation of pore pressure and heat transfer. Dissipation of thermally induced pore pressure with time is a coupled thermo-hydro-mechanical (THM) phenomenon, involving gradients of pore pressure and temperature, hydraulic and thermal flows within the mass of soil and changes in the mechanical properties with temperature. The objective of this paper is presentation of a numerical method to determine the effect of temperature on consolidation of clays. In this regard, the finite element code, PISA is used for one dimensional THM analysis of porous media. The analysis performed using both linear elastic and elastoplastic Cam clay models. Modified Cam clay model was applied in elastoplastic analysis. Variation of temperature, displacements and pore pressure determined with time and compared with numerical solutions of other researchers. Also it was indicated that implementation of coupled THM analysis yields better results for displacements compared to the hydro mechanical (HM) one. Application of elastoplastic constitutive model instead of linear elastic one indicated that preconsolidation pressure has an important effect on results of analysis.

---

In the present study, the application of Endurance Time Analysis (ETA) method is investigated on seismic analysis of a high arch dam. In this method the coupled system is excited using the predesigned intensifying acceleration functions instead of the real ground motions. Finite element model of an arch dam considering the dam-rock-water interaction effects was developed in which the concrete and rock were assumed to have linear elastic behavior. The effect of the large displacement in dam body was considered in numerical model using co-rotational approach. The coupled system was analyzed by conventional Time History Analysis (THA) method in various seismic performance levels and the results were compared with those obtained from ETA at the equivalent target time. It was found that ETA method provides the close results to THA with acceptable accuracy while it reduces the total time of the analyses considerably.