• Geomechanics and Engineering
• /
• 제5권2호
• /
• pp.119-142
• /
• 2013

Pile load tests using Osterberg cells (O-cell) were conducted on cast-in-place concrete piles founded in oil sand fill and in situ oil sand at an industrial plant site in Fort McMurray, Alberta, Canada. Interpreted pile test results show that very high pile shaft resistance (with the Bjerrum-Burland or Beta coefficient of 2.5-4.5) against oil sand could be mobilized at small relative displacements of 2-3% of shaft diameter. Finite element simulations based on linear elastic and elasto-plastic models for oil sand materials were used to analyze the pile load test measurements. Two constitutive models yield comparable top-down load versus pile head displacement curves, but very different behaviour in mobilization of pile shaft and end bearing resistances. The elasto-plastic model produces more consistent matching in both pile shaft and end bearing resistances whereas the linear elastic under- and over-predicts the shaft and end bearing resistances, respectively. The mobilization of high shaft resistance in oil sand under pile load is attributed to the very dense and interlocked structure of oil sand which results in high matrix stiffness, high friction angle, and high shear dilation.

• Structural Engineering and Mechanics
• /
• 제49권4호
• /
• pp.479-489
• /
• 2014

In the present study a parametric analysis is conducted to study the effect of pile dimension and soil properties on the nonlinear dynamic response of pile subjected to lateral sinusoidal load at the pile head. The study is conducted on soil-pile model of different pile diameter, pile length and soil modulus, and results are compared to get the effect. The soil-pile system is modelled using Finite element method. The programming is done in MATLAB. Time history analysis of model is done for varying non-dimensional frequency of load and the results are compared to get the non-dimensional frequency at which pile head displacement is maximum in each case. Maximum possible bending moment and soil-pile interacting forces for the dynamic excitation of the pile is also compared. When results are compared with the linear response, it is observed that non-dimensional frequency is reduced in nonlinear response on account of reduction in the soil stiffness due to yielding. Nonlinear response curve shows high amplitude as compared to linear response curve.

• Structural Engineering and Mechanics
• /
• 제58권6호
• /
• pp.949-966
• /
• 2016

The second order analysis taking place due to non-linear behavior of the structures under the mechanical and geometric factors through implementing exact and approximate methods is an indispensible issue in the analysis of such structures. Among the exact methods is the slope-deflection method that due to its simplicity and efficiency of its relationships has always been in consideration. By solving the differential equations of the modified slope-deflection method in which the effect of axial compressive force is considered, the stiffness matrix including trigonometric entries would be obtained. The complexity of computations with trigonometric functions causes replacement with their Maclaurin expansion. In most cases only the first two terms of this expansion are used but to obtain more accurate results, more elements are needed. In this paper, the effect of utilizing higher order terms of Maclaurin expansion on reducing the number of required elements and attaining more rapid convergence with less error is investigated for the Bernoulli beam with various boundary conditions. The results indicate that when using only one element along the beam length, utilizing higher order terms in Maclaurin expansion would reduce the relative error in determining the critical buckling load and kinematic parameters in the second order analysis.

• Steel and Composite Structures
• /
• 제32권5호
• /
• pp.657-670
• /
• 2019

This paper proposes a one-dimensional fiber beam element model taking account of materially non-linear behavior, benefiting the highly efficient elastic-plastic analysis of girders with shear-lag effects. Based on the displacement-based fiber beam-column element, two additional degrees of freedom (DOFs) are added into the proposed model to consider the shear-lag warping deformations of the slabs. The new finite element (FE) formulations of the tangent stiffness matrix and resisting force vector are deduced with the variational principle of the minimum potential energy. Then the proposed element is implemented in the OpenSees computational framework as a newly developed element, and the full Newton iteration method is adopted for an iterative solution. The typical materially non-linear behaviors, including the cracking and crushing of concrete, as well as the plasticity of the reinforcement and steel girder, are all considered in the model. The proposed model is applied to several test cases under elastic or plastic loading states and compared with the solutions of theoretical models, tests, and shell/solid refined FE models. The results of these comparisons indicate the accuracy and applicability of the proposed model for the analysis of both concrete box girders and steel-concrete composite girders, under either elastic or plastic states.

• Structural Engineering and Mechanics
• /
• 제12권6호
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.

• 대한토목학회논문집
• /
• 제29권1A호
• /
• pp.31-43
• /
• 2009

