• Computers and Concrete
• /
• v.2 no.5
• /
• pp.345-366
• /
• 2005

This work presents an assessment of the computational performance of a vector-parallel implementation of probabilistic model for concrete cracking in 3D. This paper shows the continuing efforts towards code optimization as reported in earlier works Paz, et al. (2002a,b and 2003). The probabilistic crack approach is based on the direct Monte Carlo method. Cracking is accounted by means of 3D interface elements. This approach considers that all nonlinearities are restricted to interface elements modeling cracks. The heterogeneity governs the overall cracking behavior and related size effects on concrete fracture. Computational kernels in the implementation are the inexact Newton iterative driver to solve the non-linear problem and a preconditioned conjugate gradient (PCG) driver to solve linearized equations, using an element by element (EBE) strategy to compute matrix-vector products. In particular the paper analyzes code behavior using OpenMP directives in parallel vector processors (PVP), such as the CRAY SV1 and CRAY T94. The impact of the memory architecture on code performance, and also some strategies devised to circumvent this issue are addressed by numerical experiment.

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

• Coupled systems mechanics
• /
• v.3 no.4
• /
• pp.329-344
• /
• 2014

In this work, quantum molecular dynamics simulations (QMD) are preformed to study the hydrogen molecules in three types of carbon nanostructures, $C_{60}$ fullerene, (5,5) and (9,0) carbon nanotubes and graphene layers. Interactions between hydrogen and the nanostructures is of importance to understand hydrogen storage for the development of hydrogen economy. The QMD method overcomes the difficulties with empirical interatomic potentials to model the interaction among hydrogen and carbon atoms in the confined geometry. In QMD, the interatomic forces are calculated by solving the Schrodinger's equation with the density functional theory (DFT) formulation, and the positions of the atomic nucleus are calculated with the Newton's second law in accordance with the Born-Oppenheimer approximation. It is found that the number of hydrogen atoms that is less than 58 can be stored in the $C_{60}$ fullerene. With larger carbon fullerenes, more hydrogen may be stored. For hydrogen molecules passing though the fullerene, a particular orientation is required to obtain least energy barrier. For carbon nanotubes and graphene, adsorption may adhere hydrogen atoms to carbon atoms. In addition, hydrogen molecules can also be stored inside the nanotubes or between the adjacent layers in graphite, multi-layer graphene.

• Nuclear Engineering and Technology
• /
• v.45 no.4
• /
• pp.439-458
• /
• 2013

REX-10 is a fully-passive small modular reactor in which the coolant flow is driven by natural circulation, the RCS is pressurized by a steam-gas pressurizer, and the decay heat is removed by the PRHRS. To confirm design decisions and analyze the transient responses of an integral PWR such as REX-10, a thermal-hydraulic system code named TAPINS (Thermal-hydraulic Analysis Program for INtegral reactor System) is developed in this study. Based on a one-dimensional four-equation drift-flux model, TAPINS incorporates mathematical models for the core, the helical-coil steam generator, and the steam-gas pressurizer. The system of difference equations derived from the semi-implicit finite-difference scheme is numerically solved by the Newton Block Gauss Seidel (NBGS) method. TAPINS is characterized by applicability to transients with non-equilibrium effects, better prediction of the transient behavior of a pressurizer containing non-condensable gas, and code assessment by using the experimental data from the autonomous integral effect tests in the RTF (REX-10 Test Facility). Details on the hydrodynamic models as well as a part of validation results that reveal the features of TAPINS are presented in this paper.

• Tribology and Lubricants
• /
• v.34 no.6
• /
• pp.325-331
• /
• 2018

Automobiles and systems requiring high gear ratios and high power densities generally use worm gears. In particular, as worm gears have a small volume and self-locking function, home appliances such as refrigerators and washers consist of worm gears. We can classify worm gears into cylindrical worms and rectangular worms. According to the AGMA standard, there are four types of cylindrical worms, ZA, ZN, ZK and ZI, depending on the machining of the worm shaft. It is preferable to use a ZI-type worm shaft, which is a combination of a worm wheel having an involute helical tooth surface and a conjugate tooth surface. However, in many cases, industries mostly use ZK, ZN, and ZA worm shafts because of the ease of processing. This paper presents numerical approaches to produce ZI and ZA worm surfaces and worm wheel. For the analysis of the transmission error of a worm gear system, this study (1) generates surface profile functions of ZI profile worm gear and worm shaft based on the common rack theory, (2) adopts the Newton-Raphson method for the analysis of the gear surface contact condition, and (3) presents and compares the corresponding transmission errors of ZI and ZA worm gears.

• Steel and Composite Structures
• /
• v.32 no.5
• /
• pp.657-670
• /
• 2019

