• Title/Summary/Keyword: Finite Elements in Time

Search Result 401, Processing Time 0.031 seconds

Direct Time Domain Method for Nonlinear Earthquake Response Analysis of Dam-Reservoir Systems (댐-호소계 비선형 지진응답의 직접시간영역 해석기법)

  • Lee, Jin-Ho;Kim, Jae-Kwan
    • Journal of the Earthquake Engineering Society of Korea
    • /
    • v.14 no.3
    • /
    • pp.11-22
    • /
    • 2010
  • An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

Multi-scale Analysis of Thin film Considering Surface Effects (표면효과를 고려한 박막구조의 멀티스케일 해석)

  • Cho, Maeng-Hyo;Choi, Jin-Bok;Jung, Kwang-Sub
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.20 no.3
    • /
    • pp.287-292
    • /
    • 2007
  • In general, the response of bulk material is independent of its size when it comes to considering classical elasticity theory. Because the surface to bulk ratio of the large solids is very small, the influence of surface can be negligible. But the surface effect plays important role as the surface to bulk ratio becomes larger, that is, the contribution of the surface effect must be considered in nano-size elements such as thin film or beam structure. Molecular dynamics computation has been a conventional way to analyze these ultra-thin structures but this method is limited to simulate on the order of $10^6{\sim}10^9$ atoms for a few nanoseconds, and besides, very time consuming. Analysis of structures in submicro to micro range(thin-film, wire etc.) is difficult with classical molecular dynamics due to the restriction of computing resources and time. Therefore, in this paper, the continuum-based method is considered to simulate the overall physical and mechanical properties of the structures in nano-scale, especially, for the thin-film.

A Study on the Vibration Phenomenon of 6 Bundle Boltless Spacer Damper (6도체 무볼트형 Spacer Damper의 진동현상에 관한 연구)

  • 김영달
    • Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
    • /
    • v.17 no.3
    • /
    • pp.110-118
    • /
    • 2003
  • Spacer dampers maintain the constant gaps between each conductor in a bundle conductor-transmission line, and are installed at proper intervals to keep a line from all sorts of damages derived from the vibration energy caused by mechanical or electrical external factors. It is most important to embody a technology which considers difficulties of maintenance and repair, and has optimum elements in order to prevent accidents such as destruction by fire or the snapping of a wire by the effect of vibration phenomenon coming from transmission lint In the resent thesis, therefore, the analysis of vibratory characteristics of spacer damper is set up by analytical methods such as the analysis of conductor motion's governing equation, the equation of spacer damper's motion, spacer damper-fastened wire's motion in a span and the numerical analysis of finite difference method. Furthermore, the installation distance between spacer dampers was scrutinized by simulations of various vibration phenomena which change at any time as actual conditions do, and hereafter difference method. Furthermore, the installation distance between spacer dampers was scrutinized by simulations of various vibration phenomena which change at any time as actual conditions do, and hereafter we will be able to analyze all kinds of vibration phenomena coming from a boltless spacer damper with 6 bundle conductor for 765kV transmission line based on new analytical methods.

Fracture and Hygrothermal Effects in Composite Materials (복합재의 파괴와 hygrothermal 효과에 관한 연구)

  • Kook-Chan Ahn;Nam-Kyung Kim
    • Journal of the Korean Society of Safety
    • /
    • v.11 no.4
    • /
    • pp.143-150
    • /
    • 1996
  • This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

  • PDF

