SMALL-STENCIL PADE SCHEMES TO SOLVE NONLINEAR EVOLUTION EQUATIONS

A set of small-stencil new Pade schemes with the same denominator are presented to solve high-order nonlinear evolution equations. Using this scheme, the fourth-order precision can not only be kept, but also the final three-diagonal discrete systems are solved by simple Doolittle methods, or ODE systems by Runge-Kutta technique. Numerical samples show that the schemes are very satisfactory. And the advantage of the schemes is very clear compared to other finite difference schemes.     

有限元网格自动生成及自动调整技术

本文给出了一种任意平面域内三角形网格自动生成及自动调整算法，并编制了相应的程序。通过对带Ｕ形槽的三点变曲试件的分析，证明本文算法是可行的，并且具有稳定性能好、收敛速度快的优点。     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

桥梁空间分析中预应力效应分析方法研究

预应力技术在桥梁中的广泛使用,使得桥梁结构受力状况复杂,对预应力混凝土桥梁进行空间分析是必要的;但由于预应力筋数量庞大且形状复杂,对预应力空间效应的模拟是预应力混凝土桥梁空间分析的难点。从预应力筋的实际作用机理出发,提出在预应力张拉阶段用等效节点荷载来模拟预应力对结构的作用;而在预应力管道灌浆后,将预应力筋作为结构的一部分,用实体退化单元进行结构分析,并用两算例对本文的方法进行了验证。结果表明,该方法能真正模拟预应力的空间效应,可用于预应力混凝土桥的空间分析。     

基于分离式共结点模型的钢筋混凝土结构爆破拆除数值模拟

通过大型有限元分析软件ANSYS/LS-DYNA,充分考虑钢筋和混凝土在结构倒塌时各自的受力状态和强度差异,运用分离式共节点模型模拟钢筋混凝土结构建筑物的爆破拆除倒塌过程.通过对倒塌过程中钢筋混凝土支撑立柱内侧和外侧的钢筋单元、混凝土单元的承载失效过程分析研究,对数值模拟效果与实际爆破效果对比分析得出:分离式共结点模型...     

框架剪力墙结构的定向爆破倒塌过程

利用ANSYS/LS-DYNA有限元软件,分别采用整体式模型和共用节点分离式模型对某高层框架剪力墙结构的定向爆破拆除进行数值模拟.通过与实际工程的对比发现,分离式模型模拟的效果更贴近实际,模拟结果更准确;此外,待拆除结构的爆破切口适当上移,可以减少结构的倒塌范围;最后,分析了钢筋和混凝土单元的应力时程曲线,指出共用结点...     

基于节点密度的自适应EFG法拓扑优化方法研究

针对无网格法拓扑优化中设计变量庞大、计算效率低的问题,提出了基于节点密度的自适应EFG拓扑优化方法.该方法以节点相对密度为设计变量,建立了以结构柔度最小为优化目标的自适应EFG拓扑优化模型,给出了拓扑优化中自适应节点加密准则、节点加密方法和新增节点设计变量更新规则.算例结果表明,采用论文提出的自适应EFG拓扑优化方法进行结构拓扑优化时,可以减少结构分析和优化的设计变量数目,提高优化效率.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

甲状腺癌转移性淋巴结的红外光谱研究

联合运用衰减全反射(attenuated total reflectance,ATR)探头与FTIR光谱仪,测定了新鲜离体的甲状腺癌转移性颈部淋巴结20枚和非转移性淋巴结69枚的FTIR光谱,比较两组光谱13个谱带的峰位及相对峰强等28个指标,找出甲状腺癌转移性淋巴结光谱的特征。与非转移性淋巴结相比,转移性淋巴结光谱中与蛋白、脂质、糖类、核酸相关的谱带均发生了明显的变化,具体如下：(1)与蛋白质相关的谱带的相对峰强I_{3 280}/I1 460(P=0.029),I_{1 640}/I_{1 460}(P=0.019),I_{1 546}/I_{1 460}(P=0.004)均显著升高;(2)脂类相关的相对峰强I_{1 743}/I_{1 460}显著降低(P=0.026),而I_{1 400}/I_{14 60}则明显增高(P=0.001);(3)与糖类相关的1 165与1 120 cm-1处峰位P_{1 165}(P=0.019),P_{1 120}(P=0.000)发生明显蓝移(向高波数移动),相对峰强I1 165/I_{1 460}显著降低(P=0.006);(4)与核酸相关的1 085 cm-1处峰位P_{1 085}发生明显地红移(向低波数移动)(P=0.029);(5)其他指认未明的谱带的相对峰强I_{1 303}/I_{1 460}(P=0.008),I_{1 303}/I1 240(P=0.039)均明显升高。与非转移性淋巴结光谱相比较,甲癌转移的淋巴结光谱在多个谱带的峰位或相对峰强方面存在显著的差异。FTIR可应用于术中甲状腺癌转移性淋巴结的探测。     

高压电气设备支架串联体系地震响应半解析法

针对具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系,通过应用分布参数梁振动理论与集中参数通过边界条件的引入,推导出其频率方程,然后通过数值方法求得频率及振型。通过应用Betti定律,推导出具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系振型的正交条件,应用此正交条件对集中参数与分布参数的振动方程解耦,从而得到广义质量及广义荷载,然后通过振型叠加法求解结构的地震响应。最后应用该半解析法与有限元法数值积分法对一330 KV电压互感器支架串联体系地震响应进行分析,分析结果表明,该半解析法求解的地震响应与有限元直接积分法求解的地震响应相一致,从而为具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系的地震响应分析和抗震设计提供了一种新的途径。