Optimization problems may be devided into geometry optimization problems and topology optimization problems. In this paper, a method using tile equivalent material properties prediction techniques of a particulate-reinforced composites is proposed for the topology optimization. This method makes use of penalty factor in order that regions with intermediate value of design variables can be penalized. The computational results being obtained from PLBA algorithm of some values of penalty factor are presented.

In this paper the stochastic analysis of semi-infinite domain is presented using the weighted integral method, which is expanded to include the infinite finite elements. The semi-infinite domain can be thought as to have more uncertainties than the ordinary finite domain in material constants, which shows the needs of and the importance of the stochastic finite element analysis. The Bettess's infinite element is adopted with the theoretical decomposition of the strain matrix to calculate the deviatoric stiffness of the semi-infinite domains. The calculated value of mean and the covariance of the displacement are revealed to be larger than those given by the finite domain assumptions giving the rational results which should be considered in the design of structures on semi-infinite domains.

The isoparametric spline finite strip method for degenerated shells is presented. In the formulation, both the geometry and the displacement field are represented by uniform cubic B-spline curves. Spline shell strip is degenerated stress-resultant shell with 6 dof at each node;and the penalty function method is used to incorporated the six dof,

A versatile flat shell element has been developed by combining a membrane element with drilling degree-of-freedom and a plate bending element. This element is formulated by the enhanced displacement field with the additional non-conforming displacement modes. Thus the element possesses six degrees-of-freedom (DOF) per node which permits an easy connection to other six DOF elements as well as the improvement of the element behavior. In plate bending part, this element is established by the combined use of the addition of non-conforming modes, the reduced (or selective) integration scheme, and the construction of the substitute shear strain fields. The achieved improvement may be attributable to the fact that the merits of these individual techniques are merged into the new element in a complementary manner. In membrane part, this element shows better membrane behavior as the nonconforming displacement mode is added to drilling mode.

Response spectrum analysis method(RSA) rather than time history analysis method(THA) continues to e used by the profession for evaluating maximum dynamic responses of structures subjected to earthquake excitations. Nevertheless, this simple and practical method can cause significant errors in some cases with unproper modal combination method and so on. To obtain more exact responses based n RSA many studies have been carried out considering displacement of top story, base shear and overturning moment. The purpose of this study is to verify error characteristics in RSA with respect to various responses including displacement shear force and overturning moment of each story. It's shown that RSA appears to yield underestimated responses when compared to THA calculations. Also, errors involved in RSA computations grow with an increase in total number of stories.

The current design equations for ultimate strength of tubular joints are based on a limited number of experimental results performed on simple joints with simple loading conditions and depend on value of the branch to the chord diameter- ratio $\beta$ too much. Therefore, the purpose of this study is to estimate the ultimate strength of CHS tilbular joints considering the effects of branch inclination angles $\theta$, chord length to diametel ratio $\alpha$ and chord end conditions by finite element analysis. The analyses are performed using finite element software ADINA that is capable of modeling elasto-plastic material behavior as well as geometric nonlinearities. The results show that the current use of sin $\theta$ in normalized design equations for inclined branches is reasonable, but somewhat conservative. When compared with the previous experimental database, the close numerical results are obtained from the parametric studies on the static strength of T-, Y-, DT- and X-joints. Also, a new design equation for ultimate stregth of CHS tubular joints is derived using a modified version of the ring model which can include the effects of $\alpha$ and chord end condtion.

Efficiency of design process for large scale structures highly depends on the development of efficient structural analysis and structural response control algorithms because a successful design involves a number of structural analysis based on iterative structural response control process. In this paper, distributed structural analysis model on multiple personal computers connected by ethernet network is presented. To reduce communication cost required in the process of analysis, substructuring techniques are adopted to evenly distribute computational loads on each processor. With its applications on structural analysis of plane frame structures, performance of the proposed computational model are presented in detail.

It has become important to evaluate the qualitive reliability and condition assessment of existing structural systems in order to establish a rational program for repair and maintenance. Since most of if existing structural system may suffer from defect corrosion and damage, it is necessary to account for their effects in fuzzy reliability analysis, In this paper, an attempt is made to develope a reliability analysis for damaged structural systems using failure possibility theory. Damage state is specified in terms of linguistic valiables using natural language information and numerical information, which are defined by fuzzy sets. Using a subjective condition index of failure possibility and information of the damage state is introduced into the calculation of failure probability. The subjective condition index of quantitative and qualitative analysis method is newly proposed based on the fuzzy set operations, namely logical product, drastic product, logical sum and drastic sum