Memory Organization for a Fuzzy Controller.

  • Jee, K.D.S.;Poluzzi, R.;Russo, B.
    • Proceedings of the Korean Institute of Intelligent Systems Conference
    • /
    • 1993.06a
    • /
    • pp.1041-1043
    • /
    • 1993
  • Fuzzy logic based Control Theory has gained much interest in the industrial world, thanks to its ability to formalize and solve in a very natural way many problems that are very difficult to quantify at an analytical level. This paper shows a solution for treating membership function inside hardware circuits. The proposed hardware structure optimizes the memoried size by using particular form of the vectorial representation. The process of memorizing fuzzy sets, i.e. their membership function, has always been one of the more problematic issues for the hardware implementation, due to the quite large memory space that is needed. To simplify such an implementation, it is commonly [1,2,8,9,10,11] used to limit the membership functions either to those having triangular or trapezoidal shape, or pre-definite shape. These kinds of functions are able to cover a large spectrum of applications with a limited usage of memory, since they can be memorized by specifying very few parameters ( ight, base, critical points, etc.). This however results in a loss of computational power due to computation on the medium points. A solution to this problem is obtained by discretizing the universe of discourse U, i.e. by fixing a finite number of points and memorizing the value of the membership functions on such points [3,10,14,15]. Such a solution provides a satisfying computational speed, a very high precision of definitions and gives the users the opportunity to choose membership functions of any shape. However, a significant memory waste can as well be registered. It is indeed possible that for each of the given fuzzy sets many elements of the universe of discourse have a membership value equal to zero. It has also been noticed that almost in all cases common points among fuzzy sets, i.e. points with non null membership values are very few. More specifically, in many applications, for each element u of U, there exists at most three fuzzy sets for which the membership value is ot null [3,5,6,7,12,13]. Our proposal is based on such hypotheses. Moreover, we use a technique that even though it does not restrict the shapes of membership functions, it reduces strongly the computational time for the membership values and optimizes the function memorization. In figure 1 it is represented a term set whose characteristics are common for fuzzy controllers and to which we will refer in the following. The above term set has a universe of discourse with 128 elements (so to have a good resolution), 8 fuzzy sets that describe the term set, 32 levels of discretization for the membership values. Clearly, the number of bits necessary for the given specifications are 5 for 32 truth levels, 3 for 8 membership functions and 7 for 128 levels of resolution. The memory depth is given by the dimension of the universe of the discourse (128 in our case) and it will be represented by the memory rows. The length of a world of memory is defined by: Length = nem (dm(m)+dm(fm) Where: fm is the maximum number of non null values in every element of the universe of the discourse, dm(m) is the dimension of the values of the membership function m, dm(fm) is the dimension of the word to represent the index of the highest membership function. In our case then Length=24. The memory dimension is therefore 128*24 bits. If we had chosen to memorize all values of the membership functions we would have needed to memorize on each memory row the membership value of each element. Fuzzy sets word dimension is 8*5 bits. Therefore, the dimension of the memory would have been 128*40 bits. Coherently with our hypothesis, in fig. 1 each element of universe of the discourse has a non null membership value on at most three fuzzy sets. Focusing on the elements 32,64,96 of the universe of discourse, they will be memorized as follows: The computation of the rule weights is done by comparing those bits that represent the index of the membership function, with the word of the program memor . The output bus of the Program Memory (μCOD), is given as input a comparator (Combinatory Net). If the index is equal to the bus value then one of the non null weight derives from the rule and it is produced as output, otherwise the output is zero (fig. 2). It is clear, that the memory dimension of the antecedent is in this way reduced since only non null values are memorized. Moreover, the time performance of the system is equivalent to the performance of a system using vectorial memorization of all weights. The dimensioning of the word is influenced by some parameters of the input variable. The most important parameter is the maximum number membership functions (nfm) having a non null value in each element of the universe of discourse. From our study in the field of fuzzy system, we see that typically nfm 3 and there are at most 16 membership function. At any rate, such a value can be increased up to the physical dimensional limit of the antecedent memory. A less important role n the optimization process of the word dimension is played by the number of membership functions defined for each linguistic term. The table below shows the request word dimension as a function of such parameters and compares our proposed method with the method of vectorial memorization[10]. Summing up, the characteristics of our method are: Users are not restricted to membership functions with specific shapes. The number of the fuzzy sets and the resolution of the vertical axis have a very small influence in increasing memory space. Weight computations are done by combinatorial network and therefore the time performance of the system is equivalent to the one of the vectorial method. The number of non null membership values on any element of the universe of discourse is limited. Such a constraint is usually non very restrictive since many controllers obtain a good precision with only three non null weights. The method here briefly described has been adopted by our group in the design of an optimized version of the coprocessor described in [10].

  • PDF

An analysis of horizontal deformation of a pile in soil using a beam-on-spring model for the prediction of the eigenfrequency of the offshore wind turbine (해상풍력터빈의 고유진동수 예측을 위한 지반에 인입된 파일의 탄성지지보 모델 기반 수평 거동 해석)

  • Ryue, Jungsoo;Baik, Kyungmin;Kim, Tae-Ryong
    • The Journal of the Acoustical Society of Korea
    • /
    • v.35 no.4
    • /
    • pp.261-271
    • /
    • 2016
  • In the prediction of response of a pile in soil, numerical approaches such as a finite element method are generally applied due to complicate nonlinear behaviors of soils. However, the numerical methods based on the finite elements require heavy efforts in pile and soil modelling and also take long computing time. So their usage is limited especially in the early design stage in which principal dimensions and properties are not specified and tend to vary. On the contrary, theoretical approaches adopting linear approximations for soils are relatively simple and easy to model and take short computing time. Therefore, if they are validated to be reliable, they would be applicable in predicting responses of a pile in soil, particularly in early design stage. In case of wind turbines regarded in this study, it is required to assess their natural frequencies in early stages, and in this simulation the supporting pile inserted in soil could be replaced with a simplified elastic boundary condition at the bottom end of the wind turbine tower. To do this, analysis for a pile in soil is performed in this study to extract the spring constants at the top end of the pile. The pile in soil can be modelled as a beam on elastic spring by assuming that the soils deform within an elastic range. In this study, it is attempted to predict pile deformations and influence factors for lateral loads by means of the beam-on-spring model. As two example supporting structures for wind turbines, mono pile and suction pile models with different diameters are examined by evaluating their influence factors and validated by comparing them with those reported in literature. In addition, the deflection profiles along the depth and spring constants at the top end of the piles are compared to assess their supporting features.

Observation of behavior of the Ahlat Gravestones (TURKEY) at seismic risk and their recognition by QR code

  • Isik, Ercan;Antep, Baris;Buyuksarac, Aydin;Isik, Mehmet Fatih
    • Structural Engineering and Mechanics
    • /
    • v.72 no.5
    • /
    • pp.643-652
    • /
    • 2019
  • Protection of cultural heritage and carrying it to the future are at the top of the significant topics of research and implementation in engineering in the 21st century. There are several historical structures in the district of Ahlat located in the east of Turkey on the Lake Van Basin that has harbored many civilizations. Some of such works are the gravestones that are found in the Ahlat Seljuk Cemetery, which is the oldest and largest cemetery in the district. This study firstly provides information about the Ahlat Seljuk Cemetery and the gravestones found in it. Observation-based structural analyses were carried out on these gravestones that are found in this area that are known to have belonged to different civilizations based on their physical and constructional characteristics. These stones were built out of Ahlat stone as single pieces. Information is provided on the damages that have occurred on the gravestones in time and their causes. In general, losses of mass, abrasions, separations, collapses and calcifications due to natural conditions, as well as vegetative formations, were observed in the gravestones. To provide an example of other gravestones within the context of the study, the gravestone that is known to belong to the person named Nureddin Ebu Hasan was selected. As a result of the modeling that was carried out for this gravestone by using the finite elements method, modal analyses were carried out. With these analyses, for the gravestone, period, effective mass participation rates and stress values were calculated. The stress values that were obtained in this study were compared to the material safety stress values that were obtained in previous studies. Additionally, QR code application was created for the gravestone that was selected as an example in the study, and information on this gravestone was transferred to an electronic environment. The QR code application includes different language options, visuals of the gravestone and information on the gravestone. The QR application was also supported with a video of the cemetery where the gravestone is located. With this application, access to information about gravestones will be possible by using tablets and smartphones. With a QR code to be created for each gravestone, these gravestones will obtain identity cards.

A STUDY ON THE HEAT GENERATION OF BONE DRILLING BURS ACCORDING TO IRRIGATION DURING IMPLANTATION OF DENIAL IMPLANT (치과 임플랜트 식립시 관주에 따른 골 천공 기구의 열 발생에 관한 연구)

  • Moon, Eun-Soo;Lim, Heun-Song;Cho, In-Ho
    • The Journal of Korean Academy of Prosthodontics
    • /
    • v.37 no.4
    • /
    • pp.433-453
    • /
    • 1999
  • There are a lot of elements affecting the success of implant prosthesis. The quality of surgical procedure is considered as one of the key factors. To avoid the excess heat generation is an important element of successful osseointegration and it can be achieved by using a gentle surgical technique with a sharp instrument in bone drilling. This study was performed to measure and analyze comparatively the heat transmitted to sur-rounding bone at a distance of 0.5mm from the periphery of the drill hole in each drilling stage. The results were as follows. In standard system, the temperature of surrounding bone tissue ranged from $29.2^{\circ}C\;to\;48.3^{\circ}C$ with irrigation and from $34.6^{\circ}C\;to\;84.3^{\circ}C$ without irrigation. And in wide system, the temperature of surrounding bone tissue ranged from $29.5^{\circ}C\;to\;52.5^{\circ}C$ with irrigation and from $34.8^{\circ}C\;to\;87.8^{\circ}C$ without irrigation. And the temperature ranges exceeded the threshold without irrigation, while showing less than the threshold by the cooling effect of irrigation. In comparing standard system with wide system, although there was no significant difference, ${\phi}4.3mm$ pilot and ${\phi}4.3mm$ twist drill of wide system showed high value and wide system showed slightly high elevation of temperature in all depth in fixture installation. In the finite element analysis, the calculated value by the Fourier's cooling law were applied to the bone drilling surface. And through analysis using different irrigation temperatures at $28^{\circ}C,\;15^{\circ}C\;and\;5^{\circ}C$, and according to the time. The result was that the cooling water at least below $15^{\circ}C$ was required to maintain the temperature of surrounding bone less than threshold in bone drilling, the cooling water below $5^{\circ}C$ was required to gain more sufficient cooling effect, and cooling over 5 seconds was needed after bone drilling for sufficient effect.

  • PDF

REMOVAL TORQUE AND BONE FORMATION OF ORTHODONTIC MINISCREW IMPLANT (교정용 미니스크류 임플랜트의 제거회전력 및 골형성에 관한 연굴)

  • Yun, Young-Kuk;Ryu, Jae-Jun;Suh, Kyu-Won
    • The Journal of Korean Academy of Prosthodontics
    • /
    • v.45 no.4
    • /
    • pp.492-505
    • /
    • 2007
  • Statement of problem: An orthodontic miniscrew implant has been used as a skeletal anchorage for orthodontic treatment. However, any relation among the influence of the cortical bone, morphologic differences of orthodontic miniscrew implants and new bone formation hasn't been made clear yet. Purpose: The purpose of this study was to evaluate whether the orthodontic miniscrew implant could work as an intraoral skeletal anchorage immediately and stably for orthodontic treatment after insertion of it. Material and methods: Two types of orthodontic miniscrew implants were used in this experiment; tapered type and straight type. One hundred and sixty eight orthodontic miniscrew implants were inserted into the tibiae of 21 rabbits and sacrificed on 3, 7, 11, 14, 21 and 28days later after insertion of them to study removal torque values and histologic and histomorphometric analyses. Results: The results were as follows. 1. The removal torque values of the tapered type were higher than those of the straight type in all groups(p<0.05). 2. There wasn't any distinguishing differences between the tapered type and the straight type about the new bone formation percentage. 3. The removal torque values for both the tapered type and the straight type were gradually decreased at early stages of the test but started to increase at the 7 days group of the straight type and the 11 days group of the tapered type. 4. New bone formation percentage was increased gradually for both the tapered and the straight types as time passed(p<0.05). 5. It was found that the tapered type showed lower values in the cortical bone about both the maximum equilibratory stress distribution and the maximum principal stress distribution than the straight type in linear finite elements analysis. Conclusion: According to the research, the removal torque values were decreased at 7 days group of the tapered type and 11 days group of the straight type after the insertion of the orthodontic miniscrew implants in tibiae of rabbits. Considering the human bone activity, it is better to apply the orthodontic force $3{\sim}4$ weeks later than to apply it immediately after the insertion of orthodontic miniscrew implants. Considering that general orthodontic force is about $250{\sim}500$ grams, the tapered type can be worked as a stable skeletal anchor age in an orthodontic treatment even if the orthodontic force is applied on it immediately after the insertion of it.

Explicit Transient Simulation of SH-waves Using a Spectral Element Method (스펙트럴 요소법을 이용한 SH파 전파의 외연적 시간이력해석)

  • Youn, Seungwook;Kang, Jun Won
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.31 no.2
    • /
    • pp.87-95
    • /
    • 2018
  • This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.