c语言如何输入矩阵_如何在 COMSOL 软件中调试外部材料

问题描述

我实现了一种用于“结构力学”的外部材料,但我的模型无法求解。

解决方法

外部材料一般通过 C 语言编码,然后被编译到共享库中,供 COMSOL Multiphysics 在运行时使用。如果您在编译外部材料代码时遇到问题,请参阅 COMSOL Multiphysics 安装手册中的如何编写和连接外部材料模型一节。外部材料的用途有两个方面:
  • 计算用于确定内力和外力是否达到平衡的机械应力。
  • 计算应力对应变的导数(雅可比),从而使方程求解器可以迭代以建立力平衡。
当 COMSOL Multiphysics 使用外部材料时,参数和状态(内部变量)将在程序和外部材料函数之间传递,此过程必须以一致的方式进行。例如,COMSOL Multiphysics 中定义的状态数必须与您在外部材料中使用的状态数相匹配。请参阅 COMSOL Multiphysics Reference Manual 中的 Working with External Materials 一节。此外,您也可以参考本页下方链接中提供的示例,尝试编译一个外部材料示例,并在简单的案例中使用它,以确保您能够正确用它来建模。如果模型无法求解或者不收敛,请尝试确定是应力计算还是雅可比矩阵计算出了问题。为了更容易测试您的外部材料代码,请执行以下操作:
  • 建立一个尽可能小的模型,最好使用单一的网格单元。
  • 对模型施加位移控制的载荷。这样,您就可以精确控制输入外部材料代码中的应变。此外,如果您要实现的材料模型包含软化,请务必使用位移控制的模型。
  • 使用简单的变形情况,例如单轴拉伸和纯剪切。

应力计算

您需要对一些与应力计算有关的事项进行检查:
  • 如果外部材料在单轴载荷下表现良好,但在其他变形模式下失效或产生不正确的结果,请确保在外部材料代码中使用了正确的剪切应变分量定义和顺序。
  • 您的代码中是否存在被零除的风险?在分析开始时,外部材料可能会调用零位移场(零应变),因此需要添加保护性代码,以避免可能出现的问题。
  • 许多非弹性材料模型都需要局部迭代来计算应力,这通常被称为应力更新算法。COMSOL Multiphysics 不会去控制这些局部迭代,因此,密切关注您实现的迭代方法的健壮性非常重要。是否存在数值溢出的风险?在蠕变模型中,应力通常大于一次幂关系。在局部迭代过程中,未收敛的应力可能相当高,以至于提高到幂指数关系时产生溢出。在代码中添加安全措施来识别这些类型的问题会很有帮助。
  • 如果您认为应力更新算法有问题,可以使用 C 语言中的 fprintf 语句将局部迭代中的信息打印到单独的文件中。但是要小心,这可能会产生大量的输出。

雅可比矩阵计算

雅可比矩阵不正确的一个标志是,模型在分析过程中的一部分运行没有问题,但随后会出现收敛问题。这种现象表明雅可比矩阵的非弹性部分不正确,您可以使用以下几种方法来识别并纠正与雅可比矩阵计算有关的问题:
  • 首先必须指出的是,在许多情况下,一个稍微偏离正确雅可比矩阵的近似雅可比矩阵通常足以达到收敛。这只会影响收敛速度,而不影响解的精度。
  • 导致收敛问题的最常见原因是没有正确定义雅可比矩阵。从形式上看,雅可比矩阵是第二类皮奥拉-基尔霍夫应力对格林-拉格朗日应变(或对变形梯度)的导数。必须计算导数,使其与计算应力所采用的应力更新算法一致。
  • 和应力更新算法一样,请密切关注剪切应变的定义和顺序。
  • 在求等效应力等物理量的微分时,一个常见的陷阱是忽略应力张量的对称性。在使用压缩的矢量格式处理应力和应变时,可能会出现问题。
  • 应力对应变的微分可能相当复杂(且容易出错),您可以尝试对雅可比矩阵表达式进行受控简化。
  • 如果雅可比项的表达式非常复杂,您可以使用 fprintf 语句将值打印到单独的文件中,然后与手动计算的值进行比较。但是请小心,这可能会产生大量的输出。

一般建议