Due to the drastic improvement of computer hardware and operating system, it is easy to break through the main defects of limited computer memory and processing time, etc. To keep up with this situation, this paper is focused on developing the preprocessor program with the input method based on vector graphic editor and the preprocessing technique including automatic node generation algorithm for the $\rho$-version finite element program. To develop this preprocessor program, the special data structure and the OOP(Object Oriented Programming) have been used by the Visual Basic 4.0. The Special data structure is proposed to describe the geometric data of node numberings and coordinates suitable for the $\rho$-version finite element program, which are quite different from the comvential h-version finite element program.

The most of the design engineers for construction has academic background of bachelors degree. Theories for advanced composite structures are too difficult for such engineers and some simple but accurate enough methods are necessary. The senior author has reported that some laminate orientations have decreasing values of B$_{16}$, B$_{26}$, D$_{26}$, and D$_{26}$ stiffnesses as the ply number increases. Such plates above behave as special orthotropic plates and simple formulas developed by the senior author［1, 3］ can be used. Most of the bridge and building slabs on girders have large aspect ratios. For such cases further simplification is possible by neglecting the effect of the longitudinal moment terms(M$_{x}$) on the relevant partial differential equations of equilibrium. In this paper, the result of the study on the subject problem is presented.ted.d.

A numerically stable technique to remove tile limitation in choosing a shift in the subspace iteration method with shift is presented. A major difficulty of the subspace iteration method with shift is that because of singularity problem, a shift close to an eigenvalue can not be used, resulting in slower convergence. This study selves the above singularity problem using side conditions without sacrifice of convergence. The method is always nonsingular even if a shiht is an eigenvalue itself. This is one of tile significant characteristics of the proposed method. The nonsingularity is proved analytically. The convergence of the proposed method is at least equal to that of the subspace iteration method with shift, and the operation counts of above two methods are almost the same when a large number of eigenpairs are required. To show the effectiveness of the proposed method, two numerical examples are considered

Nonlinear finite element analyses of one dimensional consolidation problem were performed using an anisotropic hardening constitutive model. For the analyses, the anisotropic hardening elasto-plastic constitutive model based on the generalized isotropic hardening(GIH) rule was implemented into a nonlinear finite element analysis program, PLASTIC. In order to preserve the accuracy of the finite element solution for nonlinear problems, an implicit stress integration algorithm was employed. A consistent tangent moduli could also ensure the quadratic convergence of Newton's method. As a result, the nonlinear solution was accurately calculated and was converged to be asymptotically quadratic. In a consolidation problem, the relationship between load and settlement and between settlement and time vertical was analyzed comparing with results using the Cam-clay type model and the final consolidation settlement and the duration of primary consolidation could be evaluated rigorously using the GIH constitutive model.

This monograph is based on the purpose of exploiting the program, which let designers to easily compile the input data of a widely used FEM analysis program, SAP90 - one of structural analysis programs, and making it easy to design slab bridges by studying more efficient analysis method with optimization method. The knowledge of position at which a maximum moment occurs is helpful to do a practical design. Seeing that more time and capital should be invested when engineers depend on manual work in the process of making input data repeatedly, the automatic process would increase the degree of reliability and productivity.

The method of vibration analysis used is the one developed by the senior author. He developed and reported, in 1974, a simple but exact method of calculating the natural frequency of beam and tower structures with irregular cross-sections and attached mass/masses. Since 1989, this method has been extended to two-dimensional problems with several types of given conditions and has been reported at several international conferences. This method uses the deflection influence surfaces. The finite difference method is used for this purpose, in this paper. In order to reduce the pivotal points required, the three simultaneous partial differential equations of equilibrium with three dependent variables, w, M$_{x}$, and $M_{y}$, are used instead of the one forth order partial differential equation. By neglecting the M$_{x}$ terms, the size of the matrices needed to solve the resulting linear equations are reduced to two thirds of the "non-modified" equations.tions.

New Thrie-Beam guardrail section has been developed. The Characteristic of its geometry, energy absorbing capability and response to impact has been studied and compared with those of conventional W-Beam guardrail system. To compare the response to impact computer simulation using Barrier Ⅶ program was made. Stretch tests and static loading tests were conducted for the performance verifications.

In this study the fire, due to overturning of oil tanker on the bridge induced heat transfer analysis and thermal stress analysis are carried out. The results of analysis for fire history of 1 hour present very large thermal gradient near the surface. However, the temperature increase of tendon & rebar that is the main resistant member of bridge is not sufficient to change material properties. The Von-Mises yield criteria is used to calculate the depth of delamination, The depth of delamination is about 4cm at center of fire and this value is close to measured value.

