THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.      

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.     

Convergence analysis of a finite element method based on different moduli in tension and compression

When analyzing materials that exhibit different mechanical behaviors in tension and compression, an iterative approach is required due to material nonlinearities. Because of this iterative strategy, numerical instabilities may occur in the computational procedure. In this paper, we analyze the reason why iterative computation sometimes does not converge. We also present a method to accelerate convergence. This method is the introduction of a new pattern of shear modulus that was strictly derived according to the constitutive model based on the bimodular elasticity theory presented by Ambartsumyan. We test this procedure with a numerical example concerning a plane stress problem. Results obtained from this example show that the proposed method reduces the cost of computation and accelerates the convergence of the solution.     

Stiffness properties of particulate composites containing debonded particles

This paper considers the problem of determining the nonlinear bimodular stiffness properties, i.e., the tensile and compressive Young’s moduli and Poisson’s ratios, and the shear modulus, of particulate composite materials with particle–matrix interfacial debonding. It treats the general case in which some of the particles are debonded while the others remain intact. The Mori–Tanaka approach is extended to formulate the method of solution for the present problem. The resulting auxiliary problem of a single debonded particle in an infinite matrix subjected to a remote stress equal to the average matrix stress, for which Eshelby’s solution does not exist, is solved by the finite element method accounting for the particle–matrix separation and contact at the debonded particle–matrix interface. Because of the nonlinear nature of the problem, an iterative process is employed in calculating the stiffness properties. The predicted stiffness properties are compared to the exact solutions of the stiffness properties of particulate composites with body-centered cubic packing arrangement.     

Parametric variational principle and quadratic programming method for van der Waals force simulation of parallel and cross nanotubes

A parametric variational principle for van der Waals force simulation between any two adjacent nonbonded atoms and the corresponding improved quadratic programming method for numerical simulation of mechanical behaviors of carbon nanotubes are developed. Carbon nanotubes are modeled and computed based on molecular structural mechanics model. van der Waals force is simulated by the network of bars (called bar network) with a special nonlinear mechanical constitutive law (called generalized parametric constitutive law) in the finite element analysis. Compared with conventional numerical methods, the proposed method does not depend on displacement and stress iteration, but on the base exchanges in the solution of a standard quadratic programming problem. Thus, the model and method developed present very good convergence behavior in computation and provide accurate predictions of the mechanical behaviors and displacement distributions in the nanotubes. Numerical results demonstrate the validity and the efficiency of the proposed method.     

参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

蠕变－弹塑性－损伤耦合响应的有限元分析

本文处理各向异性非线性材料的蠕变。弹塑性－损伤耦合响应的数值计算。建议了一个计算应力的三级向后欧拉积分算法。导出了一个利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ迭代的一般的直接应力返回映射算法。同时求解应力向量和蠕变、塑性、损伤的内状态变量。也导出了用于全局Ｎｅｗｔｏｎ－Ｒａｐｈ－ｓｏｎ迭代过程的一致性切线矩阵公式。给出的数值例题结果表明所提出的算法和公式在模拟耦合本构行为上的能力和可靠性。     

蚁群算法求解二维拉压不同模量反问题

利用光滑函数技术对二维拉压不同模量本构关系进行光滑化处理,采用初应力方法求解二维拉压不同模量正问题的有限元方程。在此基础上,建立了基于连续域蚁群算法的二维拉压不同模量反问题的数值求解模型,考虑了区域非均质的影响,实现了对拉压弹性模量和泊松比的单一/组合识别。通过两个数值算例,对所提算法进行了数值验证,分别探讨了蚁群算法相关参数、测点分布和数据噪音等对识别结果的影响。数值验证表明,所提算法可有效地求解二维拉压不同模量反问题,并具有较好的计算精度。     

Auxiliary principle and three-step iterative algorithms for generalized set-valued strongly nonlinear mixed variational-like inequalities

An auxiliary principle technique to study a class of generalized set-valued strongly nonlinear mixed variational-like inequalities is extended. The existence and uniqueness of the solution of the auxiliary problem for the generalized set-valued strongly nonlinear mixed variational-like inequalities are proved, a novel and innovative three-step iterative algorithm to compute approximate solution is constructed, and the existence of the solution of the generalized set-valued strongly nonlinear mixed variational-like inequality is shown using the auxiliary principle technique. The convergence of three-step iterative sequences generated by the algorithm is also proved.     

Combination,principal parametric and internal resonances of an accelerating beam under two frequency parametric excitation

