• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.173-179
• /
• 2000

The influence of the plastic flow rules on the elasto-plastic behaviour of a discrete joint element was investigated by performing the numerical direct shear tests under both constant normal displacement and normal displacement conditions. The finite interface elements obeying Plesha’s joint constitutive law was used to allow the relative motion of the rock blocks on the joint surface. Realistic results were obtained in the tests adopting the non-associated flow rule, while the associated flow rule overestimated the joint dilation. To overcome the computational drawbacks coming from the non-symmetric element stiffness matrix in the conventional non-associated plasticity, the symmetric formulation of the tangential stiffness matrix for a non-associated joint element was proposed. The symmetric elasto-plastic matrix it derived by assuming an imaginary equivalent joint with associated flow rule which shows the same plastic response as that of original Joint with non-associated flow rule. The validity of the formulation was confirmed through the numerical direct shear tests under constant normal stress condition.

• Composites Research
• /
• v.30 no.4
• /
• pp.247-253
• /
• 2017

In this study, the mechanical behavior and interface properties of boron nitride nanotube-poly(methyl methacrylate) nanocomposites are predicted using the molecular dynamics simulations and the double inclusion model. After modeling nanocomposite unit cell embedding single-walled nanotube and polymer, the stiffness matrix is determined from uniaxial tension and shear tests. Through the orientation average of the transversely isotropic stiffness matrix, the effective isotropic elastic constants of randomly dispersed microstructure of nanocomposites. Compared with the double inclusion model solution with a perfect interfacial condition, it is found that the interface between boron nitride nanotube and polymer matrix is weak in nature. To characterize the interphase surrounding the nanotube, the two step domain decomposition method incorporating a linear spring model at the interface is adopted. As a result, various combinations of the interfacial compliance and the interphase elastic constants are successfully determined from an inverse analysis.

• Composites Research
• /
• v.32 no.5
• /
• pp.258-264
• /
• 2019

In this paper, a progressive failure analysis method was developed using the Hashin failure criterion and crack band model. Using the failure criterion, the failure initiation was evaluated. If the failure initiation is occurred, the damage variables at each failure modes (fiber tension & compression, matrix tension & compression) was calculated according to linear softening degradation behavior and the variables are used to derive the damaged stiffness matrix. The damaged stiffness matrix is reflected to damaged material and the progressive failure analysis is continued until the damage variables to be 1 that complete failure of material. A series of processes were performed using FE commercial code ABAQUS with user defined material subroutine (UMAT). To evaluate the proposed progressive failure model, the experimental results of open hole composite laminate tests was compared with numerical result. Using digital image correlation system, the strain behavior also was compared. The proposed numerical results were coincided well with the experimental results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.6
• /
• pp.92-99
• /
• 2005

Numerical analyses of critical speed and mass unbalance response are performed for a 30 ton thrust turbopump. The stiffness and damping of ball bearings and non-contact seals are quantified under aerodynamic and hydrodynamic loads induced by a fuel pump and turbine. Critical speed margin and tip displacements of the rotating parts are evaluated using a three-dimensional finite element method. The results are used to ensure the soundness of the rotordynamic design using an one-dimensional transfer matrix method. A further study shows that sufficient resonance margin may be assured via controlling the stiffness of the rotor support by employing an additional elastic ring to the bearing support.

• Steel and Composite Structures
• /
• v.18 no.6
• /
• pp.1405-1421
• /
• 2015

Semi-rigid connections are the actual behavior of beam-to-column connections in steel frames. However, the behavior of semi-rigid connections is not taken into account for the simplicity in the conventional analysis and design of steel frames. A computer-based analysis and design has been studied for the three-dimensional steel frames with semi-rigid connections. The nonlinear analysis which includes the effects of the flexibility of connections is used for this study. It is designed according to the buckling and combined stress constraints under the present loading after the joint deformations and the member end forces of the space frame are determined by the stiffness matrix method. The semi-rigid connection type is limited to the top and bottom angles with a double web angle connection. The Frye-Morris polynomial model is used to describe the non-linear behavior of semi-rigid connections. Various design examples are presented to demonstrate the efficiency of the method. The results of design and analysis of unbraced semi-rigid frames are compared to the results of unbraced rigid frames under the same design requirements.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.2 no.1
• /
• pp.23-34
• /
• 1998