This paper proposes a one-dimensional fiber beam element model taking account of materially non-linear behavior, benefiting the highly efficient elastic-plastic analysis of girders with shear-lag effects. Based on the displacement-based fiber beam-column element, two additional degrees of freedom (DOFs) are added into the proposed model to consider the shear-lag warping deformations of the slabs. The new finite element (FE) formulations of the tangent stiffness matrix and resisting force vector are deduced with the variational principle of the minimum potential energy. Then the proposed element is implemented in the OpenSees computational framework as a newly developed element, and the full Newton iteration method is adopted for an iterative solution. The typical materially non-linear behaviors, including the cracking and crushing of concrete, as well as the plasticity of the reinforcement and steel girder, are all considered in the model. The proposed model is applied to several test cases under elastic or plastic loading states and compared with the solutions of theoretical models, tests, and shell/solid refined FE models. The results of these comparisons indicate the accuracy and applicability of the proposed model for the analysis of both concrete box girders and steel-concrete composite girders, under either elastic or plastic states.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.14 no.1
• /
• pp.147-154
• /
• 2019

In this research, a molecular dynamics system was designed and developed to calculate trajectories of molecules in nozzles and turbin blades. The Lennard-Jones potential model was used to approximate the interaction between a pair of molecules and the Verlet integration is used as a numerical method to integrate Newton's equations of motion. To compute Lennard-Jones potential functions, for the number of molecules N, the computation complexity $O(N^2)$ for interactions of all pairs of molecules is reduced to O(N) by using cutoff radius $r_c$. This was implemented to save CPU times.

• Steel and Composite Structures
• /
• v.30 no.6
• /
• pp.517-534
• /
• 2019

In this paper, the effects of inevitable out-of-plane defects on the postbuckling behavior of single-layered graphene sheets (SLGSs) under in-plane loadings are investigated based on nonlocal first order shear deformation theory (FSDT) and von-Karman nonlinear model. A generic imperfection function, which takes the form of the products of hyperbolic and trigonometric functions, is employed to model out-of-plane defects as initial geometrical imperfections of SLGSs. Nonlinear equilibrium equations are derived from the principle of virtual work and variational formulation. The postbuckling equilibrium paths of imperfect graphene sheets (GSs) are presented by solving the governing equations via isogeometric analysis (IGA) and Newton-Raphson iterative method. Finally, the sensitivity of the postbuckling behavior of GS to shape, amplitude, extension on the surface, and location of initial imperfection is studied. Results showed that the small scale and initial imperfection effects on the postbuckling behavior of defective SLGS are important and cannot be ignored.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.15 no.7
• /
• pp.2339-2355
• /
• 2021

In this paper, a joint probabilistic data association algorithm based on loss function (LJPDA) is proposed so that the computation load and accuracy of the multi-target tracking algorithm can be guaranteed simultaneously. Firstly, data association is divided in to three cases based on the relationship among validation gates and the number of measurements in the overlapping area for validation gates. Also the contribution coefficient is employed for evaluating the contribution of a measurement to a target, and the loss function, which reflects the cost of the new proposed data association algorithm, is defined. Moreover, the equation set of optimal contribution coefficient is given by minimizing the loss function, and the optimal contribution coefficient can be attained by using the Newton-Raphson method. In this way, the weighted value of each target can be achieved, and the data association among measurements and tracks can be realized. Finally, we compare performances of LJPDA proposed and joint probabilistic data association (JPDA) algorithm via numerical simulations, and much attention is paid on real-time performance and estimation error. Theoretical analysis and experimental results reveal that the LJPDA algorithm proposed exhibits small estimation error and low computation complexity.

• Journal of IKEEE
• /
• v.25 no.4
• /
• pp.650-663
• /
• 2021

Subsurface topology estimation is an important factor in the geophysical survey. Electrical impedance tomography is one of the popular methods used for subsurface imaging. The EIT inverse problem is highly nonlinear and ill-posed; therefore, reconstructed conductivity distribution suffers from low spatial resolution. The subsurface region can be approximated as piece-wise separate regions with constant conductivity in each region; therefore, the conductivity estimation problem is transformed to estimate the shape and location of the layer boundary interface. Each layer interface boundary is treated as an open boundary that is described using front points. The subsurface domain contains multi-layers with very complex configurations, and, in such situations, conventional methods such as the modified Newton Raphson method fail to provide the desired solution. Therefore, in this work, we have implemented a 7-layer artificial neural network (ANN) as an inverse problem algorithm to estimate the front points that describe the multi-layer interface boundaries. An ANN model consisting of input, output, and five fully connected hidden layers are trained for interlayer boundary reconstruction using training data that consists of pairs of voltage measurements of the subsurface domain with three-layer configuration and the corresponding front points of interface boundaries. The results from the proposed ANN model are compared with the gravitational search algorithm (GSA) for interlayer boundary estimation, and the results show that ANN is successful in estimating the layer boundaries with good accuracy.