A symmetric substructuring method for analyzing the natural frequencies of conical origami structures

Conical origami structures are characterized by their substantial out-of-plane stiffness and energy-absorptioncapacity.Previous investigations have commonly focused on the static characteristics of these lightweight struc-tures.However,the efficient analysis of the natural vibrations of these structures is pivotal for designing conicalorigami structures with programmable stiffness and mass.In this paper,we propose a novel method to analyzethe natural vibrations of such structures by combining a symmetric substructuring method(SSM)and a gener-alized eigenvalue analysis.SSM exploits the inherent symmetry of the structure to decompose it into a finiteset of repetitive substructures.In doing so,we reduce the dimensions of matrices and improve computationalefficiency by adopting the stiffness and mass matrices of the substructures in the generalized eigenvalue analysis.Finite element simulations of pin-jointed models are used to validate the computational results of the proposedapproach.Moreover,the parametric analysis of the structures demonstrates the influences of the number of seg-ments along the circumference and the radius of the cone on the structural mass and natural frequencies of thestructures.Furthermore,we present a comparison between six-fold and four-fold conical origami structures anddiscuss the influence of various geometric parameters on their natural frequencies.This study provides a strategyfor efficiently analyzing the natural vibration of symmetric origami structures and has the potential to contributeto the efficient design and customization of origami metastructures with programmable stiffness.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

Weakly Singular Symmetric Galerkin Boundary Element Method for Fracture Analysis of Three-Dimensional Structures Considering Rotational Inertia and Gravitational Forces

The Symmetric Galerkin Boundary Element Method is advantageous for the linear elastic fracture and crackgrowth analysis of solid structures,because only boundary and crack-surface elements are needed.However,for engineering structures subjected to body forces such as rotational inertia and gravitational loads,additional domain integral terms in the Galerkin boundary integral equation will necessitate meshing of the interior of the domain.In this study,weakly-singular SGBEM for fracture analysis of three-dimensional structures considering rotational inertia and gravitational forces are developed.By using divergence theorem or alternatively the radial integration method,the domain integral terms caused by body forces are transformed into boundary integrals.And due to the weak singularity of the formulated boundary integral equations,a simple Gauss-Legendre quadrature with a few integral points is sufficient for numerically evaluating the SGBEM equations.Some numerical examples are presented to verify this approach and results are compared with benchmark solutions.

THE PARALLEL MULTISPLITTING METHOD FOR CONSISTENT SYMMETRIC POSITIVE(SEMI-)DEFINITE SYSTEMS

Main resultsTheorem 1 Let A be an n×n symmetric positive semidefinite matrix and let

M-TIMES SECANT-LIKE MULTI-PROJCTION METHOD FOR SPARSE MINIMIZATION PROBLEM

In this paper, we present m time secant like multi projection algorithm for sparse unconstrained minimization problem. We prove this method are all q superlinearly convergent to the solution about m≥1 . At last, we from some numerical results, discuss how to choose the number m to determine the approximating matrix properly in practical use.

An LQP-Based Symmetric Alternating Direction Method of Multipliers with Larger Step Sizes

Symmetric alternating directionmethod of multipliers(ADMM)is an efficient method for solving a class of separable convex optimization problems.This method updates the Lagrange multiplier twice with appropriate step sizes at each iteration.However,such step sizes were conservatively shrunk to guarantee the convergence in recent studies.In this paper,we are devoted to seeking larger step sizes whenever possible.The logarithmic-quadratic proximal(LQP)terms are applied to regularize the symmetric ADMM subproblems,allowing the constrained subproblems to then be converted to easier unconstrained ones.Theoretically,we prove the global convergence of such LQP-based symmetric ADMM by specifying a larger step size domain.Moreover,the numerical results on a traffic equilibrium problem are reported to demonstrate the advantage of the method with larger step sizes.

Symmetric and symplectic methods for gyrocenter dynamics in time-independent magnetic fields

We apply a second-order symmetric Runge–Kutta method and a second-order symplectic Runge–Kutta method directly to the gyrocenter dynamics which can be expressed as a noncanonical Hamiltonian system.The numerical simulation results show the overwhelming superiorities of the two methods over a higher order nonsymmetric nonsymplectic Runge–Kutta method in long-term numerical accuracy and near energy conservation.Furthermore,they are much faster than the midpoint rule applied to the canonicalized system to reach given precision.

An implicit upwind parabolized Navier-Stokes code for chemically nonequilibrium flows

