排序方式: 共有59条查询结果,搜索用时 501 毫秒
1.
2.
建立了插值矩阵法的基本理论,用于解非线性混合阶常微分方程组多点边值问题,制作了常微分方程组求解器IVMMS,可以支持计算力学中的有限元线法。 相似文献
3.
4.
正交各向异性位势问题边界元法中几乎奇异积分的解析算法 总被引:3,自引:0,他引:3
几乎奇异积分的计算困难阻碍了边界元法的工程应用。本文针对二维正交各向异性位势问题边界元法中近边界点的几乎奇异积分,采用分部积分法,导出一种直接的解析计算公式。该解析公式可以精确计算线性单元上的几乎奇异积分。对二次单元,可将其细分为几个线性元,采用该解析公式近似计算其边界积分。当内点离当前积分单元较远时,仍保持常规高斯数值积分模式;而当内点离其较近时,因常规高斯积分结果失效,则采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元计算结果比线性元计算结果更精确。 相似文献
5.
研制了一种适用于二维正交各向异性位势问题的高阶单元(线性单元和二次单元)快速多极边界元法. 在快速多极边界元法中, 源点对于远场区域的积分采用快速多极展开式计算, 而对于近场区域的积分则直接进行计算. 高阶单元的使用使得近场积分, 尤其是奇异积分和几乎奇异积分的计算更加复杂. 通过引入复数表达对其进行简化, 若边界采用线性单元插值, 近场积分可直接解析计算; 若采用二次单元插值, 则给出一个半解析算法计算近场积分. 高阶单元奇异积分和几乎奇异积分计算难题的解决, 使得高阶单元快速多极边界元法不仅能够计算一般结构, 也能被应用于超薄体结构, 拓宽了高阶单元快速多极边界元法的适用范围. 数值算例表明, 若计算精度一定, 高阶单元快速多极边界元法较常值单元快速多极边界元法使用的单元数量显著减少, 且高阶单元快速多极边界元法计算时间与自由度数量成线性关系, 其计算效率仍处于$O(N)$量级, 因此高阶单元快速多极边界元法可更加高效求解大规模问题. 相似文献
6.
7.
针对功能梯度材料矩形板问题,基于三维弹性理论,将位移和应力分量作为基本变量,通过双三角级数将其控制微分方程转化为常微分方程组的边值问题。采用插值矩阵法直接对常微分方程组边值问题进行求解,得到了功能梯度材料矩形板三维位移、应力场的半解析解。通过算例给出了材料参数按指数形式和幂函数形式变化情况下的功能梯度板的弯曲问题。对比有限元法和状态空间法,结果表明:本文提出的状态空间与插值矩阵法相结合的半解析法能有效地分析材料参数按任意形式连续变化的功能梯度矩形板问题,且具有良好的精度,精度可达10-4量级,能够满足工程需要;与其他方法相比,本文方法具有实施便捷、计算量小等优点,根据其力学场分析结果可设计出满足各种不同需求的功能梯度材料。 相似文献
8.
提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段. 相似文献
9.
10.
基于忽略了梁截面剪切变形和转动惯量效应的Euler-Bernoulli梁理论,研究了轴向力作用下轴向功能梯度变截面梁的横向自由振动问题,将轴向功能梯度Euler-Bernoulli梁自由振动固有频率和临界荷载的计算转化为变系数常微分方程特征值问题。运用插值矩阵法可一次性计算出轴向功能梯度变截面梁各阶振动固有频率和临界荷载,分析了轴向荷载对轴向功能梯度Euler-Bernoulli梁自由振动固有频率的影响,即轴向压力使梁的第1阶固有频率降低,轴向拉力使梁的第1阶固有频率增大。在简支-简支梁(H-H)边界条件下、不同截面宽锥度系数c_b和截面高锥度系数c_h,且区间划分点数n为40时,本文计算结果与已有文献计算结果之间的最大相对误差不超过0.00768%;在简支-简支梁(H-H)、固端-自由梁(C-F)、固端-固端梁(C-C)这三种不同边界条件下,不同c_b和c_h,且n为40时,最大相对误差不超过0.101%,说明了本文方法的有效性和良好的计算精度。 相似文献