• 한국전산구조공학회논문집
• /
• 제14권4호
• /
• pp.435-444
• /
• 2001

대변형 탄소성 프레임의 해석은 운동학적, 재질적, 수치해법상으로 매우 복잡하여 정확하고 효율적인 요소나 알고리즘이 부족하다. 본 논문에서는 3차원 보요소에 대하여 봉이론으로부터 유한요소를 만들고 소성을 적용하기 위하여 Cauchy응력과 공학적변형률을 구성방정식에 적용하는 방법으로 객관화된 증분해법에 의한 보요소를 제안하였다. 특히 항복조건의 만족을 위하여 정확하고도 효율적인 소성방정식의 적분법을 개발하고 적용하였으며 연속법과 일차원 탐색 등을 고려하여 광역수렴성과 2차 수렴속도를 갖는 수치해법을 개발하고 예제를 통하여 확인하였다.

• Computers and Concrete
• /
• 제4권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.

• 대한기계학회논문집
• /
• 제14권5호
• /
• pp.1166-1180
• /
• 1990

본 연구에서는 열간 반경단조를 해석하는데 있어서 재료의 변형 해석은 강-점 소성해석법으로 해석하고 온도분포도 유한 요소법으로 해석한다.

• Structural Engineering and Mechanics
• /
• 제54권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.

• 콘크리트학회논문집
• /
• 제19권6호
• /
• pp.803-810
• /
• 2007

FRP 합성재료로 구속된 콘크리트의 응력-변형률 응답을 합리적으로 예측할 수 있는 해석 모델이 제시되었다. 제안된 모델은 하중이 증가함에 따라 점진적으로 발생하는 미세균열에 의한 부피팽창이 미세 재료 구조의 손상을 나타내는 중요한 척도이며, 이에 손상 정도에 따라 하중 지지 능력을 일관되게 산정할 수 있다는 기본 개념에 근거한다. 이를 위하여 제안 모델은 면적 변형률 및 공극의 함수로 표시된 탄성계수, 팽창 콘크리트와 구속 매체의 상호작용을 나타내는 에너지 평형식, 변화하는 구속력 및 점증 계산 논리를 포함한다. 따라서 실험으로부터 유도된 팽창비 관계식으로부터 횡방향 혹은 부피팽창변형률을 산정하는 기존의 해석 모델과는 달리 역학적 거동 및 에너지 평형식으로부터 연속적으로 변화하는 횡방향 변형률을 산정한다. 구속된 콘크리트의 전체 응답을 예측할 수 있는 기존의 여러 해석 모델에 대하여 검토하였으며, 특히 부피 팽창을 고려하는 방법에 초점을 맞추어 토의하였다. 제안된 모델 및 기존 Samaan의 2선식 모델을 사용하여 실험 결과를 예측한 결과, 만족할 만한 범위 내에서 일치를 나타냈으나, Samaan의 2선식 모델은 부피 팽창 거동을 위하여 단지 초기 포아송비와 최종 수렴 팽창변화율 만을 고려하기 때문에 횡방향 변형률 응답을 예측하는 데는 한계가 있는 것으로 판단된다. 제안된 모델은 역학적 거동에 근거하여 다양한 관련 응답을 산정하므로 다른 분야에도 쉽게 적용할 수 있다.

• 소성∙가공
• /
• 제15권6호
• /
• 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
• /
• 제2권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
• /
• 제40권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
• /
• 제89권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.

• 대한조선학회논문집
• /
• 제34권1호
• /
• pp.68-76
• /
• 1997

본 논문에서는 선박 및 해양 구조물을 구성하는 용접 판구조물의 변형 및 잔류응력 등 용접 초기 결함을 효율적으로 평가할 수 있는 열탄소성 해석방법을 정립하였다. 또한 평판의 선상가열공법(line heating process)을 열탄소성 시뮬레이션함으로써 작업의 고능률화 및 자동화를 위한 기초연구를 수행하였다. 판구조물에 대한 열전도 해석을 수행하기 위하여 해석적인 방법과 수치적인 방법을 병행하였다. 판구조물의 열탄소성 해석방법에 있어서는 초기변형도법에 근거한 유한요소법을 사용하였으며, 증분 구간중 소성 구역에서는 응력 증분 및 항복응력 중분의 2 차항까지 고려해서 시간 증분을 최적 제어함으로써 해석 불능에 빠지는 문제를 극복하였다. 특히 응력 중분에 탄성계수의 온도에 따른 증분을 2 차항까지 포함시켰다. 평판의 두께와 입열량을 변화시키면서 일련의 용접 열탄소성 실험을 수행하였는 바, 온도 및 용접 변형을 계측하여 수치 계산 결과와 비교하여 상호 부합성이 양호함을 확인하였다.