任意形状孔口双边裂纹平板的应力强度因子计算

本文采用Muskhclishvili弹性力学的复变函数和边界配位方法对不同形状孔口双边裂纹问题进行了研究,计算了圆孔、椭圆孔、矩形孔、菱形孔等不同形状孔口双边裂纹,以及Ⅰ型和复合型等不同类型断裂试件的应力强度因子,本文方法简单方便,精度较高,与某些已有计算结果的问题比较,本文方法所得的结果是令人满意的.同时,本方法可以应用于不同几何形状和加载条件下的孔口双边裂纹有限大板的计算,是解这一类问题的一致有效方法.     

带两个形状参数的四次Bézier曲线的扩展

给出了一组含有两个形状参数α,β的四次多项式基函数,是四次Bernstein基函数的扩展,分析了这组基的性质;基于这组基定义了带两个形状参数的多项式曲线,所定义的曲线不仅保留了四次Bézier曲线一些实用的几何特征,而且具有形状的可调性,在控制多边形不变的情况下,改变参数α,β的取值,可以生成不同的逼近控制多边形的曲线;通过分析该曲线与四次Bézier曲线之间的关系,给出了α和β的几何意义,并利用Bézier曲线递归分割算法给出了这种曲线的几何作图法,同时还讨论了曲线间的拼接问题.     

任意形状孔口单边裂纹问题的边界配位解法

本文提出了一组复应力函数,采用边界配位方法对不同形状孔口(包括圆、椭圆、矩形及菱形孔口)的单边裂纹平板的应力强度因子进行了计算.计算结果表明,对长度和宽度远大于孔口和裂纹几何尺寸的试件,配位法与用其他方法所得的无限大板含圆或椭圆孔边裂纹问题的解符合得很好.同时,对其他孔口问题,特别是有限大板情形,本文给出了一系列计算结果.本文所提出的函数及计算过程可以应用于任意形状孔口单边裂纹平板的计算.     

空心柱形杆的弹性扭转问题

本文证明了对具有任意形状截面的空心柱形杆的弹性扭转问题,当杆的内孔以一定的方式收缩为一点时,杆的应力在L_2空间中强收敛于相应实心杆的应力,从而证明了对柱形杆的弹性扭转问题不存在孔边应力集中现象。     

恶性肿瘤的传质问题（Ⅱ）——药物传输部分

本文提出肿瘤中药物传输的三重介质模型·基于第（Ⅰ）部分流体动力学的计算结果，采用有限元方法，在不同初始和边界条件下求解了对流—扩散方程组·根据计算结果，详细分析了药物注射方式、分子量、淋巴、药物结合和坏死区等对抗体免疫球蛋白ＩｇＧ及其片断Ｆａｂ在肿瘤中浓度分布的影响·     

考虑几何非线性的杆系结构弯曲变形样条有限点法

研究了杆系结构考虑几何非线性的大挠度弯曲变形问题,推算并验证了一种考虑几何非线性的杆系结构弯曲变形计算新方法.以三次B样条函数为基函数,采用广义参数法,构造出梁的样条基函数,通过最小势能原理,建立了杆系结构考虑几何非线性的刚度方程,对处于弹性范围内的杆系结构的大变形弯曲问题进行了计算,提出了考虑几何非线性时杆系结构弯曲变形计算的样条有限点法.结果表明:方法不用进行单元坐标变换、划分等分数少、收敛速度快且计算精度较高,是一种较传统有限元法更简单且可行的方法.     

RBF-PU方法求解二维非局部扩散问题和近场动力学问题

采用单位分解径向基函数(radial basis function partition of unity,RBF-PU)方法,数值求解了二维非局部扩散问题和近场动力学问题。主要思想是对求解区域进行局部划分,在局部子区域上分别进行函数逼近,然后加权得到未知函数的全局逼近。这种基于方程强形式的径向基函数方法在求解非局部问题时,不需要处理网格与球形邻域求交的问题,避免了额外的一层积分计算,实施简便,计算量小。数值实验显示计算结果与解析解吻合较好,RBF-PU方法可以准确有效地求解非局部扩散方程和近场动力学方程。     

基于变分辨率网格的模拟退火算法在形状优化问题上的应用研究

