• Title/Summary/Keyword: Finite-difference

Search Result 3,273, Processing Time 0.03 seconds

Behavior of lightweight aggregate concrete voided slabs

  • Adel A. Al-Azzawi;Ali O, AL-Khaleel
    • Computers and Concrete
    • /
    • v.32 no.4
    • /
    • pp.351-363
    • /
    • 2023
  • Reducing the self-weight of reinforced concrete structures problem is discussed in this paper by using two types of self-weight reduction, the first is by using lightweight coarse aggregate (crushed brick) and the second is by using styropor block. Experimental and Numerical studies are conducted on (LWAC) lightweight aggregate reinforced concrete slabs, having styropor blocks with various sizes of blocks and the ratio of shear span to the effective depth (a/d). The experimental part included testing eleven lightweight concrete one-way simply supported slabs, comprising three as reference slabs (solid slabs) and eight as styropor block slabs (SBS) with a total reduction in cross-sectional area of (43.3% and 49.7%) were considered. The holes were formed by placing styropor at the ineffective concrete zones in resisting the tensile stresses. The length, width, and thickness of specimen dimensions were 1.1 m, 0.6 m, and 0.12 m respectively, except one specimen had a depth of 85 mm (which has a cross-sectional area equal to styropor block slab with a weight reduction of 49.7%). Two shear spans to effective depth ratios (a/d) of (3.125) for load case (A) and (a/d) of (2) for load case (B), (two-line monotonic loads) are considered. The test results showed under loading cases A and B (using minimum shear reinforcement and the reduction in cross-sectional area of styropor block slab by 29.1%) caused an increase in strength capacity by 60.4% and 54.6 % compared to the lightweight reference slab. Also, the best percentage of reduction in cross-sectional area is found to be 49.7%. Numerically, the computer program named (ANSYS) was used to study the behavior of these reinforced concrete slabs by using the finite element method. The results show acceptable agreement with the experimental test results. The average difference between experimental and numerical results is found to be (11.06%) in ultimate strength and (5.33%) in ultimate deflection.

Reasonably completed state assessment of the self-anchored hybrid cable-stayed suspension bridge: An analytical algorithm

  • Kai Wang;Wen-ming Zhang;Jie Chen;Zhe-hong Zhang
    • Structural Engineering and Mechanics
    • /
    • v.90 no.2
    • /
    • pp.159-175
    • /
    • 2024
  • In order to solve the problem of calculating the reasonable completed bridge state of a self-anchored hybrid cable-stayed suspension bridge (SA-HCSB), this paper proposes an analytical method. This method simplifies the main beam into a continuous beam with multi-point rigid supports and solves the support reaction forces. According to the segmented catenary theory, it simultaneously solves the horizontal forces of the main span main cables and the stay cables and iteratively calculates the equilibrium force system on the main beam in the collaborative system bridge state while completing the shape finding of the main span main cable and stay cables. Then, the horizontal forces of the side span main cables and stay cables are obtained based on the balance of horizontal forces on the bridge towers, and the shape finding of the side spans are completed according to the segmented catenary theory. Next, the difference between the support reaction forces of the continuous beam with multiple rigid supports obtained from the initial and final iterations is used to calculate the load of ballast on the side span main beam. Finally, the axial forces and strains of each segment of the main beam and bridge tower are obtained based on the loads applied by the main cable and stay cables on the main beam and bridge tower, thereby obtaining analytical data for the bridge in the reasonable completed state. In this paper, the rationality and effectiveness of this analytical method are verified through a case study of a SA-HCSB with a main span of 720m in finite element analysis. At the same time, it is also verified that the equilibrium force of the main beam under the reasonably completed bridge state can be obtained through iterative calculation. The analytical algorithm in this paper has clear physical significance, strong applicability, and high accuracy of calculation results, enriching the shape-finding method of this bridge type.

A Study on Torsional Stress ratio and Torsional ratio of Curved Girder Bridge by Transfer Matrix Method (전달행렬법에 의한 곡선거더교의 비틀림 응력비와 비틀림 정수비에 관한 연구)

  • Lee, Wong-Hong;Lee, Yoon-Young
    • Journal of the Korea institute for structural maintenance and inspection
    • /
    • v.10 no.1
    • /
    • pp.173-182
    • /
    • 2006
  • In the design of curved girder bridges, the engineer is faced with a complex stress situatiorl. since these types of mutiple-I girder. mono-box girder and twin-box girder are subjected to both bending and torsional force. In general, the torsional forces consist of two part, St. venant's and warping. Thus the procedure for determining the induced stresses in a curved girder is difficult. The transfer matrix method is extensively used for the structural analysis because its merit in the theoretical background and applicability. The technique is attractive for implementation on a numerical solution by means of a computer program coded in Fortran language with a few elements. To demonstrate this fact. it gives good results which compare well with finite difference method. Therefore, in this paper, to clarify the range where the torsional warping stress can be approximated by pure torsional analyzed a critical value of relationships between the torsional stress ratio and torsional ratio.

