• Structural Engineering and Mechanics
• /
• 제60권4호
• /
• pp.633-653
• /
• 2016

One of the major threats to the stability of classical columns and colonnades are earthquakes. The behavior of columns under high seismic excitation loads is non-linear and complex since rocking, wobbling and sliding failure modes can occur. Therefore, three dimensional simulation approaches are essential to investigate the in-plane and out-of-plane response of such structures during harmonic and seismic loading excitations. Using a software based on the Distinct Element Method (DEM) of analysis, a three dimensional numerical study has been performed to investigate the parameters affecting the seismic behaviour of colonnades' structural systems. A typical section of the two-storey colonnade of the Forum in Pompeii has been modelled and studied parametrically, in order to identify the main factors affecting the stability and to improve our understanding of the earthquake behaviour of such structures. The model is then used to compare the results between 2D and 3D simulations emphasizing the different response for the selected earthquake records. From the results analysis, it was found that the high-frequency motion requires large base acceleration amplitude to lead to the collapse of the colonnade in a shear-slip mode between the drums. However, low-frequency harmonic excitations are more prominent to cause structural collapse of the two-storey colonnade than the high-frequency ones with predominant rocking failure mode. Finally, the 2D analysis found to be unconservative since underestimates the displacement demands of the colonnade system when compared with the 3D analysis.

• Smart Structures and Systems
• /
• 제22권3호
• /
• pp.259-276
• /
• 2018

This paper studies the energies and energy release rate (ERR) for the initially rotationally symmetric compressed (or stretched) in the inward (outward) radial direction of the PZT/Elastic/PZT sandwich circular plate with interface penny-shaped cracks. The investigations are made by utilizing the so-called three-dimensional linearized field equations and relations of electro-elasticity for piezoelectric materials. The quantities related to the initial stress state are determined within the scope of the classical linear theory of piezoelectricity. Mathematical formulation of the corresponding problem and determination of the quantities related to the stress-strain state which appear as a result of the action of the uniformly normal additional opening forces acting on the penny-shaped crack's edges are made within the scope of the aforementioned three-dimensional linearized field equations solution which is obtained with the use of the FEM modelling. Numerical results of the energies and ERR and the influence of the problem parameters on these quantities are presented and discussed for the PZT- 5H/Al/PZT-5H, PZT-4/Al/PZT-4, $BaTiO_3/Al/BaTiO_3$ and PZT-5H/StPZT-5H sandwich plates. In particular, it is established that the magnitude of the influence of the piezoelectricity and initial loading on the ERR increases with crack radius length.

• Computers and Concrete
• /
• 제3권2_3호
• /
• pp.79-90
• /
• 2006

This paper describes the results of an investigation into high mass-low velocity impact behaviour of reinforced concrete beams. Tests have been conducted on fifteen 2.7 m or 1.5 m span beams under drop-weight loads. A high-speed video camera has been used at rates of up to 4,500 frames per second in order to record the crack formation, propagation, particle spallation and scabbing. In some tests the strain in the reinforcement has been recorded using "Durham" strain gauged bars, a technique developed by Scott and Marchand (2000) in which the strain gauges are embedded in the bars, so that the strains in the reinforcement can be recorded without affecting the bond between the concrete and the reinforcement. The impact force acting on the beams has been measured using a load cell placed within the impactor. A high-speed data logging system has been used to record the impact load, strains, accelerations, etc., so that time histories can be obtained. This research has led to the development of computational techniques based on combined continuum/discontinuum methods (finite/discrete element methods) to permit the simulation of impact loaded reinforced concrete beams. The implementation has been within the software package ELFEN (2004). Beams, similar to those tested, have been analysed using ELFEN a good agreement has been obtained for both the load-time histories and the crack patterns.

• 한국강구조학회 논문집
• /
• 제12권3호통권46호
• /
• pp.303-315
• /
• 2000

본 논문에서는 일반화된 21개 강성 매트릭스의 독립 변수를 모두 사용하였고 비등방성 3차원 탄성체의 지배 방정식 및 수치 해석 근사식을 유도하였다. 일반화된 3차원 해석은 2차원 해석의 제한성을 극복하는 정밀해를 보여줄 수 있으며, 두꺼운 보나 판, 또는 쉘에서 전단 변형 효과에 의한 처짐의 증가 효과를 더욱 정밀하게 해석할 수 있다. 따라서 본 논문은 3차원 비등방성 탄성체에 대하여 다양한 경계조건에 따른 유한 차분식을 유도하였으며 이를 전산화하여 해석 프로그램을 개발하였다. 특히, 본 논문에서는 자유경계조건에 대하여 개선된 유한차분법의 적용 방식을 제시하였다. 또한 탄성체의 각 방향 자유경계면에서 경계조건을 해결할 수 있는 일반화된 방식을 제시하였다. 몇가지 수치예제를 통하여 이러한 유한차분 경계처리 방식에 의한 비등방성 3차원 탄성체 해석의 타당성 및 거동을 분석하였다.