形状优化是一类复杂的优化问题,在工业上经常作为结构优化的一个分支出现.它以几何形状作为优化对象,需要求得某种性能条件最优的几何形状,一般来说对建模和计算量的要求都比较高.模拟退火是一种用途非常广泛的优化算法,可以处理各种复杂的优化问题.但标准的模拟退火算法在处理形状优化问题时,由于搜索空间的范围太大,经常陷入局部最优解,所以需要耗费巨大的计算量,实用性有限.本文讨论了在变分辨率的离散网格中的模拟退火算法,将低分辨率网格下的最优结果过渡到高分辨率网格下作为一个良好的初始解,这样可以有效地规避局部最优点,缩小搜索范围,极大地提高模拟退火算法的效率.并以蚱蜢问题为算例,对算法的效果进行了验证.     

有理Bzier曲线

§1.前言 计算几何中曲线造型的主要工具是代数参数曲线。其中,按照Bernstein基和B样条基表示的参数曲线,称为Bezier曲线和B样条曲线,尤其应用广泛。 苏步青教授最早把代数曲线论的仿射不变量理论导进计算几何领域,用以研究仿射平面参数曲线的几何性质,特别是关于那些以实拐点和实奇点个数为特征的仿射分类,从而获得一系列具有重要应用价值的结果,推动了计算几何的理论发展。近来,这些结果被应用到CAGD的工程技术课题中去,收到了成效。     

周期压力驱动下管壁局部不规则的圆管Poiseuille流感受性问题

采用渐近分析方法,建立了在周期压力驱动下,完全发展的圆管Poiseuille流当管壁存在局部不规则几何形状时的感受性问题模型．通过特征函数的双正交系统,应用Chebyshev配点法进行数值求解．通过算例计算,获得周期压力和矩形突起激发起的流体系统中的各种空间发展模态以及相应的感受性系数．从计算和分析可以知道,在流场的不同发展阶段不同的模态起着主导作用,这与在试验中观察到的扰动流场在不同位置的特性是一致的．     

Theoretical study of diffusional release of a dispersed solute from cylindrical polymeric matrix: A novel configuration for providing zero-order release profile

In the context of controlled release drug delivery approaches, the systems providing zero-order release kinetics have special advantages. Through employing these systems, drug concentration could be maintained within the therapeutic window over release time; thus maximum effectiveness alongside minimized side effects of the drug are achieved. However, obtaining zero-order drug release is extremely challenging. One of the main obstacles is the fact that implemented devices should be designed to overcome the decreasing mass transfer driving force, especially, in polymeric systems in which diffusion mechanism is dominant. In this study, we developed a new configuration of a polymeric matrix containing dispersed solute which provides sustained zero-order release. A combination of two innovative approaches including separating baffles and a two-layer coating was proposed to be incorporated into the conventional cylindrical polymeric matrix to induce zero-order release behavior. Then, an approximate mathematical model was developed to investigate the performance of the system under different conditions. The simulated results showed the potential of proposed configuration to be used as a carrier for sustained zero-order release.     

Mathematical models for drug delivery to chronic patients

Mathematically, analysis of drug delivery kinetics involves two moving boundary problems: diffusion front and eroding front. In this paper, we have models for drug delivery for the sites which can be enclosed by spherical shaped matrices covered by membranes and these problems are helpful for designing the drug delivery devices to deliver the drug inside from outside and a corresponding device supplying drug from inside. Once the time required for treatment and rate of drug delivery is known from medical diagnosis, this analysis can design a device releasing the drug/active agent over a long period of time. The purpose of such drug delivery is to achieve more effective therapies while eliminating the effect of over dosing and maintaining drug levels within the desired levels. The device may work on optimal use of drug and increase the patient’s convenience. The proposed models provide design for eroding tumor or chemotherapy to cancerous regions. The results have been obtained for steady state release rate, zero order release time and life time of the device and discussed. It has been observed that zero order time and life time increase by introducing a membrane of uniform thickness.     

A second‐order finite difference method for fractional diffusion equation with Dirichlet and fractional boundary conditions

In this paper, we develop a practical numerical method to approximate a fractional diffusion equation with Dirichlet and fractional boundary conditions. An approach based on the classical Crank–Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second‐order accurate numerical estimates. The solvability, stability, and convergence of the proposed numerical scheme are proved via the Gershgorin theorem. Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

Linearized factorization techniques for multidimensional reaction—diffusion equations

An iterative predictor—corrector technique for the elimination of the approximate factorization errors which result from the factorization of linearized θ-methods in multidimensional reaction—diffusion equations is proposed, and its convergence and linear stability are analyzed. Four approximate factorization techniques which do not account for the approximate factorization errors are developed. The first technique uses the full Jacobian matrix of the reaction terms, requires the inversion of, in general, dense matrices, and its approximate factorization errors are second-order accurate in time. The second and third methods approximate the Jacobian matrix by diagonal or triangular ones which are easily inverted but their approximate factorization errors are, however, first-order accurate in time. The fourth approximately factorized method has approximate factorization errors which are second-order accurate in time and requires the inversion of lower and upper triangular matrices. The techniques are applied to a nonlinear, two-species, two-dimensional system of reaction—diffusion equations in order to determine the approximate factorization errors and those resulting from the approximations to the Jacobian matrix as functions of the allocation of the reaction terms, space and time.     

