变厚度簿板弯曲问题的任意网格差分解法

本文提出一种求解变厚度薄板弯曲问题的任意网格差分格式,可适应不同边界,各种荷载和复杂形状板.计算实例表明,该方法具有格式简单、通用性强、计算精度高,计算量少等特点.     

Differential quadrature method based on the highest derivative and its applications

The numerical differentiation is often used when dealing with the differential equations. Using the numerical differentiation, the differential equations can be transformed into algebraic equations. Then we can get the numerical solution from the algebraic equations. But the numerical differentiation process is very sensitive to even a small level of errors. In contrast, it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, we provide a new method using the DQ method based on the interpolation of the highest derivative (DQIHD) for the differential equations. The original function is then obtained by integration. In this paper, the DQIHD method was applied to the buckling analysis of thin isotropic plates and Winkler plates, the numerical results agree well with the analytic solutions, and the results show that our method is of high accuracy, of good convergence with little computational efforts. And it is easy to deal with the boundary conditions.     

用双向三角级数法解悬臂矩形薄板在均布荷载下的弯曲

悬臂矩形板的弯曲问题是平板理论中的一个难题.多年来,对于这种板只有能量法与数值解法的近似解.1979年以来清华大学张福范教授等用迭加法陆续得出悬臂矩形板在均布荷载和一些集中荷载作用下的解析解.对于在均布荷载作用下的悬臂矩形薄板,本文用双向三角级数法获得了其挠度函数的解析解,并将所得结果与迭加法所得的结果进行了比较.通过比较表明,两种方法计算的结果符合得十分好,因而相互印证了它们的正确性.     

A four-node discrete singular convolution for geometric transformation and its application to numerical solution of vibration problem of arbitrary straight-sided quadrilateral plates

A four-node discrete singular convolution (DSC) method is developed for free vibration analysis of arbitrary straight-sided quadrilateral plates. The straight-sided quadrilateral domain is mapped into a square domain in the computational space using a four-node element. By using the geometric transformation, the governing equations and boundary conditions of the plate are transformed from the physical domain into a square computational domain. Numerical examples illustrating the accuracy and convergence of the DSC method for skew, trapezoidal, rhombic and arbitrary quadrilateral plates are presented. The results obtained by DSC method were compared with those obtained by the other numerical methods.     

Free vibration and buckling analyses of composite plates with straight-sided quadrilateral domain based on DSC approach

Discrete singular convolution (DSC) method has been proposed to obtain the frequencies and buckling loads of composite plates. By using geometric transformation, the straight-sided quadrilateral domain is mapped into a square domain in the computational space using a four-node element. Plates having different geometries such as rectangular, skew, trapezoidal and rhombic plates are presented. The obtained results are compared with those of other numerical methods. Numerical results indicate that the DSC is a simple, accurate and reliable algorithm for vibration and buckling analyses of composite plates.     

矩形薄板的屈曲状态*

本文根据薄板的广义变分原理,用有限元法对单向受压矩形板的屈曲状态进行了讨论,并用延续算法获得了在不同边界条件下的屈曲状态.     

任意载荷作用下矩形弹性薄板横向弯曲的一个高精度近似解*

李兹法是近似求解弹性薄板横向弯曲的一种广泛使用的有效方法,其精度完全取决于基函数的选择.本文根据矩形薄板横向弯曲的特点,将基函数选择为正弦三角级数与多项式函数的叠加,不但公式简单易程序化,而且有着很高的精度.本文最后给出了两个算例,并与经典结果进行了比较.     

Bifurcations and chaos of the nonlinear viscoelastic plates subjected to subsonic flow and external loads

The subharmonic bifurcations and chaotic motions of the nonlinear viscoelastic plates subjected to subsonic flow and external loads are studied by means of Melnikov method. The critical conditions for the occurrence of chaotic motions are obtained. The chaotic features on the system parameters are discussed in detail. The conditions for subharmonic bifurcations are also obtained. For the system with no structural damping, chaotic motions can occur through infinite subharmonic bifurcations of odd orders. Furthermore, we confirm our theoretical predictions by numerical simulations. The theoretical results obtained here can help us to eliminate or suppress large nonlinear vibrations and chaotic motions of the nonlinear viscoelastic plates. Based on Melnikov method, complex dynamical behaviors of the nonlinear viscoelastic plates can be controlled by modifying the system parameters.     

Use of Adomian decomposition method in the study of parallel plate flow of a third grade fluid

In this paper, Adomian’s decomposition method is used to solve non-linear differential equations which arise in fluid dynamics. We study basic flow problems of a third grade non-Newtonian fluid between two parallel plates separated by a finite distance. The technique of Adomian decomposition is successfully applied to study the problem of a non-Newtonian plane Couette flow, fully developed plane Poiseuille flow and plane Couette–Poiseuille flow. The results obtained show the reliability and efficiency of this analytical method. Numerical solutions are also obtained by solving non-linear ordinary differential equations using Chebyshev spectral method. We present a comparative study between the analytical solutions and numerical solutions. The analytical results are found to be in good agreement with numerical solutions which reveals the effectiveness and convenience of the Adomian decomposition method.     

