• Title/Summary/Keyword: discrete-element

Search Result 686, Processing Time 0.032 seconds

CAD Interface using Topology Optimization (위상최적설계 결과를 이용한 CAD 인터페이스)

  • Kim, Seong-Hoon;Min, Seung-Jae;Lee, Sang-Hun
    • Korean Journal of Computational Design and Engineering
    • /
    • v.14 no.4
    • /
    • pp.281-289
    • /
    • 2009
  • Topology optimization has been widely used for the optimal structure design for weight reduction and high performance. Since the result of three-dimensional topology optimization is represented by the discrete material distribution in finite elements, it is hard to interpret from a design point of view. In this paper, the method for interpreting three-dimensional topology optimization resuIt into a series of cross-sectional curve representation is proposed and interfaced with the existing CAD system for the practical use. The concept of node density and virtual grid is introduced to transform element density values into grid density and material boundaries in each cross section are identified based on the element volume rate to satisfy the amount of material specified in the original design intent. Design exampIes show that three-dimensional topology result can be converted into a form of curve CAD model and the seamless interface with CAD software can be achieved.

Linear shell elements for active piezoelectric laminates

  • Rama, Gil;Marinkovic, Dragan Z.;Zehn, Manfred W.
    • Smart Structures and Systems
    • /
    • v.20 no.6
    • /
    • pp.729-737
    • /
    • 2017
  • Piezoelectric composite laminates are a powerful material system that offers vast options to improve structural behavior. Successful design of piezoelectric adaptive structures and testing of control laws call for highly accurate, reliable and numerically efficient numerical tools. This paper puts focus onto linear and geometrically nonlinear static and dynamic analysis of smart structures made of such a material system. For this purpose, highly efficient linear 3-node and 4-node finite shell elements are proposed. Both elements employ the Mindlin-Reissner kinematics. The shear locking effect is treated by the discrete shear gap (DSG) technique with the 3-node element and by the assumed natural strain (ANS) approach with the 4-node element. Geometrically nonlinear effects are considered using the co-rotational approach. Static and dynamic examples involving actuator and sensor function of piezoelectric layers are considered.

Ray Effect Analysis Using the Discrete Elements Method in X-Y Geometry (2차원 직각좌표계에서 DEM을 이용한 ray effect의 해석)

  • Choi, Ho-Sin;Kim, Jong-Kyung
    • Journal of Radiation Protection and Research
    • /
    • v.17 no.1
    • /
    • pp.43-56
    • /
    • 1992
  • As one of the methods to ameliorate the ray effects which are the nature of anomalous computational effects due to the discretization of the angular variable in discrete ordinates approximations, a computational program, named TWODET (TWO dimensional Discrete Element Transport), has developed in 2 dimensional cartesian coordinates system using the discrete elements method, in which the discrete angle quadratures are steered by the spatially dependent angular fluxes. The results of the TWODET calculation with K-2, L-3 discrete angular quadratures, in the problem of a centrally located, isotropically emitting flat source in an absorbing square, are shown to be more accurate than that of the DOT 4.3 calculation with S-10 full symmetry angular quadratures, in remedy of the ray effect at the edge flux distributions of the square. But the computing time of the TWODET is about 4 times more than that of the DOT 4.3. In the problem of vacuum boundaries just outside of the source region in an absorbing square, the results of the TWODET calculation are shown severely anomalous ray effects, due to the sudden discontinuity between the source and the vacuum, like as the results of the DOT 4.3 calculation. In the probelm of an external source in an absorbing square in which a highly absorbing medium is added, the results of the TWODET calculation with K-3, L-4 show a good ones like as, somewhat more than, that of the DOT 4.3 calculation with S-10.

  • PDF

A Boundary Element Analysis for Damage and Failure Process of Brittle Rock using ERACOD (FRACOD를 이용한 취성 암석의 손상 및 파괴에 대한 경계요소 해석)

  • ;Baotang Shen;Ove Stephansson
    • Tunnel and Underground Space
    • /
    • v.14 no.4
    • /
    • pp.248-260
    • /
    • 2004
  • Damage in brittle rock due to stress increase starts from initiation of microcracks, and then results in failure by forming macro failure planes due to propagation and coalescence of these discrete cracks. Conventionally, continuum approaches using macro-failure criteria or a number of elasto-plastic models have been major solution to implement rock damage and failure. However, actual brittle failure processes can be better described in phenomenological approach if initiation and propagation of discrete fractures are explicitly considered. This study presents damage and failure process of rock using a boundary element code, FRACOD, which has been developed to model fracturing process of rocks. Through a series of numerical uniaxial compressive tests, the feasibility of the developed model was verified, and realistic rock failure process was reproduced considering scale effects in rocks. In addition, the fracturing process and the corresponding rock damage in the vicinity of deep shaft in rock mass were presented as an application of this approach. This approach will be expected to contribute to finding better engineering solutions for the analysis of stability problems in brittle rock masses.

A Study on Estimation of the Collapse Pattern of Road Sink Using Distinct Element Method (개별요소법을 이용한 도로함몰 붕괴양상 추정에 관한 연구)

  • Ham, Myoung Soo;Park, Seon Woo;Lee, Hyun Dong
    • Journal of Korean Society of Disaster and Security
    • /
    • v.12 no.2
    • /
    • pp.57-63
    • /
    • 2019
  • The road sinks in the sewer line or subway section are affected by the ground characteristics. Therefore, it is necessary to accurately identify the relationship between the soil properties and the ground motion in the area where cavities occurred in order to establish a countermeasure against the road sink. In this paper, simulation was performed by using EDEM program, which is one of the Discrete Element Method programs, for sandy soil and clayey soil, which are most common in alluvial deposits, with different locations and sizes of cavities in the underground. As a result, it was found that the sink size occurred more in the sandy soil than in the cohesive soil. Deeper and larger cavity is more likely to occur the road sink In the sand soil model while road sink in the clay model is easy to occur when the cavity is more shallower.

