首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
在考虑横向剪切变形对层合板弹性解的影响时,本文提出一种数值计算方法。由边界条件给出边界结点位移的表达式,根据薄板的经典理论和一阶横向剪切变形理论导出位移增量所满足的平衡微分方程,引用经典理论计算的横向剪力修正了荷载列阵。致使在较粗的网格划分时、宽广的层合板长厚比范围内,仍能得到与解析解颇为一致的数值解。  相似文献   

2.
喷水推进器推力的CFD计算方法研究   总被引:5,自引:1,他引:4  
简要介绍获取喷水推进器推力的理论、试验及数值计算(CFD)方法,重点研究采用动量流量法和壁面积分法计算喷水推进器推力的CFD方法。采用多块网格技术,用六面体结构化网格和四面体非结构化网格相结合的混合网格离散计算区域,采用稳态多参考系方法求解RANS方程,对喷水推进器进水流道、叶轮、导叶体和喷口所组成的整个流场进行数值计算。计算中采用了k-ε湍流模型和标准壁面函数,对用动量流量法计算推力方法中所需的假想流管分界面和进口面的求取做了分析,将两种方法计算的推力与厂商提供的推力特性曲线进行了对比。结果表明,采用CFD计算和分析方法来研究喷水推进系统推力性能是可行、可信的。  相似文献   

3.
基于径向基函数(RBF) 的网格变形方法是一种可靠的网格变形技术,对于任意拓扑的网格都能获得高质量的变形网格. 缩减控制点的RBF 网格变形方法可以大幅提高网格变形效率,但也存在变形后物面误差较大、边界层网格交错的问题. 在缩减控制点方法的基础上,提出了一种适合于带有边界层的黏性网格变形的方法,该方法从物面中选择两组控制点,利用其中一组控制点粗略计算变形后网格位置及变形误差,利用第二组控制点与变形误差插值得到更为精确的变形网格. 利用该方法完成带襟翼的NLR 7301 二维构型和带发动机短舱的DLR F6 翼身组合体的网格变形问题,结果表明该方法可以较大幅度降低变形网格的物面误差,并且有效避免边界层网格交错问题.  相似文献   

4.
二维弹塑性有限元程序全自动重分   总被引:1,自引:1,他引:0  
对二维弹塑性流体力学有限元计算中的非结构三角网格 ,提出了一种简单、实用的全自动网格重分方法。并在LTZ 2D程序上编制了全自动重分模块 ,重分使拉格朗日流体程序能计算包含大变形的一些问题。给出两个算例 ,圆柱体撞击 (Taylor实验 )和射流形成 ,通过网格重分数值计算结果与实验和文献结果比较 ,定性、定量均有较好的符合。  相似文献   

5.
计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM)的深入发展.为了提高EFGM的计算速度,本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法.该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数;基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变,根据广义光滑伽辽金弱形式建立系统离散方程.两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成.为了提高方法的精度,由Richardson外推法确定两层光滑网格上的最优光滑应变.几个数值算例验证了该方法的精度和计算效率.数值结果表明,随着光滑积分网格数目的增加,光滑无网格伽辽金法的计算精度逐步接近EFGM的,但计算效率要远远高于EFGM的.另外,光滑无网格伽辽金法的边界条件可以像有限元那样直接施加.从计算精度和效率综合考虑,光滑无网格伽辽金法比EFGM具有更好的数值表现,具有十分广阔的发展空间.  相似文献   

6.
周蕊  李理  田保林 《力学学报》2023,(11):2675-2692
凝聚炸药爆轰驱动惰性金属材料形成的多介质流动问题广泛存在于工程应用领域, Lagrange方法由于其物质界面的高保真特性一直在相关问题的数值模拟中发挥着不可替代的作用.加密网格是提高爆轰驱动多介质问题模拟精度的常用途径之一.然而, Lagrange框架下整体密网格计算常会遇到网格畸变、计算效率低等问题.为此,针对爆轰驱动多介质流动问题,提出了一种Lagrange框架下的非结构网格多层自适应方法,在保证所关心区域局部网格分辨率的前提下,大幅缩减了整体计算规模,提升了Lagrange计算的健壮性.设计了非结构网格多层数据结构,提出了多层网格分层存储、有效网格压至一层进行Lagrange计算的AMR策略,同时还发展了自适应接触滑移耦合算法,实现了AMR计算与多分区接触滑移计算的“紧耦合”.相比于已有工作,所提出的AMR方法既保持了非结构网格多层数据结构的灵活性优势,又避免了Lagrange框架下多层网格分别计算带来的层间耦合困难,同时因实现了与接触滑移的自适应耦合,使得它能很好地适应多分区的多介质问题.在一维、二维爆轰算例验证所提出方法正确性的前提下,开展了拐角爆轰、多层炸药隔爆和有限尺寸弯...  相似文献   

