• Structural Engineering and Mechanics
• /
• v.12 no.5
• /
• pp.507-526
• /
• 2001

In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.3
• /
• pp.12-18
• /
• 2012

The Lattice Boltzmann Method has developed for solving the Boltzmann equation in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. When a immersed boundaries are sweeping the fixed fluid node, refilling the node information in a vicinity of fluid nodes is one of the important issues in Lattice Boltzmann Method. In this study, we propose a simple refill algorithm for the particle distribution function based on a proper velocity, density and strain rate to enhance accuracy and stability of the method. The refill scheme based on a asymptotic analysis of LBGK model has improved accuracy than interpolation schemes. The proposed scheme in this study is validated by the simulations of an impulsively started rotating circular cylinder to investigate adaptability for fluid-structure interaction (FSI) problem. This refill scheme has improved stability and accuracy especially at high Reynolds number region.

• Structural Engineering and Mechanics
• /
• v.32 no.5
• /
• pp.649-665
• /
• 2009

This paper reports part of a comprehensive research study conducted at the University of Queensland on the ability of CFRP web-bonded systems in strengthening an exterior beam-column joint subjected to monotonic loads. One 1/2.2 scaled plain and four CFRP repaired/retrofitted joints subjected to monotonic loads were analysed using the nonlinear finite-element program ANSYS and the results were calibrated against experiments. The ANSYS model was employed in order to account for tension stiffening in concrete after cracking and a modified version of the Hognestad's model was used to model the concrete compressive strength. The stress-strain properties of main steel bars were modelled using multilinear isotropic hardening model and the FRPs were modelled as anisotropic materials. A perfect bond was assumed as nodes were shared between adjacent elements irrespective of their type. Good agreement between the numerical predictions and the experimental observation of the failure mechanisms for all specimens were observed. Closeness of these results proved that the numerical analysis can be used by design engineers for the analysis of web-bonded FRP strengthened beam-column joints with confidence.

• Structural Engineering and Mechanics
• /
• v.66 no.5
• /
• pp.557-568
• /
• 2018

An investigation is made in the present work on the post-buckling and geometrically nonlinear behaviors of moderately thick perfect and imperfect rectangular plates made-up of functionally graded materials. Spectral collocation approach based on Legendre basis functions is developed to analyze the functionally graded plates while they are subjected to end-shortening strain. The material properties in this study are varied through the thickness according to the simple power law distribution. The fundamental equations for moderately thick rectangular plates are derived using first order shear deformation plate theory and taking into account both geometric nonlinearity and initial geometric imperfections. In the current study, the domain of interest is discretized with Legendre-Gauss-Lobatto nodes. The equilibrium equations will be obtained by discretizing the Von-Karman's equilibrium equations and also boundary conditions with finite Legendre basis functions that are substituted into the displacement fields. Due to effect of geometric nonlinearity, the final set of equilibrium equations is nonlinear and therefore the quadratic extrapolation technique is used to solve them. Since the number of equations in this approach will always be more than the number of unknown coefficients, the least squares technique will be used. Finally, the effects of boundary conditions, initial geometric imperfection and material properties are investigated and discussed to demonstrate the validity and capability of proposed method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

• Advances in materials Research
• /
• v.1 no.4
• /
• pp.245-268
• /
• 2012

Results of isothermal torsional oscillation tests are reported on melts of linear low density polyethylene and isotactic polypropylene. Prior to rheological tests, specimens were annealed at various temperatures ranging from $T_a$ = 180 to $310^{\circ}C$ for various amounts of time (from 30 to 120 min). Thermal treatment induced degradation of the melts and caused pronounced decreases in their molecular weights. With reference to the concept of transient networks, constitutive equations are developed for the viscoelastic response of polymer melts. A melt is treated as an equivalent network of strands bridged by junctions (entanglements and physical cross-links). The time-dependent response of the network is modelled as separation of active strands from and merging of dangling strands with temporary nodes. The stress-strain relations involve three adjustable parameters (the instantaneous shear modulus, the average activation energy for detachment of active strands, and the standard deviation of activation energies) that are determined by matching the dependencies of storage and loss moduli on frequency of oscillations. Good agreement is demonstrated between the experimental data and the results of numerical simulation. The study focuses on the effect of molecular weight of polymer melts on the material constants in the constitutive equations.

