• Steel and Composite Structures
• /
• v.42 no.1
• /
• pp.23-32
• /
• 2022

In this paper, a nonlinear numerical method to solve the large deflection problem is introduced. And the non-dimensional load-deflection behavior of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) plates is parametrically investigated. The large deflection problem is formulated according to the von Kármán nonlinear theory and the (1,1,0)^* hierarchical model, and it is approximated by 2-D natural element method (NEM). The shear locking phenomenon is suppressed by the selectively reduced integration method. The nonlinear matrix equations are solved by combining the incremental loading scheme and the Newton-Raphson iteration method. The proposed method is validated from the benchmark experiments, where the propose method shows an excellent agreement with the reference methods. The nonlinear behavior of FG-CNTRC plates is evaluated in terms of the non-dimensional load-deflection curve, and it is parametrically investigated with respect to the existence/non-existence and gradient pattern of CNTs, the width-to-thickness and aspect ratios of plates and the type of boundary conditions. The non-dimensional central deflection is significantly reduced when CNTs and added, and it decreases with the volume fraction of CNTs. But, it shows a uniform increase in proportion to the width-to-thickness and aspect ratios. Both the gradient pattern of CNTs and the type of boundary conditions do also show the remarkable effects.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Journal of the Korean Recycled Construction Resources Institute
• /
• v.6 no.2
• /
• pp.100-107
• /
• 2018

The material and geometrical nonlinear finite elment analysis of UHPFRC 50M composite box girder was carried out. Constitute law in tension and compressive region of UHPFRC and HPC were modeled based on specimen test. The accuracy of nonlinear FEM analysis was verified by the experimental result of UHPFRC 50M composite girder. The UHPFRC 50M segmental composite box girder which has 1.5% steel fiber of volume fraction, 135MPa compressive strength and 18MPa tensile strength was tested. The post-tensioned UHPFRC composite girder consisted of three segment UHPFRC U-girder and High Strength Concrete reinforced slab. The parts of UHPFRC girder were modeled by 8nodes hexahedron elements and reinforcement bars and tendons were built by 2nodes linear elements by Midas FEA software. The constitutive laws of concrete materials were selected Multi-linear model both of tension and compression function under total strain crack model, which was included in classifying of smeared crack model. The nonlinearity of reinforcement elements and tendon was simulated by Von Mises criteria. The nonlinear static analysis was applied by incremental-iteration method with convergence criteria of Newton-Raphson. The validation of numerical analysis was verified by comparison with experimental result and numerical analysis result of load-deflection response, neutral axis coordinate change, and cracking pattern of girder. The load-deflection response was fitted very well with comparison to the experimental result. The finite element analysis is seen to satisfactorily predict flexural behavioral responses of post-tensioned, reinforced UHPFRC composite box girder.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.15 no.6
• /
• pp.3448-3454
• /
• 2014

Generally existing structure analysis was applied to elastic analysis basically in practice. Considering the nonlinear material and the nonlinear geometric to be a more precise analysis, for this reason, The necessity for a structual analysis have been constantly required. Therefore, after optimization is performed, designed a simple model which is applied the principle of nonlinear in this study, a structural analysis of existing experienced users, have a aims at presenting theory and a method in order to perform anyone the analysis easily. In this study, the proposed model applied to die ribs, Regarding the shear load, less strain and stress was generated but strength was sufficient. The initial strain and stress was reconfigured to fit the size and shape, A hyperstudy in conjunction with Abaqus with nonlinear structural analysis, revealed an acceptable maximum and minimum range of stress and under the conditions of minimum strain, the plate made with a constant increment. In the experimental models, the plate thickness was given a power of 40 Newton, according to the thickness of the press die through an iterative process. When the stress and strain was applied to the die thickness, 7-8mm thickness could be obtained by optimizing.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.28-37
• /
• 2002

In this study nonlinear dynamic behaviors of towed tow-tension cables are numerically analysed. In the case of a taut cable analysis, a bending stiffness term is usually neglected due to its minor effect but it plays an important role in a low-tension cable analysis. A low-tension cable may experience large displacements due to relatively small restoring forces and thus the effects of fluid and geometric non-linearities become predominant. The bending stiffness and non-linearity effects are considered in this work. In order to obtain dynamic behaviors of a towed low-tension cable, three-dimensional nonlinear dynamic equation is described and discretized by employing a finite difference method. An implicit method and Newton-Raphson iteration are adopted for the time integration and nonlinear solutions. For the calculation of huge size of matrices. block tri-diagonal matrix method is applied, which is much faster than the well-known Gauss-Jordan method in two point boundary value problems. Some case studies are carried out and the results of numerical simulations are compared with those of a in-house program of WHOI Cable with good agreements.

