解黏弹性地基梁受迫振动的时域DQ法

采用时域DQ法探讨了黏弹性地基上长度为L的Euler-Bernoulli梁的横向振动.该方法直接从控制微分方程出发,在空间域和时间域均采用离散的DQ法,得到求位移场全部待定参数的可解线性方程组.针对策动力F(x,t)=Q·sinwt作用下黏弹性地基梁的动力学初—边值条件问题的计算结果表明,方法的数学原理依据充分,精度好,效率高,具有推广应用的价值.     

自由振动矩形薄板动力响应的DQ空-时半解析法

对矩形薄板自由振动的动力响应问题提出了DQ半解析法,该方法针对矩形薄板的线性振动控制微分方程,在空间域采用DQ法(differential quadrature method,DQ),在时间域取级数,采用时域配点的方法,得到求解以板各节点动力响应位移场为全部待定参数的线性方程组,只需一次求解该方程组即可得到各节点的动力响应位移场,再由高阶Lagrange插值得到全域内的动力响应位移场.结果表明,该方法具有很好的精度和计算效率.     

复杂载荷下梁弯曲问题的微分求积法应用研究

文章直接从细长梁的小挠度弯曲控制微分方程出发,利用微分求积法(DQM),沿梁的轴线把梁在空间域离散成若干个点,对每一点都可得到一关于各离散点挠度的DQ方程,从而得到关于求解全部离散点挠度的线性方程组,求解该方程组即可得到各点挠度,再由高阶Lagrange插值即可得到全域内的位移场。     

矩形薄板瞬态动力响应的DQ空-时半解析法研究

在矩形薄板的瞬态动力响应问题的研究中,提出一种新的方法——DQ空-时半解析法.该方法针对矩形薄板的振动控制微分方程,在空间域采用一种高效的数值方法,微分求积法(differential quadrature method),即DQ法,在时间域取解析形式,采用时域配点的方法,得到求解以板各节点动力响应位移场为全部待定参数的线性方程组,只需一次性求解该方程组即可得到各节点的动力响应位移场,再由高阶Lagrange插值即可得到全域内的动力响应位移场,算例结果显示,该方法是一种新的精度高、效率好的处理结构动力响应的计算方法.     

矩形薄板动力响应的DQ半解析法研究

针对矩形薄板的动力响应问题,提出了一种有效的方法:DQ半解析法,本方法针对矩形薄板的振动控制微分方程,在空间域采用DQ法,即微分求积法(differential quadrature method),在时间域取级数,采用时域配点的方法,得到求解以板各节点动力响应位移场为全部待定参数的线性方程组,只需一次求解该方程组即得到全部待定参数,进而得到各节点的动力响应位移场,再由高阶Lagrange插值得到全域内的动力响应位移场.算例结果表明,本方法具有很高的精度和极佳的计算效率,且不受边界条件约束.     

求强非线性动力系统周期解的MH法和EMH法

提出了改进的谐波平衡法即MH法和EMH法,用于求强非线性系统周期解.MH法直接应用于求解大参数参变系统:x a(t)x b(t)x=c(t).应用谐波平衡法求解强非线性系统,当谐波项取得较少时,求解结果精度低,为此引入最小二乘原理对谐波平衡法加以改进,计算结果精度高.EMH法用于求强非线性自治与非自治系统:x F(x,x)=0和x F(x,x,Ωt)=0的周期解,由于F呈非线性,不能直接应用MH法,为此,先由能量原理得出一次近似解,再引入牛顿迭代原理,得到关于修正量的周期系数方程;用MH法求此修正量,得到的结果精度较高.     

泰勒级数在数值法中的应用

图示法求解具有给定X和Y初值的常微分联立方程组: (dx)/(dt)=F(x,y,t) (dy)/(dt)=G(x,y,t)     

变厚度矩形薄板弯曲问题的GD法研究

针对矩形薄板的线性挠曲控制微分方程及边界条件,在空间域采用GD法离散,得到求解以薄板各节点弯曲挠度为全部待定参数的线性方程组,只需一次求解该方程组即可得到各节点的弯曲挠度,再用高阶La-grange插值即可得到薄板全域内任意点的挠度.     

非线性方程组的逆Broyden秩1拟Newton方法及其在MATLAB中的实现

对于非线性方程组F(x)=0,Newton迭代公式x(k+1)=x(k)-[F′(x(k))]-1F(x(k))(k=0,1,2,…)的最大优点在于其形式简单且是超线性收敛的,而最大的缺点在于对初值依赖性强且每一次迭代均需要计算Jacobi矩阵及其逆矩阵,计算量大,易导致误差累积传播.通过对Newton迭代公式的逐步改进,展现了逆Broy-den秩1拟Newton方法的形成过程,并以一具体例子,实现该方法在MATLAB7.5环境中的数值求解过程.     

利用最小作用原理研究一类非自治二阶系统周期解

主要研究以下二阶系统{u(t)-A(t)u(t)=▽F(t,u(t)),u(0)-u(T)=u(0)-u(T)=0,a.e.t∈[0,T]的周期解的存在性。当F(t,x)=F1(t,x)+F2(x)满足条件A且具有局部有界性T1lim inf x→+∞x 2α∫F(t,x)dt0T2T∫(r1(t)dt)2/0T12-T∫k(t)dt及A(t)满足条件(A(t)x,x)≥h(t)|x|β+w(t)时,通过使用最小作用原理得到了一个新的周期解的存在性定理,改进了已有结果。