7.
多维微分求积DC(differential cubature)法是近年迅速发展起来的一种处理多维微分方程的高精度数值方法。本文围绕多元函数插值适定性这一基本问题,探讨了实现DC法时插值结点组和插值空间应满足的条件,研究了离散结点组选择及其构造适定试函数的方法。提出了一种沿各坐标方向可以布置不同数量结点的广义交错网格DC方法,并给出了确定相应的插值空间及其选择试函数的具体方案。通过弹性薄板变形分析研究了广义交错网格DC法解决实际结构力学问题的可靠性及其潜力。研究结果表明,广义交错网格DC法较传统的交错网格DC格式具有更强的适用性和灵活性,采用数量少得多的离散结点即可达到传统交错网格DC法相同的数值精度,表现出精度高和计算工作量小的优点。  相似文献   

8.
马文涛 《力学学报》2018,50(5):1115-1124
计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM) 的深入发展. 为了提高EFGM 的计算速度, 本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法. 该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数; 基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变, 根据广义光滑伽辽金弱形式建立系统离散方程. 两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成. 为了提高方法的精度, 由Richardson外推法确定两层光滑网格上的最优光滑应变. 几个数值算例验证了该方法的精度和计算效率. 数值结果表明, 随着光滑积分网格数目的增加, 光滑无网格伽辽金法的计算精度逐步接近EFGM 的, 但计算效率要远远高于EFGM的. 另外, 光滑无网格伽辽金法的边界条件可以像有限元那样直接施加. 从计算精度和效率综合考虑, 光滑无网格伽辽金法比EFGM具有更好的数值表现, 具有十分广阔的发展空间.   相似文献   

9.
以映射法为基础并结合网格划分经验,提出了原理简单的六面体网格生成新办法.该方法根据物体轮廓选择初始网格,设定表面结点强制变形到目的曲面,经由有限元弹性计算确定内部节点的位置.在检查全体单元质量以后,调整畸形单元从而生成目的网格.通过为一个复杂的马头门模型构建全六面体网格,最后证明了本文所述方法的可行性.  相似文献   

10.
物质点法MPM(Material Point Method)是无网格方法之一.它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应交;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格.本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析.  相似文献   

11.
利用主元分析法对数字图像相关技术所测结构变形信息进行数据压缩,并进一步用于力学模型未知参量的反求计算。首先,为降低数字图像相关技术所测庞大数据的应用成本,提出利用主元分析法对结构表面变形场数据进行压缩,实现在保留结构表面变形信息主要特征的前提下显著降低数据量的目的;其次,针对压缩后的数据建立了基于最小二乘法的力学模型参量反求模型,并利用高斯牛顿法进行求解;最后,以具体算例从计算精度、收敛速度和抗噪性等方面验证了数据压缩对模型参量反求的效果。研究结果表明,所提方法在显著降低使用数据量的前提下,能够有效提高力学模型参量反求计算的收敛速度,特别是对于包含多个模型参数的反求问题,具有较高的精度和较好的稳定性。  相似文献   

