• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.13 no.3
• /
• pp.237-245
• /
• 2008

We developed a computational method on coupled dynamics of tracked vehicle on seafloor and long flexible pipe. The tracked vehicle is modeled as rigid-body vehicle, and the linked flexible pipe is discretized according to a lumped-parameter model. The equations of motion of the rigid-body vehicle on the soft seafloor are combined with the governing equations of flexible pipe dynamics. Four Euler parameters method is used to express the orientations of the vehicle and the flexible pipe. In order to solve the nonlinear coupled dynamics of vehicle and flexible pipe an incremental-iterative formulation is implemented. For the time-domain integration $Newmark-\beta$ method is adopted. The total Jacobean matrix has been derived based on the incremental-iterative formulation. The interactions between the dynamics of flexible pipe and the mobility of the tracked vehicle on soft seafloor are investigated through numerical simulations in time domain.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.4 s.10
• /
• pp.85-92
• /
• 2003

In this work, a finite element model is presented for geometrically non-linear analysis of shell structures. Finite element by using a three-node flat triangular shell element is formulated. The non-linear incremental equilibrium equations are formulated by using an updated Lagrangian formulation and the solutions are obtained with the incremental/iterative Newton-Raphson method and arc length method. Some of results are presented for shell structures. The obtained results are in good agreement with the results available in existing literature.

• Proceedings of the KIEE Conference
• /
• 1999.07c
• /
• pp.1506-1509
• /
• 1999

This paper presents a new economic dispatch algorithm to improve the unit commitment solution while guaranteeing the near optimal solution without reducing calculation speed. The conventional economic dispatch algorithms have the problem that it is not applicable to the unit commitment formulation due to the frequent on/off state changes of units during the unit commitment calculation. Therefore, piecewise linear iterative method have generally been used for economic dispatch algorithm for unit commitment. In that method, the approximation of the generator cost function makes it hard to obtain the optimal economic dispatch solution. In this case, the solution can be improved by introducing a inverse of the incremental cost function. The proposed method is tested with sample system. The results are compared with the conventional piecewise linear iterative method. It is shown that the proposed algorithm yields more accurate and economical solution without calculation speed reduction.

• Structural Engineering and Mechanics
• /
• v.12 no.1
• /
• pp.35-50
• /
• 2001

A numerical methodology is presented in this paper for the geometrically non-linear analysis of slender uni-dimensional structural elements under unilateral contact constraints. The finite element method together with an updated Lagrangian formulation is used to study the structural system. The unilateral constraints are imposed by tensionless supports or foundations. At each load step, in order to obtain the contact regions, the equilibrium equations are linearized and the contact problem is treated directly as a minimisation problem with inequality constraints, resulting in a linear complementarity problem (LCP). After the resulting LCP is solved by Lemke's pivoting algorithm, the contact regions are identified and the Newton-Raphson method is used together with path following methods to obtain the new contact forces and equilibrium configurations. The proposed methodology is illustrated by two examples and the results are compared with numerical and experimental results found in literature.

• Structural Engineering and Mechanics
• /
• v.26 no.6
• /
• pp.725-739
• /
• 2007

In recent years there are many plate bending elements that emerged for solving both thin and thick plates. The main features of these elements are that they are based on mix formulation interpolation with discrete collocation constraints. These elements passed the patch test for mix formulation and performed well for linear analysis of thin and thick plates. In this paper a member of this family of elements, namely, the Discrete Reissner-Mindlin (DRM) is further extended and developed to analyze both thin and thick plates with geometric nonlinearity. The Von K$\acute{a}$rm$\acute{a}$n's large displacement plate theory based on Lagrangian coordinate system is used. The Hu-Washizu variational principle is employed to formulate the stiffness matrix of the geometrically Nonlinear Discrete Reissner-Mindlin (NDRM). An iterative-incremental procedure is implemented to solve the nonlinear equations. The element is then tested for plates with simply supported and clamped edges under uniformly distributed transverse loads. The results obtained using the geometrically NDRM element is then compared with the results of available analytical solutions. It has been observed that the NDRM results agreed well with the analytical solutions results. Therefore, it is concluded that the NDRM element is both reliable and efficient in analyzing thin and thick plates with geometric non-linearity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.3
• /
• pp.237-244
• /
• 2016

This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

• Structural Engineering and Mechanics
• /
• v.49 no.6
• /
• pp.727-743
• /
• 2014

Geometrically nonlinear analysis of planar beam and frame structures made of functionally graded material (FGM) by using the finite element method is presented. The material property of the structures is assumed to be graded in the thickness direction by a power law distribution. A nonlinear beam element based on Bernoulli beam theory, taking the shift of the neutral axis position into account, is formulated in the context of the co-rotational formulation. The nonlinear equilibrium equations are solved by using the incremental/iterative procedure in a combination with the arc-length control method. Numerical examples show that the formulated element is capable to give accurate results by using just several elements. The influence of the material inhomogeneity in the geometrically nonlinear behavior of the FGM beam and frame structures is examined and highlighted.

• Structural Engineering and Mechanics
• /
• v.35 no.3
• /
• pp.349-364
• /
• 2010

The impact behaviour and the impact-induced damage in laminated composite cylindrical shell subjected to transverse impact by a foreign object are studied using three-dimensional non-linear transient dynamic finite element formulation. A layered version of 20 noded hexahedral element incorporating geometrical non-linearity is developed based on total Langragian approach. Non-linear system of equations resulting from non-linear strain displacement relation and non-linear contact loading are solved using Newton-Raphson incremental-iterative method. Some example problems of graphite/epoxy cylindrical shell panels are considered with variation of impactor and laminate parameters and influence of geometrical non-linear effect on the impact response and the resulting damage is investigated.

• Computers and Concrete
• /
• v.4 no.3
• /
• pp.187-206
• /
• 2007

A rational three-dimensional nonlinear finite element model is described and implemented for evaluating the behavior of high strength concrete slabs under transverse load. The concrete was idealized by using twenty-nodded isoparametric brick elements with embedded reinforcements. The concrete material modeling allows for normal (NSC) and high strength concrete (HSC), which was calibrated based on experimental data. The behavior of concrete in compression is simulated by an elastoplastic work-hardening model, and in tension a suitable post-cracking model based on tension stiffening and shear retention models are employed. The nonlinear equations have been solved using the incremental iterative technique based on the modified Newton-Raphson method. The FE formulation and material modeling is implemented into a finite element code in order to carry out the numerical study and to predict the behavior up to ultimate conditions of various slabs under transverse loads. The validity of the theoretical formulations and the program used was verified through comparison with available experimental data, and the agreement has proven to be very good. A parametric study has been also carried out to investigate the influence of different material and geometric properties on the behavior of HSC slabs. Influencing factors, such as concrete strength, steel ratio, aspect ratio, and support conditions on the load-deflection characteristics, concrete and steel stresses and strains were investigated.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.