A dynamic fluid-structure-soil interaction analysis method is developed to investigate the effects of translational and/or rocking motions on the seismic response of flexible rectangular liquid storage tanks founded on the deformable ground. The governing equation for the dynamics of 3-D rectangular tanks subjected to the translational and/or rocking motion is abtained by applying Rayleigh-Ritz method. The dynamic stiffness matrices of a rigid rectangular foundation resting on the surface of a stratum overlaid bedrock are calculated by hyperelement method. The seismic responses of 3-D flexible tank model founded on the deformable ground is calculated by combining the governing equation for the fluid-tank system with the dynamic stiffness matrix of th rigid surface foundation.

• Journal of Advanced Marine Engineering and Technology
• /
• v.21 no.3
• /
• pp.246-255
• /
• 1997

The wave washer springs are expected to behave non-linearly between forces and displace¬ments due to contractions of the height and due to expansions in radial direction. To find out the non -linearity of wave washer springs, the three dimensional plate analysis theory using the finite element method is adopted in this paper. The wave washer springs are considered to be three dimensional plate structures rather than frame structures, because their thickness is normally much smaller than their width. The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displacements are added to X - Y Z coordinates of nodal points. The new stiffness matrix of the system using the new coordinates of nodal points is adopted to calculate the another nodal displacements, that is, the step by step method is used in this paper. The relations between the increments of forces and displacements in each step are recorded and plotted in chart. The experimental results are compared with the calculated chart and it is shown that there are good coincidences between measured values and calculated ones.

• Journal of Advanced Marine Engineering and Technology
• /
• v.23 no.5
• /
• pp.702-710
• /
• 1999

U-shaped bellows are usually used to piping system pressure sensor and controller for refriger-ator. Bellows subjected to internal pressure are designed for the purpose of absorbing deformation. Internal pressure on the convolution sidewall and end collar will be applied to an axial load tend-ing to push the collar away from the convolutions. To find out deformation behavior of bellow sub-jected to internal pressure the axisymmetric shell theory using the finite element method is adopted in this paper. U-shaped bellows can be idealized by series of conical frustum-shaped ele-ments because it is axisymmetric shell structure. The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displace-ments are added to r-z cylindrical coordinates of nodal points. The new stiffness matrix of the sys-tem using the new coordinates of nodal points is adopted to calculate the another increments of nodal displacement that is the step by step method is used in this paper. The force required to deflect bellows axially is a function of the dimensions of the bellows and the materials from which they are made. Spring constant is analyzed according to the changing geometric factors of U-shaped bellows. The FEM results were agreed with experiment. Using developed FORTRAN PROGRAM the internal pressure vs. deflection characteristics of a particu-lar bellows can be predicted by input of a few factors.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1518-1525
• /
• 2006

Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.

• Structural Engineering and Mechanics
• /
• v.11 no.4
• /
• pp.357-372
• /
• 2001

In this paper, cantilevered tall structures are treated as cantilever bars with varying cross-section for the analysis of their free longitudinal (or axial) vibrations. Using appropriate transformations, exact analytical solutions to determine the longitudinal natural frequencies and mode shapes for a one step non-uniform bar are derived by selecting suitable expressions, such as exponential functions, for the distributions of mass and axial stiffness. The frequency equation of a multi-step bar is established using the approach that combines the transfer matrix procedure or the recurrence formula and the closed-form solutions of one step bars, leading to a single frequency equation for any number of steps. The Ritz method is also applied to determine the natural frequencies and mode shapes in the vertical direction for cantilevered tall structures with variably distributed stiffness and mass. The formulae proposed in this paper are simple and convenient for engineering applications. Numerical example shows that the fundamental longitudinal natural frequency and mode shape of a 27-storey building determined by the proposed methods are in good agreement with the corresponding measured data. It is also shown that the selected expressions are suitable for describing the distributions of axial stiffness and mass of typical tall buildings.