考虑墙趾阻力的加筋土挡墙的数值分析与筋材拉力计算方法

现行加筋土挡墙设计规范由于忽略了墙趾阻力对挡墙内部稳定性的贡献,计算的筋材拉力值远大于实测值,在设计上过于保守．本文对55座假定的模块式加筋土挡墙开展数值模拟,研究工作应力条件下筋材最大拉力随墙趾阻力的变化规律及其影响因素,在此基础上,采用响应面模型提出考虑墙趾阻力的筋材拉力计算方法,并用现场挡墙实测数据对该方法进行验证．结果表明,随着墙趾阻力增大,挡墙下部1/3高度范围内的筋材最大拉力值减小;墙趾阻力对筋材拉力的影响随着挡墙高度和填土内摩擦角的增加而减弱,随墙面仰角的增加而增强,受筋材刚度和间距影响不大;本文提出的方法可以合理预测加筋土挡墙在工作应力状态下的筋材拉力,并能反映出墙趾阻力使挡墙底部筋材拉力减小的现象．     

波动数值模拟的一种显式方法——一维波动

本文提出了波动时域数值模拟的一种新的显式方法,并通过一维非规则网格节点递推公式的建立说明此方法。为了阐明应用此法构建高精度、稳定递推公式的可行性,详细论述了一维均匀网格标量波动数值模拟的精度和稳定性;提出了构建时空精度皆为 阶( 为正整数)的稳定递推公式的技术途径,并以构建二阶(M=1)和四阶(M=2)公式为例予以说明。最后,通过算例详细说明了本文理论结果。     

一种改进的分块隐式数值方法

提出了一种改进的分块隐式数值方法，在贴体坐标和交错网格下以逆变速度分量和压力Ｕ，Ｖ，Ｗ，ｐ为基本求解变量，由此克服了原分块隐式数值方法求解复杂边界流动时的困难．９０°弯管流动数值计算初步表明，本文提出的方法合理、可行     

时间推进方法在叶轮机械内部流场计算中的进展

时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法．本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式时间推进方法和隐式时间推进方法计算叶轮机内部流场的最新进展,并简要介绍世界上一些著名的研究机构开发的时间推进方法计算程序所采用的数值计算方法      

一种结构动力时程分析的积分求微方法

传统采用微分求积(differential quadrature,DQ)法求解动力问题时都是以位移响应作为基本未知量,而将速度响应和加速度响应表示为位移响应的加权和的形式.如此做法需要处理线性方程组或者矩阵方程(Sylvester方程)才能求得动力响应,导出的算法一般为有条件稳定算法.本文利用动力响应的Duhamel积分解,逆用DQ原理,提出了一种计算卷积的高精度显式算法.该算法可以逐时段地求解出动力时程响应,当各时段内DQ节点分布完全一致时,仅须进行一次Vandermonde矩阵求逆计算即可应用于各个时段,一次性获得时段内多个时刻的位移响应值,因而具有计算效率高的优点.通过分析动力方程积分格式,证明本文动力算法传递矩阵的谱半径恒等于1,因而该算法具有无条件稳定特性,且计算过程中不会产生数值耗散. 本文算法的数值精度取决于分析时段内布置的DQ节点数量$N$,具有$N-1$阶代数精度.实际操作时可以取10个甚至更多的DQ节点数,从而获得比较高的数值精度.      

一种基于DIC测量单边缺口拉伸试件CTOD阻力曲线的简便方法

针对测量单边缺口拉伸(SENT)试样的裂纹尖端张口位移(CTOD)阻力曲线时,延性裂纹扩展量的计算过于复杂的问题,提出了一种由数字图像相关技术(DIC)计算的变形场来确定裂纹尖端位置,并快速推算延性裂纹扩展量的方法。该方法首先在缺口和裂纹上下方的弹性区域设置两条平行横线;然后通过两条横线间的位移差构建曲线,曲线变化幅度的临界点即为裂纹端部的位置;根据临界点位置计算裂纹扩展量,最后测量CTOD并构建阻力曲线。将本文方法与另一种通过DIC测量SENT试样CTOD阻力曲线的方法进行对比,结果表明本文方法经过简便的计算后拟合的阻力曲线效果更好。     

分析薄壳塑性动力响应问题的一种数值方法

探索薄壳塑性动力响应问题的分析与数值方法,是当今力学中的重要课题之一。笔者曾给出联合应用拉普拉斯变换和最小二乘法分析梁、板、壳塑性动力响应的方法。本文分别针对承受中载的短壳和长壳问题,给出只用最小二乘法分析其塑性动力响应的方法。     

基于显式有限元方法的两相介质弹塑性动力反应计算分析

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性...     

一种连续索滑移的处理方法

目的是建立一种快速简便的拉索滑移问题的有限元计算方法。采用基于弹性悬链线精确解的两节点悬链线单元模拟索结构中的拉索，利用求解非线性方程组的延拓法，再配以常用的牛顿法，建立了快速的给定索原长情况下的拉索刚度计算方法。此法对于一般的给定初值，能很快逼近真实解，从而扩大了解的收敛范围，更重要的是与牛顿法相比更加高效。进而采用二分法建立了索结构中拉索滑移问题处理的算法，编制了相应的程序段，并通过算例证明了其正确性。本文的拉索滑移处理方法可用于索结构的设计与施工。     