Local and post local buckling of stepped and perforated thin plates

The nonlinear mathematical theory for initial and post local buckling analysis of plates of abruptly varying stiffness based on the principle of virtual work is established. The method is programmed, and several numerical examples are presented to demonstrate the scope and efficacy of the procedure. Local buckling coefficients of perforated and stepped plates are obtained and the results are compared with known solutions. Post-buckling behaviour of perforated and stepped plates is studied for different geometries. The non-dimensional applied loads (P/P_cr), dimensionless lateral displacements and stress distribution of plates with varying stiffness in the post buckling region are given in several tables and graphs.     

基于厚板拉伸振动精确化方程求解动应力集中问题

过去,对拉伸平板考虑应力集中的工程设计多借鉴弹性力学平面问题分析求解结果,例如弹性力学Kirsch问题的解或弹性动力学平面问题的解．基于厚板拉伸振动精确化方程,对含圆孔平板中弹性波散射与动应力集中问题进行了研究．研究结果表明:1) 两种模型得到的开孔附近的应力是不同的;2) 当入射波波数变大或者说入射波频率变高时,动应力集中系数最大值趋于单位1．含孔平板拉伸振动的动应力集中系数最大值达到3.30,以及基于弹性动力学平面问题模型得到的结果为2.77．对数值计算结果做了分析讨论, 可以看到,当孔径厚度比是a/h=0.10,基于平板拉伸振动精确化方程得到的动应力集中系数可以达到最大值,超出基于弹性动力学平面问题所得到结果的19%．分析方法和数值计算结果可望能在工程平板结构的动力学分析和强度设计中得到应用．     

横向剪切对双模量复合材料叠层矩形板非线性弯曲的影响

本文探讨了动力松弛(DR)法在双模量复合材料叠层矩形板非线性弯曲问题中的应用。在分析中分别采用叠层板大挠度经典理论和计及大转动(在Von Karman意义上)的复合材料叠层板剪切变形理论。我们发现,对于考虑横向剪切变形的非线性弯曲问题,如何计算虚拟密度以控制数值计算的稳定性,仍然需要进一步研究。本文提出了一种虚拟密度的计算方法,从而保证了本课题数值计算的稳定性。文中介绍了用DR法求解双模量复合材料叠层板非线性弯曲的主要步骤,给出了由轻度双模量材料(Born-Epxy(B-E))和高度双模量材料(Aramid-Rubber(A-R)和Polyester-Rubber(P-R))的两层正交叠层简支矩形板在正弦分布载荷及均布载荷作用下的非线性弯曲特性的数值结果。将所得结果和小挠度分析结果及普通复合材料的结果作了比较,并分析了横向剪切变形对无量纲中心挠度的影响。     

A symplectic analytical wave propagation model for damping and steady state forced vibration of orthotropic composite plate structure

An analytical wave propagation model is proposed in this paper for damping and steady state forced vibration of orthotropic composite plate structure by using the symplectic method. By solving an eigen-problem derived in the symplectic dual system of free bending vibration of orthotropic rectangular thin plates, the wave shape of plate is obtained in symplectic analytical form for any combination of simple boundary conditions along the plate edges. And then the specific damping capacity of wave mode is obtained symplectic analytically by using the strain energy theory. The steady state forced vibration of built-up plates structure is calculated by combining the wave propagation model and the finite element method. The vibration of the uniform plate domain of the built-up plates structure is described using symplectic analytical waves and the connector with discontinuous geometry or material is modeled using finite elements. In the numerical examples, the specific damping capacity of orthotropic rectangular thin plate with three different combinations of boundary condition is first calculated and analyzed. Comparisons of the present method results with respect to the results from the finite element method and from the Rayleigh–Ritz method validate the effectiveness of the present method. The relationship between the specific damping capacity of wave mode and that of modal mode is expounded. At last, the damped steady state forced vibration of a two plates system with a connector is calculated using the hybrid solution technique. The availability of the symplectic analytical wave propagation model is further validated by comparing the forced response from the present method with the results obtained using the finite element method.     

Analytical solution of piezolaminated rectangular plates with arbitrary clamped/simply-supported boundary conditions under thermo-electro-mechanical loadings

This paper extends an analytical method for static analysis of general cross-ply piezolaminated rectangular plates with any combination of clamped/simply-supported boundary conditions under uncoupled thermo-electro-mechanical loadings. This method is based on the novel superposition method and the first-order shear deformation theory (FSDT). The FSDT enables this expanded method to consider the effect of shear deformation of the plate. The process of applying electrical and thermal resultants causes some advantages due to its simplicity and less computational process. In this analysis displacement components are written in terms of unknown force and moment resultants. Using Fourier series for displacement components, mechanical, thermal, and/or electrical stress resultants, the complex governing differential equations of the plate are reduced to a set of linear algebraic equations with non-trivial solution. The obtained equations may be solved analytically to determine the unknown stress resultants. Several examples are proposed, and their obtained numerical results are compared with those available in the literature to verify the convergence, high accuracy, and the capability of the present method to analyze the static behavior of piezolaminated plates. It is found that there is high agreement between the present results with those obtained by other investigators.     

