- 博客(48)
- 收藏
- 关注
RSS订阅原创 弹性力学的基本定理1：圣维南原理
所谓“静力等效”，是指两种不同的载荷分布如果合力相等、合力矩也相等（也就是满足力的平衡和力矩平衡条件），就可以看作是静力等效的。举个例子，一根钢轴两端传递扭矩时，无论是通过均匀分布的力偶，还是通过键槽传递的集中力，只要总扭矩大小相同，这两种载荷在静力上就是等效的——它们的总作用效果一致。圣维南原理用简洁而深刻的描述，架起了弹性力学理论复杂性和工程实用性之间的桥梁。它不追求绝对精确的理论解，而是在合理误差范围内给工程师提供简化工具，让复杂的应力分析从“不可能”变成“可操作”。
2025-07-17 03:05:31
601
原创 有限元中应力集中与应力奇异辨析
它源于模型中存在的 “数学奇点”—— 即几何或边界条件的极端简化，如绝对的尖角（内角为 0，曲率无穷大）、点约束（将面约束简化为一个点）、集中力（理论上作用面积为 0 的载荷）等。其中，应力集中与应力奇异是两种极易混淆的概念，二者看似都是局部应力升高的表现，但其本质、成因及工程意义却有着天壤之别。应力集中与应力奇异，虽同为有限元分析中的局部高应力现象，却分属物理真实与数学假象的两个范畴。在有限元分析中，网格是连接数学模型与物理现实的桥梁，而网格细化后的应力响应，是区分应力集中与应力奇异的 “试金石”。
2025-07-08 12:41:01
740
原创 金属也会 “累”？聊聊金属疲劳（一）
从工业革命时期的蒸汽机车断裂，到现代航空航天的安全标准，人类对抗金属疲劳的历史，也是材料科学进步的缩影。毕竟，再坚硬的钢铁，也经不起反复的 “折腾”—— 就像生活中的我们，懂得休息，才能走得更远。它藏在高铁的轮轨里、飞机的机翼上、桥梁的钢索中，甚至我们日常用的剪刀、自行车链条里，是工业史上许多重大事故的 “隐形杀手”。这里的 “交变应力” 很好理解：比如火车轮轨每过一次车就受一次力，飞机机翼在起飞降落时反复拉伸压缩，汽车减震弹簧在颠簸中不断伸缩 —— 这些方向或大小周期性变化的力，就是交变应力。
2025-07-08 12:40:03
113
原创 玩低空经济必看！碳纤维单向带和织物的选型避坑
在碳纤维复合材料成型领域，单向带（Unidirectional Tape, UD Tape）与织物（Woven Fabric）的成型特性存在显著差异， “铺层困难、易开裂分层” 是单向带应用中常见的工艺挑战。本文将从材料特性、成型难点、解决方案及延伸问题等方面展开介绍。欢迎读者移步到我的公众号交流学习冬生亦东生。
2025-06-24 00:11:24
607
原创 别让设计变 “事故”！柯西公式——工程师必须掌握的应力分析基础
在固体力学的世界里，柯西公式与主应力堪称描述材料内部力传递的两块重要基石。它们不仅搭建起 “面力 - 应力张量” 之间的桥梁，更通过特征值分析，精准揭示应力状态的本质。接下来，我们将从基础定义出发，结合理论推导与工程实例，带你深入理解它们的内涵与应用。请移步到公众号冬生亦东生阅读，方便交流和获取学习资料。主应力指的是应力张量在主平面上的正应力。主平面具有两个重要特征：主平面上剪应力为零，也就是说在主平面上仅存在正应力；主平面的法向量（即主方向）与面力方向共线，可以表示为，其中就是主应力。概念核心作用。
2025-06-16 23:52:33
854
原创 固体力学核心应力概念全解析：从基础到客观率
在固体力学中，应力是描述材料内部力传递的核心物理量，其衍生的主应力、剪应力、Mises应力等概念，构成了强度分析、塑性屈服、断裂准则的理论基石。本文从基础定义出发，逐层解析关键概念的理论来源、工程意义及相互关联，并深入探讨客观率对动态应力描述的重要性。请移步到公众号冬生亦东生阅读，方便交流。应力是单位面积上的内力，由法国数学家柯西（Cauchy）正应力（法向应力）：位于正面（法向为坐标轴正向）且方向为坐标轴正向为正，也可以简单的理解与截面垂直的分量（如 ，拉伸为正，压缩为负。剪应力（切向应力）
2025-06-16 23:50:05
1001
原创 张量分析2：微积分理论在连续介质力学中的应用简介
张量微积分是微分几何和连续介质力学的核心数学工具，它将向量微积分推广到高阶张量场，为描述复杂物理现象提供了坐标无关的统一语言。本文从张量微分、积分的基本理论出发，深入解析其在连续介质力学中的关键应用。欢迎读者移步到我的公众号交流学习冬生亦东生。一、张量微积分基础理论1. 张量的导数与协变微分（1）张量场的偏导数。
2025-06-13 00:04:30
940
原创 变分法与力学：从泛函极值到物理规律的统一描述
变分法是研究泛函极值问题的数学工具，其核心思想是通过比较函数的微小变化（变分）来寻找最优解。在力学领域，变分法不仅为经典力学提供了简洁的理论框架（如虚功原理、最小势能原理），更是现代计算力学（如有限元法）的理论基石。本文将系统解析变分法与力学的深层关联，并通过具体实例展现其强大应用。欢迎读者移步到我的公众号交流学习冬生亦东生。一、变分法数学基础：从函数到泛函的极值理论1. 泛函与变分的基本概念。
2025-06-13 00:03:26
733
原创 张量分析：从基础到连续介质力学应用
张量分析是连续介质力学、弹性力学、塑性力学等力学分支的核心数学工具。本指南将从张量的基本定义、分类表示、运算规则到坐标变换特性进行系统讲解，帮助学习者建立扎实的张量分析基础。张量的本质是多重线性映射，但通过指标记法和爱因斯坦求和约定，能够将复杂的多维关系简洁地表达出来。掌握张量分析不仅能够理解力学理论的数学表述，还能够为后续学习连续介质力学中的应变张量、应力张量和本构方程等核心概念奠定基础。张量分析通过指标符号和坐标变换理论，将复杂的物理关系转化为简洁的方程形式，成为连续介质力学、流体力学等领域的通用语言。
2025-05-05 00:40:02
878
原创 弹性力学基本方程及有限元编程
在工程领域，无论是设计一座桥梁、制造一台精密机械，还是研发航空航天部件，都需要弄清楚材料在受力时会发生怎样的变化。弹性力学就是专门研究这个问题的学科，它通过微分平衡方程、几何方程和物理方程，为我们揭示材料内部应力、应变和位移的规律。而有限元方法，则是帮助我们把这些理论应用到实际工程计算中的强大工具。弹性力学三大方程与边界条件构成理论基础，有限元方法通过离散化实现工程问题的数值求解。
2025-04-26 10:16:49
1065
原创 弹性力学基本方程及有限元编程
在工程领域，无论是设计一座桥梁、制造一台精密机械，还是研发航空航天部件，都需要弄清楚材料在受力时会发生怎样的变化。弹性力学就是专门研究这个问题的学科，它通过微分平衡方程、几何方程和物理方程，为我们揭示材料内部应力、应变和位移的规律。而有限元方法，则是帮助我们把这些理论应用到实际工程计算中的强大工具。弹性力学三大方程与边界条件构成理论基础，有限元方法通过离散化实现工程问题的数值求解。
2025-04-26 10:13:38
1116
原创 弹性力学中的“单元体“辨析
在弹性力学中，单元体是从物体内某一点处截取的边长为 dx,dy,dz的无限小正六面体，其核心作用是将连续介质的力学行为局部化。通过以下基本假设构建分析框架：连续性假设：单元体内物质无间隙填充，应力、应变、位移等物理量可表示为坐标的连续函数。小变形假设：变形量远小于单元体尺寸，高阶小项（如 (∂u/∂x)2) 可忽略，几何方程线性化。均匀性与各向同性假设（针对特定材料）：单元体材料性质均匀，且各方向力学性能相同（各向同性材料适用）。完全弹性假设：单元体的受力(应力)与变形（应变）成线性关系。
2025-04-24 20:34:38
833
原创 解码刚度矩阵：重新认识有限元静力分析
在有限元静力分析的核心方程 F=KU 中，刚度矩阵 K 不仅是结构刚度的数学表征，其特征值与特征向量更蕴含着对结构行为的深层描述。尽管静力分析通常聚焦于静载荷下的平衡解，但刚度矩阵的特征谱（特征值与特征向量）依然是理解结构刚度特性、变形模式及求解效率的关键。本文从数学本质、物理意义及工程应用三方面，剖析其深层次内涵。刚度矩阵的特征分解，本质上是将复杂结构的整体行为分解为一系列 “基本事件”（模态）的叠加。
2025-04-20 21:50:16
881
原创 Abaqus显示提交，但是不求解怎么办？终于知道原因了
昨天上午老师突然发消息说晚上开个组会，我这周一直在看文献也没啥实质性进展，就想着随便做个模拟应付下老师。卧槽，没想到Abaqus提交作业后，一直不求解。具体问题如下：时间紧急，加上有点慌张，一整就是两个多小时也没解决。后果可想而知。。。。。。因此中秋节我要修电脑了。
2023-09-29 12:11:19
24570
4
原创 Ansys 2023 R2软件安装教程（附安装包）
在结构力学模块方面，ANSYS Products可以模拟各种静态和动态载荷条件下的结构行为，包括线性和非线性分析、疲劳分析、接触分析等，并可以应用于机械、航空航天、汽车、建筑等行业。此外，在Ansys 2023 R2版本中，还优化了块记录功能以提高可读性、简化了光源和sensor等的定义方式、增强了GPU功能以更好地支持仿真模拟运行等。在流体力学模块方面，ANSYS Products可以模拟各种流体流动和传热问题，包括内部和外部流动、多相流动、热交换器等，并可以应用于各种行业。33、设置系统变量。
2023-09-24 10:06:20
62531
130
原创 Abaqus子程序：Uvarm子程序——自定义云图
在包含用户定义输出变量规范的材料定义中，会在所有的元素材料计算点上调用该子程序；因为Abaqus/Standard迭代到收敛解，可能会在每个增量中的每个材料点上多次调用；会在每个步骤的每个增量中调用；允许您定义作为任何可用积分点数量的函数的输出量，这些数量在“使用Abaqus/Standard输出变量标识符”表中列出；允许您将材料方向定义为输出变量；可以用于垫片元素；可以调用实用程序GETVRM以访问材料点数据；除了静态扰动过程外，不能与线性扰动程序一起使用；
2023-09-23 19:40:10
1620
原创 Mathematical13 快速入门上篇（很详细）
原文链接：https://mp.weixin.qq.com/mp/appmsgalbum?__biz=MzkxMzQ0NzY1OQ==&action=getalbum&album_id=2928367776166248448#wechat_redirect
2023-09-14 16:23:09
161
原创 一起学习C++：2.3 计算机中的常量和变量
C++语言变量必须"先定义，后使用”，定义变量的一般形式是:变量定义后，变量值是未确定的(除了静态存储情形），即变量值是随机的。直接使用此时的变量参与运算，运算结果也是随机的。
2023-09-13 16:12:36
146
原创 有限元中的约束算法——拉格朗日乘子法与罚方法
今天简单介绍下有限元中求解有类似于方程约束问题的两类常见算法：拉格朗日乘子法与罚方法分享一个杆单元在约束条件下受力F的求解案例，旨在分享求解算法思想：首先根据几何方程、物理方程等计算系统势能，然后将约束条件条件考虑进来，求势能的最小值（最小势能原理）。这一过程本质上是将有约束的优化问题转化为无约束优化问题。问题简图如下：（具体的两类算法求解程序如后）1.拉格朗日乘子法拉格朗日乘子法是一种常用的数学工具，用于处理约束条件下的优化问题。
2023-09-12 11:53:22
1215
原创 C++学习：2.3计算机中的常量和变量
由于每个字符都有一个确定的ASCII码，肯定都可以用\ooo或\xhh表示，因此，\ooo和\xhh称为通用转义字符，其中ooo表示可以用1至3位八进制数作为码值表示一个ASCII字符，hh表示可以用1至2位十六进制数作为码值表示ASCII字符。例如\t为水平制表符，实际上电脑屏幕被分为一些制表区，每8列为一个制表区，每次输入\t，光标就会跳到下一个制表区起始列，要注意的是输入一次\t，是每步跳一个空格，最多跳8步，也就是从当前制表区所在列，一步一步的跳到下一个制表区的起始列。
2023-09-06 20:19:53
158
原创 有限元分析中，怎么理解节点自由度？
对于过拟合产生原因的一个解释是离散点数据太少了，距离离散点较远处的信息随着位移插值（试）函数阶次的提高，越来越难被插值点控制。因此只能适当提高位移插值（试）函数阶次。好了现在回归正题。
2023-09-04 09:02:36
2646
原创 2.2 计算机中的数据类型
（这里如果用0.00000678就会被默认解释为双精度浮点数类型，即在二进制转换时用64位存储对应的方法来表示，也就是说十进制的数字要转化为二进制，但是同一个数字转化为单精度和双精度的结果是不同的。C++中这里的0.00000678f叫做常量或者字面常量，别的语言可能就是说是个数字，不强调啥量，这里结尾f是表示这是个float型浮点型数据，因为默认的话，编译器会把他解释为double型，然后用double型对应的二进制转化方式表示。（红色标注的是C++特有的数据类型，其余是C中就有的数据类型）
2023-09-02 11:14:14
344
原创 一起学习C++：2.1计算机中的数据存储与表示
其中十进制的0~31和127属于控制字符（就是起控制作用的），其余的是图形字符（可以显示到屏幕上的），例如”0”这个字符对应的ASCII码为十进制的48，它在内存中存储的也就是48对应的二进制值，其余的数字依次加1；首先将50.0转化为二进制形式110010.0B（B是来标识数值是二进制，十进制用D，八进制用O，十六进制用H），然后按照规范化数的形式写成(-1)0 * 25 * 1.100100B(二进制的指数记法和十进制一样，小数点左移几位，位权就对应几次方），进而得到。例1中的结果上面已经分析过。
2023-08-31 23:25:39
164
原创 Abaqus子程序：基础知识概述2
局部坐标系统最初可能是相对于零件坐标系统定义的，但根据给定的零件实例的定位数据进的坐标系统定义是用户子程序中使用的局部方向的定义。可以将实用程序包含在用户子程序中。单元的解依赖变量可以定义为这些子程序中出现的任何其他变量的函数，并且可以相应地更新。如果在分析中定义了由状态变量控制的元素删除，那么控制删除标志的状态变量的值可以由使用状态变量的任何用户子例程进行修改，前提是用户子程序在材料点处被调用。对于大多数子程序，需要在每个材料的材料定义中输入在点或节点处所需的这些变量的数量，以考虑依赖于解的状态变量。
2023-08-28 22:05:36
1903
原创 Abaqus子程序：基础知识概述1
在Abaqus/Standard 中，有时候需要建立运行时间环境并管理与外部数据文件之间的交互作用，或者与用户子程序连接中使用的并行进程之间的交互作用。或者可能为了后处理,需要将用户子程序中的COMMON 块变量所具有的在多个单元上积累得到的输出量，在一个收敛的增量末尾处写入到外部文件中。以便从堆而不是堆栈中分配内存。用户子程序必须执行它们的既有的功能，而不覆盖Abaqus 的其他部分。开发用户子程序时，在尝试将它们用于生产分析工作之前，最好在较小的示例中测试它们，其中用户子程序是模型中唯一复杂的方面。
2023-08-27 21:23:41
3217
原创 Abaqus实用技巧：自适应网格技术（ALE）
Abaqus/Standard通过周期性地创建一个新声学网格来解决计算声学域当前构型的问题，所创建的声学域使用与原来网格一样的拓扑结构，但是节点进行了调整，这样结构-声学边界的变形不会导致声学单元的严重扭曲。拉格朗日网格是一种常用的网格类型，也是有限元分析中最常见的网格类型之一。在可预见的大变形的问题中，从自适应网格划分得到的改进的网格质量可以防止分析因严重的网格扭曲而终止。Abaqus的自适应网格技术可以根据模型的特性和求解需求，自动调整网格的细化或粗化，以获得更准确的结果或更高的计算效率。
2023-08-24 19:58:37
7297
1
原创 Abaqus实用技巧：后处理云图动画多窗口同步显示
Abaqus是一种广泛使用的有限元分析软件，它提供了强大的后处理功能，用于对仿真结果进行可视化、分析和解释。Abaqus的后处理模块允许用户对仿真结果进行全面的评估和理解，并从数据中提取关键信息。Abaqus的后处理功能包括以下方面：结果可视化：Abaqus可以将仿真结果以丰富的方式进行可视化展示。用户可以查看模型的应力、应变、位移等物理量的分布情况，并通过不同的图形和动画效果来呈现结果。这有助于直观地理解仿真结果以及模型的行为和性能。
2023-08-19 12:54:28
2915
原创 Abaqus实用技巧：后处理云图动画多窗口同步显示
用户可以查看模型的应力、应变、位移等物理量的分布情况，并通过不同的图形和动画效果来呈现结果。Abaqus是一种广泛使用的有限元分析软件，它提供了强大的后处理功能，用于对仿真结果进行可视化、分析和解释。Abaqus的后处理模块允许用户对仿真结果进行全面的评估和理解，并从数据中提取关键信息。用户可以选择感兴趣的位置，比如断面、剪切线或自定义路径，并获得相应的物理量的数值和图表展示。后处理脚本编程：Abaqus提供了Python脚本接口，允许用户编写自定义的后处理脚本来实现更高级的数据处理、结果分析和报告生成。
2023-08-18 23:27:18
5187
1
原创 一起学习C++：1.计算机基础概述
考虑到数值计算对求解速度的要求较高，以及C++语言的通用性等主要优点，之前入门学习了C++高级编程语言，未来会更新C++基础语法到面向对象部分的内容。内容会讲的比较细和基础，无论是0编程基础还是有其他语言基础的朋友都可以跟着一起学习，学过C++的也可以作为复习和查漏补缺，加深记忆。C++是一种通用的高级编程语言，具有以下优势：高效性：C++是一种编译型语言，可以产生高效的机器代码。它提供了底层的内存管理和直接访问硬件的能力，使得程序的运行速度非常快，适用于对性能要求较高的应用领域。
2023-08-18 23:11:07
142
1
原创 Abaqus实用技巧：弯矩、轴力、剪力提取方法详解
最近在学习和对比一些界面本构的差异和效果，发现网上一些cohesive单元插入插件使用起来不仅方便，而且可以实现一些特殊的效果，如果大家有这方面的插件希望分享下，我也是刚开始接触，这些对我也很重要，后续我也会继续分享使用所得的技巧和收集的相关资料。3.进入report，点击free body cut,然后进行下图的设置。4.把导出的数据复制到excel中即可，分析和绘图等。1.进入后处理，点击view cut,如下图。2.选择截面位置和截面数量，具体见下图。
2023-08-07 18:20:10
6847
原创 Abaqus二次开发：Vfric摩擦子程序
（1）用于定义接触对表面之间的摩擦行为；（2）当经典库仑摩擦模型限制性太强且需要定义更复杂的接触面之间剪切传递时，可以使用；（3）必须提供接触面之间剪切相互作用的完整定义；（4）可以使用和更新与求解相关的状态变量；（5）不能与软化的切向表面行为结合使用；（6）不能与通用接触算法一起使用。主要内容1.调用方法2.用户子程序接口3.要定义的变量4.传递信息的变量5.案例讲解该数组必须更新为局部切线方向上所有接触点的摩擦力分量的当前值。重置之前此数组将为零（无摩擦力）。
2023-08-06 23:22:33
2564
1
原创 Abaqus子程序单步调试技巧
对于一些复杂的子程序报错问题，尤其是代码量较大时，往往我们需要使用一些程序调试技巧，对于abaqus子程序的单步调试往往有以下优点：1.追踪代码执行流程：通过单步调试，可以逐行查看代码的执行顺序，了解代码是如何运行的。这对于理解代码的逻辑以及发现可能的错误非常有帮助。2.观察变量数值：单步调试允许我们在每个执行步骤中观察变量的数值。这可以帮助检查变量是否按照预期进行计算，是否存在错误或异常情况。3.检查条件和循环：通过单步调试，仔细检查代码中的条件判断和循环语句是否按照预期执行。
2023-07-23 22:11:48
2349
原创 Abaqus子程序之UTEMP
节点号 NODE 处的温度值数组。如果节点未连接到梁或壳单元，则仅需返回一个温度值 (NSECPT=1)。否则，要返回的温度数取决于为梁或壳截面选择的温度模式和场变量输入。以下情况是可能出现的：梁截面的温度和场变量作为梁截面描述中显示的点处的值给出。所需值的数量 NSECPT 由指定的特定截面类型确定，如梁截面库中所述。温度和场变量作为壳截面每一层的 n 个等距点的值给出。所需值的数量 NSECPT 等于 n。
2023-07-09 19:31:11
890
1
原创 有限元非线性求解算法——Newton-Raphson
我一直都在分享笔记，是因为当前我也处于学习理论的状态，不过笔记较多，我分享的主要是之前一段时间的学习笔记，加上我当前的部分理解，重新整理呈现了出来。目前我在寻找损伤力学，断裂力学，张量分析，离散元等方面的优秀教材资料等，希望大家能分享给我，一起学习交流。。
2023-06-10 19:33:30
7458
1
原创 有限元分析—高斯积分与减缩积分
传统的高斯积分在每个单元上使用一组完整的高斯积分点，这可能会导致计算量过大，尤其是对于高阶元素或复杂几何形状的单元。如上图所示，高斯积分的结果与解析解几乎一致，减缩积分的误差略大，可见减缩积分的确会降低精度。具体而言，对于每个单元，将其内部积分区域划分为一系列高斯积分点，然后根据高斯积分点的权重和位置，对积分表达式进行数值近似。两种常见的积分方法是高斯积分和减缩积分。以平面应力问题为例，使用mathematical计算了四边形单元的刚度矩阵，对比了精确解，高斯积分解，减缩积分解的差异性。
2023-05-15 01:36:55
2200
1
原创 Abaqus子程序之UMAT
可用于定义材料的力学本构行为；建议在具有规定牵引（简单拉伸）载荷的单个元素模型上进行初始试验。将在材料定义包括用户定义的材料行为的元素的所有材料计算点调用；可用于包括力学行为的任何程序；可以使用依赖于求解的状态变量；对于力学本构模型，必须提供材料雅可比矩阵，C=(1/J)∂Δ(Jσ)/∂Δε可以与用户子程序USDFLD一起使用，以便于在任何场变量传入（uamt）之前重新定义它们；（1）所有情况都要定义的变量。
2023-05-09 00:13:33
9625
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人