A Two-Demensional Nonlinear Analysis of Precast Segmental PSC-I Girder with Dry Joint (건식접합부를 갖는 프리캐스트 세그먼트 PSC-I형 거더의 2차원 비선형해석)

  • Kim, Kwang-Soo;Kim, Tae-Wan;Park, Jun-Myung;Park, Sun-Kyu
    • Journal of the Korea institute for structural maintenance and inspection
    • /
    • v.11 no.5
    • /
    • pp.99-106
    • /
    • 2007
  • To satisfy with the increased requirements of cost reduction, labor saving, and rapid construction, the purpose of this study is to investigate the structural behavior of PSC monolithic and spliced girders. Three tests were conducted on small-scale girder specimens. This paper presents the result of experimental studies in terms of the load-deflection behavior. Different joint type and tendon amount were investigated as major variables. The monolithic girder was arranged with three tendons. The spliced girder consisted of five segments connected by three tendons. In addition, five-segmented girders connected by more than three tendons were built to examine the effect of the tendon amount. The experimental results show the difference of behavior between monolithic and spliced girders. Moreover, nonlinear finite element method analysis was utilized to verify the experimental result.

Erection Sequence Analysis of Suspension Bridge Considering to Sliding of Main Cable (주케이블의 슬라이딩을 고려한 현수교의 시공단계 해석)

  • Yhim, Sung-Soon;Kong, Min-Sik;Kim, Byeong-Cheol
    • Journal of the Korea institute for structural maintenance and inspection
    • /
    • v.13 no.3 s.55
    • /
    • pp.164-172
    • /
    • 2009
  • Anchors and saddles are used to have sufficient geometrical rigidity and make target configuration of main cable of suspension bridge. Neglecting the sliding effect at saddles, points at them have been idealized as fixed nodes in lots of former studies. In general, sliding effects are reported to show significant structural behaviors of main cable and cause to the different responses of bridges. During early erection steps of the suspension bridge, especially, the sliding effect occurs easily because there is large difference of cable tension between main and side span in removing set-back ropes or not applying set-back. This study presents the finite element analysis considering to cable sliding effect and shows the comparison of differences between sliding and non-sliding at election sequence. The analysis of sliding between main cables and saddles needs to obtain more realistic responses because the analysis result can represent unfavorably different responses of bridges. Moreover, the sliding analysis method and results in this study can be used to basic criteria in engineering design and construction steps.

Quantum transport of doped rough-edged graphene nanoribbons FET based on TB-NEGF method

  • K.L. Wong;M.W. Chuan;A. Hamzah;S. Rusli;N.E. Alias;S.M. Sultan;C.S. Lim;M.L.P. Tan
    • Advances in nano research
    • /
    • v.17 no.2
    • /
    • pp.137-147
    • /
    • 2024
  • Graphene nanoribbons (GNRs) are considered a promising alternative to graphene for future nanoelectronic applications. However, GNRs-based device modeling is still at an early stage. This research models the electronic properties of n-doped rough-edged 13-armchair graphene nanoribbons (13-AGNRs) and quantum transport properties of n-doped rough-edged 13-armchair graphene nanoribbon field-effect transistors (13-AGNRFETs) at different doping concentrations. Step-up and edge doping are used to incorporate doping within the nanostructure. The numerical real-space nearest-neighbour tight-binding (NNTB) method constructs the Hamiltonian operator matrix, which computes electronic properties, including the sub-band structure and bandgap. Quantum transport properties are subsequently computed using the self-consistent solution of the two-dimensional Poisson and Schrödinger equations within the non-equilibrium Green's function method. The finite difference method solves the Poisson equation, while the successive over-relaxation method speeds up the convergence process. Performance metrics of the device are then computed. The results show that highly doped, rough-edged 13-AGNRs exhibit a lower bandgap. Moreover, n-doped rough-edged 13-AGNRFETs with a channel of higher doping concentration have better gate control and are less affected by leakage current because they demonstrate a higher current ratio and lower off-current. Furthermore, highly n-doped rough-edged 13-AGNRFETs have better channel control and are less affected by the short channel effect due to the lower value of subthreshold swing and drain-induced barrier lowering. The inclusion of dopants enhances the on-current by introducing more charge carriers in the highly n-doped, rough-edged channel. This research highlights the importance of optimizing doping concentrations for enhancing GNRFET-based device performance, making them viable for applications in nanoelectronics.

Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit (그래픽 프로세서를 이용한 탄성파 수치모사의 계산속도 향상)

  • Nakata, Norimitsu;Tsuji, Takeshi;Matsuoka, Toshifumi
    • Geophysics and Geophysical Exploration
    • /
    • v.14 no.1
    • /
    • pp.98-104
    • /
    • 2011
  • Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

