• Computers and Concrete
• /
• v.22 no.6
• /
• pp.539-550
• /
• 2018

Concrete is highly non-linear material which is originating from the transition zone in the form of micro-cracks, governs material response under various loadings. In this paper, the constitutive models published by many researchers have been used to generate novel stiffness parameters and constitutive curves for concrete. Following such linear material formulations, where the energy is conservative during the curvature, and a nonlinear contribution to the concrete has been made and investigated. In which, nonlinear concrete elastic modulus modeling has been developed that is capable-of representing concrete elasticity for grades ranging from 10 to 140 MPa. Thus, covering the grades range of concrete up to the ultra-high strength concrete, and replacing many concrete models that are valid for narrow ranges of concrete strength grades. This has been followed by the introduction of the nonlinear Hooke's law for the concrete material through the replacement of the Young constant modulus with the nonlinear modulus. In addition, the concept of concrete elasticity index (${\varphi}$) has been proposed and this factor has been introduced to account for the degradation of concrete stiffness in compression under increased loading as well as the multi-stages micro-cracking behavior of concrete under uniaxial compression. Finally, a sub-routine artificial neural network model has been developed to capture the concrete behavior that has been introduced to facilitate the prediction of concrete properties under increased loading.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.347-356
• /
• 2006

In this paper, practical nonlinear inelastic analysis method of 3-dimensional steel structures accounting for gradual yielding with fibers on a section is developed. Geometric nonlinearities of member(p-$\delta$) and frame(p-$\Delta$) are accounted for by using stability functions. Residual stresses are considered by assigning initial stresses to the fiber on the section. The elastic core in a section is investigated at every loading step to determine the axial and bending stiffness reduction. The strain reversal effect is captured by investigating the stress change of each fiber. The proposed analysis proves to be useful in applying for practical analysis and design of three-dimensional steel frames.

• Earthquakes and Structures
• /
• v.26 no.5
• /
• pp.327-336
• /
• 2024

Evaluation of tunnel performance in seismic-prone areas demands efficient means of estimating performance at different hazard levels. The present study introduces an innovative push-over analysis approach which employs the standard earthquake spectrum to simulate the performance of a tunnel. The numerical simulation has taken into account the lining and surrounding rock to calculate the rock-tunnel interaction subjected to a static push-over displacement regime. Elastic perfectly plastic models for the lining and hardening strain rock medium were used to portray the development of plastic hinges, nonlinear deformation, and performance of the tunnel structure. Separately using a computational algorithm, the non-linear response spectrum was approximated from the average shear strain of the rock model. A NATM tunnel in Turkey was chosen for parametric study. A seismic performance curve and two performance thresholds are introduced that are based on the proposed nonlinear seismic static loading approach and the formation of plastic hinges. The tunnel model was also subjected to a harmonic excitation with a smooth response spectrum and different amplitudes in the fully-dynamic phase to assess the accuracy of the approach. The parametric study investigated the effects of the lining stiffness and capacity and soil stiffness on the seismic performance of the tunnel.

• Computers and Concrete
• /
• v.29 no.1
• /
• pp.15-29
• /
• 2022