The previously developed single-sweep parabolized Navier-Stokes （SSPNS） space marching code for ideal gas flows has been extended to compute chemically nonequilibrium flows. In the code, the strongly coupled set of gas dynamics, species conservation, and turbulence equations is integrated with the implicit lower-upper symmetric GaussSeidel （LU-SGS） method in the streamwise direction in a space marching manner. The AUSMPW＋ scheme is used to calculate the inviscid fluxes in the crossflow direction, while the conventional central scheme for the viscous fluxes. The k-g two-equation turbulence model is used. The revised SSPNS code is validated by computing the Burrows-Kurkov non-premixed H2/air supersonic combustion flows, premixed H2/air hypersonic combustion flows in a three-dimensional duct with a 15° compression ramp, as well as the hypersonic laminar chemically nonequilibrium air flows around two 10° half-angle cones. The results of these calculations are in good agreement with those of experiments, NASA UPS or Prabhu＇s PNS codes. It can be concluded that the SSPNS code is highly efficient for steady supersonic/ hypersonic chemically reaction flows when there is no large streamwise separation.

Application of SSOR-PCGM with improved iteration format in FEM simulation of massive concrete

In this study, for the purpose of improving the efficiency and accuracy of numerical simulation of massive concrete, the symmetric successive over relaxation-preconditioned conjugate gradient method （SSOR-PCGM） with an improved iteration format was derived and applied to solution of large sparse symmetric positive definite linear equations in the computational process of the finite element analysis. A three-dimensional simulation program for massive concrete was developed based on SSOR-PCGM with an improved iteration format. Then, the programs based on the direct method and SSOR-PCGM with an improved iteration format were used for computation of the Guandi roller compacted concrete （RCC） gravity dam and an elastic cube under free expansion. The comparison and analysis of the computational results show that SSOR-PCGM with the improved iteration format occupies much less physical memory and can solve larger-scale problems with much less computing time and flexible control of accuracy.

HIGH-ORDER SYMPLECTIC AND SYMMETRIC COMPOSITION METHODS FOR MULTI-FREQUENCY AND MULTI-DIMENSIONAL OSCILLATORY HAMILTONIAN SYSTEMS

The multi-frequency and multi-dimensional adapted Runge-Kutta^NystrSm （ARKN） integrators, and multi-frequency and multi-dimensional extended Runge-Kutta-NystrSm （ERKN） integrators have been developed to efficiently solve multi-frequency oscillatory Hamiltonian systems. The aim of this paper is to analyze and derive high-order sym- plectic and symmetric composition methods based on the ARKN integrators and ERKN integrators. We first consider the symplecticity conditions for the multi-frequency and multi-dimensional ARKN integrators. We then analyze the symplecticity of the adjoint in- tegrators of the multi-frequency and multi^dimensional symplectic ARKN integrators and ERKN integrators, respectively. On the basis of the theoretical analysis and by using the idea of composition methods, we derive and propose four new high-order symplectic and symmetric methods for the multi-frequency oscillatory Hamiltonian systems. The numer- ical results accompanied in this paper quantitatively show the advantage and efficiency of the proposed high-order symplectic and symmetric methods.

Fault diagnosis and protection of motors in coal shearer

Symmetrical components method is employed in analysis of the characteristic motor faults.Motor protection method is put forward based on detecting positive sequence,negative sequence and zero sequence current.And problems of lack of motor overload capacity in existing mining motor protection system,impact of dynamic current on stage and definite-time delay operation and inaccuracy of criterion phase failure protection are analyzed.The unbalanced faults protection and inverse-time overload protection,which can make protection time change with the current movement,are proposed.The above problems can be solved,and the reliability and intelligent of coal shearer are improved.

A Symmetric Inertial Alternating Direction Method of Multipliers for Elliptic Equation Constrained Optimization Problem

A new algorithm,called symmetric inertial alternating direction method of multipliers(SIADMM),is designed for separable convex optimization problems with linear constraints in this paper.The convergence rate of the SIADMM is proved to be O(1/√k).Two kinds of elliptic equation constrained optimization problems,the un-constrained cases as well as the box-constrained cases of the distributed control and the Robin boundary control,are analyzed theoretically and solved numerically.First,the existence and uniqueness of the solutions to these problems are proved.Second,these continuous optimization problems are transformed into discrete optimization problems by thefinite element method,and then the discrete optimization problems are solved by the proposed SIADMM.Numerical experiments with different problems are investigated to demonstrate the efficiency of the SIADMM.And the numerical per-formance of the SIADMM is better than the performance of the ADMM.Moreover,the numerical results show that the convergence rate of the SIADMM tends to be faster than O(1/√k)in calculation process.

