• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Advances in aircraft and spacecraft science
• /
• v.10 no.2
• /
• pp.179-202
• /
• 2023

This article presents a numerical analysis to investigate the natural frequencies and harmonic response of a perforated cantilever beam attached to two layers of piezoelectric materials by using the finite element method for the first time. The bimorph piezoelectric is composed of 3 layers; two of them at the outer are piezoelectric, and the inner isotropic material. A higher order 3-D 20-node solid element that exhibits quadratic displacement behavior is exploited to discretize the isotropic layer, and coupled piezoelectric 3D element with twenty nodes is used to mesh the top and bottom layers. CIRCU94 element is added to act as a resistor part of the model. The proposed model is validated with previous works. The numerical parametric studies are presented to illustrate the effects of perforation geometry, the number of rows, the resistance on the natural frequencies, frequency response, and power. It is found that the thickness has a positive relationship with the natural frequency. Perforations help in producing higher voltage, and the best shape is rectangular perforations, and to produce higher voltage, two rows of rectangular perforations should be applied.

• Journal of Mechanical Science and Technology
• /
• v.14 no.3
• /
• pp.306-313
• /
• 2000

The homogenization method and the density function method are common approaches to evaluate the equivalent material properties for design cells composed of matter and void. In this research, using a new topology optimization method based on the homogenized material with a penalty factor and the chessboard prevention strategy, we obtain the optimal layout of a structure for the natural frequency of a designated mode. The volume fraction of nodes of each finite element is chosen as the design variable and a total material usage constraint is imposed. In this paper, the subspace method is used to evaluate the eigenvalue and its corresponding eigenvector of the structure for the designated mode and the recursive quadratic programming algorithm, PLBA algorithm, is used to solve the topology optimization problem.

• Transactions of Materials Processing
• /
• v.29 no.6
• /
• pp.347-355
• /
• 2020

The mechanical responses of multilayered aluminum sheet fabricated by roll bonding, i.e., A1050/A3004 (65% A1050, 35% A3004 by thickness), were investigated via combined experiment and finite element (FE) analysis. The mechanical properties were measured using uniaxial tensile tests in various loading directions for the multilayered sheet. The corresponding tests for individual layers were also conducted. The testing samples were prepared by wire electro discharge machining (EDM). Stress-strain curves and Lankford coefficients of the multilayered sheet were then predicted by FE simulations. The measured mechanical properties of the individual layers were utilized as inputs for the simulation. Two yield functions, i.e., isotropic von-Mises and anisotropic non-quadratic Hill1948, were employed. Predicted results were compared with the experimental data and further discussed.

• East Asian mathematical journal
• /
• v.34 no.5
• /
• pp.661-681
• /
• 2018

In this paper, we discuss a reduced-order modeling for the Benjamin-Bona-Mahony-Burgers (BBMB) equation and its application to a distributed feedback control problem through the centroidal Voronoi tessellation (CVT). Spatial distcritization to the BBMB equation is based on the finite element method (FEM) using B-spline functions. To determine the basis elements for the approximating subspaces, we elucidate the CVT approaches to reduced-order bases with snapshots. For the purpose of comparison, a brief review of the proper orthogonal decomposition (POD) is provided and some numerical experiments implemented including full-order approximation, CVT based model, and POD based model. In the end, we apply CVT reduced-order modeling technique to a feedback control problem for the BBMB equation.

• Structural Engineering and Mechanics
• /
• v.18 no.5
• /
• pp.691-707
• /
• 2004

This paper presents the micro-mechanical modeling for predicting concrete behavior under compressive loading. The model is able to represent the heterogeneities in the microstructure up to three phases, i.e., aggregate particles, matrix and interfaces. The smeared crack concept based on non-linear fracture mechanics is implemented in order to formulate the constitutive relation for each component. The splitting tensile strength is considered as a fracture criterion for cracking in micro-level. The finite element method is employed to simulate the model based on plane stress condition by using quadratic triangular elements. The validation of the model is verified by comparing with the experimental results. The influence of tensile strength from both aggregate and matrix phases on the concrete compressive strength is demonstrated. In addition, a guideline on selecting appropriate tensile strength for each phase to obtain specified concrete compressive strength is also presented.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2003.05a
• /
• pp.427-432
• /
• 2003

An analytical method to predict the flexural behavior of composite girder is presented in which the early-age properties of concrete are specified including maturing of elastic modulus, creep and shrinkage. The time dependent constitutive relation accounting for the early-age concrete properties is derived in an incremental format by expanding the total form of stress-strain relation by the first order Taylor series with respect to the reference time. The sectional analysis calculates the axial and curvature strains based on the force and moment equilibriums. The deflection curve of the box girder approximated by the quadratic polynomial function is calculated by applying to the proper boundary conditions in the consecutive segments. Numerical applications are made for the 3-span double composite steel box girders which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The one dimensional finite element analysis results are compared with those of the three dimensional finite element analysis and the analytical method based on the sectional analysis. Close agreement is observed among the three methods.

• Transactions of Materials Processing
• /
• v.8 no.3
• /
• pp.252-261
• /
• 1999

The 3-dimensional finite element program is developed to analyze the non-isothermal forming processes of aluminum-alloy sheet metals. Bishop's method is introduced to solve the heat balance and force equilibrium equations. Also, Barlat's non-quadratic anisotropic yield function depicts the planar anisotropy of the aluminum-alloy sheet. To find an appropriate constitutive equation, four different forms are reviewed. For the verification of the reliability of the developed program, the computational try-outs of the non-isothermal cylindrical cupping processes of AL5052-H32 and Al1050-H16 are carried out. As results, the constitutive equation relating to strain and strain-rate, in which the constants are represented by the 5th-degree polynomials of temperature, is in good agreement with measurement. The computational try-outs can predict optimal forming conditions in non-isothermal forming processes.

• Transactions of Materials Processing
• /
• v.26 no.3
• /
• pp.144-149
• /
• 2017

In order to characterize the macroscopic mechanical response of ultra-thin (0.1 mm thick) ferritic stainless steel sheet at various loading paths, a crystal plasticity finite element method (CP-FEM) was introduced. The accuracy of the prediction results was validated by comparing with the experimental data. Based on the results, the forming limit diagram (FLD) was predicted using a modified Marchinicak-Kuczinski model coupled to a non-quadratic anisotropic yield function, namely, Yld2000-2d. The predicted FLD was found to be in good agreement with the experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.274-282
• /
• 2014

Time series of the dynamic response of a slender marine structure was predicted using quadratic Volterra series. The wave-structure interaction system was identified using the NARX(Nonlinear Autoregressive with Exogenous Input) technique, and the network parameters were determined through the supervised training with the prepared datasets. The dataset used for the network training was obtained by carrying out the nonlinear finite element analysis on the freely standing riser under random ocean waves of white noise. The nonlinearities involved in the analysis were both large deformation of the structure under consideration and the quadratic term of relative velocity between the water particle and structure in Morison formula. The linear and quadratic frequency response functions of the given system were extracted using the multi-tone harmonic probing method and the time series of response of the structure was predicted using the quadratic Volterra series. In order to check the applicability of the method, the response of structure under the realistic ocean wave environment with given significant wave height and modal period was predicted and compared with the nonlinear time domain simulation results. It turned out that the predicted time series of the response of structure with quadratic Volterra series successfully captures the slowly varying response with reasonably good accuracy. It is expected that the method can be used in predicting the response of the slender offshore structure exposed to the Morison type load without relying on the computationally expensive time domain analysis, especially for the screening purpose.