由于应力计算或雅可比矩阵计算存在问题,有时很难确定外部材料代码是否失败。基于应变场的数值扰动来实现雅可比矩阵可能是一种有效的做法。Miehe 于 1996 年提出了一种扰动技术,用六个额外的应力计算来代替雅可比矩阵计算。这项技术在计算上可能非常昂贵,但在外部材料开发过程中非常有用。此外,还请注意,本软件针对本构建模中经常出现的许多张量和矢量运算提供了一个实用程序函数库,可以帮助您减少自己编写相应功能代码所花费的时间。请参阅 Structural Mechanics Module User's Guide 中的 Using External Materials 一节。关于如何在 COMSOL Multiphysics 中实现外部材料的演示,我们提供了许多资料供您参考,例如:外部材料示例,结构力学 演示模型,COMSOL Multiphysics 软件提供一个丰富的材料模型集合,在每个版本中,我们都会不断添加更多的内置材料模型,当前版本中现在可能已经包含您正在作为外部材料编码的材料特性。
参考文献

Miehe, C., "Numerical computation of algorithmic (consistent) tangent moduli in large-strain computational inelasticity", Comp. Meth. Appl. Mech. Engng., Vol 134, No. 3-4, pp. 223-240, Aug 1996.

本文内容来自 COMSOL 知识库,点击“”,获取实现外部材料的演示的案例模型。

237e327bff39167111341270db995d73.gif

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

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

相关文章

lwip网络通信socket_lwIP在Socket模式下接口:BSD Socket API

http://home.eeworld.com.cn/my/space-uid-162102-blogid-52270.html2011原文:关于BSD Socket API在网上找到的两个网站,是关于BSD Sockets API的,这是与lwIP在Socket模式下兼容的。里面对API函数做了较为详细的介绍,先记下来&…

idal 创建springboot 项目_SpringBoot教程:Maven方式创建SpringBoot项目

一.Maven方式创建Spring Boot项目1.配置Maven环境在以maven方式创建Spring Boot项目之前,请先确保自己的Maven并配置好环境变量,并且在ideal中关联好maven环境。2.创建一个新的maven项目3.设置项目坐标信息到此为止,就已经成功的…

mysql xa_Mysql对XA的支持

mysql8.0文档:https://dev.mysql.com/doc/refman/8.0/en/xa-statements.html。13.3.8.1 XA Transaction SQL Syntax章节讲述了Mysql对于XA事务的语法。XA {START|BEGIN} xid [JOIN|RESUME] XA END xid [SUSPEND [FOR MIGRATE]] XA PREPARE xid XA COMMIT xid [ONE P…

mysql创建新的表_Mysql 创建表

Mysql 创建表1.首先打开终端输入如下命令:/usr/local/MySQL/bin/mysql -u root -p其中root为用户名。2.这时会出现如下命令:Enter password:此时如果你没有改密码,直接敲回车。否则,输入你的密码。这样就可以访问你的数据库服务器…

mysql gfs2_Mysql_HA+IP_SAN+Clvmd+GFS2

Mysql_HAIP_SANClvmdGFS2GFS:全局文件系统,支持DLM锁管理,支持并发写,Centos6不超过16个节点适用于HA集群使用的共享存储适用于LB集群RealServer使用的共享存储mysql高可用集群架构设计1、节点主机系统:Centos6.52、高…

mysql 使用sum limit_mysql踩坑记录之limit和sum函数混合使用问题

