• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.2
• /
• pp.215-227
• /
• 2020

A new multiscale finite element method for elliptic problems with highly oscillating coefficients are introduced. A hybridization yields a locally flux-conserving numerical scheme for multiscale problems. Our approach naturally induces a homogenized equation which facilitates error analysis. Complete convergence analysis is given and numerical examples are presented to validate our analysis.

• Journal of Industrial Technology
• /
• v.25 no.A
• /
• pp.191-201
• /
• 2005

A numerical analysis was done for a tunnel fire in a 1000m road tunnel. A cartesian coordinate was adopted to make a computational grid sytem which has 448,000 computational cells. A transient flow phenomena in the tunnel was simulated by the commercial code of PHEONICS from the ignition of fire to 600 seconds by the interval of 100 seconds. Total computational time of about 44 hours was required to get a convered solution in each time step. The purpose of this research is to analyze of the backlayering pheonomena and recirculation flow in a tunnel. The compuational results say that the backlayering does not happens near the fire of vehicle in this case because the vehicle fire is located at the outside of recirculation zone of flow ocuured near the jet fan. In this research, onset of backlayering pheonomena could be escaped if jet fan is set 95m in front of the the fire of vehicle.

• Computational Structural Engineering
• /
• v.9 no.2
• /
• pp.115-120
• /
• 1996

The main concern of numerical dynamic analysis of large structures is to find an acceptable solution with fewer mode shapes and less computational efforts. Component mode method utilizes substructure technique to reduce the degree of freedom but have a disadvantage to not consider the dynamic characteristics of loads. Ritz Vector method consider the load characteristics but requires many integrations and errors are accumulated. In this study, to improve the effectiveness of component mode method, Lanczos algorithm is introduced. To prove the effectiveness of this method, example structure are analyzed and the results are compared with SAP90.

• Computers and Concrete
• /
• v.10 no.3
• /
• pp.307-330
• /
• 2012

This work intends to contribute for the improvement of the procedure suggested by Brazilian Technical Code that takes into account the cracked concrete stiffness in the estimative of the displacement of reinforced concrete beams submitted to service loads. A damage constitutive model accounting for induced anisotropy, plastic deformations and bimodular elastic response is used in order to simulate the concrete behaviour, while an elastoplastic behaviour is admitted for the reinforcement. The constitutive models were implemented in a program for bars structures analysis with layered finite elements. Initially, the damage model is briefly presented as well as the parametric identification of the materials that have been used in the reinforced concrete beams. After that, beams with different geometries and reinforcement area are analyzed and a statistical method (ANOVA) is employed in order to identify the main variables in the problem. Soon after, the same procedure is used with another resistance of concrete, where the compression strength is changed. The numerical responses are compared with the ones obtained by Brazilian Technical Code and experimental tests in order to validate the use of the damage model. Finally, some remarks are discussed based on responses presented in this work.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.395-409
• /
• 2022

The presence of infill generally neglected in design despite the fact that infill contribution significantly increase the lateral stiffness and strength of the reinforced concrete frame structure. Several experimental studies and computational models have been proposed to capture the rational response of infill-frame interaction at global level. However, limited studies are available on explicit finite element modelling to study the local behavior due to high computation and convergence issues in numerical modelling. In the current study, the computational modelling of RC frames is done with various configurations of infill masonry in terms of types of blocks, lateral loading and reinforcement detailing employed with material nonlinearities, interface contact issues and bond-slip phenomenon particularly near the beam-column joints. To this end, extensive computational modelling of five variant characteristics test specimens extracted from the detailed experimental program available in literature and process through nonlinear static analysis in FEM code, ATENA generally used to capture the nonlinear response of reinforced concrete structures. Results are presented in terms of damage patterns and capacity curves by employing the finest possible detail provided in the experimental program. Comparative analysis shows that good correlation amongst the experimental and numerical simulated results both in terms of capacity and crack patterns.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.182-187
• /
• 2003

Numerical analysis of three-dimensional vicous flow is used to compute the design speed operating line of a transonic axial-flow compressor. The Navier-Stokes equation was solved by an explicit finite-difference numerical scheme and the Baldwin-Lomax turbulence model was applied. A spatially-varying time-step and an implicit residual smoothing were used to improve convergence. Two-stage axial compressor of a turboshaft engine developed KARI was chosen for the analysis. Numerical results show reasonably good agreements with experimental measurements made by KARI. Numerical solutions indicate that there exist a strong shock-boundary layer interaction and a subsequent large flow separation. It is also observed that the shock is moved ahead of the blade passage at near-stall condition.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.596-601
• /
• 2007

A simplified finite element analysis method is proposed to calculate elasto-plastic responses of general hardening materials. The method provides an effective tool to calculate structural elasto-plastic responses. Numerical examples have demonstrated that its computational efficiency is very much higher than that of the incremental elasto-plastic finite element analysis, and computational results are accurate enough to meet the need of engineering practice. Compared with the general elasto-plastic incremental finite element analysis, the proposed method can avoid the incremental iteration of nodal displacements and the constitutive equation integration at each Gauss integral point, and computational results are accurate enough to meet the need of engineering practice.

• Korean Journal of Computational Design and Engineering
• /
• v.11 no.5
• /
• pp.384-392
• /
• 2006

In this research, the feasibility of distributed finite element analysis system with STEP and CORBA has been investigated. The enabling technologies such as CORBA and Java play key roles in the development of integrated and geographically distributed application software. In addition to the distribution of analysis modules, numerical solution process itself is again divided into parallel processes using multi-frontal method for computational efficiency. In contrast to the specially designed parallel process for specific hardware, CORBA-based parallel process is well suited for heterogeneous platforms over the network. The idea of Web-based distributed analysis system may be applied to the engineering ASP for design and analysis in the product development processes. We believe that the proposed approach for the analysis can be extended to the entire product development process for sharing and utilizing common product data in the distributed engineering environment, thus eventually provide basis for virtual enterprise.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.20-26
• /
• 2008

A numerical study is carried out to analyze the steady three-dimensional turbulent flow and convective heat transfer in a staggered pin-fin array with diamond shaped elements at various geometrical configurations. Steady Reynolds-averaged Navier-Stokes equations and energy equation are solved using a finite volume based solver. Shear stress transport (SST) model is used as turbulence closure. The computational domain is composed of one pitch of pin-fin displacement with periodic boundary conditions on the surfaces normal to the streamwise direction and the cross-streamwise direction. The numerical results for Nusselt number and friction factor are validated with experimental results. The effects of pin angle, pin height and pitch on Nusselt number, friction factor and efficiency index are investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.17-22
• /
• 1999

Convergence characteristics and efficiency of three implicit approximate factorization schemes(ADI, DDADI and MAF) are examined using 2-Dimensional compressible upwind Navier-Stokes code. Second-order CSCM(Conservative Supra Characteristic Method) upwind flux difference splitting method with Fromm scheme is used for the right-hand side residual evaluation, while generally first-order upwind differencing is used for the implicit operator on the left-hand side. Convergence studies are performed using an example of the flow past a NACA0012 airfoil at steady transonic flow condition, i. e. Mach number 0.8 at $1.25^{\circ}$ angle of attack. The results were compared with other computational results in order to validate the current numerical analysis. The results from the implicit AF algorithms were compared well in low surface with the other computational results; however, not well in upper surface. It might be due to lack of the grid around the shock position. Because the algorithm minimizes the errors of the approximate decomposition, the improved convergence rate with MAF were observed.