模型用于待建工程计算前,应进行天然河道率定。最大输沙率计算值为0.16615kg/s,发生在7920s,沙峰试验值较计算值滞后。最大输沙率试验值与计算值之比1.69。验证计算结果见图9-25。......
2023-06-22
1985年,Seaman 和Curran 采用统计平均的方法,建立了模拟准脆性材料动态拉伸断裂与破碎的微裂纹细观损伤模型,得到了岩石材料在加载过程中的裂纹尺寸及分布的演化[14]。后来,Seaman 在此基础上建立了黏性内损伤本构模型,用于模拟固体推进剂和炸药在冲击作用下的黏弹塑性变形、拉伸断裂和破碎。在Seaman 工作的基础上,Olsen 等提出了新的黏性内损伤模型(Viscous Internal Damage,VID),可考虑黏弹塑性力学行为,以及微裂纹的成核、成长、聚合、强度损失、破碎等动态拉伸损伤。在此基础上,Mathe⁃son[20~23]提出了损伤与化学反应耦合模型 (Coupled Damage and Reaction,CDAR)。CDAR 模型包括一个黏弹塑性模块和一个拉伸损伤膨胀模块(Ten⁃sile Damage and Distension,TDD),可以描述含能材料内部脱黏颗粒附近孔洞变化以及黏结剂开裂产生的微裂纹的张开和闭合。
1.黏弹塑性模块(VEP)
VEP 模块由一个多参数的黏弹性元件和一个黏塑性元件组成,用来计算偏应力,其中的黏弹性元件由一个弹性元件和N 个Maxwell 体并联而成,可以表示为
式中:Sij为偏应力;为长时剪切模量;Gi和Fi分别为第i 个Maxwell 体的剪切模量和黏性系数;Dt 为微分算子d/dt;为总的偏应变;为黏塑性偏应变,描述非弹性变形,,evp为标量黏塑性偏应变率,evp可由下式表示:
式中:Seff为von Mises 有效偏应力;A 为缺陷尺寸;N0为初始缺陷密度;M 为塑性应变系数;σDr为缺陷拉应力;H 为硬化系数;n1为率相关程度指数;n2为应变硬化指数。
2.拉伸损伤膨胀模块(TDD)
TDD 模块主要考虑了含能材料中的颗粒与黏结剂的界面脱黏和黏结剂的开裂两类损伤。损伤的演化根据状态方程计算得到的压力和由VEP 计算的应力偏张量来计算。假设材料内部的裂纹尺寸符合指数函数分布,平均裂纹尺寸可以表示为式中:R0为初始夹杂半径;Rd为脱黏引起的裂纹变化;Rs为黏结剂开裂引起的裂纹变化。
R=R0+Rd+Rs
假设脱黏先发生,脱黏引起的裂纹增长速率表示为
式中:c0为体声速;fc为系数(0~1);σeff为有效主应力;为脱黏应力阈值;为控制裂纹增长速率的脱黏特征应力;εeff为有效主应变;为脱黏应变阈值;为控制裂纹增长速率的脱黏特征应变。
脱黏完全发展后,发生黏结剂开裂损伤,其产生的裂纹增长速率表示为
式中:为黏结剂开裂应力阈值;为黏结剂开裂特征应力;为黏结剂开裂应变阈值;为黏结剂开裂特征应变。
对于黏结剂开裂损伤,裂纹增长的应力阈值为
式中:R 为平均裂纹半径;Ks为静态断裂韧性;Kd为动态断裂韧性;τK为控制静态向动态转变速率的参数。
假设损伤度与裂纹尺寸的立方成正比,其损伤为
式中:为固体参考密度;n0为单位体积的夹杂数量。TDD 模型中对膨胀的描述是在Baer 等和Drumheller 等的多相连续混合理论基础上进行了扩展,可以描述脱黏颗粒附近孔洞以及黏结剂开裂产生的微裂纹的张开和闭合。体积膨胀分为弹性和塑性两部分,弹性体积膨胀应变率可以表示为
式中:μβ为体积膨胀黏度;γs为材料的真实密度;ps为固体压力;βs为构形应力。
塑性体积膨胀应变率模型比较复杂,包括两个黏塑性过程:一个是高度体积膨胀材料在很低应力下发生的孔隙不可逆闭合;另一个是低体积膨胀材料在应力松弛条件下孔隙的闭合和张开。塑性体积膨胀应变率表示为
有关高聚物黏结炸药损伤断裂的数值模拟的文章
模型用于待建工程计算前,应进行天然河道率定。最大输沙率计算值为0.16615kg/s,发生在7920s,沙峰试验值较计算值滞后。最大输沙率试验值与计算值之比1.69。验证计算结果见图9-25。......
2023-06-22
1978年,美国著名运筹学家A.Charnes与W.W.Copper等人开创的效率评价新方法——数据包络分析方法。由于农户可以控制投入,但难以控制产出,因此,本研究采用VRS条件下的DEA估计小型水利设施的技术效率。......
2023-06-28
回归模型分析预测法是根据事物内部因素变化因果关系来预测事物未来发展趋势。4)利用模型进行预测。图6-11 社会总产值预测值置信区域示意图6.3.3.2 多元线性回归预测法当预测对象Y受到多个因素X1、X2、…、Xm的影响,如果各个因素Xi与Y相关关系可以同时近似用函数关系来表示,则可以用多元线性回归预测模型进行预测。......
2023-07-15
炮控系统是一个较为复杂的系统,在研究过程中,可以先以炮塔电动机和其功放部件晶体驱动器为被控对象,进行研究。用示波器观测测速电动机的输出得到控制结果,检验自抗扰控制器的控制性能。图5.4系统原理2.高低向稳定器建模以图5.5所示高低向电液炮控伺服系统为对象讨论并建立其数学模型[106-110],它由陀螺仪、放大器、液压系统、火炮等部件组成。......
2023-06-24
炸药颗粒采用各向同性的弹塑性本构模型。式是根据松弛条件本构方程,通过将一点的应变分解为应变球张量和应变斜张量两部分推导得到的。用Prony 级数表示黏弹性属性的基本形式为式中:和Gi为剪切模量;和Ki为体积模量;和为各Prony 级数分量的松弛时间。表7-7剪切松弛模量与松弛时间的关系......
2023-06-27
结构模型解析法是美国Bottele研究所华费尔特教授于1973年在分析复杂社会经济系统问题时提出的。4)根据要素明细表作构思模型,并建立邻接矩阵和可达矩阵。5)对可达矩阵进行分解后建立结构分析模型。下一步利用上面获得的可达矩阵迸行系统结构的层次分析。......
2023-07-15
通过ADAMS软件提供的图形接口模块,将三维实体造型软件UG中建立的机构三维模型导入到ADAMS软件中的。考虑到所关心的问题和仿真的方便,对操作机构作适当简化后,用ADAMS软件包建模。约束关系的添加要充分考虑实际操作机构中可能存在的约束,遗漏任何一个约束关系都会严重影响到输出结果。3)仿真过程的控制 在ADAMS软件中,为了实现对仿真过程的控制,需要添加传感器来控制仿真过程的自动终止。......
2023-06-15
相关推荐