用摄动法和有限条法分析矩形板的大挠度

本文将摄动法和有限条法结合起来进行矩形板的大挠度弯曲分析.用摄动的概念,将非线性微分方程组化为一系列线性微分方程组,然后用有限条法解这些线性微分方程组.     

在材料非线性问题中的摄动有限元法

摄动法是解决非线性连续介质力学问题的一种有效方法.这种方法是建立在该问题的线性解析解的基础上的,因此,若得不到一个简单的解析解,应用这种方法去解决一些复杂的非线性问题将遇到困难.有限元法对解非线性问题也是一种十分有用的工具,然而一般来说,它需要相当长的计算时间. 本文介绍摄动有限元法.这种方法吸取上述两种方法的优点,能够解决更复杂的非线性问题,而且也能大量节省计算机的计算时间. 本文讨论了比例加载下的弹塑性力学问题,并提出一个带孔拉板的数值解.     

非线性阻尼作用下标准线性固体粘弹性Ⅲ型破裂的解析解

把非线性Rayleigh阻尼引入标准线性固体粘弹性介质的Ⅲ型破裂的控制方程中,此方程是一个偏微分积分方程;首先设法消去积分项,得到一个三阶非线性偏微分方程,然后用小参数摄动法,得出线性化的各阶渐近控制方程;把每一个具有变系数的三阶线性控制方程分解为弹性部分及剩余部份,而前者的解析解为已知,后者是一个二阶变系数线性偏微分方程;它化不成Mathieu方程,也化不成Hill方程,故采用WKBJ的方法得出其渐近的解析解。     

分析土壤半无限域的一种有限层方法(Ⅰ)——一般理论

本文研究了一种用于横观各向同性体动力学的有限层法.将土壤介质半无限域视为一个横观各向同性半空间体,介质的材料函数沿深度变化,将介质沿深度方向分成若干层,在每一层内材料函数用一个指数函数来模拟.这样,使求解问题的方程得到简化,利用Fourier变换我们得到了各层处“节点”力与位移的关系.这种有限层法实际上是一种半解析法.具有一般半解析法所具备的数据前后处理量少及计算量小等优点.本文所研究的这种有限层法为分析土壤与结构相互作用问题提供了一条新途径.     

解半定规划的二次摄动方法

半定规划在系统论,控制论,组合优化,和特征值优化等领域有着广泛的应用。本文将半定规划摄动成二次半定规划,它的唯一解恰为原问题的解,并且对其偶问题等价于一个线性对称的投影方程,可方便地用投影收缩方法求解,从而获得原半定规划问题的解。文章给出了算法及其收敛性分析,数值试验结果表明摄动方法是解半定规划的一种有效的方法。     

介质与结构三维动力相互作用分析的半解析边界元—半解析有限元结合法

本文将半解析边界元一半解析有限无结合法用于介质与结构的动力相互作用研究:用半解析边界元法分析具有复杂地表面的半无限介质,用半解析有限元法分析具有任意截面形状的柱体结构,利用介质与结构交界面上的位移相容条件和力平衡条件,将介质与结构联系起来。联立京解上述半解析边界元方程和半解析有限元方程,对应每一时间步进,可同时求出介质与结构交界面上的位移、速度、加速度和相互作用力以及地表面的运动情况.与目前广泛研究的边界元—有限元结合法相比,本方法在介质与结构二个个区域各降低了一维空间,因而离散单元数和计算工作量大幅度减少,人工输入数据非常简单.文中还考虑了地下结构的长跨比效应、厚度效应和介质效应.     

半线性抛物型方程奇异摄动问题的数值解

本文讨论具有抛物边界层的半线性抛物型方程奇异摄动问题的数值解法,在非均匀网格上构造了两层非线性差分格式,证明了差分格式是一致收敛的,给出了一些数值例子.     

扁球壳和扁锥壳的轴对称非线性热弹耦合振动

假设温度场与应变场相互耦合,研究了旋转扁薄球壳和锥壳的轴对称非线性热弹振动问题．基于von Krmn理论和热弹性理论,导出了本问题的全部控制方程及其简化形式．应用Galerkin技术进行时空变量分离后,得到了一个关于时间的非线性常微分方程组．根据方程的特点,分别用多尺度法和正则摄动法求得了壳体振动的频率与振幅间特征关系和振幅衰减规律的一次近似解析解,并讨论了壳体几何参数、热弹耦合参数以及边界条件等因素对其非线性热弹耦合振动特性的影响．     

三阶奇异摄动问题及其在数值分析中的应用

§1.引言本文给出一类三阶非线性奇异摄动问题的解析渐近解法,此类问题人们经常在非线性数值稳定性分析中遇到,其解析渐近解有助于我们对若干非线性数值不稳定性问题的认识,也有助于我们对一些数值分析中奇怪的现.象(例如[1][2][3]中讨论的怪影解(GhostSolution))有个全局性认识.在§2中讨论一类三阶奇异摄动问题,这类问题和F.Howes     

均布载荷下矩形板大挠度问题的摄动变分解

本文以中心挠度为摄动参数,将矩形板大挠度问题的非线性偏微分方程组转化成几个线性偏微分方程,然后用变分法求解,得出了具有任意长宽比的板的解答,给出了位移、挠度及各内力的解析表达式;并给出中心点和边界中心的应力数值计算公式。本文还以长宽比λ为参数,作出了最大挠度——载荷曲线及最大应力曲线。其结果与实验进行了比较,表明二者是一致的。     

基于时域精细积分算法的瞬态传热多宗量反演