Concrete-filled steel tubes (CFSTs) are increasingly used as composite sections in structures owing to their excellent load bearing capacity. Therefore, predicting the mechanical behavior of CFST sections under axial compression loading is vital for design purposes. This paper presents the first study on the nonlinear analysis of heated CFSTs with high-strength concrete core containing steel fiber and waste tire rubber under axial compression loading. CFSTs had steel fibers with 0, 1, and 1.5% volume fractions and 0, 5, and 10% rubber particles as sand alternative material. They were subjected to 20, 250, 500, and 750℃ temperatures. Using flow rule and analytical analysis, a model is developed to predict the load bearing capacity of steel tube, and hoop strain-axial strain relationship, and axial stress-volumetric strain relationship of CFSTs. An elastic-plastic analysis method is applied to determine the axial and hoop stresses of the steel tube, considering elastic, yield, and strain hardening stages of steel in its stress-strain curve. The axial stress in the concrete core is determined as the difference between the total experimental axial stress and the axial stress of steel tube obtained from modeling. The results show that steel tube in CFSTs under 750℃ exhibits a higher load bearing contribution compared to those under 20, 250, and 500℃. It is also found that the ratio of load bearing capacity of steel tube at peak point to the load bearing capacity of CFST at peak load is noticeable such that this ratio is in the ranges of 0.21-0.33 and 0.31-0.38 for the CFST specimens with a steel tube thickness of 2 and 3.5 mm, respectively. In addition, after the steel tube yielding, the load bearing capacity of the tube decreases due to the reduction of its axial stiffness and the increase of hoop strain rate, which is in the range of about 20 to 40%.

• Earthquakes and Structures
• /
• v.15 no.2
• /
• pp.123-132
• /
• 2018

In this paper the strain energy density (SED) model is used to analyze the seismic behavior of suspen-domes and a new criterion is established for judging the seismic failure based on a characteristic point in the SED model. Firstly, a nonlinear time-history response analysis was carried out using the finite-element package ANSYS for typical suspen-domes subjected to different ground motions. The seismic responses including nodal displacements, ratios of yielding members, strain energy density and structural maximum deformation energy were extracted corresponding to the increasing peak ground acceleration (A). Secondly, the SED sum ($I_d$) was calculated which revealed that the $I_d-A$ curve exhibited a relatively large change (called a characteristic point) at a certain value of A with a very small load increment after the structures entered the elastic-plastic state. Thirdly, a SED criterion is proposed to judge the seismic failure load based on the characteristic point. Subsequently, the case study verifies the characteristic point and the proposed SED criterion. Finally, this paper describes the unity and application of the SED criterion. The SED method may open a new way for structural appraisal and the SED criterion might give a unified criterion for predicting the failure loads of various structures subjected to dynamic loads.

• Proceedings of the KWS Conference
• /
• 2009.11a
• /
• pp.97-97
• /
• 2009

The use of thin plate increases due to the need for light weight in large ship. Thin plate is easily distorted and has residual stress by welding heat. Therefore, the thin plate should be carefully joined to minimize the welding deformation which costs time and money for repair. For one effort to reduce welding deformation, it is very useful to predict welding deformation before welding execution. There are two methods to analyze welding deformation. One is simple linear analysis. The other is nonlinear analysis. The simple linear analysis is elastic analysis using the equivalent load method or inherent strain method from welding experiments. The nonlinear analysis is thermo-elastic analysis which gives consideration to the nonlinearity of material dependent on temperature and time, welding current, voltage, speed, sequence and constraint. In this study, the welding deformation is analyzed by using thermo-elastic method for PCTC(Pure Car and Truck Carrier) which carries cars and trucks. PCTC uses thin plates of 6mm thickness which is susceptible to welding heat. The analysis dimension is 19,200mm(length) * 13,825mm(width) * 376mm(height). MARC and MENTAT are used as pre and post processor and solver. The boundary conditions are based on the real situation in shipyard. The simulations contain convection and gravity. The material of the thin block is mild steel with $235N/mm^2$ yield strength. Its nonlinearity of conductivity, specific heat, Young's modulus and yield strength is applied in simulations. Welding is done in two pass. First pass lasts 2,100 second, then it rests for 900 second, then second pass lasts 2,100 second and then it rests for 20,000 second. The displacement at 0 sec is caused by its own weight. It is maximum 19mm at the free side. The welding line expands, shrinks during welding and finally experiences shrinkage. It results in angular distortion of thin block. Final maximum displacement, 17mm occurs around welding line. The maximum residual stress happens at the welding line, where the stress is above the yield strength. Also, the maximum equivalent plastic strain occurs at the welding line. The plastic strain of first pass is more than that of second pass. The flatness of plate in longitudinal direction is calculated in parallel with the direction of girder and compared with deformation standard of ${\pm}15mm$. Calculated value is within the standard range. The flatness of plate in transverse direction is calculated in perpendicular to the direction of girder and compared with deformation standard of ${\pm}6mm$. It satisfies the standard. Buckle of plate is calculated between each longitudinal and compared with the deformation standard. All buckle value is within the standard range of ${\pm}6mm$.