• Pediatric Infection and Vaccine
• /
• v.21 no.2
• /
• pp.139-143
• /
• 2014

Bacille Calmette-$Gu\acute{e}rin$ (BCG) lymphadenitis is the most common complication of BCG vaccination. It commonly occurs in infants aged <6 months involving ipsilateral axillary lymph nodes. We described BCG lymphadenitis in a 22-month-old boy presenting swelling of left supraclavicular lymph node that was confirmed by real-time polymerase chain reaction (PCR) and the multiplex PCR targeting the region of difference (RD).

• Nuclear Engineering and Technology
• /
• v.52 no.11
• /
• pp.2565-2571
• /
• 2020

Background: Understanding the behaviour of nuclear fuel claddings by conducting burst test on single cladding tube under simulated loss-of-coolant accident conditions and developing theoretical cum empirical predictive computer codes have been the focus of several investigations. The developed burst criterion (a) assumes symmetrical deformation of cladding tube in contrast to experimental observation (b) interpolates the properties of Zircaloy-4 cladding in mixed α+β phase (c) does not account for azimuthal temperature variations. In order to overcome all these drawbacks of burst criterion, it is reasoned that artificial intelligence technique may be a better option to predict the burst parameters. Methods: Artificial neural network models based on feedforward backpropagation algorithm with logsig transfer function are developed. Results: Neural network architecture of 2-4-4-3, that is model with two hidden layers having four nodes in each layer is found to be the most suitable. The mean, maximum, and minimum prediction errors for this optimised model are 0.82%, 19.62%, and 0.004%, respectively. Conclusion: The burst stress, burst temperature, and burst strain obtained from burst criterion have average deviation of 19%, 12%, and 53% respectively whereas the developed neural network model predicted these parameters with average deviation of 6%, 2%, and 8%, respectively.

• Korean Journal of Materials Research
• /
• v.19 no.12
• /
• pp.674-679
• /
• 2009

Yttria stabilized zirconia single crystals show plastic deformation at high temperatures by activating dislocations. The plastic deformation is highly dependent on crystallographic orientation. When the samples were deformed at different orientations, stress-strain curves changed by operating different slip systems. The strength of samples was also highly dependent on crystallographic orientation, i.e., samples without yield drop showed higher strength than that of samples exhibiting yield drop. The slip systems in the sample deformed along <112>, <111> and <001> agreed with the theoretical values of the plastic deformation, following Schmid's Law. Dislocations play a major role in the plastic deformation of this crystal. At the early stages of plastic deformation, all samples exhibited dislocation dipoles and, in the later stages, dislocation interactions occurred by forming nodes, tangles and networks. In this study, three different orientations, [11-2], [111] and [001] were employed to explain the plastic deformation behavior. A microstructural analysis was performed to elucidate the mechanism of the plastic behavior of this crystal.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.81-94
• /
• 2010

This paper deals with the fracture simulation of UHPFRC girder with the interface type model. Based on the existing numerical simulation of quasi-brittle fracture in normal strength concrete, constitutive modeling for UHPFRC I-girder has been improved by including a tensile hardening at the failure surface. The finite element formulation is based on a triangular unit, constructed from constant strain triangles, with nodes along its sides and neither at the vertex nor the center of the unit. Fracture is simulated through a hardening/softening fracture constitutive law in tension, a softening fracture constitutive law in shear as well as in compression at the boundary nodes, with the material within the triangular unit remaining linear elastic. LCP is used to formulate the path-dependent hardening-softening behavior in non-holonomic rate form and a mathematical programming algorithm is employed to solve the LCP. The piece-wise linear inelastic yielding-failure/failure surface is modeled with two compressive caps, two Mohr-Coulomb failure surfaces, a tensile yielding surface and a tensile failure surface. The comparison between test results and numerical results indicates this method effectively simulates the deformation and failure of specimen.