• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.435-444
• /
• 2001

Even todays, accurate and efficient algorithms for the large deformation analysis of elastoplastic frame structures lack due to the complexities of kinematics, material nonlinearities and numerical methods to cater for. The author suggests appropriate beam element based upon the incremental formulation from the 3D rod theory where Cauchy stress and engineering strain are variables to incorporate plasticity equations so that objectivity may be satisfied. A rectum mapping methods which can integrate and satisfy yield criteria efficiently is suggested and a continuation method which has global convergency and quadratic speed is developed as well. leading-unloading example problems are tested and the ideas are proved to be valuable.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.5
• /
• pp.1166-1180
• /
• 1990

The study is concerned with the two-dimensional thermo-viscoplastic finite element analysis for radial forging as an incremental forging process. The deformation and temperature distribution of the workpiece during radial forging are studied. The analysis of deformation and the analysis of heat transfer are carried out for simple upsetting of cylinder by decoupling the above two analyses. A method of treatment for heat transfer through the contact region between the die and the workpiece is suggested, in which remeshing of the die elements is not necessary. Radial forging of a mild steel cylinder at the elevated temperature is subjected to the decoupled finite element analysis as well as to the experiment. The computed results in deformation, load and temperature distribution are found to be in good agreement with the experimental observations. As an example of viscoplastic decoupled analysis of hot radial forging, forging of a square section into a circular section is treated. The stresses, strains, strain rates and temperature distribution are computed by superposing material properties as the workpiece is rotated and forged incrementally. It was been thus shown that proposed method of analysis can be effectively applied to the hot radial forging processes.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Journal of the Korea Concrete Institute
• /
• v.19 no.6
• /
• pp.803-810
• /
• 2007

An analytical method capable of predicting various stress-strain responses in axially loaded concrete confined with FRP (fiber reinforced polymers) composites in a rational manner is presented. Its underlying idea is that the volumetric expansion due to progressive microcracking in mechanically loaded concrete is an important measure of the extent of damage in the material microstructure, and can be utilized to estimate the load-carrying capacity of concrete by considering the corresponding accumulated damage. Following from this, an elastic modulus expressed as a function of area strain and concrete porosity, the energy-balance equation relating the dilating concrete to the confining device interactively, the varying confining pressure, and an incremental calculation algorithm are included in the solution procedure. The proposed method enables the evaluation of lateral strains consecutively according to the related mechanical model and the energy-balance equation, rather than using an empirically derived equation for Poisson's ratio or dilation rate as in other analytical methods. Several existing analytical methods that can predict the overall response were also examined and discussed, particularly focusing on the way of considering the volumetric expansion. The results predicted by the proposed and Samaan's bilinear equation models correlated with observed results with a reasonable degree, however it can be judged that the latter is not capable of predicting the response of lateral strains correctly due to incorporating the initial Poisson's ratio and the final converged dilation rate only. Further, the proposed method seems to have greater benefits in other applications by the use of the fundamental principles of mechanics.

• Transactions of Materials Processing
• /
• v.15 no.6 s.87
• /
• pp.445-452
• /
• 2006

In most of automobile body panel manufacturing, trimming process is generally performed before flanging. To find feasible trimming line is crucial in obtaining accurate edge profile after flanging. Section-based method develops blank along manually chosen section planes and find trimming line by generating loop of end points. This method suffers from inaccurate results of edge profile. On the other hand, simulation-based method can produce more accurate trimming line by iterative strategy. In this study, new fast simulation-based method to find feasible trimming line is proposed. Finite element inverse method is used to analyze the flanging process because final shape after flanging can be explicitly defined and most of strain paths are simple in flanging. In utilizing finite element inverse method, the main obstacle is the initial guess generation for general mesh. Robust initial guess generation method is developed to handle genera] mesh with very different size and undercut. The new method develops final triangular mesh incrementally onto the drawing tool surface. Also in order to remedy mesh distortion during development, energy minimization technique is utilized. Trimming line is extracted from the outer boundary after finite element inverse method simulation. This method has many advantages since trimming line can be obtained in the early design stage. The developed method is verified by shrink/stretch flange forming and successfully applied to the complex industrial applications such as door outer flanging process.

• 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.

• 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.

• Structural Engineering and Mechanics
• /
• v.89 no.1
• /
• pp.23-31
• /
• 2024

Ratcheting is the cyclic buildup of inelastic strain on a structure resulting from a combination of primary and secondary cyclic stress. It can lead to excessive plastic deformation, incremental collapse, or fatigue. Ratcheting has been numerically investigated on a cantilever beam, considering the current study's primary and secondary bending loads. In addition, the effect of input frequency on the onset of ratcheting has been investigated. The non-linear dynamic elastic-plastic approach has been utilized. Analogous to Yamashita's bending-bending ratchet diagram, a non-dimensional ratchet diagram with a frequency effect is proposed. The result presents that the secondary stress values fall sequentially with the increase of primary stress values. Moreover, a displacement amplification factor graph is also established to explain the effect of frequency on ratchet occurrence conditions. In terms of frequency effect, it has been observed that the lower frequency (0.25 times the natural frequency) was more detrimental for ratchet occurrence conditions than the higher frequency (2 times the natural frequency) due to the effect of dynamic displacement. Finally, the effect of material modeling of ratcheting behavior on a beam is shown using different hardening coefficients of kinematic hardening material modeling.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.68-76
• /
• 1997

The welding-induced initial imperfections such as residual stresses and initial strains in plated structures of ships and offshore structures can be effectively evaluated by the thermal elasto-plastic analysis method proposed in this paper. In the analysis of heat conduction of plate structures, both the analytical method and the numerical method are used. For the thermal elasto-plastic analysis of plates, the finite element analysis is performed, based on the initial strain method. In the plastic domain during incremental process, the 2nd order terms of stress increments and yield stress increments were considered, so that time increment could be controlled for more stable solution. To measure temperature distribution and angular distortion of plates during welding, bead-on-plate experiment are perform with various heat input and plate thickness. Measured data show good agreement with the calculated results.