• Geomechanics and Engineering
• /
• v.8 no.4
• /
• pp.559-594
• /
• 2015

A nonlinear creep function incorporated into the elastic visco-plastic model may describe the long-term soil deformation more accurately. However, by applying the conventional procedure, there are challenges to determine the model parameters due to limitation of suitable data points. This paper presents a numerical solution to obtain several parameters simultaneously for a nonlinear elastic visco-plastic (EVP) model using the available consolidation data. The finite difference scheme using the Crank-Nicolson procedure is applied to solve a set of coupled partial differential equations of the time dependent strain and pore water pressure dissipation. The model parameters are determined by applying the algorithm of trust-region reflective optimisation in conjunction with the finite difference solution. The proposed method utilises all available consolidation data during dissipation of the excess pore water pressure to determine the required model parameters. Moreover, the reference time in the elastic visco-plastic model can readily be adopted as a unit of time; denoting creep is included in the numerical predictions explicitly from the very first time steps. In this paper, the settlement predictions of thick soft clay layers are presented and discussed to evaluate and compare the accuracy and reliability of the proposed method against the graphical procedure to obtain the model parameters. In addition, comparison of the available experimental results to the numerical predictions confirms the accuracy of the numerical procedure.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.51-68
• /
• 2004

Since the linear elastic fracture analysis has been proved to be insufficient in predicting the failure of strain hardening materials, a number of fracture concepts have been studied which remain applicable in the presence of plasticity near a crack tip. This work thereby presents a new finite element model to predict the elastic-plastic crack-tip field and fatigue life of center-cracked panels(CCP) with ductile fracture under large-scale yielding conditions. Also, this study has been carried out to investigate the path-dependence of J-integral within the plastic zone for elastic-perfectly plastic, bilinear elastic-plastic, and nonlinear elastic-plastic materials. Based on the incremental theory of plasticity, the p-version finite element is employed to account for the accurate values of J-integral, the most dominant fracture parameter, and the shape of plastic zone near a crack tip by using the J-integral method. To predict the fatigue life, the conventional Paris law has been modified by substituting the range of J-value denoted by ${\Delta}J$ for ${\Delta}K$. The experimental fatigue test is conducted with five CCP specimens to validate the accuracy of the proposed model. It is noted that the relationship between the crack length a and ${\Delta}K$ in LEFM analysis shows a strong linearity, on the other hand, the nonlinear relationship between a and ${\Delta}J$ is detected in EPFM analysis. Therefore, this trend will be depended especially in the case of large scale yielding. The numerical results by the proposed model are compared with the theoretical solutions in literatures, experimental results, and the numerical solutions by the conventional h-version of the finite element method.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Proceedings of the Korean Society for Agricultural Machinery Conference
• /
• 1993.10a
• /
• pp.1329-1339
• /
• 1993

Creep tests were performed to determine the nonlinear viscoelastic properties of apples and pears with the creep experiment apparatus designed in this study. Compressive creep characteristics of fruits were tested at two kinds of storage conditions, four periods of storage and three levels of initial stress. Ten replications were made at each treatment combination. The creep behavior of the fruits could be well described by the nonlinear viscoelastic model as a function of initial stress and time. however, for each level of initial stress applied, the compressive behavior of the samples was satisfactorily represented by Burger's model. For all sample fruits, the longer the samples was stored, the higher the instantaneous elastic strain was observed, and the creep progressed at a high rate. These phenomena were even more remarkable on the fruit stored at the normal temperature storage rather than at the low temperature storage.