动、静摩擦对可接触型裂纹动态影响的一种算法

本文提出了一种分析考虑裂纹面摩擦影响的可接触型裂纹动力反应的时域显式计算方法，这一方法不但可以较好描述裂纹面之间碰撞过程，也能反映当裂纹面发生接触－滑移时，裂纹面间动、静摩擦力的影响．文中给出了使用这一方法的基本步骤，计算了无限大裂纹的动接触问题和直线型裂纹对入射Ｐ波的散射问题，初步验证了本文提出的方法．     

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

数值模型内力影响面分析的机动法

本文提出并证明了普遍适用于以广义位移为独立变量的各种数值模型的内力影响面分析的机动算法原理。这一原理把内力影响面分析归结为求解在一个等价荷载作用下的结构位移场,并给出了得到等价荷载的方法。用本文方法得到的是结构某一点的内力的影响面,而且对数值模型来说,是精确解。作为本文原理的应用,用机动法分析了采用有限条模型的薄板的内力影响面,并给出了若干数值例子。     

板料三维成形过程的有限元数值模拟

基于有限变形理论建立了三维金属板料成形过程的刚塑性有限元数学模拟,该数学模型采用物质坐标系中的Ｕｐｄａｔｅ－Ｌａｇｒａｎｇｅ描述,等向强化假设,考虑了板料的厚向异性,对于金属板料与模具有摩擦采用近似的库仑摩擦定律以改善计算的收敛性。为简经计算采用薄膜单元,并根据此模型编制程序。     

土体固结分析的一种有效算法

对土体固结分析,提出一种有效算法,避免了文[1]求解的迭代过程。文末附有算例。     

含能材料三轴动态压缩实验的一个计算机模拟方法

在含能材料单轴、受限三轴动态压缩实验的基础上，提出一种将计算机数值模拟和实验数据相结合，估算实验系统中钢筒与试件间的库伦摩擦力，从而正确确定含能材料的动态三轴力学性能（弹性模量Ｅ，泊桑比ν，屈服限Ｙ等）的方法。本方法用于ＴＮＴ炸药受限三轴动态压缩实验，通过建立钢筒－试件的库伦摩擦有限元模型，采用ＡＤＩＮＡ通用程序进行动力分析，估算了摩擦力对实验数据的影响，获得了考虑摩擦力修正后的ＴＮＴ三轴动态压缩响应，为含能材料三轴动态力学性能的确定提供了一种有效的方法。     

A finite element analysis of the flange earrings of strong anisotropic sheet metals in deep-drawing processes

Flange earrings of strong anisotropic sheet metals in deep-drawing process are numerically analyzed by the elastic-plastic large deformation finite element formulation based on a discrete Kirchhoff triangle plate shell element model. A Barlat-Lian anisotropic yield function and a quasi-flow corner theory are used in the present formulation. The numerical results are compared with the experimental ones of cylindrical cup drawing process. The focus of the present researches is on the numerical analysis and the constraining scheme of the flange earring of circular sheets with strong anisotropy in square cup drawing process. The project supported by the National Natural Science Foundation of China (19832020) and Provincial Natural Science Foundation of Jilin, China (200000519)     

索穹顶结构的找力分析方法--不平衡力迭代法

索穹顶的刚度和承载能力都是通过施加预应力来获取的,预应力设计是对索穹顶进行任何其它分析的基础,所以它的找力分析(Force finding)是一项很重要的工作。本文提出一种找力分析方法———不平衡力迭代法,它可以克服整体可行预应力法的一些不足:即使在索穹顶构件分组发生错误时仍可找到满足已知形状的预应力分布,并且该方法具有良好的稳定性和运算效率,算例表明不平衡力迭代法非常适合大型复杂索穹顶结构的找力分析。     

一种新的金属压缩试验方法的有限元分析

对一种新的金属压缩试验方法进行了有限元数值分析，论证了这种方法的可行性和有性，并且对不同尺寸的试件进行了有限元计算．     

力法在索-桁结构静力分析中的应用

索-梁(桁)结构是一种新型的空间杂交结构体系，它首先在中国获得工程应用和发展。通过下压平行布置在单层悬索上的实腹梁或桁架等横向加劲构件，强迫其产生向下位移后并固定其位置，以便在整个索与横向构件组成的体系之间建立起预应力，就形成具有足够刚度和形状稳定性的横向加劲单层索屋盖结构。本文根据力法原理，提出了一种新型迭代计算方法，采用Newtpm Raphson Method(N-R)方法计算索-桁结构体系的静力响应。算例计算结果表明：本文方法是可行的，计算精度满足工程需求。     

结构谐振动力响应的有限元、动凝聚、迁移子结构法

本文扩展了FETM法思想,结合子结构、动凝聚及Riccati变换理论,直接求解谐和激励下结构的受迫动力响应,从而节省了运行时间,严格控制了计算精度。本文编制了程序,并以振型叠加法运行结果为参照,对不同算例的运算效率作了比较。结果表明:本文方法在理论上、实践上是切实可行的。