• Advances in nano research
• /
• v.10 no.4
• /
• pp.385-395
• /
• 2021

This research concentrates on the effects of distributions and volume fractions of carbon nanotubes (CNT) on the nonlinear bending behavior of deep cylindrical panels reinforced by functionally graded carbon nanotubes under thermo-mechanical loading, hitherto not reported in the literature. Assuming the effects of shear deformation and moderately high value of the radius-to-side ratio (R/a), based on the first-order shear deformation theory (FSDT) and von Karman type of geometric nonlinearity, the governing system of equations is obtained. The analytical solution of field equations is carried out using the Ritz method together with the Newton-Raphson iterative scheme. The effects of radius-to-side ratio, temperature change, and boundary conditions on the nonlinear response of the functionally graded carbon nanotubes reinforced composite deep cylindrical panel (FG-CNTRC) are investigated. It is concluded that, among the five possible distribution patterns of CNT, FG-V CNTRC deep cylindrical panel is strongest with the highest bending moment and followed by UD, X, O, and Ʌ-ones. Also, considering the present deep cylindrical panel formulation increases the accuracy of the results. Hence, according to the noticeable amount of R/a in FG-CNTRC cylindrical panels, it is mandatory to apply strain-displacement relations of deep cylindrical panels for bending analysis of FG-CNTRC which certainly is desirable for industrial application.

• Advances in nano research
• /
• v.13 no.4
• /
• pp.393-406
• /
• 2022

In the present article, functionally graded small-scaled plates based on modified strain gradient theory (MSGT) are studied for analyzing the nonlinear bending and post-buckling responses. Von-Karman's assumptions are applied to incorporate geometric nonlinearity and the first-order shear deformation theory is used to model the plates. Modified strain gradient theory includes three length scale parameters and is reduced to the modified couple stress theory (MCST) and the classical theory (CT) if two or all three length scale parameters become zero, respectively. The Ritz method with Legendre polynomials are used to approximate the unknown displacement fields. The solution is found by the minimization of the total potential energy and the well-known Newton-Raphson technique is used to solve the nonlinear system of equations. In addition, numerical results for the functionally graded small-scaled plates are obtained and the effects of different boundary conditions, material gradient index, thickness to length scale parameter and length to thickness ratio of the plates on nonlinear bending and post-buckling responses are investigated and discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Journal of the Korean Recycled Construction Resources Institute
• /
• v.6 no.2
• /
• pp.100-107
• /
• 2018

The material and geometrical nonlinear finite elment analysis of UHPFRC 50M composite box girder was carried out. Constitute law in tension and compressive region of UHPFRC and HPC were modeled based on specimen test. The accuracy of nonlinear FEM analysis was verified by the experimental result of UHPFRC 50M composite girder. The UHPFRC 50M segmental composite box girder which has 1.5% steel fiber of volume fraction, 135MPa compressive strength and 18MPa tensile strength was tested. The post-tensioned UHPFRC composite girder consisted of three segment UHPFRC U-girder and High Strength Concrete reinforced slab. The parts of UHPFRC girder were modeled by 8nodes hexahedron elements and reinforcement bars and tendons were built by 2nodes linear elements by Midas FEA software. The constitutive laws of concrete materials were selected Multi-linear model both of tension and compression function under total strain crack model, which was included in classifying of smeared crack model. The nonlinearity of reinforcement elements and tendon was simulated by Von Mises criteria. The nonlinear static analysis was applied by incremental-iteration method with convergence criteria of Newton-Raphson. The validation of numerical analysis was verified by comparison with experimental result and numerical analysis result of load-deflection response, neutral axis coordinate change, and cracking pattern of girder. The load-deflection response was fitted very well with comparison to the experimental result. The finite element analysis is seen to satisfactorily predict flexural behavioral responses of post-tensioned, reinforced UHPFRC composite box girder.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.1
• /
• pp.27-48
• /
• 1984

This study is conserned with the optimization of reinforced concrete frames using limit state design theory. Formulations of the optimal design for reinforced concrete frame based on the limit state theory turn out to be the nonlinear programming problems which have to deal with the required steel area, the width and effective height of the beam and column section and the moment reduction factor as the design variables. The objective function is formulated as the total construction cost which considers the costs of steel, concrete and forming for the reinforced concrete frames, and the basic constraints are imposed upon both ultimate and serviciability limit state concepts. Also, the stress blocks assumpted in CP110 and Hognestad et al. theory are applied to analysis an ultimate resistant section force for the ultimate limit state and only the criteria of CP110 are used for serviciability limit state. The optimized technique which is applied to solve the nonlinear programming problems for the optimization of reinforced concrete frames is SUMT utilizing the modified Newton-Raphson method. This algorithm is used to test for the two reinforced concrete frames, and then is compared and analysized with the numerical results of reference(10) to examine its convergence, applicability and stability under the same conditions. The results of this study are discussed about the economy comparision of the optimal values for CP110 and Hognestad et al., and the applicability, stability, convergence and validity of this algorithm used herein through the numerical analyses.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.15 no.6
• /
• pp.3448-3454
• /
• 2014

