• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.495-504
• /
• 2017

This paper proposes a new direct numerical integration algorithm for solving equation of motion in structural dynamics problems with nonlinear stiffness. The new implicit method's degree of accuracy is higher than that of existing methods due to the higher order of the acceleration. Two parameters are defined, leading to a new family of unconditionally stable methods, which helps to take greater time steps in integration and eliminate concerns about the duration of solving. The method developed can be utilized for a number of solid plane finite elements, examples of which are given to compare the proposed method with existing ones. The results indicate the superiority of the proposed method.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.107-125
• /
• 2013

A strain smoothing meshfree formulation with stabilized conforming nodal integration is presented for modeling the consolidation process in saturated porous media. In the present method, nodal strain smoothing is consistently introduced into the meshfree approximation of strain and pore pressure gradient variables associated with the saturated porous media. Meanwhile, in order to achieve a consistent numerical implementation, a smoothing approximation of the meshfree shape function within a nodal representative domain is also proposed in the stiffness construction. The resulting discrete system of equations is all expressed in smoothed nodal measures that are very efficient for numerical evaluation. Subsequently the space-time fully discrete equations are further established by the generalized trapezoidal rule for time integration. The effectiveness of the proposed meshfree consolidation analysis method is systematically illustrated by several benchmark problems.

• Kyungpook Mathematical Journal
• /
• v.64 no.2
• /
• pp.353-369
• /
• 2024

In this study, we investigate the unknown time-dependent heat source function in inverse problems. We consider three general nonlocal conditions; two classical boundary conditions and one nonlocal over-determination, condition, these genereate six different cases. The finite integration method (FIM), based on numerical integration, has been adapted to solve PDEs, and we use it to discretize the spatial domain; we use backward differences for the time variable. Since the inverse problem is ill-posed with instability, we apply regularization to reduce the instability. We use the first-order Tikhonov's regularization together with the minimization process to solve the inverse source problem. Test examples in all six cases are presented in order to illustrate the accuracy and stability of the numerical solutions.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.6 no.4
• /
• pp.147-158
• /
• 2003

The prediction of wind field is very important fact in the radioactive and chemical warfare. In spite of advanced numerical weather prediction modelling and computing technology, the high resolution prediction of wind field is limited by the very high integration costs. In this study we coupled the mesoscale numerical model and microscale diagnostic numerical model with minimized integration costs. This coupled model has not only the ability of prediction of high resolution wind field including complex building but also microscale pollutant diffusion fields. For military operation this system can help making a practical and cost-effective decision in a battle field.

• Journal of information and communication convergence engineering
• /
• v.12 no.1
• /
• pp.39-45
• /
• 2014

A numerical method for solving Hamiltonian equations is said to be symplectic if it preserves the symplectic structure associated with the equations. Various symplectic methods are widely used in many fields of science and technology. A symplectic method preserves an approximate Hamiltonian perturbed from the original Hamiltonian. It theoretically supports the effectiveness of symplectic methods for long-term integration. Although it is also related to long-term integration, numerical stability of symplectic methods have received little attention. In this paper, we consider explicit symplectic methods defined for Hamiltonian equations with Hamiltonians of the special form, and study their numerical stability using the harmonic oscillator as a test equation. We propose a new stability criterion and clarify the stability of some existing methods that are visually based on the criterion. We also derive a new method that is better than the existing methods with respect to a Courant-Friedrichs-Lewy condition for hyperbolic equations; this new method is tested through a numerical experiment with a nonlinear wave equation.

• Journal of Digital Convergence
• /
• v.19 no.12
• /
• pp.339-345
• /
• 2021

In this paper, estimation angle performance analysis of amplitude-comparison monopulse radar under additive noise effect is dealt with. When uncorrelated white noises are added to the squinted beams, the angle estimation performance is analyzed through the mean square error(MSE). The numerical integration-based mean square error result completely overlaps the Monte Carlo-based mean square error result, which corresponds to 99.8% of the Monte Carlo-based mean square error result. In addition, the mean square error analysis method based on numerical integration has a much faster operation time than the mean square error method based on Monte Carlo. the angle estimation performance of the amplitude comparison monopulse radar can be efficiently analyzed in various noise environments through the proposed numerical integration-based mean square error method.

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

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1131-1150
• /
• 2014

One of the issues in extending the range of applicable problems of real-time hybrid simulation is the computation speed of the simulator when large-scale computational models with a large number of DOF are used. In this study, functionality of real-time dynamic simulation of MDOF systems is achieved by creating a logic circuit that performs the step-by-step numerical time integration of the equations of motion of the system. The designed logic circuit can be implemented to an FPGA-based system; FPGA (Field Programmable Gate Array) allows large-scale parallel computing by implementing a number of arithmetic operators within the device. The operator splitting method is used as the numerical time integration scheme. The logic circuit consists of blocks of circuits that perform numerical arithmetic operations that appear in the integration scheme, including addition and multiplication of floating-point numbers, registers to store the intermediate data, and data busses connecting these elements to transmit various information including the floating-point numerical data among them. Case study on several types of linear and nonlinear MDOF system models shows that use of resource sharing in logic synthesis is crucial for effective application of FPGA to real-time dynamic simulation of structural response with time step interval of 1 ms.

• Steel and Composite Structures
• /
• v.20 no.3
• /
• pp.501-512
• /
• 2016

High-strength structural bolts have been utilized for beam-to-column connections in steel-framed structural buildings. Failure of these components may be caused by the bolt shank fracture or threads stripping-off, documented in the literature. Furthermore, these structural bolts are galvanized for corrosion resistance or quenched-and-tempered in the manufacturing process. This paper adopted the finite element simulation to demonstrate discrete mechanical performance for these bolts under tensile loading conditions, the coated and uncoated numerical model has been built up for two numerical integration methods: explicit and implicit. Experimental testing and numerical methods can fully approach the failure mechanism of these bolts and their ultimate load capacities. Comparison has also been conducted for two numerical integration methods, demonstrating that the explicit integration procedure is also suitable for solving quasi-static problems. Furthermore, by using precise bolt models in T-Stub, more accurately simulate the mechanical behavior of T-Stub, which will lay the foundation of the mechanical properties of steel bolted joints.

• Computers and Concrete
• /
• v.25 no.2
• /
• pp.181-192
• /
• 2020

Nowadays, steel fiber reinforced concretes (SFRCs) are widely used in practical applications. Significant experimental research has thus been carried out to determine the constitutive equations that represent the behavior of SFRCs under multiaxial loadings. However, numerical modelling of SFRCs via FEM has been challenging due to the complexities of the implementation of these constitutive equations. In this study, following the literature, a plasticity model is constructed for the behavior of SFRCs that involves the Willam-Warnke failure surface with the relevant evolution laws and a non-associated flow rule for determining the plastic deformations. For the precise (yet rapid) integration of the constitutive equations, an explicit substepping scheme consisting of yield intersection and drift correction algorithms is employed and thus implemented in ABAQUS via UMAT. The FEM model includes various material parameters that are determined from the experimental data. Three sets of parameters are used in the numerical simulations. While the first set is from the experiments that are conducted in this study on SFRC specimens with various contents of steel fibers, the other two sets are from the experiments reported in the literature. The response of SFRCs under multiaxial compression obtained from various numerical simulations are compared with the experimental data. The good agreement between numerical results and the experimental data indicates that not only the adopted plasticity model represents the behavior of SFRCs very well but also the implemented integration scheme can be employed in practical applications of SFRCs.