Material nonlinear analyses of RC(Reinforced Concrete) beams considering the tension stiffening effect and plastic hinge length have been conducted. Instead of taking the sophisticated layer approach which has some limitations in application to the large structures with many degrees of freedom. the moment-curvature relationships of RC sections previously constructed through the section analysis have been used. To reduce the numerical instability in nonlinear analysis and to remove the imprecision in calculation of ultimate resisting capacity, according to the used finite element mesh size, the tension stiffening effect and plastic hinge length have been taken into consideration. Finally, correlation studies between analytical and experimental results have been conducted with the objective to establish the validity of the proposed algorithms.

An analytical model which can simulate the post-cracking behavior of reinforced concrete structures subjected to in-plane shear and normal stresses is presented. Based on the force equilibriums, compatibility conditions, and bond stress-slip relationship between steel and concrete, a criterion to simulate consider the tension-stiffening effect is proposed. The material behavior of concrete is described by an orthotropic constitutive model, and focused on the tension-compression region with tension-stiffening and compression softening effects defining equivalent uniaxial relations in the axes of orthotropy. Correlation studies between analytical results and available experimental data are conducted with the objective to establish the validity of the proposed model.

This paper presents an optimal decision model for minimizing the life-cycle cost of steel box girder bridges. The point is that it takes into account service life process as a whole, and the life-cycle costs include initial (design, testing, and construction) costs, maintenance costs and expected failure costs. The problem is formulated as that of minimization of expected total life-cycle cost with respect to the design variables. The optimal solution identifies those values of the decision variables that result in minimum expected total cost. The performance constraints in the form of flexural failure and shear failure are those specified in the design code. Based on extensive numerical investigations, it may be positively stated that the optimum design of steel box girder bridges based on life-cycle cost approach proposed in this study provides a lot more rational and economical design, and thus the proposed approach will propose the development of new concepts and design methodologies that may have important implications in the next generation performance-based design codes and standards.

To study temperature effect of P.C. box girder bridge, field measurement was performed for six months, along with thermal analysis for the newly constructed viaduct of Gangbyun Highway in Seoul. Thermocouples were installed inside and surface of the flange and web of the box and temperature of box section md ambient temperature was measured. Measured environmental data are incorporated in finite element thermal analysis and computed temperature of the section was compared with measured one. Temperature gradient from thermal analysis was compared with Korean Highway Specification(1996)and the New Zealand Ministry of Works and Development code(1976). Thermal stress distribution across the box section was also compared.

It is important to measure e dynamic behaviors of the structures and determine the safety and serviceability of those structures by analyzing the gathered data. It is very much easier and more economical to measure the accelerations than displacements in order to obtain the dynamic responses of a structure, but the physical meanings of the displacements are more clear and definite than those of accelerations. In this study, the algorithm which can evaluate the displacements from the measured acceleration data by integration are developed. The calculated displacements are compared with measured data through the laboratory test and the results represent good agreements. This algorithm is applied to the data of real structures.

An analytical study to trace the nonlinear structural behavior of composite beams is undertaken to include the nonlinear material properties of steel sheeting, reinforcing steel bar and concrete. To trace Moment-curvature relations, sectioning analysis method and two simple formulas are developed. A simple power model which has been originally used to expect the flexural capacity of the beam to column connections is proposed and the second formula is composed of two experimental functions to express the Moment-curvature relation in the elastic and plastic range separately. The load-deflection behavior of the beams has been simulated by the step-by-step numerical integration method and is compared with the test results available.

In this study, I-type girders, main members of a steel composite bridge, are designed by Load and Resistance Factor Design method as well as Allowable Stress Design method. The width, height web thickness and flange thickness of main girders are set as design variables. The design program connects optimization program ADS, which is coded with FORTRAN, and a main program coded with $C^{++}$. In this study, it is shown that in this particular steel composite bridge, the design by The Load and Resistance Factor Design method is more economical than that by The Allowable Stress Design method.d.

The precast prestressed concrete girders of I-type section are frequently employed to design the short-to-medium span bridge. However, its beam depth is greatly increased as its span length is increased over than about 30m. Therefore, the economic and aesthetic effectiveness are rapidly decreased in case of the span length over 30m. The purpose of this paper is to verify the structural safety on the new spliced two span bridge and analyze the variation of member forces and stress distribution according to the construction stages and time. The new spliced technique is performed by partial post tensioning and release in the U-type girders. The structural characteristics of this technique is the introduction of secondary moment to reduce the bending moment by self weight of precast U-type girders constructed in simply supported beam type. So, it is expected that the structural efficiency of this spliced bridge may be improved more than other techniques.

In this study, a geometric nonlinear analysis formulation for spatial frames is developed using the 3D stability functions. For the formulation, the relationships of local and global coordinate systems in force, deformation, and the initial and current configurations of a frame are derived. The force-deformation relationship in global coordinate system is derived as well. The developed formulation is verified in each derivation by reducing the derived equations into 2D equations. The gradual plastification of connections and critical sections can be implemented effectively to this formulation for the complete second order inelastic advanced analysis of spatial frames.