• Interaction and multiscale mechanics
• /
• 제5권3호
• /
• pp.267-284
• /
• 2012

In most of the design offices, analysis of the frame is carried out without considering the effect of the rigidity of mat. The analysis of the superstructure without modelling the foundation properly and conversely analysing the foundation system without considering the stiffness of the superstructure may mislead the estimation of the forces. This paper examines the parameters, which affect the interaction and they are grouped into relative stiffness factors ${\kappa}_{rs}$ and ${\kappa}_{sb}$. An interaction analysis is performed for the five storeyed space frame of 3 bays ${\times}$ 5 bays, using ANSYS finite element code. The soil was treated as an isotropic, homogenous and elastic half space medium and the following conclusions were drawn from the analyses. The differential settlement is reduced due to interaction and the performance of the mat depends on ${\kappa}_{sb}$ values. The moments $M_x$ and $M_y$ in the corner column at all the storey levels are higher in the case of the interaction analysis than in the conventional analysis. The axial forces in the peripheral columns increased and to that extent, the inner column axial loads are reduced. In the beam, more variation is seen in the support moments than in the span moments.

• Computers and Concrete
• /
• 제32권6호
• /
• pp.567-576
• /
• 2023

In this study, a 3D microstructure-based model is established to simulate the effective thermal conductivity of cement paste, covering varying influencing factors associated with microstructure and thermal transfer mechanisms. The virtual cement paste divided into colloidal C-S-H and heterogeneous paste are reconstructed based on its structural attributes. Using the two-level hierarchical cement pastes as inputs, a lattice Boltzmann model for heat conduction is presented to predict the thermal conductivity. The results suggest that due to the Knudsen effect induced by the nanoscale pore, the thermal conductivity of air in C-S-H gel pore is significantly decreased, maximumly accounting for 3.3% thermal conductivity of air at the macroscale. In the cement paste, the thermal conductivities of dried and saturated cement pastes are stable at the curing age larger than 100 h. The high water-to-cement ratio can decrease the thermal conductivity of cement paste.

• Advances in Computational Design
• /
• 제3권3호
• /
• pp.233-267
• /
• 2018

Conventional seismic design of concentrically braced frame (CBF) structures suggests that the gusset plate connecting a steel brace to beams and/or columns should be designed as non-dissipative in earthquakes, while the steel brace members should be designed as dissipative elements. These design intentions lead to thicker and larger gusset plates in design on one hand and a potentially under-rated contribution of gusset plates in design, on the other hand. In contrast, research has shown that compact and thinner gusset plates designed in accordance with the elliptical clearance method rather than the conventional standard linear clearance method can enhance system ductility and energy dissipation capacity in concentrically braced steel frames. In order to assess the two design methods, six cyclic push-over tests on full scale models of concentric braced steel frame structures were conducted. Furthermore, a 3D finite element (FE) shell model, incorporating state-of-the-art tools and techniques in numerical simulation, was developed that successfully replicates the response of gusset plate and bracing members under fully reversed cyclic axial loading. Direct measurements from strain gauges applied to the physical models were used primarily to validate FE models, while comparisons of hysteresis load-displacement loops from physical and numerical models were used to highlight the overall performance of the FE models. The study shows the two design methods attain structural response as per the design intentions; however, the elliptical clearance method has a superiority over the standard linear method as a fact of improving detailing of the gusset plates, enhancing resisting capacity and improving deformability of a CBF structure. Considerations were proposed for improvement of guidelines for detailing gusset plates and bracing members in CBF structures.

• Advances in aircraft and spacecraft science
• /
• 제5권3호
• /
• pp.363-383
• /
• 2018

One-dimensional (1D) models of incompressible flows, can be of interest for many applications in which fast resolution times are demanded, such as fluid-structure interaction of flows in compliant pipes and hemodynamics. This work proposes a higher-order 1D theory for the flow-field analysis of incompressible, laminar, and viscous fluids in rigid pipes. This methodology is developed in the domain of the Carrera Unified Formulation (CUF), which was first employed in structural mechanics. In the framework of 1D modelling, CUF allows to express the primary variables (i.e., velocity and pressure fields in the case of incompressible flows) as arbitrary expansions of the generalized unknowns, which are functions of the 1D computational domain coordinate. As a consequence, the governing equations can be expressed in terms of fundamental nuclei, which are invariant of the theory approximation order. Several numerical examples are considered for validating this novel methodology, including simple Poiseuille flows in circular pipes and more complex velocity/pressure profiles of Stokes fluids into non-conventional computational domains. The attention is mainly focused on the use of hierarchical McLaurin polynomials as well as piece-wise nonlocal Lagrange expansions of the generalized unknowns across the pipe section. The preliminary results show the great advantages in terms of computational costs of the proposed method. Furthermore, they provide enough confidence for future extensions to more complex fluid-dynamics problems and fluid-structure interaction analysis.

