• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.717-733
• /
• 2003

A time domain method is presented for soil-structure interaction analysis under seismic excitations. It is based on the finite element formulation incorporating infinite elements for the far field soil region. Equivalent earthquake input forces are calculated based on the free field responses along the interface between the near and far field soil regions utilizing the fixed exterior boundary method in the frequency domain. Then, the input forces are transformed into the time domain by using inverse Fourier transform. The dynamic stiffness matrices of the far field soil region formulated using the analytical frequency-dependent infinite elements in the frequency domain can be easily transformed into the corresponding matrices in the time domain. Hence, the response can be analytically computed in the time domain. A recursive procedure is proposed to compute the interaction forces along the interface and the responses of the soil-structure system in the time domain. Earthquake response analyses have been carried out on a multi-layered half-space and a tunnel embedded in a layered half-space with the assumption of the linearity of the near and far field soil region, and results are compared with those obtained by the conventional method in the frequency domain.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.5
• /
• pp.483-492
• /
• 2009

In this study, nonlinear analysis of steel plane frame was performed using the refined plastic hinge method of advanced analysis techniques. In deterioration of stiffness in plastic zone, influences by flexural bending, residual stress, geometrical non-linearity, and semi-rigid connection are considered. And also, further reduced tangent modulus was used for geometrical non-linearity, top and seat angle were chosen for semi-rigid connection. Furthermore, 3 parameter power model was used for moment-rotation behaviour of beam to column connection. The loading conditions are combined with axial and lateral force and the inverse triangle distribution of lateral and eight type of analytical models were used in analysis. The results of analyses were compared with semi-rigid and rigid connection on behaviour of numerical analysis models. And also, the behaviors of frame with changes of semi-rigidity were analyzed by using the results obtained from MIIDAS-GENw.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.298-309
• /
• 1990

In this paper, a new method is suggested to analyze impulsive stresses at loading poing of concentrated impact load under certain impact conditions determined by impact velocity, stiffness of plate and mass of impact body, etc. The impulsive stresses are analyzed by using the three dimensional dynamic theory of elasticity so as to analytically clarify the generation phenomenon of cone crack at the impact fracture of fragile materials (to be discussed if the second paper). The Lagrange's plate theory and Hertz's law of contact theory are used for the analysis of impact load, and the approximate equation of impact load is suggested to analyze the impulsive stresses at the impact point to decide the ranage of impact load factor. When impact load factors are over and under 0.263, approximate equations are suggested to be F(t)=Aexp(-Bt)sinCt and F(t)=Aexp(-bt) {1-exp(Ct)} respectively. Also, the inverse Laplace transformation is done by using the F.F.T.(fast fourier transform) algorithm. And in order to clarity the validity of stress analysis method, experiments on strain fluctuation at impact point are performed on a supported square glass plate. Finally, these analytical results are shown to be in close agreement with experimental results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.249-257
• /
• 2022

This paper presents an improved computational framework for the direct-solution-based finite element tearing and interconnecting (FETI) algorithm. The FETI-local algorithm is further improved herein, and localized Lagrange multipliers are used to define the interface among its subdomains. Selective inverse entry computation, using a property of the Boolean matrix, is employed for the computation of the subdomain interface stiffness and load, in which the original FETI-local algorithm requires a full matrix inverse computation of a high computational cost. In the global interface computation step, the original serial computation is replaced by a parallel multi-frontal method. The performance of the improved FETI-local algorithm was evaluated using a numerical example with 64 million degrees of freedom (DOFs). The computational time was reduced by up to 97.8% compared to that of the original algorithm. In addition, further stable and improved scalability was obtained in terms of a speed-up indicator. Furthermore, a performance comparison was conducted to evaluate the differences between the proposed algorithm and commercial software ANSYS using a large-scale computation with 432 million DOFs. Although ANSYS is superior in terms of computational time, the proposed algorithm has an advantage in terms of the speed-up increase per processor increase.

• Composites Research
• /
• v.30 no.4
• /
• pp.247-253
• /
• 2017

In this study, the mechanical behavior and interface properties of boron nitride nanotube-poly(methyl methacrylate) nanocomposites are predicted using the molecular dynamics simulations and the double inclusion model. After modeling nanocomposite unit cell embedding single-walled nanotube and polymer, the stiffness matrix is determined from uniaxial tension and shear tests. Through the orientation average of the transversely isotropic stiffness matrix, the effective isotropic elastic constants of randomly dispersed microstructure of nanocomposites. Compared with the double inclusion model solution with a perfect interfacial condition, it is found that the interface between boron nitride nanotube and polymer matrix is weak in nature. To characterize the interphase surrounding the nanotube, the two step domain decomposition method incorporating a linear spring model at the interface is adopted. As a result, various combinations of the interfacial compliance and the interphase elastic constants are successfully determined from an inverse analysis.

• Journal of Biomedical Engineering Research
• /
• v.24 no.6 s.81
• /
• pp.515-522
• /
• 2003