A stability problems of isotropic shells under pure bending is investigated based on the classical shells theory. The governing equations of stability problem presented by Donnell and Love, are developed and the solutions for the cylindrical shells are obtained by using Galerkin method. Bending moment is applied at the ends of the cylindrical shell as a from of distributed load in the shape of sine curve. For the isotropic materials, the result of the general purpose structural analysis program based on the finite element method are compared with the critical moment obtained from the classical shell theories. The critical loads for the cylindrical shells with various geometry can not be evaluated with a simple equation. However, accurate solutions for the stability problems of cylindrical shells can be obtained through the equilibrium equation developed in the study.

Effects of translational and rotational motions of the foundation on the dynamic behaviors of a bridge under seismic excitations are examined by utilizing a simplified 3 degree-of-freedom of system. To consider the nonlinear characteristics of the RC pier, a hysteresis model is adapted, which can simulate the inelastic motion of the pier with the stiffness degradation. From results, the portion of the total displacement due to rotational motion of the foundation becomes larger as applied seismic excitation increases.

The most of the design engineers for construction has academic background of bachelors degree. Theories for advanced composite structures are too difficult for such engineers and some simple but accurate enough methods are necessary. The senior author has reported that some laminate orientations have decreasing values of D$_{16}$, B$_{16}$, D$_{26}$ and B$_{26}$ stiffnesses as the ply number increases. Such plates behave as special orthotropic plates and simple formulas developed by the author can be used. Most of the bridge and building slabs on girders have large aspect ratios. For such cases further simplification is possible by neglecting the effect of the longitudinal moment terms(M$_{x}$) on the relevant partial differential equations of equilibrium. In this paper, the result of the study on the subject problem is presented.ented.d.

Many papers which deal with the dynamic instability for shell-like structures under the step load have been published, but there are few papers which treat the essential phenomenon of the dynamic buckling using the phase plane for investigating occurrence of chaos. Dynamic buckling process in the phase plane is a very important thing for understanding why unstable phenomena are sensitively originated in nonlinear dynamics by various initial conditions. In this study, the direct and the indirect snap-buckling of shallow arches considering geometrical nonlinearity are investigated numerically and compared with the static critical load.

In this paper, for the fuzzy constraints not only fuzziness of the constraints relation but also uncertainties of the response of the structures, allowable limits of the constraints and structural design variables, etc. are considered,. so that the fuzzy optimization of the structures can involve more wide scope of the problem and the fuzzy optimal problem is more generalized. In the decision making of the structural design scheme, every possible cases of the fuzzy variables, random variables and fuzzy-random variables, etc. for the uncertainties of the optimization problem are all considered, so the most general method of the decision making is presented. And a numerical example for the three bar truss is offered to demonstrate the reliability and execution possibility proposed method in this paper.

Spherical space frame structure with triangular network pattern, which has the various characteristics for the mechanic property, a funtional property, an aesthetic property and so on, has often been used as one of the most efficient space structures. It is expected that this type will be used widely in large-span structural roofs. But because this structure is made of network by combination of line elements there me many nodes therefore, the structure behavior is very complicated and there can be an overall collapse of structure by buckling phenomenon if the external force reaches a limitation. This kind of buckling is due to geometric shape, network pattern, the number of layer and so on, of structure. Therefore spherical space frame with triangle network pattern have attracted many designers and researchers attention all over the world. The number of layer of space frame is divided in to the simgle, double, multi layer. That is important element which is considered deeply in the beginning of structural design. The buckling characteristics of single-layer model and double-layer model for the spherical space frame structure with triangular network pattern are evaluated and the buckling loads of these types are compared with investigation their structural efficiency in this study.

Application of the flat shell element with drilling D.O.F to linear buckling analysis of thin-walled structures is presented in this paper. The shell element has been developed basically by combining a membrane element with drilling D.O.F. and Mindlin plate bending element. Thus, the shell element possesses six degrees-of-freedom per node which, in addition to improvement of the element behavior, permits an easy connection to other six degrees-of-freedom per node elements(CLS, Choi and Lee, 1995). Accordingly, structures like folded plate and stiffened shell structure, for which it is hard to find the analytical solutions, can be analyzed using these developed flat shell elements. In this paper, linear buckling analysis of thin-walled structures like folded plate structures using the shell elements(CLS) with drilling D.O.F. to be formulated and then fulfilled. Subsequently, buckling modes and the critical loads can be output. Finally. finite element solutions for linear buckling analysis of folded plate structures are compared with available analytic solutions and other researcher's results.

