• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.1032-1040
• /
• 2001

It is widely accepted that the pressure variation of interstitial fluid is one of the most important factors in bone physiology. In order to understand the role of interstitial fluid on porous bony structure, a consideration for the biomechanical interactions between fluid and solid constituents within bone is required. In this study, a poroelastic theory was applied to investigate the elastic behavior of calf vertebral trabecular bone composed of the porous solid trabeculae and the viscous bone marrow. The poroelastic behavior of trabecular bone in a uniaxial stress condition was simulated using a commercial finite difference analysis software (FLAC, Itasca Consulting Group, USA), and tested for 5 different strain rates, i.e., 0.001, 0.01, 0.1, and 10 per second. The material properties of the calf vertebral trabecular bone were utilized from the previous experimental study. Two asymptotic poroelastic responses, the drained and undrained deformations, were predicted. From the predicted results for the simulated five strain rates, it was found that the pore pressure generation has a linearly increasing behavior when the strain rate is the highest at 10 per second, otherwise it showed a nonlinear behavior. The pore pressure generation with respect to the strain was found to be increased as the strain rate increased. The elastic moduli predicted at each strain were 208.3, 212.2, 337.6, 593.1, and 602.2 MPa, respectively. Based on the results of the present study, it was suggested that the calf vertebral trabecular bone could be modeled as a poroelastic material and its strain rate dependent material behavior could be predicted.

• Computers and Concrete
• /
• v.12 no.3
• /
• pp.303-318
• /
• 2013

The purpose of this study is to present different composed material models for reinforced concrete structures (RC). For this aim a nonlinear finite element analysis program is coded in MATLAB. This program contains several yield criteria and stress-strain relationships for compression and tension behavior of concrete. In this study, the well-known criteria, Drucker-Prager, von Mises, Mohr Coulomb, Tresca, and two new criteria, Hsieh-Ting-Chen and Bresler-Pister, are taken into account. It is concluded that the coded program, the new yield criteria, and the models considered can be effectively used in the nonlinear analysis of reinforced concrete beams.

• Steel and Composite Structures
• /
• v.21 no.5
• /
• pp.1121-1144
• /
• 2016

This paper presents a displacement-based finite element procedure for second-order distributed plasticity analysis of planar steel frames with semi-rigid beam-to-column connections under static loadings. A partially strain-hardening elastic-plastic beam-column element, which directly takes into account geometric nonlinearity, gradual yielding of material, and flexibility of semi-rigid connections, is proposed. The second-order effects and distributed plasticity are considered by dividing the member into several sub-elements and meshing the cross-section into several fibers. A new nonlinear solution procedure based on the combination of the Newton-Raphson equilibrium iterative algorithm and the constant work method for adjusting the incremental load factor is proposed for solving nonlinear equilibrium equations. The nonlinear inelastic behavior predicted by the proposed program compares well with previous studies. Coupling effects of three primary sources of nonlinearity, geometric imperfections, and residual stress are investigated and discussed in this paper.

• Steel and Composite Structures
• /
• v.28 no.3
• /
• pp.389-401
• /
• 2018

An isoparametric six-node triangular element is utilized for geometrically nonlinear analysis of functionally graded (FG) shells. To overcome the shear and membrane locking, the element is improved by using strain interpolation functions. The Total Lagrangian formulation is employed to include the large displacements and rotations. Finding the nonlinear behavior of FG shells via laminated modeling is also the goal. A power function is employed to formulate the variation of elastic modulus through the thickness of shells. The results are presented in two ways, including the general FGM formulation and the laminated modeling. The equilibrium path is obtained by using the Generalized Displacement Control Method. Some popular benchmarks, including hyperbolical shell structures are solved to declare the correctness and accuracy of proposed formulations.

• 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.

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

• Computers and Concrete
• /
• v.33 no.3
• /
• pp.275-284
• /
• 2024

This work presents a comprehensive investigation of buckling behavior of bidirectional functionally graded imperfect beams exposed to several thermal loading with general boundary conditions. The nonlinear governing equations are derived based on 2D shear deformation theory together with Von Karman strain-displacement relation. The beams are composed of two different materials. Its properties are porosity-dependent and are continuously distributed over the length and thickness of the beams following a defined law. The resulting equations are solved analytically in order to determine the thermal buckling characteristics of BDFG porous beams. The precision of the current solution and its accuracy have been proven by comparison with works previously published. Numerical examples are presented to explore the effects of the thermal loading, the elastic foundation parameters, the porosity distribution, the grading indexes and others factors on the nonlinear thermal buckling of bidirectional FG beam rested on elastic foundation.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.26 no.3
• /
• pp.92-102
• /
• 2022

Various numerical analysis models can be used to evaluate the behavior characteristics of tunnel facilities which are representative underground structures. In general, the Mohr-Coulomb model, which is most often used for numerical analysis, is an elastic-perfect plastic behavior model. And the deformation characteristics are the same during the load increase-load reduction phase. So there is a problem that the displacement may appear different from the field situation in the case of excavation analysis. In contrast, the HS-small strain stability model has a wide range of applications for each ground. And it is known that soil deformation characteristics can be analyzed according to field conditions by enabling input of initial elastic modulus and nonlinear curve parameter and so on. However, civil engineers are having difficulty using nonlinear models that can apply material nonlinear properties due to difficulties in estimating ground property coefficients. In this study, the necessity of rational model selection was reviewed by comparing the results of seismic performance evaluation using the Mohr-Coulomb model, which civil engineers generally apply for numerical analysis of tunnels, and the HS Small strain Stiffness model, which can consider ground nonlinearity.

• Korea-Australia Rheology Journal
• /
• v.15 no.3
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.3
• /
• pp.199-209
• /
• 2009

This paper deals with the three-dimensional finite element analysis of failure behavior of UHPFRC I-beam under monotonic load. Different from the constitutive law of normal and high strength concrete, an elastic-plastic fracture model that considers the tensile strain hardening is proposed to describe the material properties of UHPFRC. A multi-directional fixed crack criterion with tensile strain hardening is defined in the tensile region, and Drucker-Prager criterion with an associated flow rule is adopted in the compressive region. The influence of span, prestressing force and section on the behavior of UHPFRC I-beam are investigated. The comparison of the numerical results with the test results indicates a good agreement.