应用DQE增量迭代法分析斜直井内管柱的非线性屈曲

基于斜直井内管柱屈曲的平衡微分方程,构建DQE (Differential Quadrature Element)增量迭代法,对斜直井内管柱的非线性屈曲进行计算.通过与有限元的计算结果对比,验证方法的正确性.DQE法方法简单、易于实施,计算量少、精度较高,所得到的螺旋屈曲计算结果与实验结果吻合,最大井壁约束力随上端载荷的增加而增大.对于工程中比较长的受压段管柱,其屈曲是一个局部非线性稳定性问题,屈曲首先从下端开始发生,随着上端载荷的增加,逐渐向上扩展;上端边界条件对下端局部屈曲无明显影响.     

QL-Born迭代法电磁散射计算

积分方程法是多维电磁响应的数值模拟和反演计算的有效方法之一.在线性化散射电磁场积分方程的基础上,采用迭代Born近似法进行了电导率成像反演.针对迭代Born近似法反演依赖初始模型,而QL近似法无需给定初值的特点,提出了两步法(QL-Born迭代法),用QL近似法计算的结果作为Born迭代法的初始模型,避免人为设置初值,并进行了数值实验.     

基于概率迭代的NDP反演方法

针对中子深度分析(NDP)技术,基于概率迭代思想,推导一种迭代反演方法,对其和线性正则化方法在NDP反演问题中的应用进行比较.未引入计数随机误差时,两者都能得到较好的反演结果,迭代法的结果是非负的,线性正则化方法在源强度发生阶跃处得到的结果更好.在考虑随机误差存在的情况下,迭代法的效果更佳.同时,研究了迭代矩阵对结果的影响,矩阵元的指数放大可以实现迭代过程的加速.     

水面目标复合电磁散射的并行迭代快速计算

采用快速远场近似的降秩电磁流迭代法计算目标与粗糙水面的复合散射, 并进一步采用粗糙面分区的多进程并行迭代算法提高其计算速度, 分析了降秩电磁流迭代法对散射结果的影响, 表明该方法可用于近场复合散射问题的计算; 对平板目标与水面间的复合散射进行了计算, 计算结果与快速多极子方法(MLFMM)结果符合较好, 验证了该算法的有效性. 同时, 利用造波水池开展了目标与水面近场复合散射的模拟实测, 测试与计算结果趋势较为一致.     

修正的变分迭代法在四阶Cahn-Hilliard方程和BBM-Burgers方程中的应用

变分迭代法是一种基于变分原理,具有高数值精度的数值格式,目前已广泛应用于各类强非线性孤立波方程的数值求解中.本文利用修正的变分迭代法对两类非线性方程进行研究.该格式是对原数值方法的一种改进,即在变分项前引入了参数h.通过定义误差函数的离散二范数并在定义域内绘出h-曲线,从而确定出使误差达到最小的h,再返回原迭代过程进行求解.同时,参数的引入也扩大了原数值解的收敛域,在迭代次数一定的情况下达到了数值最优.在数值实验中,将上述结果应用于四阶的Cahn-Hilliard方程和BenjaminBona-Mahoney-Burgers方程.对于四阶的Cahn-Hilliard方程,普通的变分迭代法绝对误差在10~(-1)左右,经过修正后,绝对误差降为10~(-4),而且修正后的方法扩大了原数值解的收敛域.对于Benjamin-Bona-MahonyBurgers方程,利用带有辅助参数的变分迭代法将数值解的精度提高到10~(-3),对真解的逼近效果优于原始的变分迭代法.此数值方法也为其他强非线性孤立波微分方程的数值求解提供了方法和参考.     

斜直圆筒内基于Euler-Rodrigues参数的细长杆螺旋屈曲研究

基于Euler-Rodrigues参数描述受约束细杆的后屈曲变形状态,建立相应的能量泛函,导出非线性平衡微分方程及接触力表达式.采用有限元法对能量泛函进行分析计算,通过与现有文献中数值、理论结果比较,验证了公式和求解过程的正确性.结合算例,揭示了正弦屈曲模态与螺旋屈曲模态之间的转化过程,证实了在受约束细长杆的后屈曲响应中存在不同屈曲模态之间的互相转变.     

