• Proceedings of the KSR Conference
• /
• 2004.06a
• /
• pp.955-960
• /
• 2004

In this study, three dimensional FE analysis of concrete slab track has been performed in order to develop the realistic design of precast concrete slab track. The precast slab track system including the precast concrete slab panel and the grout layer is modeled using the three dimensional solid element with crack softening effect. The input load is computed from the one dimensional beam element model constituting the rail and several discrete springs. To investigate the effect of the longitudinal connection of slab panels, two different systems-continuous and discrete systems - are modeled. The analytical results show that the stresses of both the slab panel and the grout layer are in the range of linear elastic, and, at the interface between two adjacent panels, the primary stresses of the grout layer of the discrete system are higher than those of the continuous system. However, The overall stress levels of the grout layer are very low relative to the strength of th grout.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.399-404
• /
• 2000

Functionally graded materials(FGMs) involve dual-phase graded layers in which two different constituents are mixed continuously and functionally according to a given volume fraction. For the analysis of their thermo-mechanical response, conventional homogenized methods have been widely employed in order to estimate equivalent material properties of the graded layer. However, such overall estimations are insufficient to accurately predict the local behavior. In this paper, we compare the thermo-elastic behaviors predicted by several overall material-property estimation techniques with those obtained by discrete analysis models utilizing the finite element method, for various volume fractions and loading conditions.

• Journal of the Korean Mathematical Society
• /
• v.60 no.2
• /
• pp.273-302
• /
• 2023

In this paper, a fully discrete numerical scheme for the viscoelastic Oldroyd flow is considered with an introduced auxiliary variable. Our scheme is based on the finite element approximation for the spatial discretization and the backward Euler scheme for the time discretization. The integral term is discretized by the right trapezoidal rule. Firstly, we present the corresponding equivalent form of the considered model, and show the relationship between the origin problem and its equivalent system in finite element discretization. Secondly, unconditional stability and optimal error estimates of fully discrete numerical solutions in various norms are established. Finally, some numerical results are provided to confirm the established theoretical analysis and show the performances of the considered numerical scheme.

• Structural Engineering and Mechanics
• /
• v.15 no.4
• /
• pp.477-493
• /
• 2003

Masonry is not a simple material, the influence of mortar joints as a plane of weakness is a significant feature and this makes the numerical modelling of masonry very difficult especially when dynamic (seismic) analysis is involved. In order to develop a simple numerical model for masonry under earthquake load, an analytical model based on Distinct Element Method (DEM) is being developed. At the first stage, the model is applied to simulate the in-plane shear behaviour of an unreinforced masonry wall with and without opening where the testing results are available for comparison. In DEM, a solid is represented as an assembly of discrete blocks. Joints are modelled as interface between distinct bodies. It is a dynamic process and specially designed to model the behaviour of discontinuities. The numerical solutions obtained from the distinct element analysis are validated by comparing the results with those obtained from existing experiments and finite element modelling.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.19 no.5
• /
• pp.27-37
• /
• 2020

Structural optimization problems with discrete design variables require more function calculations (or finite element analyses) than those in the continuous design space. In this study, a method to find an optimal solution in the discrete design of the truss structure is presented, reducing the number of function calculations. Because a continuous optimal solution is the Karush-Kuhn-Tucker point that satisfies the optimality condition, it is assumed that the discrete optimal solution is around the continuous optimum. Then, response values such as weight, displacement, and stress are predicted using approximate models-referred to as hybrid metamodels-within specified design ranges. The discrete design method using the hybrid metamodels is used as a post-process of the continuous optimization process. Standard truss design problems of 10-bar, 25-bar, 15-bar, and 52-bar are solved to show the usefulness of this method. The results are compared with those of existing methods.

• Computers and Concrete
• /
• v.21 no.2
• /
• pp.189-197
• /
• 2018

A discrete element approach is used to investigate the effects of confining stress on the shear behaviour of joint's bridge area. A punch-through shear test is used to model the concrete cracks under different shear and confining stresses. Assuming a plane strain condition, special rectangular models are prepared with dimension of $75mm{\times}100mm$. Within the specimen model and near its four corners, four equally spaced vertical notches of the same depths are provided so that the central portion of the model remains intact. The lengths of notches are 35 mm. and these models are sequentially subjected to different confining pressures ranging from 2.5 to 15 MPa. The axial load is applied to the punch through the central portion of the model. This testing and models show that the failure process is mostly governed by the confining pressure. The shear strengths of the specimens are related to the fracture pattern and failure mechanism of the discontinuities. The shear behaviour of discontinuities is related to the number of induced shear bands which are increased by increasing the confining pressure while the cracks propagation lengths are decreased. The failure stress and the crack initiation stress both are increased due to confining pressure increase. As a whole, the mechanisms of brittle shear failure changes to that of the progressive failure by increasing the confining pressure.

