• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.6
• /
• pp.349-357
• /
• 2019

Finite element method (FEM) is utilized in the development of products to realistically analyze and predict the mechanical behavior of materials in various fields. However, the approach based on the numerical analysis of glass fiber reinforced plastic (GFRP) composites, for which the fiber orientation and strain rate affect the mechanical properties, has proven to be challenging. The purpose of this study is to define and evaluate the mechanical properties of glass fiber reinforced plastic composites using the numerical analysis models of Digimat, a linear, nonlinear multi-scale modeling program for various composite materials such as polymers, rubber, metal, etc. In addition, the aim is to predict the behavior of realistic polymeric composites. In this regard, the tensile properties according to the fiber orientation and strain rate of polybutylene terephthalate (PBT) with short fiber weight fractions of 30wt% among various polymers were investigated using references. Information on the fiber orientation was calculated based on injection analysis using Moldflow software, and was utilized in the finite element model for tensile specimens via a mapping process. LS-Dyna, an explicit commercial finite element code, was used for coupled analysis using Digimat to study the tensile properties of composites according to the fiber orientation and strain rate of glass fibers. In addition, the drawbacks and advantages of LS-DYNA's various anisotropic material models were compared and evaluated for the analysis of glass fiber reinforced plastic composites.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.37-43
• /
• 2006

Dynamic numerical analysis of geotechnical problems requires a way to simulate the decrease of energy as the domain of interest gets larger. This phenomenon is usually referred to as radiation damping or geometric attenuation and it is distinguished from material damping in which elastic energy is actually dissipated by viscous, hysteretic, or other mechanism. The fact that the domain of analysis in numerical modeling must be chosen, however, causes a need for special attention at the boundary. This observation leads directly to the idea of determining the dynamic response of the interior region from a finite model consisting of the interior region subjected to a boundary condition which ensures that all energy arriving at the boundary is absorbed. This paper presents a simple methodology to simulate transmitting boundaries condition using viscoelastic infinite elements within the recently developed "OpenSees" finite element code. The methodology used here provides that the level of absorption for traveling waves is efficient enough for practical purposes, but unsatisfactory for the case of sharp incident angles. The effectiveness of the infinite elements for the absorption of incident waves at boundaries is evaluated via example analysis.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.269-279
• /
• 1994

This paper concentrates on the finite element analysis of concrete structures considering the material nonlinearity and time-dependent structural behavior. Using the rotating crack model among the smeared cracking model, the structural behavior up to ultimate load is simulated, and concrete is assumed to be an orthotropic material. Especially to include the tension stiffening effect in bending behavior, a criterion based on the fracture mechanics concept is introduced and the numerical error according to the finite element mesh size can be minimized through the application of the proposed criterion. Besides, the governing equation for steel is systematized by embeded model to cope with the difficulty in modeling of complex geometry. Finally, to trace the structural behavior with time under cracked and/or uncracked section, an algorithm for the purpose of time-dependent analysis is formulated in plane stress-strain condition by the age-adjusted effective modulus method.

• Journal of Korean Society of Steel Construction
• /
• v.25 no.6
• /
• pp.659-669
• /
• 2013

The purpose of this study is to investigate the ultimate strength and curling influence using finite element analysis based on the previous test results of ferritic stainless steel bolted connections. Results obtained from finite element analysis (FEA) were compared with those of test results and the validity of numerical modeling was verified. The conditions of curing occurrence for edge/end distance were investigated and ultimate strength reduction ratio caused by curling was estimated quantitatively. Moreover, the ultimate behaviors such as fracture mode and ultimate strength by FEA were compared with those predicted by current design specifications such as AISC, KBC2009, AIJ and AISI.

• Structural Engineering and Mechanics
• /
• v.61 no.4
• /
• pp.497-510
• /
• 2017

A framed structure may be composed of two sub-structures, which are linked by a hinged joint. One sub-structure is the primary system and the other is the secondary system. The primary system, which is subjected to the periodic external load, can give rise to an auto-parametric resonance of the second system. Considering the geometric-stiffness effect produced by the axially internal force, the element equation of motion is derived by the extended Hamilton's principle. The element equations are then assembled into the global non-homogeneous Mathieu-Hill equations. The Newmark's method is introduced to solve the time-history responses of the non-homogeneous Mathieu-Hill equations. The energy-growth exponent/coefficient (EGE/EGC) and a finite-time Lyapunov exponent (FLE) are proposed for determining the auto-parametric instability boundaries of the structural system. The auto-parametric instabilities are numerically analyzed for the two frames. The influence of relative stiffness between the primary and secondary systems on the auto-parametric instability boundaries is investigated. A phenomenon of the "auto-parametric internal resonance" (the auto-parametric resonance of the second system induced by a normal resonance of the primary system) is predicted through the two numerical examples. The risk of auto-parametric internal resonance is emphasized. An auto-parametric resonance experiment of a ${\Gamma}$-shaped frame is conducted for verifying the theoretical predictions and present calculation method.

