• Proceedings of the KSME Conference
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.66-71
• /
• 2004

Estimation of the CANDU pressure tube deflection is important since the deflection may cause significant structural failure due to hydrogen diffusion and blister. However, there is no appropriate engineering model to estimate it exactly. The purpose of this paper is to propose a new analysis method and program to resolve this issue. For development of proper analysis method, a series of finite element analyses has been carried under elastic-creep condition. In addition, for effective estimation of the creep deflection, an analysis program named PC-DAS was developed based on the proposed method. Comparison of simple case study results with corresponding reference ones showed good agreement. Therefore, the proposed method and program can be utilized as one of valuable toolkit for integrity assessment of CANDU pressure tube.

• Structural Engineering and Mechanics
• /
• 제35권6호
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• Structural Engineering and Mechanics
• /
• 제73권1호
• /
• pp.97-108
• /
• 2020

This study presents a 4 node, 11 DOF/node plate element based on higher order shear deformation theory for lamina composite plates. The theory accounts for parabolic distribution of the transverse shear strain through the thickness of the plate. Differential field equations of composite plates are obtained from energy methods using virtual work principle. Differential field equations of composite plates are obtained from energy methods using virtual work principle. These equations were transformed into the operator form and then transformed into functions with geometric and dynamic boundary conditions with the help of the Gâteaux differential method, after determining that they provide the potential condition. Boundary conditions were determined by performing variational operations. By using the mixed finite element method, plate element named HOPLT44 was developed. After coding in FORTRAN computer program, finite element matrices were transformed into system matrices and various analyzes were performed. The current results are verified with those results obtained in the previous work and the new results are presented in tables and graphs.

• Journal of Advanced Marine Engineering and Technology
• /
• 제32권2호
• /
• pp.250-255
• /
• 2008

Springs are widely utilized in machine element. To find out stiffness of frustum-shaped coil spring, the space beam theory using the finite element method is adopted in this paper. In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. To find out load vector of coil spring subjected to distributed compression. principle of virtual work is adapted. And this theory was programming using MATLAB software. To compare FEM using MATLAB software was applied MSC. Nastran software. The geometry model for MSC. Patran was produced by 3-D design modeling software. Finite element model was produced by MSC. Patran. Finite element was applied tetra (CTETRA) having 10 node. The analysis results of the MATLAB and MSC. Nastran are fairly well agreed with those of various experiments. Using MATLAB program proposed in this paper and MSC. Nastran, spring constants and stresses can be predicted by input of few factors.

• Interaction and multiscale mechanics
• /
• 제1권1호
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Earthquakes and Structures
• /
• 제16권3호
• /
• pp.359-368
• /
• 2019

This paper focuses on the study of dynamic analysis of thick plates resting on Winkler foundation. The governing equation is derived from Mindlin's theory. This study is a parametric analysis of the reflections of the thickness / span ratio, the aspect ratio and the boundary conditions on the earthquake excitations are studied. In the analysis, finite element method is used for spatial integration and the Newmark-${\beta}$ method is used for the time integration. While using finite element method, a new element is used. This element is 17-noded and it's formulation is derived from using higher order displacement shape functions. C++ program is used for the analyses. Graphs are presented to help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that the 17-noded finite element is used in the earthquake analysis of thick plates. It is shown that the changes in the aspect ratio are more effective than the changes in the aspect ratio. The center displacements of the reinforced concrete thick clamped plates for b/a=1, and t/a=0.2, and for b/a=2, and t/a=0.2, reached their absolute maximum values of 0.00244 mm at 3.48 s, and of 0.00444 mm at 3.48 s, respectively.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 한국소음진동공학회 2005년도 춘계학술대회논문집
• /
• pp.846-851
• /
• 2005

The purpose of this paper is that investigates the dynamic behavior characteristic of W.T.S(Wind turbine System) and carries out the evaluation analysis during operating W.T.S. To investigate the dynamic behavior characteristic of W.T.S, the experiments to measure vibration of the blade from the attached accelerometer on the flap and edge section of the blade that is one of the most important elements of dynamic characteristic of W.T.S are performed. Natural frequency and mode shape are calculated with commercial program (ANSYS) using the measured vibration acceleration that receives the signal with F.F.T Analyzer from the accelerometer. For validation of these experiments, the finite element analysis is performed with commercial F.E.M program (ANSYS) on the basis of the natural frequency and mode shape. The results indicate that experimental values have good agreements with the finite element analysis.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• 제15권11호
• /
• pp.1232-1240
• /
• 2005

The purpose of this paper is that investigates the dynamic behavior characteristic of W.T.S (wind turbine system) and carries out the evaluation analysis during operating W.T.S. To investigate the dynamic behavior characteristic of W.T.S, the experiments to measure vibration of the blade from the attached accelerometer on the flap and edge section of the blade that is one of the most important elements of dynamic characteristic of W.T.S are performed. Natural frequency and mode shape are calculated with commercial program ( ANSYS) using the measured vibration acceleration that receives the signal with F.F.T Analyzer from the accelerometer For validation of these experiments, the finite element analysis is performed with commercial F.E.M program (ANSYS) on the basis of the natural frequency and mode shape. The results indicate that experimental values have good agreements with the finite element analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• 제35권3호
• /
• pp.196-203
• /
• 2007

The postbuckling behavior and failure of composite cylinders subjected to external hydrostatic pressure are investigated by a finite element method and test. A nonlinear finite element program, ACOS, is used for the postbuckling progressive failure analysis of composite cylinders. A total of 5 carbon/epoxy composite cylinders were fabricated and tested to verify the finite element results. For comparison, analyses by MSC/NASTRAN and MSC/MARC are additionally conducted. Among the softwares, the finite element program, ACOS, predicts the buckling loads the best with about 11 to 26% deviation from experimental results except for one specimen. While the finite element analysis shows global buckling modes with 4 waves in hoop direction, in the experiments the local buckling appears first and results in the final failure without global buckling.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• 제24권2호
• /
• pp.526-534
• /
• 2000

Near-net-shape forming of zirconia powder was investigated under the combination of cold die and isostatic pressing and pressureless sintering. A novel combination pressing technique, i.e., die com paction under cold isostatic pressing, allowed to produce a complex shaped ceramic powder compact with the controlled dimensions and relatively uniform density distributions. The constitutive models proposed by Kim and co-workers for densification of ceramic powder under cold compaction and high temperature were implemented into a finite element program (ABAQUS). Experimental data for relative density distributions and deformations of zirconia powder compacts produced by cold combination pressing and pressureless sintering were compared with finite element results. Finite element results agreed well with experimental data.