This study analyses the nonlinear transverse vibration of an axially moving beam subject to two frequency excitation. Focus has been made on simultaneous resonant cases i.e. principal parametric resonance of first mode and combination parametric resonance of additive type involving first two modes in presence of internal resonance. By adopting the direct method of multiple scales, the governing nonlinear integro-partial differential equation for transverse motion is reduced to a set of nonlinear first order ordinary partial differential equations which are solved either by means of continuation algorithm or via direct time integration. Specifically, the frequency response plots and amplitude curves, their stability and bifurcation are obtained using continuation algorithm. Numerical results reveal the rich and interesting nonlinear phenomena that have not been presented in the existent literature on the nonlinear dynamics of axially moving systems.     

正交各向异性弹塑性摩擦接触问题的数值求解

采用正交各向异性摩擦定律对三维弹塑性摩擦接触问题进行分析,基于参变量变分原理,经过有限元离散,将问题化为线性互补问题模型,之后给出一个求解互补问题的非内点光滑化算法．对三维接触问题,滑动方向的确定一直是个难点,为此,该文采用作者提出的组合规划法和迭代法对各向异性摩擦本构模型进行分析,数值结果说明了模型与算法的正确性。     

Widening basins of attraction of optimal iterative methods

In this work, we analyze the dynamical behavior on quadratic polynomials of a class of derivative-free optimal parametric iterative methods, designed by Khattri and Steihaug. By using their parameter as an accelerator, we develop different methods with memory of orders three, six and twelve, without adding new functional evaluations. Then a dynamical approach is made, comparing each of the proposed methods with the original ones without memory, with the following empiric conclusion: Basins of attraction of iterative schemes with memory are wider and the behavior is more stable. This has been numerically checked by estimating the solution of a practical problem, as the friction factor of a pipe and also of other nonlinear academic problems.     

Numerical simulation of contact problems in vane machinery by a parametric quadratic programming method

Interference fits are widely used for connecting impeller and shaft assembly that are forced together slowly by pressing. The interference fit design ensures stable balance behavior and allows for positive contact between the impeller and shaft assembly throughout the range of operating speeds. In addition to maintaining radial contact, sufficient net radial interface pressure must remain in order to transmit torque when the rotational speed is very high. Therefore, the interference fit between the impeller and the shaft assembly is one of the most important factors influencing the performance of the turbo unit in the design of turbocharger compressor. A suitable fit tolerance needs to be considered in the structural design. A locomotive-type turbocharger compressor with 24 blades under combined centrifugal and interference fit loading is used for the analysis. The finite-element (FE) parametric quadratic programming (PQP) method developed based on the parametric variational principle (PVP) is used for the analysis of the stress distribution in the three-dimensional (3D) contact problem of impeller. The advantages of the parametric programming method compared with conventional approaches are that the penalty factors can be canceled and that solutions can be obtained directly without tedious iterative procedures such as the general incremental iterative method. To save time in the computation, a~multi-substructure technique is adopted for structural modeling. This not only simplifies the calculation, but also provides a convenient service for process computer-aided design (CAD) by means of FE simulation. The effects of the fit tolerance, coefficient of friction and rotational speed (centrifugal force), wall thickness of the shaft sleeve and the contact stress on the interference-fitting surfaces are studied in detail in the numerical computation. It is found that a nonuniform initial amount of interference in the structural design avoids the relative displacement generated and ensures uniformity of the contact stress. To assure quality of press-fitting, the amount of interference between the shaft sleeve and shaft should be strictly controlled to avoid the rapid increase of the contact stress. The numerical results demonstrate the high accuracy and good convergence of the algorithm presented here, which provides an effective approach that achieves more-reliable interference-fitted connections and more-precise assembly accuracy with lower manufacturing cost in the structural design.     

Nonlinear finite element analysis of inflatable beams made from orthotropic woven fabric

This paper was devoted to the three-dimensional nonlinear finite element analysis of inflatable beams. The beams under consideration are made of modern textile materials and can be used as a load-bearing beams or arches when inflated. A 3D Timoshenko beam with a homogeneous orthotropic woven fabric (OWF) was proposed. The model took into account the geometric nonlinearities and the follower force resulting from the inflation pressure. The use was made of the usual total Lagrangian form of the virtual work principle to perform the nonlinear equilibrium equations which were discretized by the finite element method. Two kinds of solutions were then investigated: finite elements solutions for linearized problems which were obtained by the means of the linearization around the prestressed reference configuration of the nonlinear equations and nonlinear finite element solutions which were performed by the use of an optimization algorithm based on the Quasi-Newton method. As an example, the bending problem of a cantilever inflated beam under concentrated load was considered and the deflection results improve the existing theoretical models. As these beams are made from fabric, the beam models were validated through their comparison with a 3D thin-shell finite element model. The influence of the material effective properties and the inflation pressure on the beam response was also investigated through a parametric study. The finite elements solutions for linearized problems were found to be close to the theoretical results existing in the literature. On the other hand, the results for the nonlinear finite element model were shown to be close to the results for the linearized finite elements model in the case of high mechanical properties and the nonlinear finite element model was used to improve the linearized model when the mechanical properties of the fabric are low.     

