首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 140 毫秒
1.
METHOD OF GREEN'S FUNCTION OF CORRUGATED SHELLS   总被引:1,自引:1,他引:0  
By using the fundamental equations of axisymmetric shallow shells of revolution, the nonlinear bending of a shallow corrugated shell with taper under arbitrary load has been investigated. The nonlinear boundary value problem of the corrugated shell was reduced to the nonlinear integral equations by using the method of Green's function. To solve the integral equations, expansion method was used to obtain Green's function. Then the integral equations were reduced to the form with degenerate core by expanding Green's function as series of characteristic function. Therefore, the integral equations become nonlinear algebraic equations. Newton' s iterative method was utilized to solve the nonlinear algebraic equations. To guarantee the convergence of the iterative method, deflection at center was taken as control parameter. Corresponding loads were obtained by increasing deflection one by one. As a numerical example,elastic characteristic of shallow corrugated shells with spherical taper was studied.Calculation results show that characteristic of corrugated shells changes remarkably. The snapping instability which is analogous to shallow spherical shells occurs with increasing load if the taper is relatively large. The solution is close to the experimental results.  相似文献   

2.
The truncated expansion method for finding explicit and exact soliton-like solution of variable coefficient nonlinear evolution equation was described. The crucial idea of the method was first the assumption that coefficients of the truncated expansion formal solution are functions of time satisfying a set of algebraic equations, and then a set of ordinary different equations of undetermined functions that can be easily integrated were obtained. The simplicity and effectiveness of the method by application to a general variable coefficient KdV-MKdV equation with three arbitrary functions of time is illustrated.  相似文献   

3.
Stability for basic system of equations of atmospheric motion   总被引:1,自引:1,他引:0  
The topological characteristics for the basic system of equations of atmospheric motion were analyzed with the help of method provided by stratification theory. It was proved that in the local rectangular coordinate system the basic system of equations of atmospheric motion is stable in infinitely differentiable function class. In the sense of local solution, the necessary and sufficient conditions by which the typical problem for determining solution is well posed were also given. Such problems as something about "speculating future from past" in atmospheric dynamics and how to amend the conditions for determining solution as well as the choice of underlying surface when involving the practical application were further discussed. It is also pointed out that under the usual conditions, three motion equations and continuity equation in the basic system of equations determine entirely the property of this system of equations.  相似文献   

4.
The topological characteristics for the basic system of equations of atmo- spheric motion were analyzed with the help of method provided by stratification theory.It was proved that in the local rectangular coordinate system the basic system of equations of atmospheric motion is stable in infinitely differentiable function class.In the sense of local solution,the necessary and sufficient conditions by which the typical problem for determining solution is well posed were also given.Such problems as something about“speculating future from past”in atmospheric dynamics and how to amend the condi- tions for determining solution as well as the choice of underlying surface when involving the practical application were further discussed.It is also pointed out that under the usual conditions,three motion equations and continuity equation in the basic system of equations determine entirely the property of this system of equations.  相似文献   

5.
NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS   总被引:3,自引:0,他引:3  
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.  相似文献   

6.
A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.  相似文献   

7.
In this paper we consider the singular perturbation boundary-value problem of thefollowing coupling type system of convection-diffusion equationsWe advance two methods:the first one is the initial value solving method,by which theoriginal boundary-value problem is changed into a series of unperturbed initial-valueproblems of the first order ordinary differential equation or system so that an asymptoticexpansion is obtained;the second one is the boundary-value solving method,by which theoriginal problem is changed into a few boundary-value problems having no phenomenon ofboundary-layer so that the exact solution can be obtained and any classical numericalmethods can be used to obtain the numerical solution of consismethods can be used to obtainthe numerical solution of consistant high accuracy with respect to the perturbationparameterε  相似文献   