问题复盘本次复盘会用一个很简单的订单表作为示例。数据准备订单表建表语句如下(这里偷懒了,使用了自增ID,实际开发中不建议使用自增ID作为订单ID)CREATE TABLE order (id int(11) NOT NULL AUTO_INCREMENT COMMENT ‘订单ID‘,amount decimal(10,2) NOT…

mysql 视图 过程 函数_MySQL视图,函数,触发器,存储过程

1. 视图视图是一个虚拟表,它的本质是根据SQL语句获取动态的数据集,并为其命名,用户使用时只需使用【名称】即可获取结果集,可以将该结果集当做表来使用。使用视图我们可以把查询过程中的临时表摘出来,用视图去实现&…

ubuntu安装mysql5.7.17_ubuntu 16.04安装mysql-server_5.7.17

在MySQL官网上下载最新版的Ubuntu Linux专用的MySQL。我这里下载的是:mysql-server_5.7.17-1ubuntu16.04_amd64.deb-bundle.tar也可从以下mirrors下载:解压文件 命令为:rootubuntu:/# tar -zxvf mysql-server_5.7.17-1ubuntu16.04_amd64.deb-…

socket抓包_64、抓包分析tcp与udp

从前面的两个案例,我们了解到了如何通过原生socket函数分别创建tcp和udp服务,以及通过相应的客户端进行连接测试。在本文中,我们将继续深入地去了解tcp和udp的差别,和思考它们两者如何应该多个客户端请求1、抓包分析这里采用tcpdu…

c#水晶报表连接mysql,如何将数据库绑定到水晶报表?

i have C# program, i have DataSet that i want to bind to crystal report.how to do it ?解决方案Create a Typed DataSet from the "Add New Item > DataSet and give the dataseta useful name.In Server Explorer, drag a view with the fields you want in the…

fw313r手机登录_迅捷(FAST)fw313r路由器手机设置教程

本文中,鸿哥主要给大家介绍,迅捷(FAST)fw313r路由器用手机设置的方法。一台新买回来的迅捷(FAST)fw313r路由器,要用手机来设置它连接Internet上网,需要以下几个步骤:1、正确连接FW313R路由器2、手机连接FW313R的信号3、…

mysql 5.6.21不能选择安装路径_mysql5.6.21服务器安装图解(选择安装和数据目录)

mysql安装图解(服务器上安装mysql,只安装mysql服务)1.3 mysql安装及配置版本:MySQL5.6.21 64位1.3.1 安装.net framework.exe 4.0注意:mysql5.6版本需要安装.net framwork.exe 4.0环境1.3.2 安装mysql数据库(1)打开安装程序后,…

mysql替换sql中rank函数_MySQL sql Rank()函数实现

一字符串类 Concat函数:连接字符串 Instr函数:返回字符串在某一个字段的内容中的位置, 没有找到字符串返回0,否则返回位置(从1开始) 字符串大小写转换[upper()、ucase()和lower()、lcase()] upper()和ucase():把字符串中的小…

python决策树生成规则_ID3决策树(python实现)

决策树的生成(该函数是一个递归的过程)CreateTree输入:数据集、特征输出:字典型数据——决策树a、判断是否满足停止划分的条件若当前数据集的属性值为空,则投票表决当前样本中最多的类别若当前所有的样本类别相同,则返回当前数据的…

while循环python的范围_python-无法在while循环中从列表中排除一定范围内的项目

所以我发布了一个问题before,但是它过于简化了,正确地被标记为重复.我现在将更详细地发布我的问题,以便希望可以解决我的问题.简而言之如下:我有两个列表:a [10.0,20.0,25.0,40.0]和b [1.0,10.0,15.0,20.0,30.0,100.0]使用列表推导,我想从b中排除a中指…

怎么使用mysql打表_MySQL的表使用

-- 创建表CREATE TABLE teacher(id INT,NAME VARCHAR(20))-- 查看所有表SHOW TABLES;DESC student;DROP TABLE student;CREATE TABLE student(id INT,NAME VARCHAR(20),gender VARCHAR(2),age INT)-- ********一、增删改数据********* ----- 1.1 增加数据-- 插入所有字段。一定…

表级锁的mysql读写_Mysql的表级锁

我们首先需要知道的一个大前提是:mysql的锁是由具体的存储引擎实现的。所以像Mysql的默认引擎MyISAM和第三方插件引擎 InnoDB的锁实现机制是有区别的。可根据不同的场景选用不同的锁定机制。Mysql有三种级别的锁定:表级锁定、页级锁定、行级锁定一、定义…

account表里有什么 银行_模拟一个银行账户类Account,账户类中包括所有者、账号、余额、账户总数、存款、取款等信息。_学小易找答案...

【单选题】廉价磁盘冗余阵列RAID利用冗余技术实现高可靠性,其中RAID1的磁盘利用率为() 。【简答题】与直流调速相比,交流调速有何优点?【论述题】方法【单选题】以《老子注》一书驰名汉学界的学者是下列哪一位 ______【填空题】本任务中,所需的实训设备有:( )、( )、( )【单选…

创建mysql视图语法正确的是_MySQL创建视图的语法格式

视图,具有简化查询语句、安全性和保证逻辑数据独立性等作用创建视图的语法格式视图中,包含SELECT查询的结果,因此,视图的创建基于SELECT语句,和已经存在的数据表,视图可以建立在一张表上,也可以…

爬空气质量MySQL_爬虫:利用selenium采集某某环境网站的空气质量数据

前言:在上一篇文章中,我们介绍了在http://PM2.5.in这个网站采集空气质量的数据,本篇文章是对其产生的一些问题的另一种解决方案,提供更加权威的数据采集。技术框架:selenium、json、etree这里的selenium是一种自动化测…