任意厚度具有自由边叠层板的精确解析解

自由边问题一直是三维弹性力学中的难题,通常很难满足自由边上一个正应力和两个剪应力都等于0．基于三维弹性力学基本方程和状态空间方法,引入自由边界位移函数并考虑全部弹性常数,建立了正交异性具有自由边单层和叠层板的状态方程．对状态方程中的变量以级数形式展开,通过边界条件的满足精确求解任意厚度具有自由边叠层板的位移和应力,此解满足层间应力和位移的连续条件．算例计算表明,采用引入的位移函数形式,简化了计算过程并且采用较少的级数项可以获得收敛解．与有限元方法计算结果进行了对比,可以得到较高精度的数值结果．其解可以作为其它数值方法和半解析方法的参考解．     

轴对称圆板(含叠层板)的三维非线性分析

本文提出了轴对称固支圆板(含叠层板)受均布横向载荷作用下的三维非线性摄动解答.文中所考虑的是一种中等大挠度的几何非线性,并采用一种发展的摄动方法对复杂的三维非线性平衡微分方程进行求解.该方法的基本思想是以二维解答为基础,对板的厚度参数进行摄动而求得相应的三维解答.文中给出了一般板及叠层板的三维非线性理论结果及数值结果,并图示出了各个应力的分布情况.而且,该三维非线性结果能退化为完全一致的相应的二维板理论非线性结果.结果表明,该方法对板的三维非线性分析是一种行之有效的方法.     

New exact series solutions for transverse vibration of rotationally-restrained orthotropic plates

The exact series solutions of plates with general boundary conditions have been derived by using various methods such as Fourier series expansion, improved Fourier series method, improved superposition method and finite integral transform method. Although the procedures of the methods are different, they are all Fourier-series based analytical methods. In present study, the foregoing analytical methods are reviewed first. Then, an exact series solution of vibration of orthotropic thin plate with rotationally restrained edges is obtained by applying the method of finite integral transform. Although the method of finite integral transform has been applied for vibration analysis of orthotropic plates, the existing formulation requires of solving a highly non-linear equation and the accuracy of the corresponding numerical results can be questionable. For that reason, an alternative formulation was proposed to resolve the issue. The accuracy and convergence of the proposed method were studied by comparing the results with other exact solutions as well as approximate solutions. Discussions were made for the application of the method of finite integral transform for vibration analysis of orthotropic thin plates.     

Large amplitude vibration of heated corrugated circular plates with shallow sinusoidal corrugations

This paper presents a nonlinear free vibration analysis of corrugated circular plates with shallow sinusoidal corrugations under uniformly static ambient temperature. Based on the nonlinear bending theory of thin shallow shells, the governing equations for corrugated plates are established from Hamilton’s principle. These partial differential equations are reduced to corresponding ordinary ones by elimination of the time variable with Kantorovich method following an assumed harmonic time mode. The resulting equations, which form a nonlinear two-point boundary value problem in spatial variable, are then solved numerically by shooting method, and the temperature-dependent characteristic relations of frequency vs. amplitude for nonlinear vibration of heated corrugated plates are obtained successfully. The comparison with available published results shows that the numerical analysis here is of good reliability. A detailed parametric study is conducted involving the dependency of nonlinear frequency on the depth and density of corrugations along with the temperature change. Effects of these variables on the trend of nonlinearity are plotted and discussed.     

Application of the Haar wavelet-based discretization technique to problems of orthotropic plates and shells

The Haar wavelet discretization technique for solving the elastic bending problems of orthotropic plates and shells is proposed. Free transverse vibrations of orthotropic rectangular plates with a variable thickness in one direction are considered as a model problem. In the case of constant plate thickness, the numerical results are validated by comparing them with an exact solution. The results obtained are found to be in good agreement with those available in the literature.     

Secondary Buckling Analysis of Thin Rectangular Plates Based on the Wavelet Galerkin Method北大核心CSCD

Application of the wavelet Galerkin method (WGM) to numerical solution of nonlinear buckling problems was studied with classical elastic thin rectangular plates. First, the discretized scheme of the von Kármán equation were introduced, then a simple calculation approach to the Jacobian and Hessian matrices based on the WGM was proposed, and the wavelet discretized scheme-based eigenvalue equation method, the extended equation method and the pseudo arc-length method for nonlinear buckling analysis were discussed. Second, the secondary post-buckling equilibrium paths of elastic thin rectangular plates and the effects of aspect ratios, boundary conditions and bi-directional compression on the mode jumping behaviors, were discussed in detail. Numerical results show that, the WGM possesses good convergence for solving buckling loads on rectangular plates, and the obtained equilibrium paths are in good agreement with those from the stability experiments, the 2-step perturbation method and the nonlinear finite element method. Given the feasibility of combination with different bifurcation computation methods, the WGM makes an efficient spatial discretization method for complex nonlinear stability problems of typical plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.