基于有限元法和精细积分算法,提出了一种求解瞬态热传导多宗量反演问题的新方法．采用有限元法和精细积分算法分别对空间、时间变量进行离散,可以得到正演问题高精度的半解析数值模型,由此建立了多宗量反演的计算模式,并给出敏度分析的计算公式．对一维和二维的热物性参数、热源项、边界条件等进行了单宗量和多宗量的反演求解,初步考虑了初值和噪音等对反演结果的影响,数值算例验证了该方法的有效性．     

摄动有限元法解一般轴对称壳几何非线性问题

本文在处理几何非线性问题时,利用在变分方程中引入振动过程,得到各级变分摄动方程,并通过有限元法求解.由于有限元法能成功地处理各种复杂边界条件、几何形状的力学问题,摄动法又可将非线性问题转化为线性问题求解.若结合这两种方法的优点,将能够解决大量复杂的非线性力学问题.并能够消除单独使用有限元法或摄动法求解复杂非线性问题所出现的困难. 本文应用摄动有限元法求解了一般轴对称壳的几何非线性问题.     

时域自适应精细算法求解二维非线性湿热耦合问题

采用时域自适应精细算法求解二维非线性湿热耦合问题．通过在离散时间段上将变量展开,更精确地描述变量的变化过程,将非线性时空耦合的初边值问题转化为一系列递推的线性边值问题,并利用有限元法进行求解．计算中不需任何附加假设和迭代过程,并提出了自适应计算的收敛判则,以保持时间步长选择不同时的计算精度．数值比较中,分别考虑了材料属性随温度变化、随湿度变化、以及随温度和湿度变化的情形,结果令人满意,体现了这种算法处理复杂非线性问题的能力．     

Asymptotics of eigenvalues of the nonlinear eigenvalue problem arising from the near mixed-mode crack-tip stress-strain field problems

In the present paper, approximate analytical and numerical solutions to nonlinear eigenvalue problems arising in nonlinear fracture mechanics in studying stress-strain fields near a crack tip under mixed-mode loading are presented. Asymptotic solutions are obtained by the perturbation method (the artificial small parameter method). The artificial small parameter is the difference between the eigenvalue corresponding to the nonlinear eigenvalue problem and the eigenvalue related to the linear “undisturbed” problem. It is shown that the perturbation technique is an effective method of solving nonlinear eigenvalue problems in nonlinear fracture mechanics. A comparison of numerical and asymptotic results for different values of the mixity parameter and hardening exponent shows good agreement. Thus, the perturbation theory technique for studying nonlinear eigenvalue problems is offered and applied to eigenvalue problems arising in fracture mechanics analysis in the case of mixed-mode loading.     

Homotopy perturbation-reproducing kernel method for nonlinear systems of second order boundary value problems

In this paper, based on homotopy perturbation method (HPM) and reproducing kernel method (RKM), a new method is presented for solving nonlinear systems of second order boundary value problems (BVPs). HPM is based on the use of traditional perturbation method and homotopy technique. The HPM can reduce a nonlinear problem to a sequence of linear problems and generate a rapid convergent series solution in most cases. RKM is also an analytical technique, which can solve powerfully linear BVPs. Homotopy perturbation-reproducing kernel method (HP-RKM) combines advantages of these two methods and therefore can be used to solve efficiently systems of nonlinear BVPs. Three numerical examples are presented to illustrate the strength of the method.     

B-Spline Collocation Method for Nonlinear Singularly-Perturbed Two-Point Boundary-Value Problems

A B-spline collocation method is presented for nonlinear singularly-perturbed boundary-value problems with mixed boundary conditions. The quasilinearization technique is used to linearize the original nonlinear singular perturbation problem into a sequence of linear singular perturbation problems. The B-spline collocation method on piecewise uniform mesh is derived for the linear case and is used to solve each linear singular perturbation problem obtained through quasilinearization. The fitted mesh technique is employed to generate a piecewise uniform mesh, condensed in the neighborhood of the boundary layers. The convergence analysis is given and the method is shown to have second-order uniform convergence. The stability of the B-spline collocation system is discussed. Numerical experiments are conducted to demonstrate the efficiency of the method.     

摄动初参数法解轴对称壳几何非线性问题

作者在文[7]中提出轴对称壳任意大挠度问题的一阶微分方程组和以载荷变量为尺度的变特征无量纲化方法。在此基础上,本文选取挠角非线性偏差的加权均方根作摄动参数,给出该问题的无量纲摄动微分方程组。从而把非线性问题转化为n个线性问题来解决。本文采用数值积分的初参数方法对摄动后的各阶线性问题进行了计算。摄动结果与实验[4]相符合。     

A finite differences MATLAB code for the numerical solution of second order singular perturbation problems

We show the main features of the MATLAB code HOFiD_UP for solving second order singular perturbation problems. The code is based on high order finite differences, in particular on the generalized upwind method. Within its simplicity, it uses order variation and continuation for solving any difficult nonlinear scalar problem. Several numerical tests on linear and nonlinear problems are considered. The best performances are reported on problems with perturbation parameters near the machine precision, where most of the codes for two-point BVPs fail.     

Method for solving nonlinear initial value problems by combining homotopy perturbation and reproducing kernel Hilbert space methods

The nonlinear singular initial value problems including generalized Lane–Emden-type equations are investigated by combining homotopy perturbation method (HPM) and reproducing kernel Hilbert space method (RKHSM). He’s HPM is based on the use of traditional perturbation method and homotopy technique and can reduce a nonlinear problem to some linear problems and generate a rapid convergent series solution in most cases. RKHSM is also an analytical technique, which can overcome the difficulty at the singular point of non-homogeneous, linear singular initial value problems; especially when the singularity appears on the right-hand side of this type of equations, so it can solve powerfully linear singular initial value problems. Therefore, using advantages of these two methods, more general nonlinear singular initial value problems can be solved powerfully. Some numerical examples are presented to illustrate the strength of the method.