一种研究敞口容器非稳态导热问题的新方法

在变导热系数条件下,把分段迭代法应用于解决非稳态导热物理问题的集总参数法,为敞口容器恒温热源加热过程中容器内的平均温度变化提供了一个简便的计算方法.利用该方法只需三次分段迭代,即可使得理论计算与实验数据结果吻合.     

面向纯相位型全息显示的自适应混合约束迭代算法

纯相位型计算全息图舍弃全息面上的振幅信息,降低了全息显示的图像质量.迭代法通过多次正、逆向的波前传播计算,优化了纯相位型计算全息图的相位轮廓,提升了全息再现的图像质量.然而,传统迭代法可能存在迭代发散、收敛速度偏慢以及再现质量受限等问题.本文提出了一种基于角谱传播模型的自适应混合约束迭代算法,通过在迭代中加入自适应的反馈机制和频带约束策略,增加了迭代优化的自由度,提升了迭代的收敛性和全息再现的图像质量.仿真和实验结果表明,提出的算法能够在较短的计算时间内获得更高重建质量的纯相位型计算全息图,对于高质量的全息显示具有重要应用价值.     

高超声速可压缩流中粗糙壁热流研究

用计算流体力学（CFD）数值模拟和理论方法对高超声速可压缩湍流中粗糙壁面热增量进行研究. 着重考虑粗糙单元密度和粗糙单元形状对粗糙表面热流的影响. 结果表明： ①粗糙单元密度变化时, CFD数值方法计算所得的粗糙单元等效热流随着粗糙单元密度的降低而增加, 理论方法预测结果的规律随方法不同而不同;②在相同粗糙单元密度和高度时,如果粗糙单元形状改变,CFD计算结果也随之发生改变. 理论预测方法所得结果不发生变化.     

基于低温共烧陶瓷工艺的一种新型层叠式多层结构的波概念迭代方法研究

提出了基于波概念的一种有效迭代方法,详细给出了波概念迭代法（wave concept iterative process）理论推导公式,并用来处理基于低温共烧陶瓷工艺的新型多层结构计算问题,涉及到了混合磁场与电场方程,包括基于一种简单形式的匹配负载模拟上的S参数提取问题.利用matlab语言编写程序, 从而确定S参数的频率响应,并将数值结果与仿真结果进行比较,两者呈现出很好的一致性,同时验证了波概念迭代法相对于其他算法的优势所在. 关键词： 波概念迭代法 快速傅里叶模式变换 低温共烧陶瓷工艺     

Surface effect on buckling configuration of nanobeams containing internal flowing fluid: A nonlinear analysis

In this paper, the post-buckling behavior of supported nanobeams containing internal flowing fluid with two surface layers is studied based on a nonlinear theoretical model. The nonlinear governing equation, in which the surface effect and stretching-related nonlinearity are accounted for, is analytically solved for both clamped-clamped and pinned-pinned systems. The effects of nanobeam length, bulk thickness and several dimensionless parameters on the post-buckling behavior are analyzed. It is found that, the nanobeam with low flow velocity is stable at its original static equilibrium position and then undergoes a buckling instability at a critical flow velocity, which depends on the system parameters. When buckled, in all cases, the amplitude of the resultant buckling increases with the increasing flow velocity. Typically, the surface effect is explored by considering different nanobeam lengths and bulk thicknesses. The buckling amplitude is found to be length-dependent and thickness-dependent, showing that the effect of surface layers is considerably strong.     

Orientation dependent size effects in thermal buckling and post-buckling of nanoplates with cubic anisotropy

In this work, a continuum model is presented for size and orientation dependent thermal buckling and post-buckling of anisotropic nanoplates considering surface and bulk residual stresses. The model with von-Karman nonlinear strains and material cubic anisotropy of single crystals contains two parameters that reflect the orientation effects. Using Ritz method, closed form solutions are given for buckling temperature and post-buckling deflections. Regarding self-instability states of nanoplates and their recovering at higher temperatures, an experiment is discussed based on low pressurized membranes to verify the predictions. For simply supported nanoplates, the size effects are lowest when they are aligned in [100] direction. When the edges get clamped, the orientation dependence is ignorable and the behavior becomes symmetric about [510] axis. The surface residual stress makes drastic increase in buckling temperature of thinner nanoplates for which a minimum thickness is pointed to stay far from material softening at higher temperatures. Deflection of [100]-oriented buckled nanoplates is higher than [110] ones but this reverses at higher temperatures. The results for long nanoplates show that the buckling mode numbers are changed by orientation which is verified by FEM.     