Prediction of Heat-treatment Time of Black Pine Log Damaged by Pine Wilt Disease (소나무재선충병 피해를 받은 곰솔 원목의 열처리 소요시간 예측)

  • Han, Yeonjung;Seo, Yeon-Ok;Jung, Sung-Cheol;Eom, Chang-Deuk
    • Journal of the Korean Wood Science and Technology
    • /
    • v.44 no.3
    • /
    • pp.370-380
    • /
    • 2016
  • The black pine logs damaged by pine wilt disease in Jeju-do were heat-treated to extend the utilization of domestic trees damaged by pine wilt disease. The heat-treatment of wood requires wood to be heated to $56^{\circ}C$ for 30 min at the core. The average moisture content and top-diameter of the black pine logs were ranged from 46% to 141% and from 180 mm to 500 mm, respectively. And the basic specific gravity and oven-dry specific gravity of the black pine logs were 0.47 and 0.52, respectively. The time required for heat-treatment at $105^{\circ}C$ temperature was ranged from 7.7 h to 44.2 h, depending on moisture content and top-diameter. The temperature distribution was used to predict the time required for heat-treatment of black pine log with various moisture contents and top-diameters using finite difference method. The thermal properties of wood including the thermal conductivity and specific heat in accordance with moisture content were calculated. Heat transfer coefficient for mixed convection in form of adding natural convection and forced convection was used for heat transfer analysis. The error between the measured and predicted values ranged from 3% to 45%. The predicted times required for heat-treatment of black pine log with 50% moisture content and 200 mm, 300 mm, and 400 mm top-diameter were 10.9 h, 18.3 h, and 27.0 h, respectively. If the initial moisture content of black pine log is 75%, heat treatment times of 13.6 h, 22.5 h, and 32.8 h were predicted in accordance with top-diameter. And if the initial moisture content of black pine log is 100%, heat treatment times of 16.2 h, 26.5 h, and 38.2 h were predicted in accordance with top-diameter. When the physical properties of logs damaged by pine wilt disease are presented, these results can be applicable to the heat-treatment of red pine and Korean pine logs as well.

Simulation Model of Two Dimensional Heat Transfer in Grain Bin (저장곡물(貯藏穀物) Bin내(內)의 삼차원(三次元) 열전달모형(熱傳達模型))

  • Han, Kwang Jin;Kim, Man Soo
    • Korean Journal of Agricultural Science
    • /
    • v.12 no.1
    • /
    • pp.118-127
    • /
    • 1985
  • The grain temperature is a crucial factor determining the deterioration rate of stored grain. Therefore, it is used to be predicted in order to evaluate the various stored methods rapidly and inexpensively. In this study, a mathematical model was developed to simulate the temperatures of grain stored in a cylinderical bin. It was formulated for the two dimensional heat transfer by the finite difference method. Then, it was verified statistically using the actual test deta and the predicted. The changes of grain temperature were analyed using the simulated data of one year for a safe stoarge and the following results were obtained: 1. Simulation model developed by the finite difference method was validated with the actual and the predicted grain temperatures and it's result showed that it could predict the grain temperature of storage bin reasonably well. 2. Grain temperature near the wall of storage bin were changed with $6-7^{\circ}C$ higher then average atmospheric temperature from June to September. Therefore, the parts of stored grain near the wall is supposed to be deteriorated fast. 3. When the dimension of bin diameter is about the same as the bed height, the changes of grain temperature of radial direction was higher than the verticals. 4. The predicted temperature showed that the grain temperature of which were from the end of April to mid October were higher than the safe storage limit at Yusung, Korea.

  • PDF

Numerical Analysis of Pile Foundation Considering the Thawing and Freezing Effects (융해-동결작용을 고려한 말뚝 기초에 관한 수치해석 연구)

  • Park, Woo-Jin ;Park, Dong-Su;Shin, Mun-Beom;Seo, Young-Kyo
    • Journal of the Korean Geotechnical Society
    • /
    • v.39 no.5
    • /
    • pp.51-63
    • /
    • 2023
  • Numerical analysis was conducted to determine the effect of soil behavior by thawing and freezing of seasonal frozen soil on pile foundations. The analysis was performed using the finite element method (FEM) to simulate soil-pile interaction based on the atmosphere temperature change. Thermomechanical coupled modeling using FEM was applied with the temperature-dependent nonlinear properties of the frozen soil. The analysis model cases were applied to the MCR and HDP models to simulate the elastoplastic behavior of soil. The numerical analysis results were analyzed and compared with various conditions having different length and width sizes of the pile. The results of the numerical analysis showed t hat t he HDP model was relat ively passive, and t he aspect and magnit ude of t he bearing capacit y and displacement of the pile head were similar depending on the length and width of the pile conditions. The vertical displacement of the pile head by thawing and freezing of the ground showed a large variation in displacement for shorter length conditions. In the MCR model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0387 and 0.0277 m, respectively. In the HDP model, the vertical displacement appeared in the maximum thaw settlement and frost heaving of 0.0367 and 0.0264 m, respectively. The results of the pile bearing capacity for the two elastoplastic models showed a larger difference in the width condition than the length condition of the pile, with a maximum of about 14.7% for the width L condition, a maximum of about 5.4% for M condition, and a maximum of about 5.3% for S condition. The significance of the effect on the displacement of the pile head and the bearing capacity depended on the pile-soil contact area, and the difference depended on the presence or absence of an active layer in the soil and its thickness.