Modelling of Large Triaxial Test with Rockfill Materials by Distinct Element Method (개별요소법에 의한 락필재료의 대형삼축압축시험 모델링)

  • Jeon, Je-Sung;Kim, Ki-Young;Shin, Dong-Hoon
    • Journal of the Korean Geotechnical Society
    • /
    • v.22 no.10
    • /
    • pp.111-120
    • /
    • 2006
  • In this research, numerical simulations by PFC considering discrete element method are conducted to predict experimental results of large triaxial compression test with rockfill material for dam construction. For generation of compacted assembly with specific grain size distribution and initial material porosity, the clump logic method and expansion of generated particles are adapted. To predict stress-stain behavior of large triaxial test, discrete particle modelling is applied with micro parameters which are chosen by calibration process. It is expected that distinct particle modelling method could be used as a useful tool to investigate micro and macro behavior associated with geotechnical problems and develop a numerical laboratory.

KRULL RING WITH UNIQUE REGULAR MAXIMAL IDEAL

  • Chang, Gyu Whan
    • Korean Journal of Mathematics
    • /
    • v.15 no.2
    • /
    • pp.115-119
    • /
    • 2007
  • Let R be a Krull ring with the unique regular maximal ideal M. We show that R has a regular prime element and reg-$dimR=1{\Leftrightarrow}R$ is a factorial ring and reg-$dim(R)=1{\Rightarrow}M$ is invertible ${\Leftrightarrow}R{\varsubsetneq}[R:M]{\Leftrightarrow}M$ is divisorial ${\Leftrightarrow}$ reg-$htM=1{\Rightarrow}R$ is a rank one discrete valuation ring. We also show that if M is generated by regular elements, then R is a rank one discrete valuation ring ${\Rightarrow}$ R is a factorial ring and reg-dim(R)=1.

  • PDF

MINIMAL LOCALLY STABILIZED Q1-Q0 SCHEMES FOR THE GENERALIZED STOKES PROBLEM

  • Chibani, Alima;Kechkar, Nasserdine
    • Journal of the Korean Mathematical Society
    • /
    • v.57 no.5
    • /
    • pp.1239-1266
    • /
    • 2020
  • In this paper, some novel discrete formulations for stabilizing the mixed finite element method Q1-Q0 (bilinear velocity and constant pressure approximations) are introduced and discussed for the generalized Stokes problem. These are based on stabilizing discontinuous pressure approximations via local jump operators. The developing idea consists in a reduction of terms in the local jump formulation, introduced earlier, in such a way that stability and convergence properties are preserved. The computer implementation aspects and numerical evaluation of these stabilized discrete formulations are also considered. For illustrating the numerical performance of the proposed approaches and comparing the three versions of the local jump methods alongside with the global jump setting, some obtained results for two test generalized Stokes problems are presented. Numerical tests confirm the stability and accuracy characteristics of the resulting approximations.

Aeroelastic analysis of bridges using FEM and moving grids

  • Selvam, R. Panneer;Govindaswamy, S.;Bosch, Harold
    • Wind and Structures
    • /
    • v.5 no.2_3_4
    • /
    • pp.257-266
    • /
    • 2002
  • In the recent years flow around bridges are investigated using computer modeling. Selvam (1998), Selvam and Bosch (1999), Frandsen and McRobie (1999) used finite element procedures. Larsen and Walther (1997) used discrete vorticity procedure. The aeroelastic instability is a major criterion to be checked for long span bridges. If the wind speed experienced by a bridge is greater than the critical wind speed for flutter, then the bridge fails due to aeroelastic instability. Larsen and Walther (1997) computed the critical velocity for flutter using discrete vortex method similar to wind tunnel procedures. In this work, the critical velocity for flutter will be calculated directly (free oscillation procedure) similar to the approaches reported by Selvam et al. (1998). It is expected that the computational time required to compute the critical velocity using this approach may be much shorter than the traditional approach. The computed critical flutter velocity of 69 m/s is in reasonable comparison with wind tunnel measurement. The no flutter and flutter conditions are illustrated using the bridge response in time.

Efficacy of CFRP configurations for shear of RC beams: experimental and NLFE

  • Shuraim, Ahmed B.
    • Structural Engineering and Mechanics
    • /
    • v.39 no.3
    • /
    • pp.361-382
    • /
    • 2011
  • This paper presents the results of an investigation on shear strengthening of RC beams externally reinforced with CFRP composite. A total of six full-scale beams of four CFRP strengthened and two unstrengthened were tested in the absence of internal stirrups in the shear span. The strengthening configurations contained two styles: discrete uniformly spaced strips and customized wide strips over B-regions. The composite systems provided an increase in ultimate strength as compared to the unstrengthened beams. Among the three layouts that had the same area of CFRP, the highest contribution was provided by the customized layout that targeted the B-regions. A comparative study of the experimental results with published empirical equations was conducted in order to evaluate the assumed effective strains. The empirical equations were found to be unconservative. Nonlinear finite element (NLFE) models were developed for the beams. The models agreed with test results that targeting the B-region was more effective than distributing the same CFRP area in a discrete strip style over shear spans. Moreover, the numerical models predicted the contribution of different configurations better than the empirical equations.