Generally existing structure analysis was applied to elastic analysis basically in practice. Considering the nonlinear material and the nonlinear geometric to be a more precise analysis, for this reason, The necessity for a structual analysis have been constantly required. Therefore, after optimization is performed, designed a simple model which is applied the principle of nonlinear in this study, a structural analysis of existing experienced users, have a aims at presenting theory and a method in order to perform anyone the analysis easily. In this study, the proposed model applied to die ribs, Regarding the shear load, less strain and stress was generated but strength was sufficient. The initial strain and stress was reconfigured to fit the size and shape, A hyperstudy in conjunction with Abaqus with nonlinear structural analysis, revealed an acceptable maximum and minimum range of stress and under the conditions of minimum strain, the plate made with a constant increment. In the experimental models, the plate thickness was given a power of 40 Newton, according to the thickness of the press die through an iterative process. When the stress and strain was applied to the die thickness, 7-8mm thickness could be obtained by optimizing.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.773-784
• /
• 2008

The present study was carried out to examine flow properties in linear drainage channels such as road surface drainage facilities. The finite difference formulation for the varied flow analysis was solved for flow profiles in the channels. Starting the first step at the control section, the Newton-Raphson method was applied for producing numerical solutions of the equation. We considered two types of linear drainage channels, a channel with one outlet at downstream end and a channel with two outlets at both ends. Moreover, the flow analysis for various channel slopes was performed. However, we considered channels with the two outlets of slopes satisfying the condition that the both ends are the control section. The maximum of those slopes was decided from the relation between the channel slope and the location of control section. The flow of a channel with one outlet was calculated upward and downward from the control section existing in channel or upward from the control section at downstream end. The flow of a channel with two outlets at both ends were calculated for upstream and downstream channel segments divided by the water dividend, respectively and the flow analysis was completed when the water depth at the water dividend calculated from upstream end was equal to that calculated from downstream end. If the slope was larger than the critical slope, the channel with two outlets was likely to behave like the channel with one outlet. The maximum water depth was investigated and compared with that calculated additionally from the uniform flow analysis. The uniform flow analysis was likely to lead a excessive design of a drainage channel with mild slope.

• Geotechnical Engineering
• /
• v.8 no.2
• /
• pp.5-20
• /
• 1992

The physical-based and lumped-parameter hydrologic groundwater flow model for predicting the rainfall-triggered rise of groundwater levels in hillside slopes is developed in this paper to assess the risk of landslides. The developed model consists of a vertical infiltration model for unsaturated zone linked to a linear storage reservoir model(LSRM) for saturated zone. The groundwater flow model has uncertain constants like soil depttL slope angle, saturated permeability, and potential evapotranspiration and four free model parameters like a, b, c, and K. The free model parameters could be estimated from known input-output records. The BARD algorithm is uses as the parameter estimation technique which is based on a linearization of the proposed model by Gauss -Newton method and Taylor series expansion. The application to examine the capacity of prediction shows that the developed model has a potential of use in forecast systems of predicting landslides and that the optimal estimate of potential 'a' in infiltration model is the most important in the global optimum analysis because small variation of it results in the large change of the objective function, the sum of squares of deviations of the observed and computed groundwater levels. 본 논문에서는 가파른 산사면에서 산사태의 발생을 예측하기 위한 수문학적 인 지하수 흐름 모델을 개발하였다. 이 모델은 물리적인 개념에 기본하였으며, Lumped-parameter를 이용하였다. 개발된 지하수 흐름 모델은 두 모델을 조합하여 구성되어 있으며, 비포화대 흐름을 위해서는 수정된 abcd 모델을, 포화대 흐름에 대해서는 시간 지체 효과를 고려할 수 있는 선형 저수지 모델을 이용하였다. 지하수 흐름 모델은 토층의 두께, 산사면의 경사각, 포화투수계수, 잠재 증발산 량과 같은 불확실한 상수들과 a, b, c, 그리고 K와 같은 자유모델변수들을 가진다. 자유모델변수들은 유입-유출 자료들로부터 평가할 수 있으며, 이를 위해서 본 논문에서는 Gauss-Newton 방법을 이용한 Bard 알고리즘을 사용하였다. 서울 구로구 시흥동 산사태 발생 지역의 산사면에 대하여 개발된 모델을 적용하여 예제 해석을 수행함으로써, 지하수 흐름 모델이 산사태 발생 예측을 위하여 이용할 수 있음을 입증하였다. 또한, 매개변수분석 연구를 통하여, 변수 a값은 작은 변화에 대하여 목적함수값에 큰 변화를 일으키므로 a의 값에 대한 최적값을 구하는 것이 가장 중요한 요소라는 결론을 얻었다.

• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.4
• /
• pp.297-305
• /
• 2018

In a shipyard, block lifting is an important process in the production of ships and offshore structures. Block lifting is a sensitive process because lifting blocks have to be erected with exact positions and orientations. If we use a numerical method for the process, it is important to find tensions of wires and positions of equalizers to maintain the initial equilibrium position of the block. At this time, equations of motion of the block should be solved to calculate the initial equilibrium position of the block. Because the solving technique changes according to the number of equalizers, a suitable equation for the corresponding problem is required. In this study, three types of equations are proposed to find the initial equilibrium position of the block according to the number of equalizers. The Newton-Raphson's method is used to solve nonlinear simultaneous equations and the optimization method is used to determine the appropriate solution to the undetermined problem. To evaluate the applicability of the proposed methods, the dynamic simulations are performed using the tensions calculated from the proposed methods, and the results are discussed. The results show that the proposed methods can be effectively used to determine initial equilibrium position of the block for the block lifting.