Pharynx is a system transporting foods by peristaltic motion(contraction and expansion movement! into the esophagus and functioning as airway passages. In this study, structural changes of pharyngeal dysfunction are analyzed by biomechanical model using CT and FEM(finite clement method). Loading condition was assumed that equal pressure was loaded sequentially to inside of pharyngeal tissue. In order to analyze the pharyngeal muscular dysfunction by biomechanical model. the pharyngeal dysfunctions was classified into 3 cases. Taking into account the clinical complication by neuromuscular symptoms such as pharyngeal dysfunction after stroke. we assumed that a change of material property is caused by muscular tissue stiffness. A deformation of cross sectional area of the pharynx is analyzed increasing the stiffness $25\%,\;50\%,\;75\%$ in each case on the basis of stress-strain relationship. Based on three-dimensional reconstruction of pharyngeal structure using limited factor - techniques and the optimization procedure by means of inverse dynamic approach. the biomechanical model of the human pharynx is implemented. The results may be used as clinical index illustrating the degree of pharyngeal muscular dysfunction. This study may be used as useful diagnostic model in discovering early deglutitory impediment caused by physiological or pathological pharyngeal dysfunction.

• Computers and Concrete
• /
• v.25 no.3
• /
• pp.215-224
• /
• 2020

In the present study, according to the important of porosity in low specific weight in comparison of high stiffness of carbon nanotubes reinforced composite, buckling and free vibration analysis of sandwich composite beam in two configurations, of laminates using differential quadrature method (DQM) is studied. Also, the effects of porosity coefficient and three types of porosity distribution on critical buckling load and natural frequency are discussed. It is shown the buckling loads and natural frequencies of laminate 1 are significantly larger than the results of laminate 2. When configuration 2 (the core is made of FRC) and laminate 1 ([0/90/0/45/90]_s) are used, the first natural frequency rises noticeably. It is also demonstrated that the influence of the core height in the case of lower carbon volume fractions is negligible. Even though, when volume fraction of fiber increases, the critical buckling load enhances smoothly. It should be noticed the amount of decline has inverse relationship with the beam aspect ratio. Investigating three porosity patterns, beam with the distribution of porosity Type 2 has the maximum critical buckling load and first natural frequency. Among three elastic foundations (constant, linear and parabolic), buckling load and natural frequency in linear variation has the least amount. For all kind of elastic foundations, when the porosity coefficient increases, critical buckling load and natural frequency decline significantly.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.5 s.98
• /
• pp.527-535
• /
• 2005

This paper explains a general coupling system in terms of the system parameters. impedance of a cavity or mobility of a structure. To easily access the mechanism of the structural-acoustic coupled system, a simple expression is derived. A general coupled equation is also derived of a general coupled problem constituted a flexible structure and an opening boundary in terms of vector and matrix notation, and is analyzed the coupling phenomena using the understanding acquired simple coupled system. The paper shows that the general coupled equation is expanded version of the simple coupled equation by some limiting checks. The paper also shows that the degree of coupling is proportioned to a stiffness of the acoustic system and a modal coupling coefficient, but is in inverse proportion to a mass of the structural system and the difference of the excitation frequency and resonant frequency of the acoustic or structural system.

• Geomechanics and Engineering
• /
• v.13 no.3
• /
• pp.475-488
• /
• 2017

Geotechnical data contributes substantially to the cost of engineering projects due to increasing cost of site investigations. Existing information in the form of soil maps can save considerable time and expenses while deciding the scope and extent of site exploration for a proposed project site. This paper presents spatial interpolation of data obtained from soil investigation reports of different construction sites and development of soil maps for geotechnical characterization of Multan area using ArcGIS. The subsurface conditions of the study area have been examined in terms of soil type and standard penetration resistance. The Inverse Distance Weighting method in the Spatial Analyst extension of ArcMap10 has been employed to develop zonation maps at different depths of the study area. Each depth level has been interpolated as a surface to create zonation maps for soil type and standard penetration resistance. Correlations have been presented based on linear regression of standard penetration resistance values with depth for quick estimation of strength and stiffness of soil during preliminary planning and design stage of a proposed project in the study area. Such information helps engineers to use data derived from nearby sites or sites of similar subsoils subjected to similar geological process to build a preliminary ground model for a new site. Moreover, reliable information on geometry and engineering properties of underground layers would make projects safer and economical.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.8
• /
• pp.41-46
• /
• 2006

The capability for reconstructing impact forces of an aircraft wing using impact response functions and regularization methods were examined. The impact response function which expresses the relation between the structure response and the impact force was derived using the information on mass and stiffness data of a finite element model for the wing. Iterative Tikhonov regularization method and generalized singular value decomposition method were used to inverse the impact response function that was generally ill-posed. For the numerical verification, a fighter aircraft wing was used. Strain and deflection histories obtained from finite element analysis were compared with the results calculated using impact response functions. And the impact forces were reconstructed with the strain histories obtained from finite element analysis. The numerical verification results showed that this method can be used to monitor impact forces on aircraft structures.