Convergence of the matrix transformation method for the finite difference approximation of fractional order diffusion problems

Numerical solution of fractional order diffusion problems with homogeneous Dirichlet boundary conditions is investigated on a square domain. An appropriate extension is applied to have a well-posed problem on R² and the solution on the square is regarded as a localization. For the numerical approximation a finite difference method is applied combined with the matrix transformation method. Here the discrete fractional Laplacian is approximated with a matrix power instead of computing the complicated approximations of fractional order derivatives. The spatial convergence of this method is proved and demonstrated by some numerical experiments.     

On the applicability of Genocchi wavelet method for different kinds of fractional‐order differential equations with delay

A novel collocation method based on Genocchi wavelet is presented for the numerical solution of fractional differential equations and time‐fractional partial differential equations with delay. In this work, to achieve the approximate solution with height accuracy, we employed the operational matrix of integer derivative and the pseudo‐operational matrix of fractional derivative in Caputo sense. Also, based on Genocchi function properties, we presented delay and pantograph operational matrices of Genocchi wavelet functions (GWFs). Due to operational and pseudo‐operational matrices, the equations under this study can be turned into nonlinear algebraic equations with the unknown GWF coefficients. For illustrating the upper bound of error for the proposed method, we estimate the error in the sense of Sobolev space. In addition, to demonstrate the efficacy of the pseudo‐operational matrix of fractional derivative, we investigate the upper bound of error for the mentioned matrix. Finally, the algorithm based on the proposed approach is implemented for some numerical experiments to confirm accuracy and applicability.     

Fractional‐order orthogonal Bernstein polynomials for numerical solution of nonlinear fractional partial Volterra integro‐differential equations

In this paper, a new two‐dimensional fractional polynomials based on the orthonormal Bernstein polynomials has been introduced to provide an approximate solution of nonlinear fractional partial Volterra integro‐differential equations. For this aim, the fractional‐order orthogonal Bernstein polynomials (FOBPs) are constructed, and its operational matrices of integration, fractional‐order integration, and derivative in the Caputo sense and product operational matrix are derived. These operational matrices are utilized to reduce the under study problem to a nonlinear system of algebraic equations. Using the approximation of FOBPs, the convergence analysis and error estimate associated to the proposed problem have been investigated. Finally, several examples are included to clarify the validity, efficiency, and applicability of the proposed technique via FOBPs approximation.     

几乎各向同性的高维空间分数阶扩散方程的分块快速正则Hermite分裂预处理方法

一类空间分数阶扩散方程经过有限差分离散后所得到的离散线性方程组的系数矩阵是两个对角矩阵与Toeplitz型矩阵的乘积之和.在本文中,对于几乎各向同性的二维或三维空间分数阶扩散方程的离散线性方程组,采用预处理Krylov子空间迭代方法,我们利用其系数矩阵的特殊结构和具体性质构造了一类分块快速正则Hermite分裂预处理子.通过理论分析,我们证明了所对应的预处理矩阵的特征值大部分都聚集于1的附近.数值实验也表明,这类分块快速正则Hermite分裂预处理子可以明显地加快广义极小残量(GMRES)方法和稳定化的双共轭梯度(BiCGSTAB)方法等Krylov子空间迭代方法的收敛速度.     

Fractional Boundary Value Problems with Integral Boundary Conditions via Topological Degree Method

This paper examines the existence and uniqueness of solutions for the fractional boundary value problems with integral boundary conditions. Banach''s contraction mapping principle and Schaefer''s fixed point theorem have been used besides topological technique of approximate solutions. An example is propounded to uphold our results.     

Solving 2D time‐fractional diffusion equations by a pseudospectral method and Mittag‐Leffler function evaluation

Two‐dimensional time‐fractional diffusion equations with given initial condition and homogeneous Dirichlet boundary conditions in a bounded domain are considered. A semidiscrete approximation scheme based on the pseudospectral method to the time‐fractional diffusion equation leads to a system of ordinary fractional differential equations. To preserve the high accuracy of the spectral approximation, an approach based on the evaluation of the Mittag‐Leffler function on matrix arguments is used for the integration along the time variable. Some examples along with numerical experiments illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.