不等式约束的拉格朗日乘数法_Abaqus血管支架仿真|接触约束执行方式


根据支架的建模基础、几何和网格划分、单元选择、材料和截面属性、分析步设置、载荷接触和约束、收敛调整,以及后处理、参数优化等内容。接着Abaqus血管支架仿真|建模基础、Abaqus血管支架仿真|几何和网格划分、Abaqus血管支架仿真|单元选择、Abaqus血管支架仿真|材料模型和参数、Workshop: Abaqus单胞支架的材料赋值、Abaqus血管支架仿真|分析步设置(上)、Abaqus血管支架仿真|分析步设置(中):静态分析、Abaqus血管支架仿真|分析步设置(下):动态分析、Abaqus血管支架仿真|载荷和边界条件、Abaqus血管支架仿真|通用接触、Abaqus血管支架仿真|接触对和表面行为,继续推送,今天讲解接触的实施方法。

Abaqus/Standard有三种数值方法,用以实现或近似“硬Hard”接触条件:

  • 直接执行方法:用拉格朗日乘子法Lagrangemultiplier严格执行压贯Pressure-penetration关系

  • 罚函数法:使用罚刚度近似执行

  • 增广拉格朗日法:采用增进迭代的罚函数法,近似执行

8033f209d2ac82ddd69440b647147a02.png

▲ 接触约束的实施方式

直接执行 Direct enforcement拉格朗日乘数法,在方程组中加入约束方程和拉格朗日乘数

a53bceacdfdc0670e0dc4c6cd0eda71b.png

▲ 直接法

优点: 准确—满足精度约束缺点: 在每个接触约束上增加变量,扩大了待解方程组,同时,限制了稀疏求解器的消去顺序,降低了算法的性能,故而增加了求解方程的成本;从零接触刚度(非活动接触时)突变为无限接触刚度(活动接触时),具有潜在的收敛困难;接触约束与MPC的重叠等,过约束的收敛困难。罚函数法 Penalty method罚函数法是硬接触的一种严厉逼近法。

56f8e1b4a196c19ab32287c75c87c7ca.png

▲ 罚函数法

优点: 显著提高收敛速度;具有更好的方程求解性能,没有拉格朗日乘子自由度,除非接触刚度非常高;对重叠约束的良好处理。缺点:少量渗透,通常微不足道;某些情况,可能需要调整惩罚刚度。

罚刚度的接触刚度行为有线性和非线性。其中,线性更容易收敛,更适合解决涉及严格接触的问题;非线性具有较低的初始刚度,更适合涉及抖动的问题,较高的最终刚度有助于减少穿透,当然整体收敛性可能更加困难。

80552b43c8d715f9f873e896fb150852.png

▲接触刚度行为

如支架,弯曲为主的问题分析: 默认的惩罚性刚度通常可以缩小两个数量级,而没有任何明显的精度损失;减小弯曲占主问题的惩罚刚度,有时会显著提高收敛速度。766bb776cded986cd8a62f430ba118d1.png

3213ab0834dfce7bd14bf303407e2e62.png

▲ 三点折弯试件

举例:支架扩张分析

42c636b5f02f73379108a5e8538766dd.png

▲ 支架扩张分析

ac18ebc01a7f8952db4c36469c00df3c.png

推荐使用接触实施方法的设置,对有收敛困难的支架进行分析。

e590798d05df7bb03768ca387caeb3d3.png

▲ 支架分析推荐使用的接触属性设置

总结采用罚函数法:减小默认惩罚刚度,以提高收敛速度,同时,非线性惩罚刚度,更加容易建立接触。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/350521.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

OpenCV形态学运算

