Dynamics of a Deployable Mesh Reflector of Satellite Antenna: Parallel Computation and Deployment Simulation1

Author(s):  
Pei Li ◽  
Cheng Liu ◽  
Qiang Tian ◽  
Haiyan Hu ◽  
Yanping Song

The finite-element approach of absolute nodal coordinate formulation (ANCF) is a possible way to simulate the deployment dynamics of a large-scale mesh reflector of satellite antenna. However, the large number of finite elements of ANCF significantly increases the dimension of the dynamic equations for the deployable mesh reflector and leads to a great challenge for the efficient dynamic simulation. A new parallel computation methodology is proposed to solve the differential algebraic equations for the mesh reflector multibody system. The mesh reflector system is first decomposed into several independent subsystems by cutting its joints or finite-element grids. Then, the Schur complement method is used to eliminate the internal generalized coordinates of each subsystem and the Lagrange multipliers for joint constraint equations associated with the internal variables. With an increase of the number of subsystems, the dimension of simultaneous linear equations generated in the numerical solution process will inevitably increase. By using the multilevel decomposition approach, the dimension of the simultaneous linear equations is further reduced. Two numerical examples are used to validate the efficiency and accuracy of the proposed parallel computation methodology. Finally, the dynamic simulation for a 500 s deployment process of a complex AstroMesh reflector with over 190,000 generalized coordinates is efficiently completed within 78 hrs.

Author(s):  
Yasuhito Takahashi ◽  
Koji Fujiwara ◽  
Takeshi Iwashita ◽  
Hiroshi Nakashima

Purpose This paper aims to propose a parallel-in-space-time finite-element method (FEM) for transient motor starting analyses. Although the domain decomposition method (DDM) is suitable for solving large-scale problems and the parallel-in-time (PinT) integration method such as Parareal and time domain parallel FEM (TDPFEM) is effective for problems with a large number of time steps, their parallel performances get saturated as the number of processes increases. To overcome the difficulty, the hybrid approach in which both the DDM and PinT integration methods are used is investigated in a highly parallel computing environment. Design/methodology/approach First, the parallel performances of the DDM, Parareal and TDPFEM were compared because the scalability of these methods in highly parallel computation has not been deeply discussed. Then, the combination of the DDM and Parareal was investigated as a parallel-in-space-time FEM. The effectiveness of the developed method was demonstrated in transient starting analyses of induction motors. Findings The combination of Parareal with the DDM can improve the parallel performance in the case where the parallel performance of the DDM, TDPFEM or Parareal is saturated in highly parallel computation. In the case where the number of unknowns is large and the number of available processes is limited, the use of DDM is the most effective from the standpoint of computational cost. Originality/value This paper newly develops the parallel-in-space-time FEM and demonstrates its effectiveness in nonlinear magnetoquasistatic field analyses of electric machines. This finding is significantly important because a new direction of parallel computing techniques and great potential for its further development are clarified.


1996 ◽  
Vol 86 (2) ◽  
pp. 379-388 ◽  
Author(s):  
H. Takenaka ◽  
M. Ohori ◽  
K. Koketsu ◽  
B. L. N. Kennett

Abstract The Aki-Larner method is one of the cheapest methods for synthetic seismograms in irregularly layered media. In this article, we propose a new approach for a two-dimensional SH problem, solved originally by Aki and Larner (1970). This new approach is not only based on the Rayleigh ansatz used in the original Aki-Larner method but also uses further information on wave fields, i.e., the propagation invariants. We reduce two coupled integral equations formulated in the original Aki-Larner method to a single integral equation. Applying the trapezoidal rule for numerical integration and collocation matching, this integral equation is discretized to yield a set of simultaneous linear equations. Throughout the derivation of these linear equations, we do not assume the periodicity of the interface, unlike the original Aki-Larner method. But the final solution in the space domain implicitly includes it due to use of the same discretization of the horizontal wavenumber as the discrete wavenumber technique for the inverse Fourier transform from the wavenumber domain to the space domain. The scheme presented in this article is more efficient than the original Aki-Larner method. The computation time and memory required for our scheme are nearly half and one-fourth of those for the original Aki-Larner method. We demonstrate that the band-reduction technique, approximation by considering only coupling between nearby wavenumbers, can accelerate the efficiency of our scheme, although it may degrade the accuracy.