Electric Power Quality Analysis 6-10/0.4 kV Distribution Networks

The paper reports quality analysis and evaluation at 6 - 10/0.4 kV low-voltage distribution grids in Uzbekistan. Power quality frequently does not correspond to the rated value which is largely due to unbalanced phase loading in grids and which also results in increased power loss. The study of the asymmetrical operating modes of the rural distribution networks of 0.4 kV was conducted in three steps: measurement, calculations and analysis of relevant data;providing practical guidelines and finally, implementing instruments to normalize grid operation. Measuring was conducted using certified instrumentation analyzer “MALIKA” designed by authors. The study and analysis of additional power losses as the function of indicators of asymmetrical features of voltage and current in operating 0.4 kV grids reveals that, quality of electric power at grids under investigation, merely does not meet the requirements of the Interstate Standard.

Consistency evaluation of the current distribution of short-circuit tests in AT feeding systems for high-speed railways

To improve the accuracy of fault location system, several short-circuit tests need to be conducted before being brought into service in autotransformer （AT） feeding systems for high-speed railways in China. However, no systematic algorithm yet exists to evaluate the consistency of the current distribution of short-circuit tests. A methodology is proposed in this paper to address this problem. Based on Kirchhoff’s current law and the generalized method of symmetrical components, the current deviations of the AT feeding systems are analysed and then normalized with the short-circuit current as they vary greatly with systems and short-circuit sites. It is also found that the short-circuit current varies with the calculation methods, and its unbiased standard deviation also reflects the consistency of the short-circuit test. The mean and maximum of the current deviations, as well as the unbiased standard deviation of the short-circuit current, show the consistency of the short-circuit test from different aspects,although the last two items are highly relevant. Therefore, a unified evaluation index is defined as the sum of the three items, and then applied in two case studies to test its performance. The results show that, the proposed index canclearly distinguish the consistency of the short-circuit tests and may be used to sort the short-circuit tests for fault location systems. Besides, some short-circuit tests may have very poor consistency indices, and thus are not applicable to the tuning of fault location systems. In the authors’ opinion, the determination of the threshold of the proposed index needs further investigation.

Kernel Function-Based Primal-Dual Interior-Point Methods for Symmetric Cones Optimization

In this paper, we present a large-update primal-dual interior-point method for symmetric cone optimization（SCO） based on a new kernel function, which determines both search directions and the proximity measure between the iterate and the center path. The kernel function is neither a self-regular function nor the usual logarithmic kernel function. Besides, by using Euclidean Jordan algebraic techniques, we achieve the favorable iteration complexity O（ √r（1/2）（log r）^2 log（r/ ε））, which is as good as the convex quadratic semi-definite optimization analogue.

A MATRIX-FREE IMPLICIT METHOD FOR INCOMPRESSIBLE FLOW AROUND 3D COMPLEX UNDERWATER BODIES

A fast, matrix-free implicit method based on the Lower-Upper Symmetric Gauss-Seidel （LU-SGS） method was developed to solve the three-dimensional incompressible Reynolds-averaged Navier-Stokes equations on multi-block structured grids. The method was applied to the simulations of a variety of flowfields around 3D complex underwater bodies with different appendages. The numerical procedure and results show that the method is efficient, reliable and robust for steady viscous flow simulations.

Optimal L_∞ Estimates for Galerkin Methods for Nonlinear Singular Two-point Boundary Value Problems

In this paper, we apply the symmetric Galerkin methods to the numerical solutions of a kind of singular linear two-point boundary value problems. We estimate the error in the maximum norm. For the sake of obtaining full superconvergence uniformly at all nodal points, we introduce local mesh refinements. Then we extend these results to a class of nonlinear problems. Finally, we present some numerical results which confirm our theoretical conclusions.

TWO NOVEL CLASSES OF ARBITRARY HIGH-ORDER STRUCTURE-PRESERVING ALGORITHMS FOR CANONICAL HAMILTONIAN SYSTEMS

In this paper,we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems.The one class is the symplectic scheme,which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method,respectively.Each member in these schemes is symplectic for any fixed parameter.A more general form of generating functions is introduced,which generalizes the three classical generating functions that are widely used to construct symplectic algorithms.The other class is a novel family of energy and quadratic invariants preserving schemes,which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step.The existence of the solutions of these schemes is verified.Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.