• Computers and Concrete
• /
• v.28 no.5
• /
• pp.521-532
• /
• 2021

Present paper focuses on the modeling of size effect on the compressive strength of normal concrete with the application of Discrete Element Method (DEM). Test specimens with different size and shape were cast and uniaxial compressive strength test was performed on each sample. Five different concrete mixes were used, all belonging to a different normal strength concrete class (C20/25, C30/37, C35/45, C45/55, and C50/60). The numerical simulations were carried out by using the PFC 5 software, which applies rigid spheres and contacts between them to model the material. DEM modeling of size effect could be advantageous because the development of micro-cracks in the material can be observed and the failure mode can be visualized. The series of experiments were repeated with the model after calibration. The relationship of the parallel bond strength of the contacts and the laboratory compressive strength test was analyzed by aiming to determine a relation between the compressive strength and the bond strength of different sized models. An equation was derived based on Bazant's size effect law to estimate the parallel bond strength of differently sized specimens. The parameters of the equation were optimized based on measurement data using nonlinear least-squares method with SSE (sum of squared errors) objective function. The laboratory test results showed a good agreement with the literature data (compressive strength is decreasing with the increase of the size of the specimen regardless of the shape). The derived estimation models showed strong correlation with the measurement data. The results indicated that the size effect is stronger on concretes with lower strength class due to the higher level of inhomogeneity of the material. It was observed that size effect is more significant on cube specimens than on cylinder samples, which can be caused by the side ratios of the specimens and the size of the purely compressed zone. A limit value for the minimum size of DE model for cubes and cylinder was determined, above which the size effect on compressive strength can be neglected within the investigated size range. The relationship of model size (particle number) and computational time was analyzed and a method to decrease the computational time (number of iterations) of material genesis is proposed.

• Journal of Drive and Control
• /
• v.21 no.2
• /
• pp.44-52
• /
• 2024

In this study, we predicted PTO power requirements based on torque predicted by the discrete element method and the multi-body dynamics coupling method. Six different scenarios were simulated to predict PTO power requirements in different soil conditions. The first scenario was a tillage operation on cohesionless soil, and the field was modeled using the Hertz-Mindlin contact model. In the second through sixth scenarios, tillage operations were performed on viscous soils, and the field was represented by the Hertz-Mindlin + JKR model for cohesion. To check the influence of surface energy, a parameter to reproduce cohesion, on the power requirement, a simple regression analysis was performed. The significance and appropriateness of the regression model were checked and found to be acceptable. The study findings are expected to be used in design optimization studies of agricultural machinery by predicting power requirements using the discrete element method and the multi-body dynamics coupling method and analyzing the effect of soil cohesion on the power requirement.

• Journal of the Korea Concrete Institute
• /
• v.17 no.6 s.90
• /
• pp.1053-1064
• /
• 2005

The mechanical damage of concrete is normally attributed to the formation of microcracks and their propagation and coalescence into macroscopic cracks. This physical degradation is caused from progressive and hierarchical damage of the microstructure due to debonding and slip along bimaterial interfaces at the mesoscale. Their growth and coalescence leads to initiation of hairline discrete cracks at the mesoscale. Eventually, single or multiple major discrete cracks develop at the macroscale. In this paper, from this conceptual model of mechanical damage in concrete, the computational efforts were made in order to characterize physical cracks and how to quantify the damage of concrete materials within the laws of thermodynamics with the aid of interface element in traditional finite element methodology. One dimensional effective traction/jump constitutive interface law is introduced in order to accommodate the normal opening and tangential slips on the interfaces between different materials(adhesion) or similar materials(cohesion) in two and three dimensional problems. Mode I failure and mixed mode failure of various geometries and boundary conditions are discussed in the sense of crack propagation and their spent of fracture energy under monotonic displacement control.

• Proceedings of the Korean Society for Agricultural Machinery Conference
• /
• 2000.11b
• /
• pp.117-123
• /
• 2000

Using the modified DEM (Discrete Element Method), which we proposed in order to improve the accuracy of the simulation, soil behavior and reaction by lugs of rotating wheel and a soil cutting process by a high speed blade were calculated and compared with experimental data. The DEM is one of computational mechanics, where the object body is supposed as an assembly of small particles called elements and not a continuum as in the case of FEM. We can easily treat some discrete phenomena such as cracking, separating and sliding by the DEM. We had to modify the original mechanical model, which induced too free movement of elements, adding a tension spring, which would display the role of soil adhesion. The results of DEM simulations were successful from both the soil behavior and reaction points of view.