1.OpenCV形态学运算morphologyEx(开运算 、闭运算 、形态学梯度 、顶帽运算 、黑帽运算 、腐蚀运算 、膨胀运算 、击中击不中运算) void cv::morphologyEx (InputArray src,OutputArray dst,int op,InputArray kernel,Point anchor Point(-1,-1),int …

python判断成绩等级_Python计算山东新高考选考科目卷面原始成绩为等级成绩

原标题:Python计算山东新高考选考科目卷面原始成绩为等级成绩图书详情:配套资源:用书教师可以免费获取教学大纲、教案、课件、源码、习题答案、课堂管理与考试系统。山东省新高考政策33中,考生必考科目有语文、数学、英语&#xf…

java web ee_Java EE 6 Web配置文件。 在云上。 简单。

java web eeJava SE还可以。 Java EE是邪恶的。 这就是我一直想的。 好吧,现在不再了。 让我分享我的经验。 几周前,我开始考虑将旧版spring hibernate tomcat应用程序移植到新平台上: SAP NetWeaver云 。 我知道您在极客那里的想法&…

帆软报表(finereport) 饼图联动

饼图联动:点击饼图1,饼图2和饼图3显示饼图1的关联数据,接着点击饼图2,饼图3显示饼图2的关联数据,点击上方清除级联,饼图则恢复默认展示状态 下面以上图示例效果为例,说明制作过程。 1.为每个饼图…

中值滤波medianBlur

C: void medianBlur(InputArray src, OutputArray dst, int ksize) 参数解释: . InputArray src: 输入图像,图像为1、3、4通道的图像,当模板尺寸为3或5时,图像深度只能为CV_8U、CV_16U、CV_32F中的一个,如而对于较大…

wxpython下载缓慢_我可以在wxPython的wx.grid.Grid中加速优化GridCellAttr的使用吗?

设置单元格属性将新GridCellAttr添加到GridCellAttrProvider名单。 随着列表的增长,查找单元格的特定属性(通过遍历列表和比较坐标)变得越来越慢。您可以尝试通过PyGridTableBase.SetAttr和GetAttr(例如使用字典)实现自己的加快步伐:编辑:更新…

特性和混入不是面向对象的

让我立刻说,我们将在这里讨论的功能是那些迫切需要进行放线手术的人带给面向对象编程的纯粹的毒药 ,就像David West在他的《 Objecting Thought》一书中所建议的那样。 这些功能具有不同的名称,但最常见的是traits和mixins 。 我很不明白&…

模板匹配matchTemplate

C: void matchTemplate(InputArray image, InputArray templ, OutputArray result, int method); image:输入图像。必须为8位或者32位的浮点型。 templ:用于搜索的模板图像。必须小于输入图像并且是一样的数据类型。 result:匹配…

pandas处理csv

import pandas as pd from pandas import Series,DataFrame 打开csv文件 deli pd.read_csv(../得力.csv) 转换为DataFrame数据类型 df DataFrame(datadeli) 更改列名 df.rename(columns{"title":"商品名称","price":"价格","pr…

php 随机在文章中添加锚文本_SEO站长布局锚文本时的7大注意事项

通过对什么是锚文本链接、锚文本作用的讲解,以及怎样布局网站站内锚文本的介绍。相信大家都有了比较清晰的认识,也会重视网站锚文本的建设。但是菜鸟菌还是不得不再啰嗦点,锚文本虽然好处大大的,在布局锚文本的时候还是注意以下几…

OpenCV rectangle

void rectangle(Mat& img, Point pt1,Point pt2,const Scalar& color, int thickness1, int lineType8, int shift0) img 图像. pt1 矩形的一个顶点。 pt2 矩形对角线上的另一个顶点 color 线条颜色 (RGB) 或亮度(灰度图像 )(grayscale image…

【Java】NIO中Selector的select方法源码分析

该篇博客的有些内容和在之前介绍过了,在这里再次涉及到的就不详细说了,如果有不理解请看【Java】NIO中Channel的注册源码分析, 【Java】NIO中Selector的创建源码分析 Selector的创建在Windows下默认生成WindowsSelectorImpl对象,那…

简单小程序代码_小程序该如何运营

小型程序升级有三种最重要的方法:1.借助微信附近小程序功能实现客户扩展大家可能都知道,微信平台具有附近小程序的功能。此功能使用户可以查看周围五公里范围内的所有小程序。正是由于这个原因,一旦商人完成了这个小程序并加入了附近的小程序…

在javafx中界面主题_最小的JavaFX演示文稿(在JavaFX中)

在javafx中界面主题如果您要进行有关JavaFX的演示,那么使用JavaFX本身进行演示非常方便。 这样,您无需离开演示文稿即可轻松显示示例。 这是一个非常简单的示例。 在NetBeans中,设置一个新的JavaFX项目“ New Project”->“ JavaFX”->…

OpenCV cvLine

cvLine 函数原型: void cvLine( CvArr* img, CvPoint pt1, CvPoint pt2, CvScalar color, int thickness1, int line_type8, int shift0 ); img 图像。pt1 线段的第一个端点。pt2 线段的第二个端点。color 线段的颜色。thickness 线段的粗细程度。line_type 线段…

lengthOfLongestSubstring

lengthOflongestSubstring 试例: 输入: "abcabcbb" 输出: 3 解释: 因为无重复字符的最长子串是 "abc",所以其长度为 3。 /**** param {string} s* 无重复字符的最长子串*/ function lengthOfLongestSubstring(s) {if (typeof s ! s…

编译原理--词法分析C++

一、实验项目要求 1.实验目的 通过设计编制调试一个具体的词法分析程序,加深对词法分析原理的理解。并掌握在对程序设计语言源程序进行扫描过程中将其分解为各类单词的词法分析方法。 编制一个读单词过程,从输入的源程序中,识别出各个具有…

python list转map_Python 进阶之术 Map Filter Reduce

" 本文字数:763 字 || 阅读时间:3分钟"Map Map 会将⼀个函数映射到⼀个输⼊列表的所有元素上。这是它的规范:规范 map(function_to_apply, list_of_inputs) ⼤多数时候,我们要把列表中所有元素⼀个个地传递给⼀个函数&a…

JavaFX动画工具

好的,我想是时候让您讲一个小秘密了。 在过去的大约三个月中,我从事一个私有项目,目标是创建一个工具,使我可以轻松地为Java桌面应用程序创建动画。 JavaFX在API级别上提供了出色的动画制作支持,但对于初学者甚至中级程…

ImageWatch的使用

1.能够实时显示图像和矩阵Mat的内容,跟Matlab很像,方便程序调试,相当好用。 下载地址: https://marketplace.visualstudio.com/items?itemNameVisualCPPTeam.ImageWatch 2.下载好插件ImageWatch,双击ImageWatch.vs…