12.
宓思恩  刘小明  魏悦广 《力学学报》2021,53(11):3080-3096
提出了一种从离散分子动力学模拟(MD)到连续介质弹性有限元计算分析(FEA)的过渡方法, 简称MD-FEA方法. 首先通过MD计算获得晶体材料原子的移动位置, 然后根据晶体结构的周期性特征构造连续介质假设下的有限单元变形模型, 进一步结合材料的力学行为本构关系获得应变和应力场. 为了检验MD-FEA方法的有效性, 将该方法应用于详细分析Al-Ni软硬组合两相材料纳米柱体的拉伸变形问题和基底材料为Al球形压头材料为金刚石的纳米压痕问题. 采用MD-FEA方法获得了上述两种问题的应力?应变场, 并将计算结果分别与传统MD方法中通过变形梯度计算的原子应变以及原子的位力应力进行了比较, 详细讨论了用MD-FEA方法计算的应力?应变场与传统MD原子应变和位力应力的区别, 并对MD-FEA方法的有效性及其相较于传统MD方法所具有的优势进行了探讨. 结论显示, MD-FEA方法与传统MD方法在应力?应变变化平缓的区域得到的结果接近, 但在变化剧烈的区域以及材料的表/界面区域, MD-FEA方法能够得到更加精确的结果. 同时, MD-FEA方法避免了传统MD方法中, 需要人为选取截断半径以及加权函数所导致的误差. 另外, 当应变较大时, MD-FEA方法计算的小应变与传统MD方法计算的格林应变存在一定差异, 因此, MD-FEA方法更适合应变较小的情形.   相似文献   

13.
薄层弹性体功能结构表面在工程领域有着广泛的应用前景,薄层弹性体的接触特性分析对功能结构表面设计与应用非常重要.采用随机粗糙表面模拟生成技术和有限元分析技术建立了粗糙薄层弹性体表面的接触刚度确定性分析模型,研究了薄层弹性体的接触特性,对基于半数值确定性分析方法的粗糙层-基体层串联模型用于粗糙薄层弹性体接触刚度计算的适用性进行了讨论,进一步分析了粗糙薄层弹性体串联模型接触刚度的误差来源.研究结果表明:当薄层弹性体的基体层厚度小于10倍表面均方根粗糙度时,由传统粗糙层-基体层串联模型分析获得的粗糙薄层弹性体刚度误差将超过15%,误差主要来源于粗糙峰的大变形和基体层的局部变形不均匀的共同作用.  相似文献   

14.
海上箱筒型基础结构下沉到位后,舱内土体的受力变形可以看作是侧限压缩模型问题,上部荷载引起的侧壁摩擦力与土体应力之间的关系较为复杂。本文通过简化模型,列出了侧限压缩模型的平衡微分方程,分析了方程解析计算时存在的难点,结合摩擦应力边界条件的特点,构造函数迭代法对该方程进行数值求解,展示了计算的收敛过程,得到接触侧面上水平应力和土体表面位移的非线性分布规律,并将其结果与ABAQUS数值模拟结果进行了对比。结果表明,函数迭代法收敛速度较快,计算过程稳定,对其他存在变量耦合边界条件的微分方程数值求解有一定的参考价值。  相似文献   

15.
晶体塑性变形离散滑移模型及有限元分析   总被引:1,自引:0,他引:1  
基于韧性单晶体实验现象,建立了描述晶体塑性变形的离散滑移模型.该模型的主要特点是:晶体滑移变形在宏观上是不均匀的,滑移带的分布是离散的.利用晶体塑性理论对模型进行了有限变形有限元分析,计算结果揭示了晶体滑移的离散行为,模拟的应力 应变曲线与实验曲线相吻合  相似文献   

16.
The motion and deformation of soft particles are commonly encountered and important in many applications. A discrete element-embedded finite element model (DEFEM) is proposed to solve soft particle motion and deformation, which combines discrete element and finite element methods. The collisional surface of soft particles is covered by several dynamical embedded discrete elements (EDEs) to model the collisional external forces of the particles. The particle deformation, motion, and rotation are independent of each other in the DEFEM. The deformation and internal forces are simulated using the finite element model, whereas the particle rotation and motion calculations are based on the discrete element model. By inheriting the advantages of existing coupling methods, the contact force and contact search between soft particles are improved with the aid of the EDE. Soft particle packing is simulated using the DEFEM for two cases: particle accumulation along a rectangular straight wall and a wall with an inclined angle. The large particle deformation in the lower layers can be simulated using current methods, where the deformed particle shape is either irregular in the marginal region or nearly hexagonal in the tightly packed central region. This method can also be used to simulate the deformation, motion, and heat transfer of non-spherical soft particles.  相似文献   