The object of this study is shape finding and cutting pattern generation of membrane structures under the following assumptions: (1) material is linearly elastic (2) stress state is plane stress. Cable and membrane structures should introduce the nonlinear analysis considering geometric nonlinearity because these structures deform largely under the external loads. The analysis procedure is consisted of three steps considering geometric nonlinearity unlike any other structures. First step is the shape finding analysis to determine the initial equilibrium shape. Second step is the stress-deformation analysis to investigate the behaviors of structures under various external loads. Once a satisfactory shape has been found, a cutting pattern based on the shape finding analysis may be generated from the view point of construction. In this paper, after shape finding analysis, cutting pattern determination procedure using weighted least-square minimization flattening method and some results are presented.

The structural behaviors of a arch composed of flexible membrane are investigated. The membrane is considered as thin shell with internal pressure during FEM analysis by using ABAQUS. In the paper, a wind load and uniformly distributed vertical load are considered. As a vertical load, snow loads including applied over all and half of the structure are introduced. The ends of arch are fixed to the ground. Load-Deflections relationship, buckling mode of the structure are presented.

In the latticed domes which is a set of space frame, buckling is derived if the external force reaches a limitation by the lightness of the material and the minimization of the member section area. these are concerned with a geometric shape, network pattern, the number of layer, and so on. Most of all, the number of layer of the lattice dome is a important factor from the viewpoint of initial and structure design. Therefore this study compared buckling characteristics of single-layer with double-layer latticed domes and investigated the relativity of buckling-stress-ratio and member-density-ratio according to rise ratio to improve that designers could extend the range of .design selection

In this paper, we propose to a method to estimate the elasto-plastic buckling for single layer latticed domes. First, we assume that each member consists of the rigid zone and elastic spring at both end joint, the elastic element and three elasto-plastic spring to judge for yeilding the member. Next, the member which has most influence on buckling for structures is determined by a distributed pattern of the strain energy which is calculated through linear eigenvalue analysis. And then, normalized slenderness ratio of the element is derived considering the axial force at elastic buckling load. Later, we execute elasto-plastic nonlinear analysis that based on loading increasement method and displacement increasement method. From this results, we discusses the effect of the joint rigidity and the half open angle $\theta$$_{0}$ on the buckling strength of single layer lattice domes ; (1) how the joint rigidity contributes to the reduction of buckling loads, (2) how the reduction can be interrelated to compressive strength curves in terms of the generalized slenderness for the member most relevant to the overall buckling of domes.s.

The purpose of this study is to propose the method of determining the initial fabric membrane structures surface and membrane patterning procedures. Tension structure, such as, fabric membrane structures and cable-net, is stabilized by their initial prestress and boundary condition. The process to find initial structural overall shape of tension structures produced by initial prestress called Shape Finding or Shape Analysis. One of the most important factor for the design of membrane structures is to search initial smooth surface, because unlike steel or concrete building elements which resist loads in bending, all tension structure forces are carried within the surface by membrane stress or cable tension. To obtain initial surface of fabric membrane element in large deformation analysis, the membrane element is idealized as cable using a technique with Force-density method. and that result is compared with well-known nonlinear numerical method, such as Newton-raphson method and Dynamic relaxation method. The shape resulting from Force-density method has been dealt with as the initial membrane shape and used patterning procedures.

Contrary to reinforced concrete frames, the beams or girders found in actual steel frames are neither hinged nor fixed at their ends. Instead, they are usually restrained by the columns to which they are attached. Here in this paper, common finite element program for plane frame stress analysis is improved by including the effect of partial fixities of the girder ends. To include the effect of partial fixities of the girder ends, Kim's modified slope-deflection equation is adopted. The stress analysis results obtained by this proposed method are then compared with Kim's example. Finally, method of choosing the most economic girder section for a multi-story frame is suggested through the examples.

This paper presents motion based design methodology for structures. Current design methodologies are primarily strength-based. Such methods are adequate when strength is expected to govern the design. But as the slenderness of structures increases, motion such as displacement and acceleration becomes the dominant criterion. In this paper, a preliminary design approach for beam-type buildings, where motion dominates the design, is discussed by effectively distributing the magnitude of structural stiffness to control the distribution of displacement under service load. This analytic development is illustrated using a cantilever beam as the structure under static loads, free vibration, and forced vibration.

