qq_46417931的博客

昨天上午老师突然发消息说晚上开个组会，我这周一直在看文献也没啥实质性进展，就想着随便做个模拟应付下老师。卧槽，没想到Abaqus提交作业后，一直不求解。具体问题如下​：时间紧急，加上有点慌张，一整就是两个多小时也没解决。后果可想而知。。。。。。因此中秋节我要修电脑了。

在结构力学模块方面，ANSYS Products可以模拟各种静态和动态载荷条件下的结构行为，包括线性和非线性分析、疲劳分析、接触分析等，并可以应用于机械、航空航天、汽车、建筑等行业。此外，在Ansys 2023 R2版本中，还优化了块记录功能以提高可读性、简化了光源和sensor等的定义方式、增强了GPU功能以更好地支持仿真模拟运行等。在流体力学模块方面，ANSYS Products可以模拟各种流体流动和传热问题，包括内部和外部流动、多相流动、热交换器等，并可以应用于各种行业。33、设置系统变量。

在包含用户定义输出变量规范的材料定义中，会在所有的元素材料计算点上调用该子程序；因为Abaqus/Standard迭代到收敛解，可能会在每个增量中的每个材料点上多次调用；会在每个步骤的每个增量中调用；允许您定义作为任何可用积分点数量的函数的输出量，这些数量在“使用Abaqus/Standard输出变量标识符”表中列出；允许您将材料方向定义为输出变量；可以用于垫片元素；可以调用实用程序GETVRM以访问材料点数据；除了静态扰动过程外，不能与线性扰动程序一起使用；

原文链接：https://mp.weixin.qq.com/mp/appmsgalbum?__biz=MzkxMzQ0NzY1OQ==&action=getalbum&album_id=2928367776166248448#wechat_redirect

C++语言变量必须"先定义，后使用”，定义变量的一般形式是:变量定义后，变量值是未确定的(除了静态存储情形），即变量值是随机的。直接使用此时的变量参与运算，运算结果也是随机的。

今天简单介绍下有限元中求解有类似于方程约束问题的两类常见算法：拉格朗日乘子法与罚方法分享一个杆单元在约束条件下受力F的求解案例，旨在分享求解算法思想：首先根据几何方程、物理方程等计算系统势能，然后将约束条件条件考虑进来，求势能的最小值（最小势能原理）。这一过程本质上是将有约束的优化问题转化为无约束优化问题。问题简图如下：（具体的两类算法求解程序如后）​1.拉格朗日乘子法拉格朗日乘子法是一种常用的数学工具，用于处理约束条件下的优化问题。

由于每个字符都有一个确定的ASCII码，肯定都可以用\ooo或\xhh表示，因此，\ooo和\xhh称为通用转义字符，其中ooo表示可以用1至3位八进制数作为码值表示一个ASCII字符，hh表示可以用1至2位十六进制数作为码值表示ASCII字符。例如\t为水平制表符，实际上电脑屏幕被分为一些制表区，每8列为一个制表区，每次输入\t，光标就会跳到下一个制表区起始列，要注意的是输入一次\t，是每步跳一个空格，最多跳8步，也就是从当前制表区所在列，一步一步的跳到下一个制表区的起始列。

对于过拟合产生原因的一个解释是离散点数据太少了，距离离散点较远处的信息随着位移插值（试）函数阶次的提高，越来越难被插值点控制。因此只能适当提高位移插值（试）函数阶次。好了现在回归正题。

（这里如果用0.00000678就会被默认解释为双精度浮点数类型，即在二进制转换时用64位存储对应的方法来表示，也就是说十进制的数字要转化为二进制，但是同一个数字转化为单精度和双精度的结果是不同的。C++中这里的0.00000678f叫做常量或者字面常量，别的语言可能就是说是个数字，不强调啥量，这里结尾f是表示这是个float型浮点型数据，因为默认的话，编译器会把他解释为double型，然后用double型对应的二进制转化方式表示。（红色标注的是C++特有的数据类型，其余是C中就有的数据类型）

