• Structural Engineering and Mechanics
• /
• v.55 no.5
• /
• pp.1055-1086
• /
• 2015

The dynamic response of two-dimensional unbounded domain on the rigid bedrock in the time domain is numerically obtained. It is realized by the modified scaled boundary finite element method (SBFEM) in which the original scaling center is replaced by a scaling line. The formulation bases on expanding dynamic stiffness by using the continued fraction approach. The solution converges rapidly over the whole time range along with the order of the continued fraction increases. In addition, the method is suitable for large scale systems. The numerical method is employed which is a combination of the time domain SBFEM for far field and the finite element method used for near field. By using the continued fraction solution and introducing auxiliary variables, the equation of motion of unbounded domain is built. Applying the spectral shifting technique, the virtual modes of motion equation are eliminated. Standard procedure in structural dynamic is directly applicable for time domain problem. Since the coefficient matrixes of equation are banded and symmetric, the equation can be solved efficiently by using the direct time domain integration method. Numerical examples demonstrate the increased robustness, accuracy and superiority of the proposed method. The suitability of proposed method for time domain simulations of complex systems is also demonstrated.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Advances in nano research
• /
• v.12 no.5
• /
• pp.441-455
• /
• 2022

This study explores the linear and nonlinear solutions of sigmoid functionally graded material (S-FGM) nanoplate with porous effects. A size-dependent numerical solution is established using the strain gradient theory and isogeometric finite element formulation. The nonlinear nonlocal strain gradient is developed based on the Reissner-Mindlin plate theory and the Von-Karman strain assumption. The sigmoid function is utilized to modify the classical functionally graded material to ensure the constituent volume distribution. Two different patterns of porosity distribution are investigated, viz. pattern A and pattern B, in which the porosities are symmetric and asymmetric varied across the plate's thickness, respectively. The nonlinear finite element governing equations are established for bending analysis of S-FGM nanoplates, and the Newton-Raphson iteration technique is derived from the nonlinear responses. The isogeometric finite element method is the most suitable numerical method because it can satisfy a higher-order derivative requirement of the nonlocal strain gradient theory. Several numerical results are presented to investigate the influences of porosity distributions, power indexes, aspect ratios, nonlocal and strain gradient parameters on the porous S-FGM nanoplate's linear and nonlinear bending responses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.3
• /
• pp.127-132
• /
• 2018

In this paper, a scheme for the evaluation of variability in the eigen-modes of functionally graded material(FGM) beams is proposed within the framework of perturbation-based stochastic analysis. As a random parameter, the spatially varying elastic modulus of FGM along the axial direction at the mid-surface of the beam is chosen, and the thru-thickness variation of the elastic modulus is assumed to follow the original form of exponential variation. In deriving the formulation, the first order Taylor expansion on the eigen-modes is employed. As an example, a simply supported FGM beam having symmetric elastic modulus with respect to the mid-surface is chosen. Monte Carlo analysis is also performed to check if the proposed scheme gives reasonable outcomes. From the analyses it is found that the two schemes give almost identical results of the mean and standard deviation of eigen-modes. With the propose scheme, the standard deviation shape of respective eigen-modes can be evaluated easily. The deviated mode shape is found to have one more zero-slope points than the mother modes shapes, irrespective of order of modes. The amount of deviation from the mean is found to have larger values for the higher modes than the lower modes.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.8
• /
• pp.1195-1208
• /
• 1997

Bellows is a component in piping systems which absorbs mechanical deformation with flexibility. Its geometry is an axial symmetric shell which consists of two toroidal shells and one annular plate or conical shell. In order to analyze bellows, this study presents the finite element analysis using a conical frustum shell element. A finite element analysis is developed to analyze various bellows. The validity of the developed program is verified by the experimental results for axial and lateral stiffness. The formula for calculating the natural frequency of bellows is made by the simple beam theory. The formula for fatigue life is also derived by experiments. The shape optimal design problem is formulated using multiple objective optimization. The multiple objective functions are transformed to a scalar function by weighting factors. The stiffness, strength and specified stiffness are considered as the multiple objective function. The formulation has inequality constraints imposed on the fatigue limit, the natural frequencies, and the manufacturing conditions. Geometric parameters of bellows are the design variables. The recursive quadratic programming algorithm is selected to solve the problem. The results are compared to existing bellows, and the characteristics of bellows is investigated through optimal design process. The optimized shape of bellows is expected to give quite a good guideline to practical design.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.