2016 ◽  
Vol 35 (23) ◽  
pp. 1696-1711 ◽  
Author(s):  
Danilo S Victorazzo ◽  
Andre De Jesus

In this paper we extend Kollár and Pluzsik’s thin-walled anisotropic composite beam theory to include multiple cells with open branches and booms, and present a finite element formulation utilizing the stiffness matrix obtained from this theory. To recover the 4 × 4 compliance matrix of a beam containing N closed cells, we solve an asymmetric system of 2N + 4 linear equations four times with unitary section loads and extract influence coefficients from the calculated strains. Finally, we compare 4 × 4 stiffness matrices of a multicelled beam using this method against matrices obtained using the finite element method to demonstrate accuracy. Similarly to its originating theory, the effects of shear deformation and restrained warping are assumed negligible.


Author(s):  
H. Ren ◽  
W. D. Zhu

A spatial discretization and substructure method is developed to calculate the dynamic responses of one-dimensional systems, which consist of length-variant distributed-parameter components such as strings, rods, and beams, and lumped-parameter components such as point masses and rigid bodies. The dependent variable, such as the displacement, of a distributed-parameter component is decomposed into boundary-induced terms and internal terms. The boundary-induced terms are interpolated from the boundary motions, and the internal terms are approximated by an expansion of trial functions that satisfy the corresponding homogeneous boundary conditions. All the matching conditions at the interfaces of the components are satisfied, and the expansions of the dependent variables of the distributed-parameter components absolutely and uniformly converge. The spatial derivatives of the dependent variables, which are related to the internal forces/moments, such as the axial forces, bending moments, and shear forces, can be accurately calculated. Assembling the component equations and the geometric matching conditions that arise from the continuity relations leads to a system of differential algebraic equations (DAEs). When some matching conditions are linear algebraic equations, some generalized coordinates can be represented by others so that the number of the generalized coordinates can be reduced. The methodology is applied to moving elevator cable-car systems in Part II of this work.


2012 ◽  
Vol 23 (10) ◽  
pp. 1143-1160 ◽  
Author(s):  
Walid Khalil ◽  
Alain Mikolajczak ◽  
Céline Bouby ◽  
Tarak Ben Zineb

In this article, we propose a finite element numerical tool adapted to a Fe-based shape memory alloy structural analysis, based on a developed constitutive model that describes the effect of phase transformation, plastic sliding, and their interactions on the thermomechanical behavior. This model was derived from an assumed expression of the Gibbs free energy taking into account nonlinear interaction quantities related to inter- and intragranular incompatibilities as well as mechanical and chemical quantities. Two scalar internal variables were considered to describe the phase transformation and plastic sliding effects. The hysteretic and specific behavior patterns of Fe-based shape memory alloy during reverse transformation were studied by assuming a dissipation expression. The proposed model effectively describes the complex thermomechanical loading paths. The numerical tool derived from the implicit resolution of the nonlinear partial derivative constitutive equations was implemented into the Abaqus® finite element code via the User MATerial (UMAT) subroutine. After tests to verify the model for homogeneous and heterogeneous thermomechanical loadings, an example of Fe-based shape memory alloy application was studied, which corresponds to a tightening system made up of fishplates for crane rails. The results we obtained were compared to experimental ones.


Author(s):  
Christophe Geuzaine ◽  
Laurent Stainier ◽  
Francois Henrotte

In this article we propose a macroscopic model for ferromagnetic hysteresis that is well-suited for finite element implementation. The model is readily vectorial and relies on a consistent thermodynamic formulation. In particular, the stored magnetic energy and the dissipated energy are known at all times, and not solely after the completion of closed hysteresis loops as is usually the case. The obtained incremental formulation is variationally consistent, i.e., all internal variables follow from the minimization of a thermodynamic potential. This variational approach is directly inspired from the kinematic hardening theory of plasticity, which opens the door for novel energy-consistent coupled mechanical/electromagnetic models.


Sign in / Sign up

Export Citation Format

Share Document