본 논문에서는 다양한 차종의 영향을 반영할 수 있고, 차량과 교량의 연성 운동방정식을 구성하여 시간 단계별 직접해를 산정할 수 있는 수치해석기법을 제시하였다. 운동방정식의 해는 직접적분법인 Newmark ${\beta}$을 이용하여 해석 단계별로 구성된 유효강성행렬과 유효하중벡터를 바탕으로 정적평형방정식의 해를 구하는 원리와 동일하게 산정하였다. 또한 해석의 효율성을 증진시키기 위하여 유효강성행렬은 Skyline 법에 의해 재구성하였으며, Cholesky의 행렬 분해기법을 동시에 적용하여 직접적인 역행렬 계산에서 야기되는 오차의 발생을 최소화 하였다. 또한 기존선 철도차량인 새마을 PMC 열차와 디젤 견인 무궁화 열차에 대한 3차원 정밀수치해석 모델을 개발하였고, 각 차량은 차체와 전 후방 대차에 각각 6자유도씩 고려하여 총 18자유도로 수치모델을 작성하였다. 교량은 3차원 공간뼈대 요소를 이용하여 모델링하였고, 차륜과 레일 접촉면의 불규칙성은 미국의 FRA에서 규정하고 있는 연직방향 및 횡방향틀림에 대한 PSD 함수를 이용하여 궤도틀림을 수치적으로 구현하였다. 제시된 수치해석 기법은 12 m, 18 m형 판형교의 실측결과를 이용하여 타당성을 검증하였으며, 실측 및 수치해석결과는 교량의 1차 휨 고유진동수의 2.0배를 기준으로 Low pass filtering 하였다.

• Ocean Systems Engineering
• /
• 제4권4호
• /
• pp.327-342
• /
• 2014

The tension leg platform (TLP) is one of the compliant structures which are generally used for deep water oil exploration. With respect to the horizontal degrees of freedom, it behaves like a floating structure moored by vertical tethers which are pretension due to the excess buoyancy of the platform, whereas with respect to the vertical degrees of freedom, it is stiff and resembles a fixed structure and is not allowed to float freely. In the current study, a numerical study for square TLP using modified Morison equation was carried out in the time domain with water particle kinematics using Airy's linear wave theory to investigate the effect of changing the tether tension force on the stiffness matrix of TLP's, the dynamic behavior of TLP's; and on the fatigue stresses in the cables. The effect was investigated for different parameters of the hydrodynamic forces such as wave periods, and wave heights. The numerical study takes into consideration the effect of coupling between various degrees of freedom. The stiffness of the TLP was derived from a combination of hydrostatic restoring forces and restoring forces due to cables. Nonlinear equation was solved using Newmark's beta integration method. Only uni-directional waves in the surge direction was considered in the analysis. It was found that for short wave periods (i.e., 10 sec.), the surge response consisted of small amplitude oscillations about a displaced position that is significantly dependent on tether tension force, wave height; whereas for longer wave periods, the surge response showed high amplitude oscillations that is significantly dependent on wave height, and that special attention should be given to tethers fatigue because of their high tensile static and dynamic stress.

• 한국지진공학회논문집
• /
• 제5권3호
• /
• pp.1-8
• /
• 2001

교량 구성요소의 설계지진력은 현행 국내 도로교설계기준에 의하면 설계지진을 가하여 얻어진 탄성지진력을 구조형식에 따른 응답수정계수로 나눔으로써 결정되어진다. 말뚝기초가 채택된 교량시스템의 탄성지진력의 크기는 말뚝기초의 모형화 방법에 따라 크게 달라질 수 있다. 이 논문에서는 근사적이고 실용적인 말뚝기초의 모형화 기법을 제시하였다. 이 모형화 기법에서는 말뚝기초의 강도를 횡방향으로 반복하중을 가진 현장시험으로 얻은 말뚝-지반의 상호작용이 고려된 지반반력-변위 곡선을 이용한 말뚝의 수평방향 강도와 탄성 축변형은 물론 선단지지력 및 주변마찰력을 고려한 말뚝의 수직방향 강도로 나타내는 것이다. 예제 교량의 해석을 수행하여 제시된 절차가 타당성있고 적용 가능한 교량의 지진응답해석용 말뚝기초의 모형화 기법임을 검증하였다.

• 대한기계학회논문집
• /
• 제15권1호
• /
• pp.1-10
• /
• 1991

본 연구에서는 고속 대변형 탄소성 변형 과정을 해석할 수 있는 프로그램(NET )을 개발하고 이것을 실린더 및 링 성형 문제에 적용하여 마찰 및 관성 효과가 변형 거동에 미치는 영향을 규명하여 보았다.

• 한국전산구조공학회논문집
• /
• 제17권1호
• /
• pp.11-20
• /
• 2004

본 연구에서는 직선 강상자형 거더의 단면변형에 의한 변형 및 응력계산을 위한 Matlab 해석프로그램을 개발하고자 한다. 이를 위하여 단면변형이론을 요약하고 빔유사이론을 제시한다. 이후 탄성지반위의 보-기둥부재의 지배방정식을 제시하고, 일반화된 고유치해석을 통하여 집중 및 분포하중을 받는 보요소의 엄밀한 강성행렬을 계산한다. 본 연구의 효율성과 정확성을 입증하기 위하여 격벽을 갖는 상자형 거더의 뒤틀림응력을 계산하고 유한요소해와 비교한다.