• Steel and Composite Structures
• /
• v.35 no.4
• /
• pp.463-474
• /
• 2020

In this paper, the influence of adding multi-walled carbon nanotube (MWCNT) on the pull behavior of steel and GFRP bars in ultra-high-performance concrete (UHPC) was examined experimentally and numerically. For numerical analysis, 3D nonlinear finite element modeling (FEM) with the help of ABAQUS software was used. Mechanical properties of the specimens, including Young's modulus, tensile strength and compressive strength, were extracted from the experimental results of the tests performed on standard cube specimens and for different values of weight percent of MWCNTs. In order to consider more realistic assumptions, the bond between concrete and bar was simulated using adhesive surfaces and Cohesive Zone Model (CZM), whose parameters were obtained by calibrating the results of the finite element model with the experimental results of pullout tests. The accuracy of the results of the finite element model was proved with conducting the pullout experimental test which showed high accuracy of the proposed model. Then, the effect of different parameters such as the material of bar, the diameter of the bar, as well as the weight percent of MWCNT on the bond behavior of bar and UHPC were studied. The results suggest that modifying UHPC with MWCNT improves bond strength between concrete and bar. In MWCNT per 0.01 and 0.3 wt% of MWCNT, the maximum pullout strength of steel bar with a diameter of 16 mm increased by 52.5% and 58.7% compared to the control specimen (UHPC without nanoparticle). Also, this increase in GFRP bars with a diameter of 16 mm was 34.3% and 45%.

• Steel and Composite Structures
• /
• v.12 no.1
• /
• pp.53-72
• /
• 2012

In this study, the physically-based failure models for matrix and fibers in compression and tension loading are introduced. For the 3D stress based fiber kinking model a modification is proposed for calculation of the fiber misalignment angle. All of these models are implemented into the finite element code by using the advantage of damage variable and the numerical results are discussed. To investigate the matrix failure model, purely in-plane transverse compression experiments are carried out on the specimens made by Glass/Epoxy to obtain the fracture surface angle and then a comparison is made with the calculated numerical results. Furthermore, shear failure of $({\pm}45)_s$ model is investigated and the obtained numerical results are discussed and compared with available experimental results. Some experiments are also carried out on the woven laminated composites to investigate the fracture pattern in the matrix failure mode and shown that the presented matrix failure model can be used for the woven composites. Finally, the obtained numerical results for stress based fiber kinking model and improved ones (strain based model) are discussed and compared with each other and with the available results. The results show that these models can predict the kink band angle approximately.

• Structural Engineering and Mechanics
• /
• v.85 no.1
• /
• pp.89-104
• /
• 2023

This work presents a proposal for employing reduced numerical integration in the formulation of the 4-node quadrilateral solid finite element. The use of these low-order integration rules leads to numerical instabilities such as those producing the hourglass effect. The proposed procedure allows evaluating a given constitutive model only in one integration point, achieving an attractive computational cost reduction and, also, successfully controls the hourglass effect. A validation of the proposal is included and discussed throughout the paper. To show the efficiency of the proposal, several application examples of masonry structures are studied and discussed. To represent the non-linear mechanical behaviour of masonry a plastic-damage model is implemented within the application of this sub-integration scheme. Also, in order to have a full and computationally efficient strategy to determine the behaviour of masonry structures, involving its evolution to collapse, a homogenization technique with a macro-modeling approach is used. The methodology discussed throughout this paper demonstrates a substantial computational cost reduction and an improved approximation of the non-linear problem evidenced by a reduction of up to 85% of the computational time for some cases.

• Wind and Structures
• /
• v.28 no.3
• /
• pp.181-190
• /
• 2019

Light-frame wood structures have the ability to carry gravity loads. However, their performance during severe wind storms has indicated weakness with respect to resisting uplift wind loads exerted on the roofs of residential houses. A common failure mode observed during almost all main hurricane events initiates at the roof-to-wall connections (RTWCs). The toe-nail connections typically used at these locations are weak with regard to resisting uplift loading. This issue has been investigated at the Insurance Research Lab for Better Homes, where full-scale testing was conducted of a house under appropriate simulated uplift wind loads. This paper describes the detailed and sophisticated numerical simulation performed for this full-scale test, following which the numerical predictions were compared with the experimental results. In the numerical model, the nonlinear behavior is concentrated at the RTWCs, which is simulated with the use of a multi-linear plastic element. The analysis was conducted on four sets of uplift loads applied during the physical testing: 30 m/sincreased by 5 m/sincrements to 45 m/s. At this level of uplift loading, the connections exhibited inelastic behavior. A comparison with the experimental results revealed the ability of the sophisticated numerical model to predict the nonlinear response of the roof under wind uplift loads that vary both in time and space. A further component of the study was an evaluation of the load sharing among the trusses under realistic, uniform, and code pressures. Both the numerical model and the tributary area method were used for the load-sharing calculations.