其中十进制的0~31和127属于控制字符（就是起控制作用的），其余的是图形字符（可以显示到屏幕上的），例如”0”这个字符对应的ASCII码为十进制的48，它在内存中存储的也就是48对应的二进制值，其余的数字依次加1；首先将50.0转化为二进制形式110010.0B（B是来标识数值是二进制，十进制用D，八进制用O，十六进制用H），然后按照规范化数的形式写成(-1)0 * 25 * 1.100100B(二进制的指数记法和十进制一样，小数点左移几位，位权就对应几次方），进而得到。例1中的结果上面已经分析过。

局部坐标系统最初可能是相对于零件坐标系统定义的，但根据给定的零件实例的定位数据进的坐标系统定义是用户子程序中使用的局部方向的定义。可以将实用程序包含在用户子程序中。单元的解依赖变量可以定义为这些子程序中出现的任何其他变量的函数，并且可以相应地更新。如果在分析中定义了由状态变量控制的元素删除，那么控制删除标志的状态变量的值可以由使用状态变量的任何用户子例程进行修改，前提是用户子程序在材料点处被调用。对于大多数子程序，需要在每个材料的材料定义中输入在点或节点处所需的这些变量的数量，以考虑依赖于解的状态变量。

在Abaqus/Standard 中，有时候需要建立运行时间环境并管理与外部数据文件之间的交互作用，或者与用户子程序连接中使用的并行进程之间的交互作用。或者可能为了后处理,需要将用户子程序中的COMMON 块变量所具有的在多个单元上积累得到的输出量，在一个收敛的增量末尾处写入到外部文件中。以便从堆而不是堆栈中分配内存。用户子程序必须执行它们的既有的功能，而不覆盖Abaqus 的其他部分。开发用户子程序时，在尝试将它们用于生产分析工作之前，最好在较小的示例中测试它们，其中用户子程序是模型中唯一复杂的方面。

Abaqus/Standard通过周期性地创建一个新声学网格来解决计算声学域当前构型的问题，所创建的声学域使用与原来网格一样的拓扑结构，但是节点进行了调整，这样结构-声学边界的变形不会导致声学单元的严重扭曲。拉格朗日网格是一种常用的网格类型，也是有限元分析中最常见的网格类型之一。在可预见的大变形的问题中，从自适应网格划分得到的改进的网格质量可以防止分析因严重的网格扭曲而终止。Abaqus的自适应网格技术可以根据模型的特性和求解需求，自动调整网格的细化或粗化，以获得更准确的结果或更高的计算效率。

在 C++ 中，

Abaqus是一种广泛使用的有限元分析软件，它提供了强大的后处理功能，用于对仿真结果进行可视化、分析和解释。Abaqus的后处理模块允许用户对仿真结果进行全面的评估和理解，并从数据中提取关键信息。Abaqus的后处理功能包括以下方面：结果可视化：Abaqus可以将仿真结果以丰富的方式进行可视化展示。用户可以查看模型的应力、应变、位移等物理量的分布情况，并通过不同的图形和动画效果来呈现结果。这有助于直观地理解仿真结果以及模型的行为和性能。

用户可以查看模型的应力、应变、位移等物理量的分布情况，并通过不同的图形和动画效果来呈现结果。Abaqus是一种广泛使用的有限元分析软件，它提供了强大的后处理功能，用于对仿真结果进行可视化、分析和解释。Abaqus的后处理模块允许用户对仿真结果进行全面的评估和理解，并从数据中提取关键信息。用户可以选择感兴趣的位置，比如断面、剪切线或自定义路径，并获得相应的物理量的数值和图表展示。后处理脚本编程：Abaqus提供了Python脚本接口，允许用户编写自定义的后处理脚本来实现更高级的数据处理、结果分析和报告生成。

考虑到数值计算对求解速度的要求较高，以及C++语言的通用性等主要优点，之前入门学习了C++高级编程语言，未来会更新C++基础语法到面向对象部分的内容。内容会讲的比较细和基础，无论是0编程基础还是有其他语言基础的朋友都可以跟着一起学习，学过C++的也可以作为复习和查漏补缺，加深记忆。C++是一种通用的高级编程语言，具有以下优势：高效性：C++是一种编译型语言，可以产生高效的机器代码。它提供了底层的内存管理和直接访问硬件的能力，使得程序的运行速度非常快，适用于对性能要求较高的应用领域。

最近在学习和对比一些界面本构的差异和效果，发现网上一些cohesive单元插入插件使用起来不仅方便，而且可以实现一些特殊的效果，如果大家有这方面的插件希望分享下，我也是刚开始接触，这些对我也很重要，后续我也会继续分享使用所得的技巧和收集的相关资料。3.进入report，点击free body cut,然后进行下图的设置。4.把导出的数据复制到excel中即可，分析和绘图等。1.进入后处理，点击view cut,如下图。2.选择截面位置和截面数量，具体见下图。

（1）用于定义接触对表面之间的摩擦行为；（2）当经典库仑摩擦模型限制性太强且需要定义更复杂的接触面之间剪切传递时，可以使用；（3）必须提供接触面之间剪切相互作用的完整定义；（4）可以使用和更新与求解相关的状态变量；（5）不能与软化的切向表面行为结合使用；（6）不能与通用接触算法一起使用。主要内容1.调用方法2.用户子程序接口3.要定义的变量4.传递信息的变量5.案例讲解该数组必须更新为局部切线方向上所有接触点的摩擦力分量的当前值。重置之前此数组将为零（无摩擦力）。

本次为介绍Vfric摩擦子程序做下铺垫，下周末将会介绍一个稍微不太简单的教程。