• Journal of Korea Water Resources Association
• /
• v.41 no.9
• /
• pp.895-905
• /
• 2008

To analyze hydrologic processes in a watershed requires both various geographical data and hydrological time series data. Recently, not only geographical data such as DEM(Digital Elevation Model) and hydrologic thematic map but also hydrological time series from numerical weather prediction and rainfall radar have been provided as grid data, and there are studies on hydrologic analysis using these grid data. In this study, GRM(Grid based Rainfall-runoff Model) which is physically-based distributed rainfall-runoff model has been developed to simulate short term rainfall-runoff process effectively using these grid data. Kinematic wave equation is used to simulate overland flow and channel flow, and Green-Ampt model is used to simulate infiltration process. Governing equation is discretized by finite volume method. TDMA(TriDiagonal Matrix Algorithm) is applied to solve systems of linear equations, and Newton-Raphson iteration method is applied to solve non-linear term. Developed model was applied to simplified hypothetical watersheds to examine model reasonability with the results from $Vflo^{TM}$. It was applied to Wicheon watershed for verification, and the applicability to real site was examined, and simulation results showed good agreement with measured hydrographs.

• Journal of the Korean earth science society
• /
• v.37 no.5
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• Journal of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.59-66
• /
• 2005

A new ion transport code for planetary ionospheric studies has been developed with consideration of velocity differences among ion species involving ion-ion collision. Most of previous planetary ionosphere models assumed that ions diffuse through non-moving ion and neutral background in order to consolidate continuity and momentum equations for ions into a simple set of diffusion equations. The simplification may result in unreliable density profiles of ions at high altitudes where ion velocities are fast and their velocity differences are significant enough to cause inaccuracy when computing ion-ion collision. A new code solves explicitly one-dimensional continuity and momentum equations for ion densities and velocities by utilizing divided Jacobian matrices in matrix inversion necessary to the Newton iteration procedure. The code has been applied to Martian nightside ionosphere models, as an example computation. The computed density profiles of $O^+,\;OH^+$, and $HCO^+$ differ by more than a factor of 2 at altitudes higher than 200 km from a simple diffusion model, whereas the density profile of the dominant ion, $O_2^+$, changes little. Especially, the density profile of $HCO^+$ is reduced by a factor of about 10 and its peak altitude is lowered by about 40 km relative to a simple diffusion model in which $HCO^+$ ions are assumed to diffuse through non-moving ion background, $O_2^+$. The computed effects of the new code on the Martian nightside models are explained readily in terms of ion velocities that were solved together with ion densities, which were not available from diffusion models. The new code should thus be expected as a significantly improved tool for planetary ionosphere modelling.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.2 s.173
• /
• pp.329-337
• /
• 2000

An investigation has been performed to develop an analysis tool based on a nonlinear beam theory, which can be used to predict the long-term behavior of an artificial hip joint. The nonlinear behav ior of the femur arise from the coupled dependence of the bone density and the mechanical properties on each other. The beam theory together with its numerical algorithm is developed to take into account the nonlinear bone remodeling process of the femur that is long enough to be assumed as a beam. A piecewise linear curve for the bone remodeling rate is used in the bone remodeling theory and the surface area density of bone is modeled as the third order polynomial function of bone density. At each section of the beam, a constant curvature is assumed and the longitudinal strains are also assumed to vary linearly across the section. The Newton-Rhapson iteration method is used to solve the nonlinear equations for each cross section of the bone and a backward method is used to march along the time. The density and the remodeling signal ar, calculated along with time for the various time steps, and the developed beam theory has been verified by comparing with the results of finite element analysis of a remodeling bone with an artificial hip joint of titanium prosthesis subjected to uni-axial loads and pure bending moment. It is concluded that the developed beam theory can be used to predict the long-term behavior of the femur and thus to design the artificial hip prosthesis.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.185-191
• /
• 2004

A mooring system can be applied to keep the position of a floating structures. In this study, the structural analysis is carried out to analyze the dynamic characteristics of a mooring line for a floating breakwater. A three-dimensional equations of motion for a submerged chain are derived. Bending stiffness is considered for the necessary restoring force in the regions of zero tension. A fortran program is to be developed by employing finite difference method. In the algorithm, an implicit time integration and Newton-Raphson iteration are adopted. The results of simulation show good agreement in tension response pattern with the experimental results of a reference. The results of this study can contribute for the design of mooring system for a floating breakwater.