8.
Based on the Timoshenko beam theory, the finite-deflection and the axial inertia are taken into account, and the nonlinear partial differential equations for flexural waves in a beam are derived. Using the traveling wave method and integration skills, the nonlinear partial differential equations can be converted into an ordinary differential equation. The qualitative analysis indicates that the corresponding dynamic system has a heteroclinic orbit under a certain condition. An exact periodic solution of the nonlinear wave equation is obtained using the Jacobi elliptic function expansion. When the modulus of the Jacobi elliptic function tends to one in the degenerate case, a shock wave solution is given. The small perturbations are further introduced, arising from the damping and the external load to an original Hamilton system, and the threshold condition of the existence of the transverse heteroclinic point is obtained using Melnikov's method. It is shown that the perturbed system has a chaotic property under the Smale horseshoe transform.  相似文献   

9.
By the separation of singularity, a special Fourier series solution of the boundary valueproblem for plane is obtained, which can satisfy all boundary conditions and converges rapidly. It isproved that the solution is equal to the result of separation of variables. As a result, the non-linearcharacteristic equations resulting from the method of separation of variables are transformed into poly-nomial equations that can provide a foundation for approximate computation and asymptotic analysis.  相似文献   

10.
To the optimal control problem of affine nonlinear system, based on differential geometry theory, feedback precise linearization was used. Then starting from the simulative relationship between computational structural mechanics and optimal control, multiple-substructure method was inducted to solve the optimal control problem which was linearized. And finally the solution to the original nonlinear system was found. Compared with the classical linearizational method of Taylor expansion, this one diminishes the abuse of error expansion with the enlargement of used region.  相似文献   

11.
We present a robust numerical method for solving incompressible, immiscible two-phase flows. The method extends both a monolithic phase conservative level set method with embedded redistancing and a semi-implicit high-order projection scheme for variable-density flows. The level set method can be initialized conveniently via a simple phase indicator field instead of a signed distance function (SDF). To process the indicator field into a SDF, we propose a new partial differential equation-based redistancing method. We also improve the monolithic level set scheme to provide more accuracy and robustness in full two-phase flow simulations. Specifically, we perform an extra step to ensure convergence to the signed distance level set function and simplify other aspects of the original scheme. Lastly, we introduce consistent artificial viscosity to stabilize the momentum equations in the context of the projection scheme. This stabilization is algebraic, has no tunable parameters and is suitable for unstructured meshes and arbitrary refinement levels. The overall methodology includes few numerical tuning parameters; however, for the wide range of problems that we solve, we identify only one parameter that strongly affects performance of the computational model and provide a value that provides accurate results across all the benchmarks presented. This methodology results in a robust, accurate, and efficient two-phase flow model, which is mass- and volume-conserving on unstructured meshes and has low user input requirements, making it attractive for real-world applications.  相似文献   

12.
A three zone unsteady state mass transfer problem is considered for flow between parallel plates. Solute is transferred from the moving zone to the other two, in one of which a first order reaction consuming the solute occurs. (An industrially important example of this situation would be the absorption of a gas from a lean mixture into a liquid-saturated ion exchange membrane.) The double Laplace transformation is applied to the system equations. In obtaining the inversion of the transformed equations the first inversion (with respect to the transformed dimensionless axial distance) is performed by use of the residue method, and then the second inversion (with respect to the transformed dimensionless time) is performed by use of the numerical technique developed by Bellman et al. Some numerical results are presented and discussed.  相似文献   

13.
This paper presents a numerical method that couples the incompressible Navier–Stokes equations with the level set method in a curvilinear co‐ordinate system for study of free surface flows. The finite volume method is used to discretize the governing equations on a non‐staggered grid with a four‐step fractional step method. The free surface flow problem is converted into a two‐phase flow system on a fixed grid in which the free surface is implicitly captured by the zero level set. We compare different numerical schemes for advection of the level set function in a generalized curvilinear format, including the third order quadratic upwind interpolation for convective kinematics (QUICK) scheme, and the second and third order essentially non‐oscillatory (ENO) schemes. The level set equations of evolution and reinitialization are validated with benchmark cases, e.g. a stationary circle, a rotating slotted disk and stretching of a circular fluid element. The coupled system is then applied to a travelling solitary wave, and two‐ and three‐dimensional dam breaking problems. Some interesting free surface phenomena are revealed by the computational results, such as, the large free surface vortices, air entrapment and splashing of the water surge front. The computational results are in excellent agreement with theoretical predictions and experimental data, where they are available. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