• 지구물리와물리탐사
• /
• 제14권1호
• /
• pp.98-104
• /
• 2011

탐사 지구물리학에서 수치 모사는 지하매질에서의 탄성파 전파 현상을 이해하는데 중요한 통찰력을 제공한다. 탄성파 모사는 음향파 근사에 의한 수치 모사보다 계산시간이 많이 소요되지만 전단응력 성분을 포함하여 보다 현실적인 파동의 모사를 가능하게 한다. 그러므로 탄성파 모사는 탄성체의 반응을 탐사하는데 적합하다고 할 수 있다. 계산 시간이 길다는 단점을 극복하기 위해 본 논문에서는 그래픽 프로세서(GPU)를 이용하여 탄성파 수치 모사 시간을 단축하고자 하였다. GPU는 많은 수의 프로세서와 광대역 메모리를 갖고 있기 때문에 병렬화된 계산 아카텍쳐에서 사용할 수 있는 장점이 있다. 본 연구에서 사용한 GPU 하드웨어는 NVIDIA Tesla C1060으로 240개의 프로세서로 구성되어 있으며 102 GB/s의 메모리 대역폭을 갖고 있다. NVIDIA에서 개발된 병렬계산 아카텍쳐인 CUDA를 사용할 수 있음에도 불구하고 계산효율을 상당히 향상시키기 위해서는 GPU 장치의 여러 가지 다양한 메모리의 사용과 계산 순서를 최적화해야만 한다. 본 연구에서는 GPU 시스템에서 시간영역 유한차분법을 이용하여 2차원과 3차원 탄성과 전파를 수치 모사하였다. 파동전파 모사에 가장 널리 사용되는 유한차분법 중의 하나인 엇갈린 격자기법을 채택하였다. 엇갈린 격자법은 지구물리학 분야에서 수치 모델링을 위해 사용하기에 충분한 정확도를 갖고 있는 것으로 알려져 있다. 본 논문에서 제안한 모델링기법은 자료 접근 시간을 단축하기 위해 GPU 장치를 메모리 사용을 최적화하여 가능한 더 빠른 메모리를 사용한다. 이점이 GPU를 이용한 계산의 핵심 요소이다. 하나의 GPU 장치를 사용하고 메모리 사용을 최적화함으로써 단일 CPU를 이용할 경우보다 2차원 모사에서는 14배 이상, 3차원에서는 6배 이상 계산시간을 단축할 수 있었다. 세 개의 GPU를 사용한 경우에는 3차원 모사에서 계산효율을 10배 향상시킬 수 있었다.

• 한국지반환경공학회 논문집
• /
• 제17권7호
• /
• pp.15-25
• /
• 2016

본 연구에서는 압밀을 고려한 고등 3차원 유한요소해석을 통하여 압밀이 진행 중인 지반에 근입된 단독말뚝 및 군말뚝의 거동을 연구하였다. 수치해석에서는 단독말뚝 및 말뚝 중심간의 간격이 2.5D인 $4{\times}4$ 및 $6{\times}6$ 군말뚝을 고려하였다. 여기서 D는 말뚝의 직경을 의미한다. 부마찰에 의한 말뚝의 침하 및 부마찰력은 압밀초기 단계에서 비교적 빠르게 발생하는 것으로 분석되었다. 그러나 압밀도가 50~75%를 초과하는 경우 말뚝의 침하 및 부마찰력의 증가량은 상대적으로 크지 않은 것으로 분석되었다. 성토층에서의 부마찰은 신속하게 발현되며 이후의 압밀 단계에서는 일정하게 유지된다. 말뚝에 작용하는 부마찰은 상대변위 및 유효지중응력에 좌우되는 것으로 분석되었다. 압밀 초기단계에서는 상대변위가 큰 영향을 미치는데 비해, 압밀 후반기에서는 유효지중응력이 큰 영향을 주는 것으로 분석되었다. 말뚝-인접지반에서의 전단응력 전이로 인해 말뚝과 인접한 흙의 유효수직응력이 감소하며 이러한 현상은 군말뚝에서 특히 현저하다. 부마찰이 영향을 미치는 영역의 범위는 흙의 유효수직응력의 분포를 고려할 경우 말뚝으로부터 수평으로 20D 정도 되는 것으로 분석되었다. 이에 비해 말뚝에 작용하는 유효수평응력은 far field 조건의 응력과 거의 유사한 것으로 분석되었다.