Prediction of compressive post-buckling behavior of single-walled carbon nanotubes in thermal environments

In the present investigation, the axial buckling and post-buckling configurations of single-walled carbon nanotubes (SWCNTs) are studied including the thermal environment effect. For this purpose, Eringen’s nonlocal elasticity continuum theory is implemented into the classical Euler–Bernoulli beam theory to represent the SWCNTs as a nonlocal elastic beam model. A closed-form analytical solution is carried out to analyze the static response of SWCNTs in their post-buckling state in which the axial buckling load is assumed to be beyond the critical axial buckling load. Common sets of boundary conditions, named simply supported–simply supported (SS–SS), clamped–clamped (C–C), and clamped–simply supported (C–SS), are considered in the investigation. Selected numerical results are given to represent the variation of the carbon nanotube’s mid-span deflection with the applied axial load corresponding to various nonlocal parameters, length-to-diameter aspect ratios, temperature changes, and end supports. Moreover, a comparison between the post-buckling behaviors of SWCNTs at low- and high-temperature environments is presented. It is found that the size effect leads to a decrease of the axial buckling load especially for SWCNTs with C–C boundary conditions. Also, it is revealed that the value of the temperature change plays different roles in the post-buckling response of SWCNTs at low- and high-temperature environments.     

电磁内爆套筒屈曲的数值模拟

 采用瞬态非线性有限元方法对电磁内爆套筒屈曲进行了数值模拟研究。介绍了瞬态非线性有限元方法的基本原理、电磁内爆套筒模拟的基本方法及对实验的模拟结果。通过分析比较认为,建立的数值模拟方法较好地反映了实验结果,可以作为研究分析电磁内爆套筒屈曲的基本方法之一。     

轴向冲击载荷作用下双壁碳纳米管的动力屈曲

应用改进的有限元方法,建立考虑层间范德华力作用的壳-弹簧非线性有限元模型,基于B-R运动准则,系统地研究了双壁碳纳米管的动力屈曲问题,得到了轴向冲击载荷作用下双壁碳纳米管的临界动力屈曲载荷和临界动力失效载荷. 研究结果表明,在动力屈曲过程中,双壁碳纳米管层间距的变化非常小,各管的变形相互协调;碳纳米管中应力波的传播导致碳纳米管出现非对称屈曲模态,可明显观测到四个环向波瓣,沿着碳纳米管的轴线方向,四个波瓣的波峰和波谷交替变化. 对碳纳米管动力屈曲问题的研究表明,冲击载荷的大小和持续时间对碳纳米管的动力屈曲有 关键词： 碳纳米管 动力屈曲 冲击载荷     

Nonlinear analysis of imperfect squarely-reticulated shallow spherical shells

Nonlinear behavior of single-layer squarely-reticulated shallow spherical shells with geometrical imperfections subjected to a central concentrated (joint) load has been studied in this paper. Using the asymptotic iteration method, an analytical characteristic relationship between the non-dimensional load and central deflection is obtained. The resulting asymptotic solution can be used readily to perform the analysis of parameters and predict the buckling critical load. Meanwhile, numerical examples are presented and effects of imperfection factor and boundary conditions on buckling of the structures are discussed. Comparisons with data based on the finite element method show good exactness of the resulting solution.     

空间弹性细杆在扭矩作用下的屈曲

用有限元方法研究了三维弹性细杆在扭矩作用下的屈曲.利用自然坐标形式的细长空间曲杆的能量方程和2节点12个自由度的自然坐标形式的三维曲梁单元,采用特征值分析方法,研究分析了同时受有轴力和扭矩作用时的空间弹性细杆的屈曲问题.数值结果与存在的理论解极为吻合.具有一定曲率和挠率的空间细长曲杆,其临界扭矩值与扭矩的指向有着极其明显的关系.     

非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧，提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步，在粗网格上采用Oseen迭代法求解非线性问题，在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散，时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.