The use of air-supported membrane has considerable attention in recent years. There are, however, a number of problems in the behavior of these structures that have not been fully investigated. For example, the problem of local loafing on such membranes has not been analyzed yet. The paper presents an analysis of internal and external potential energy of a spherical air-supported membrane with vertical longitudinal axis, subjected to local loadings. An internal pressure value of the applied load is established at which tile potential of these structures change positive. During such change the loading portion of the membrane comes into work beyond stable state. The mathematical method is used throughout the paper in obtaining solution. For the mathematical modeling, two assumption are used. One is the theta's elimination and the other is the infinite condition. The paper is illustrating the examples of spherical air membrane dome subjected to local load.

The general method of structural analysis for building structures has been based upon the assumption that all dead loads are imposed on a building simultaneously throughout the entire structure. In reality, buildings are built floor by floor or a few floors at a time. The construction dead load is applied gradually onto the structure as the structure is being erected. The prevailing commercial software for structural analysis used to date have resulted in the representation of inaccurate structural behaviors. The actual construction sequence and the loading of the structure ere not properly represented in the analysis. This paper identifies the source of the errors and develops the algorithm to account for the differential column shortening due to construction dead load based upon a given construction sequence

The purpose of this study is to develop a central database and a database management system to store and manage information systematically from each module of an integrated structural design system. In order to efficiently express structural design process related to the data which is very complex, we used an object-oriented modeling methodology to propose the possibility to apply a database schema for application programs in an integrated system for reinforced concrete structural design. Based on this model, we developed an interface between each module and central database. After modeled by using object modeling technique, the database was mapped by the relational database table. Then the central database and the interface were programed by using Visual C/syo ++/, a windows environmental development tool.

The post-processing modules are parts of an integrated system for reinforced concrete structures. This modules are composed of two modules: member design module and calculation report module. The purpose of this paper is to develope modules that increase efficiency and usefulness of an integrated system used reinforced concrete structures design. The development of post-processing modules is necessary for user to design reinforced concrete structures conveniently and quickly. This modules are connected with central database for the benefit of storing amount of input/output data and being used system with little effort. Post-processing modules used Object-Oriented concepts and techniques include identity, classification, polymorphism, and inheritance. Member design module automatically converts no good members into satisfied members by changing section size or reinforcement bar arrangement. This module can be operated both independent member design modules with user input and a part of integrated system with database input. If user operates member design module, calculation report module is created automatically.

The computational environment in which engineers perform their designs has been rapidly evolved from coarse serial machines to massively parallel machines. Although the recent development of high-performance computers are available for a number of years, only limited successful applications of the new computational environments in computational structural engineering field has been reported due to its limited availability and large cost associated with high-performance computing. As a new computational model for high-performance engineering computing without cost and availability problems, parallel structural analysis models for large scale structures on a network of personal computers (PCs) are presented in this paper. In structural analysis solving routine for the linear system of equations is the most time consuming part. Thus, the focus is on the development of efficient preconditioned conjugate gradient (PCG) solvers on the proposed computational model. Two parallel PCG solvers, PPCG-I and PPCG-II, are developed and applied to analysis of large scale space truss structures.

3D finite element analyses of movable bearing shoe and hinged bearing shoe are performed. The finite element models are built using MSC/PATRAN and analyses are carried out using MSC/NASTRAN. Results are again completely processed using MSC/PATRAN. From the results of the analyses, trends of deformation and stress distribution are reviewed and important factors to consider in the design of bearing shoes are discussed. Furthermore optimum crowning amount for the roller of movable bearing shoe was determined according to the results of the analyses.

This paper describes the development of an expert system designed to help engineers select proper foundation types and construction methods for structures under various subsoil conditions. The system includes geographical and geological data in certain areas of Ulsan and a knowledge base for the selection of foundation types. Geological data, such as boring logs, in the areas were collected and arranged to form the subsoil database in the areas concerned Test borings at 30 holes were carried out for reference and confirmation purposes. The use of this system is twofold. It provides the users with a quick view on the geological situations of the areas concerned, and suggests proper foundation types for the specific spot together with some explanations on the selected foundation types and methods.

One solution to prevent deterioration due to expansion joint and to extend lifetime of short span bridges, is jointless integral abutment bridge. To understand behavior of pile foundation of skewed plate girder bridge with integral abutment, finite element analysis was performed for the model of different skew angle from 90。 to 50。. Comparison of stresses at pile and abutment was made for each case. It is found that effect of temperature change is major factor to influence the behavior of skewed integral abutment bridge.

A plate on the elastic half-space may generally be analyzed by the finite element method. However, there are some difficulties to obtain the flexibility matrix of the foundation based on the Boussinesq's theory. In this study, an efficient numerical procedure which uses the analysis results of the vertical displacements due to the uniformly distributed loading in a circular area is presented. Some numerical examples represent better results than those of numerical integration technique or subsection method especially in the case of irregular mesh pattern.