• Journal of the Korea Computer Graphics Society
• /
• v.4 no.1
• /
• pp.79-88
• /
• 1998

Till now there are several rendering models for water and simulating other fluids and their dynamics. Especially in order to generate a curved surface of flexible objects such as liquid and snow, the implicit metaball formulation is widely used in favor of its simplicity and flexibility. This paper proposes one excellent method for generating water droplets, which would be deformed in gravitation field. In previous works, a water droplet was simply represented by approximated curved surfaces of a symmetric metaball. Thus the final result of the rendered water droplet was far from a realistic droplet, because they do not consider the gravitational effect in droplets. We propose a new metaball model for rendering water droplets placed on an arbitrary surface considering the gravitation and friction between droplet and plate. Our new metaball model uses a new vector field isosurface function to control the basic scalar metaball with respect to the norm of gravitational force. In several experiments, we could render a photo-realistic water droplets with natural-looking shadows by applying ray-tracing.

• Steel and Composite Structures
• /
• v.28 no.3
• /
• pp.309-318
• /
• 2018

Buckling Restrained Braced (BRB) frames have been widely used as an efficient seismic load resisting system in recent years mostly due to their symmetric and stable hysteretic behavior and significant energy dissipation capacity. In this study, to provide a better understanding of the behavior of BRB frames with various beam-column connections, a numerical study using non-linear finite element (FE) analysis is conducted. All models are implemented in the Abaqus software package following an explicit formulation. Initially, the results of the FE model are verified with experimental data. Then, diverse beam-column connections are modeled for the sake of comparison from the shear capacity, energy dissipation and frame hysteresis behavior points of view until appropriate performance is assessed. The considered connections are divided into three different categories: (1) simple beam-column connections including connection by web angle and connection by seat angle; (2) semi-rigid connection including connection by web and seat angles; and (3) rigid beam-column connections by upper-lower beam plates and beam connections with web and flange splices. Results of the non-linear FE analyses show that these types of beam-column connections have little effect on the maximum story drift and shear capacity of BRB frames. However, the connection type has a significant effect on the amount of energy dissipation and hysteresis behavior of BRB frames. Also, changes in length and thickness of the angles in simple and semi-rigid connections and changes in length and thickness of plates in rigid connections have slight effects (less than 4%) on the overall frame behavior.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.4
• /
• pp.15-26
• /
• 1990

This paper describes a potential-baoed panel method formulated for the analysis cf a supercavitating two-dimensional symmetri strut. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring that the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type, With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lifting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged.

• Structural Engineering and Mechanics
• /
• v.85 no.3
• /
• pp.289-304
• /
• 2023

The main objective of this paper is to study the effect of porosity on the buckling behavior of thick functionally graded sandwich plate resting on various boundary conditions under different in-plane loads. The formulation is made for a newly developed sandwich plate using a functional gradient material based on a modified power law function of symmetric and asymmetric configuration. Four different porosity distribution are considered and varied in accordance with material propriety variation in the thickness direction of the face sheets of sandwich plate, metal foam also is considered in this study on the second model of sandwich which containing metal foam core and FGM face sheets. New quasi-3D high shear deformation theory is used here for this investigate; the present kinematic model introduces only six variables with stretching effect by adopting a new indeterminate integral variable in the displacement field. The stability equations are obtained by Hamilton's principle then solved by generalized solution. The effect of Pasternak and Winkler elastic foundations also including here. the present model validated with those found in the open literature, then the impact of different parameters: porosities index, foam cells distribution, boundary conditions, elastic foundation, power law index, ratio aspect, side-to-thickness ratio and different in-plane axial loads on the variation of the buckling behavior are demonstrated.