14.
Helmholtz方程的微分容积解法   总被引:1,自引:0,他引:1  
用一种新型的数值技术--微分容积法(Differential Cubature Method)求解二维Helmholtz方程的边值问题,几个数值算例表明,该方法稳定收敛,并具有较好的数值精度,本文方法适用于求解具有较小波数的Helmholtz方程。  相似文献   

15.
本文基于大变形的理论,采用弧坐标首先建立了具有初始位移的桩基的非线性数学模型,一组强非线性的微分-积分方程,其中,地基的抗力采用了Winkeler模型;其次,引入变数变换将微分-积分方程转化为一组非线性微分方程,并用微分求积方法离散了方程组,得到一组离散化的非线性代数方程;最后用Newton-Raphson迭代方法对离散化方程进行了求解,得到了桩基变形前后的构形、弯矩和剪力.计算中选取了两种不同类型的初始位移,并考察了它们对桩基大变形力学行为的影响.  相似文献   

16.
This paper presents a numerical method, a transmission matrix method, for the wave propagation in viscoelastic stratified saturated porous media. The wave propagation in saturated media, based on Biot theory, is a coupled problem. In this stratified three-dimensional model we do the Laplace transform for the time variable and the Fourier transform for the horizontal space coordinate. The original problem is transformed into ordinary differential equations with six independent unknown variables, which are only the function of the coordinate of depth. Thus, we get a transmission matrix of the wave problem for each layer. In the process of solution we use numerical method to calculate the eigenvalues and the eigenvectors of the transmission matrices. In the first step of the solution process we can obtain the wave field in the transformed space. The fast Fourier transform (FFT) method is used to do the inverse Laplace and the inverse Fourier transforms to get the solution in the time space. The detailed formulae are derived and some numerical examples are given.  相似文献   

17.
A hybrid numerical-analytical solution for steady-state natural convection in a porous cavity is proposed, based on application of the ideas in the generalized integral transform technique. The integral transformation process reduces the original coupled partial differential equations, for temperature and stream function, into an infinite system of non-linear ordinary differential equations for the transformed potentials, which is adaptively truncated and numerically solved through well-established algorithms. The approach is applied to a vertical rectangular enclosure subjected to uniform internal heat generation. The convergence characteristics of the explicit inversion formulae are illustrated and critical comparisons with previously reported purely numerical solutions are performed.  相似文献   

18.
This paper deals with the stress concentration in plane with swveral arbitrarily distributed elliptic holes. By using the functions of complex variables, the stress functions in which the interactions of neighbouring holes are taken into consideration can be constructed. By applying the conformed mapping method to satisfy the boundary conditions of each hole, the governing equations can then be transformed into a set of simultaneous equations through boundary integrals. Moreover, the problems with crack can be derived by changing the elliptical rates of the ellipses, thereby an approximate solution of cracking problem may be obtained. Some computing examples are given in the paper.  相似文献   

19.
Awrejcewicz  J.  Kudra  G.  Lamarque  C.-H. 《Meccanica》2003,38(6):687-698
This report is a part of the larger project of non-linear dynamics investigation of three coupled physical pendulums with damping and with arbitrary situated barriers, and externally driven. The set of differential equations and the set of algebraic inequalities (representing a barrier) governing the motion of three coupled rods are presented in the non-dimensional form. The system of governing equations is integrated between two successive impacts, and the discontinuity points are detected (by halving time step until a required precision is obtained). In each impact time, the state of the system is transformed using the extended restitution coefficient rule. The theory of Aizerman and Gantmakher is used to calculate the fundamental solution matrices in the analyzed system exhibiting discontinuities. The fundamental matrices are used during calculation of Lyapunov exponents, during stability analysis of periodic solutions (Floquet multipliers) and in shooting method applied to detect and trace periodic orbits. Some examples for three coupled identical rods with horizontal barrier are reported.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号