Coulomb traction on a penny-shaped crack in a three dimensional piezoelectric body

The axisymmetric problem of a penny-shaped crack embedded in an infinite three-dimensional (3D) piezoelectric body is considered. A general formulation of Coulomb traction on the crack surfaces can be obtained based on thermodynamical considerations of electromechanical systems. Three-dimensional electroelastic solutions are derived by the classical complex potential theory when Coulomb traction is taken into account and the poling direction of piezoelectric body is perpendicular to the crack surfaces. Numerical results show that the magnitude of Coulomb tractions can be large, especially when a large electric field in connection with a small mechanical load is applied. Unlike the traditional traction-free crack model, Coulomb tractions induced by an applied electric field influence the Mode I stress intensity factor for a penny-shaped crack in 3D piezoelectric body. Moreover, compared to the current model, the traditional traction-free crack model always overestimates the effect of the applied electric load on the field intensity factors and energy release rates, which has consequences for 3D piezoelectric fracture mechanics.     

Karman型正交异性矩形薄板非线性弯曲的统一求解方法

本文对Karman型四边支承正交异性薄板在5种不同边界条件下的几何非线性弯曲进行了统一分析。所设的位移函数均为梁振动函数。它们精确地满足边界条件，利用Galerkin方法和位移函数的正交属性，转换控制方程为非线性代数方程。用“稳定化双共轭梯度法”求解稀疏矩阵线性方程组以及“可调节参数的修正迭代法”求解非线性代数方程组，最后给出了相应的数值结果。     

Effects of transverse shear on the nonlinear bending of rectangular plates laminated of bimodular composite materials

This paper investigates the application of Dynamic-Relaxation (DR) method to the problems of nonlinear bending of rectangular plates laminated of bimodular composite materials. The classical lamination theory and a shear deformation theory of layered composite plates, taking account of large rotations (in the von Karman sense) are employed separately to analyze the subject. It has been found here that the estimation of the fictitious densities which control the convergence and numerical stability of nonlinear DR solution considering transverse shear effect still needs to be further investigated. In this paper, a procedure to calculate fictitious densities has been presented; hence the numerical stability of this topic has been ensured. In this paper the main steps of solving the nonlinear bending of bimodular composite laminates by means of DR method are outlined. The numerical results are given for simply supported, two-layer cross-ply rectangular plates made of mildly bimodular material (Boron-Epoxy (B-E)) and highly bimodular materials (Aramid-Rubber (A-R) and Polyester-Rubber (P-R)) under sinusoidally distributed and uniformly distributed transverse loads. The results obtained have been compared with linear results and those obtained for laminates fabricated from conventional composite materials, the elastic moduli of which are identical with the tensile moduli of the bimodular materials. In addition, the effect of transverse shear deformation on the nondimensionalized center deflection has been studied.The main contents of this paper were presented at the International Symposium of Composite Materials and Structures (June 1986, Beijing).The authors thank Prof. Zhou Li for his guidance.     

Complicated nonlinear responses of a simply supported FGM rectangular plate under combined parametric and external excitations

In this paper, we use the asymptotic perturbation method based on the Fourier expansion and the temporal rescaling to investigate the nonlinear oscillations and chaotic dynamics of a simply supported rectangular plate made of functionally graded materials (FGMs) subjected to a through-thickness temperature field together with parametric and external excitations. Material properties are assumed to be temperature-dependent. Based on the Reddy’s third-order plate theory, the governing equations of motion for the plate are derived using the Hamilton’s principle. The Galerkin procedure is employed to obtain a two-degree-of-freedom nonlinear system including the quadratic and cubic nonlinear terms. The resonant case considered here is 1:2 internal resonance, principal parametric resonance-1/2 subharmonic resonance. Based on the averaged equation in polar coordinate form, the stability of steady state solutions is analyzed. The phase portrait, waveform and Poincaré map are used to analyze the periodic and chaotic motions of the FGM rectangular plate. It is found that the FGM rectangular plate exhibits the chaotic motions under certain circumstances. It is seen that the nonlinear dynamic responses of the FGM rectangular plate are more sensitive to transverse excitation. The excitation force can be used as a controlling factor which can change the response of the FGM rectangular plate from periodic motion to the chaotic motion.