According to the relative stiffness between the half-space and plate or loading condition, some parts of the plate can be separated from the half-space. The finite element procedure to determine the contact area by considering the distribution of contact pressure between plate and the elastic half-space is developed. The vertical surface displacements of the elastic half-space can be obtained through the integrations of the Boussinesq's solution for a point load. The rectangular plate on the elastic half-space is modeled by the 8-node rectangular and 6-node triangular elements and the Mindlin plate theory is used in oder to consider the transverse shear effect. In this study, the contact area may be determined approximately by the analysis with rectangular elements. From this results, the mesh pattern is modified by using triangular and rectangular elements. The contact area can be determined by the new mesh pattern with a relatively sufficient accuracy.

In this study, an expert system is developed to evaluate a safety of tunnel structures. Using a dynamic finite element analysis module, this expert system predicts dynamic responses of a concrete lining surface which a transient force is applied on and estimates the condition between the concrete lining and surrounding ground. The evaluation parameter values of the module are multi-reflected wave frequency and amplitude of the dynamic responses. The multi-reflected wave frequency represents the depth of concrete lining, and the other parameter, the amplitude of the frequency, is utilized for detecting the internal flaws. A comparison of the dynamic responses between numerical and experimental model test verifies an effectiveness of this system. By this expert system, the safety of tunnel structures and the detection of internal flaws of concrete linings are estimated quantitatively.

Elastomeric bearings are frequently used as a means to isolate structural systems from earthquake loadings. The combination of rubber layers and reinforcing steel shims makes the bearings stiff axially but soft laterally The shear flexibility of these short columns can lead to relatively low buckling loads which may be further reduced when high shear strains are simultaneously imposed. The area reduction formula has been proposed to account for the reduction in buckling load due to shear. The result obtained from the formula is presumed to be conservative but the degree of conservatism is unknown. This paper describes a numerical study which aims at determining the effect of high shear strain on the critical load of rubber bearings. The results from the finite element analysis which accounts for both the material and geometric non-linearities are compared against the theoretical results in order to examine the validity of the theoretical formulas.

An efficient solution method is presented to solve the eigenvalue problem arising in tile dynamic analysis of non-proportionally damped structural systems with multiple or close eigenvalues. The proposed method is obtained by applying the modified Newton-Raphson technique and the orthonormal condition of the eigenvectors to the quadratic eigenvalue problem. Even if the shift value is an eigenvalue of the system, the proposed method guarantees nonsingularity, which is analytically proved. The initial values of the proposed method can be taken as the intermediate results of iteration methods or results of approximate methods. Two numerical examples are also presented to demonstrate the effectiveness of the proposed method and the results are compared with those of the well-known subspace iteration method and the Lanczos method.

One of the most important roles in the nonlinear dynamic structural analysis is to select a proper ground excitation, which dominates the response of a structure. Because of the lack of recorded accelerograms in Korea, a stochastic model of ground excitation with various dynamic properties rather than recorded accelerograms is necessarily required. If all information is not available at site, the information from other sites with similar features can be used by the procedure of seismic hazard analysis. Eliopoulos and Wen identified the parameters of the ground motion model by the empirical relations or expressions developed by Trifunac and Lee. Because the relations used in the parameter identification are largely empirical, it is required to apply the artificial neural networks instead of the empirical model. Additionally, neural networks have the advantage of the empirical model that it can continuously re-train the new recorded data, so that it can adapt to the change of the enormous data. Based on the redefined traditional processes, three neural-networks-based models (FAS_NN, PSD_NN and INT_NN) are proposed to individually substitute the Fourier amplitude spectrum, the parameter identification of power spectral density function and intensity function. The paper describes the first half of the research for the development of Neural-Networks-based model for the generation of an Artificial earthquake and a Response Spectrum(NNARS).

The paper describes the second half of the research for the development of Neural-Networks-based model for the generation of an Artificial earthquake and a Response Spectrum(NNARS). Based on the redefined traditional processes related to the generation of an earthquake acceleration response spectrum and design spectrum, four neural-networks-based models are proposed to substitute the traditional processes. RS_NN tries to directly generate acceleration response spectrum with basic data that are magnitude, epicentral distance, site conditions and focal depth. The test results of RS_NN are not good because of the characteristics of white noise, which is randomly generated. ARS_NN solve this problem by the introduction of the average concept. IARS_NN has a role to inverse the ARS_NN, so that is applied to generate a ground motion accelerogram compatible with the shape of a response spectrum. Additionally, DS_NN directly produces design spectrum with basic data. As these four neural networks are simulated as a step by step, the paper describes the methods to generate a response spectrum and a design spectrum using the neural networks.