17.
The size dependent deformation of Cu single crystal micropillars with thickness ranging from 0.2 to 2.5 μm subjected to uniaxial compression is investigated using a Multi-scale Dislocation Dynamics Plasticity (MDDP) approach. MDDP is a hybrid elasto-viscoplastic simulation model which couples discrete dislocation dynamics at the micro-scale (software micro3d) with the macroscopic plastic deformation. Our results show that the deformation field in these micropillars is heterogeneous from the onset of plastic flow and is confined to few deformation bands, leading to the formation of ledges and stress concentrations at the surface of the specimen. Furthermore, the simulation yields a serrated stress–strain behavior consisting of discrete strain bursts that correlates well with experimental observations. The intermittent operation and stagnation of discrete dislocation arms is identified as the prominent mechanism that causes heterogeneous deformation and results in the observed macroscopic strain bursts. We show that the critical stress to bow an average maximum dislocation arm, whose length changes during deformation due to pinning events, is responsible for the observed size dependent response of the single crystals. We also reveal that hardening rates, similar to that shown experimentally, occur under relatively constant dislocation densities and are linked to dislocation stagnation due to the formation of entangled dislocation configuration and pinning sites.  相似文献   

18.
大面积金属板材304L/Q235B的爆炸焊接过程涉及炸药爆轰、金属板材的高速碰撞和塑性变形等。采用有限元法计算模拟这个问题时,网格单元会发生扭曲畸变现象,导致计算精度下降,甚至出现单元负体积而使计算终止,并且炸药爆轰形成气体产物飞散过程也很难模拟。为了能模拟大面积金属板材的爆炸焊接整个过程并获得合理的技术工艺参数,采用物质点法进行三维数值模拟分析。物质点法作为一种无网格法,在模拟冲击动力学问题中主要采用显式积分算法。通过将拉格朗日质点单元与固定的欧拉背景网格相结合,可以实现爆炸焊接的复板与基板的高速碰撞、炸药滑移爆轰、金属板面的塑性变形过程的数值模拟,并给出爆炸复合板材的形变、有效塑性应变和复板与基板的碰撞速度的计算结果。采用物质点法模拟的复合板材变形与爆炸焊接实验结果基本一致。计算复板与基板的碰撞速度这个重要的物理参数时,物质点法与Richter理论公式的相对误差不超过13%。数值计算和实验结果表明,物质点法在数值精度和计算效率方面具有优势,物质点法是研究金属焊接爆炸问题的一种有效数值方法。  相似文献   

19.
澜沧江某水电站处于青藏高原东部边缘地带,属于高山峡谷地形地貌,高地应力环境,岩体卸荷裂隙很发育,使得倾倒变形和岩体质量、断裂活动及地震构造一样成为影响工程边坡岩体稳定的主要因素。针对工程边坡的大变形问题可采用离散元的数值模拟分析方法。通过建立理论开挖和工程边坡开挖离散元模型,可分别得出倾倒变形破坏机理发展过程为初期弱倾倒变形岩体的层内剪切错动、强倾倒变形岩体的层内拉张变形、强倾倒变形岩体的切层张-剪破裂及极强倾倒破裂岩体的折断张裂(坠覆)破裂和工程边坡的变形范围、确定开挖面及加固方式等。通过工程边坡模型的计算结果和现场地质调查成果的比较表明,计算结果和实际情况基本吻合。  相似文献   

20.
过量开采地下承压水会引起地基土体的变形,进一步导致地面沉降和地裂缝等地质灾害,抽水引起地层的变形过程是土层中孔隙水和土颗粒相互作用的复杂过程.两步法和流固耦合法是计算地面沉降常用的两种思路,它们的适用性各有不同.本文基于理想的地层模型,分别建立两步法计算模型和流固耦合计算模型,利用FLAC3D6.0有限差分软件对模型进行求解,对比分析两类方法计算结果的差异.结果表明流固耦合法符合力学机制,能够更好地反映土体变形的发展过程,但计算时间更长.而两步法适用于稳定流,在渗流达到稳定流时两步法计算的孔压分布和地面沉降能够满足工程精度要求.本文的研究成果为合理选取地面沉降计算方法提供依据.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号