Analysis and design of vibrating machine foundations subjected to dynamic loads is a very complex problem. Thus it is difficult to set up an accurate analytical modeling. Generally, the design of foundations for vibrating machines has been performed by the equivalent static analysis which is generally based on engineer's experience and various assumptions The purpose of this study is to develop an integrated system which enables structural engineers to produce results of high quality within a short time in works related to structural analysis and design of foundation for vibrating machines. As the result of this study, level-up of application software is expected as well as improvement of quality in structural engineering and reduction of engineers' effort.

The boundary/finite element analysis for the seismic wave amplifications due to nonhomogeneous alluvial deposits was performed in this study. For numerical analysis, the homogeneous linear elastic soil half-space was modeled by using the 3-node isoparametric boundary elements and the inhomogeneous alluvial deposit was modeled by using the 8-node isoparametric finite elements. The two elements at interface were coupled together by the equilibrium condition of the tractions and the compatibility condition of the displacements. As a prarmetric variable, the incident angle and the dimensionsless frequency of the SH, P and SV-waves and the shear wave velocity ratio and the mass density ratio between the half-space and the alluvial deposit were selected.

The substructuring pseudodynamic test is a hybrid testing method consisting of a numerical simulation of the earthquake response of an analytical model and a loading test of a specimen. The substructuring pseudodynamic testing technique has been applied to various seismic experiments since it has advantages over the shaking table test to study dynamic behaviors of relatively large scale structures. However, experimental errors are inevitable in substructuring pseudodynamic testing. Some of these errors can be monitored during the test, but, due to limitations in control system, they cannot be eliminated. For example, one cannot control exactly the displacements that are actually imposed on the structures at each time step. This paper focuses on a technique to minimize the cumulative effect of such control errors for MDOF system. For this purpose, the modified posterior adjustment of the time increment from a target value $\Delta$t$_{n}$ to an adjusted value is performed to minimize the effect of the control errors for MDOF system.for MDOF system.

The drift and inter-story drift control method for steel structures subjected to seismic forces is formulated into a structural optimization problem in this paper. The formulated optimization problem with constraints on drift, inter-story drifts, and member strengthes are transformed into an unconstrained optimization problem. For the solution of the tranformed optimization problem an searching algorithm based on the gradient projection method utilizing gradient information on eigenvalues and eigenvectors are developed and presented in detail. The performance of the proposed algorithm is demonstrated by application to drift control of a verifying example.

The single layer latticed domes have various behaviors with each geometrical shape and scale, and they are affected by vertical component as well as horizontal component of the dynamic load. And they represent ye different earthquake responses under each ground acceleration compared with another structural systems. Generally, all of the members of latticed domes undergo three dimensional deflections if they are subjected to arbitrary one dimensional horizontal load under earthquake motions. And their response characteristics are very different to their shapes, rise/span ratios, and damping mechanisms. In this study, the earthquake response behavior is verified according to the factor of each shape, rise/span ratio, and damping ratio of latticed domes, which undergo horizontal and vertical earthquake motions by numerical approaches.

In the nonlinear dynamic structural analysis, the given ground excitation as an input should be well defined. Because of the lack of recorded accelerograms in Korea, it is required to generate an artificial earthquake by a stochastic model of ground excitation with various dynamic properties rather than recorded accelerograms. It is well own that earthquake motions are generally non-stationary with time-varying intensity and frequency content. Many researchers have proposed non-stationary random process models. Yeh and Wen (1990) proposed a non-stationary stochastic process model which can be modeled as components with an intensity function, a frequency modulation function and a power spectral density function to describe such non-stationary characteristics. This model is based on the simulation for the strong-motion earthquakes with magnitude greater than approximately 5.0~6.0, because it will be not only expected to cause structural damage but also involved the characteristics of earthquake motions. Also, the recorded earthquake motion within this range are still very scarce in Korea. Thus, it is necessary to verify the model by the application of it to the mid-magnitude (approximately 4.0~6.0) earthquakes actually recorded in domestic or foreign area. The purpose of the paper is to generate an artificial earthquake using the model of Yeh and Wen in the area with low seismicity.

지하철과 같은 지하구조물 설계시 기존 구조물 지하를 관통하는 경우가 있다. 본 논문은 기존 역사의 바닥슬래브 밑에 새로운 지하철용 박스 구조물 상부 슬래브가 위치하게 되는 특별한 경우의 굴착 방법을 제시하고 제시된 굴착방법이 상부 구조물에 어떤 영향을 미치는지를 FLAC을 이용하여 검토해 본 것이다. 기존 구조물이 위치하고 있는 지반이 경암인 경우 발파가 주의 깊게 설계되고 시행된다면 발파로 인한 상부 구조의 진동 영향은 충분히 제어할 수 있는 수준으로 예상되었다.