背景说明
最近在通过仿真秀焦晓娟老师的课堂学习 True-load载荷识别软件,焦老师在讲解完软件的使用说明之后,以一个悬臂梁的案例帮助学生入门。作为初学者,很想将这个过程中的思考和理解记录一下以帮助自己捋清楚思路。读者想要深入学习,推荐报名焦老师的课程。
案例是一个悬壁梁,一端固定,需要识别另一端的载荷(分解为x和y方向),也就是说需要识别两个载荷。
已有的前提条件:
1.有限元单元载荷的计算结果文件:.odb或.rst,前者为Abaqus的结果文件格式;后者为Ansys的结果文件格式。这个文件是在有限元模型上在两个方向上施加单位载荷所得出来的应力应变的仿真结果。在本例中由焦老师案例直接给出。
2.应变试验设备:因为载荷识别的输入是应变片的测量数值,所以需要实测的应变片数值,从而对实际载荷进行识别。
载荷识别的基本流程
提出问题
Q1:对图中的缩写词进行解释:FEA, TFU,QSE?
A1: FEA: Finite element analysis; 后面两个暂时未知,后面补上;
Q2:请简单描述一下上图工作流程?
A2:描述如下:
- 首先,在有限元仿真软件环境中,通过在载荷未知量的位置施加单位载荷得到单位载荷激励下的FEA文件;
- 然后,在软件的Pre-Test模块载入FEA文件,创建Load Table, 选择候选区域并且确定应变片的贴方位;在Pre-Test模块最后会生成一个 .tld文件,该文件包含了应变片的位置信息。
- 其次,经过实际测试之后得到应变片csv格式的应变数据, 与Pre-Test模块创建的.tld文件一同载入到Post-Test模块进行载荷的识别。Post-Test模块会生成一个 .tfu文件,该文件包含了载荷识别之后的信息,与最开始的.odb文件一同导入QSE模块可以映射全场应变。
- 最后,在TFU Manager 中可以实现曲线绘制,并且比较测试应变和识别到的载荷所计算的应变。
PS: 单位载荷激励下的FEA文件需要实操获取一下;csv格式应变数据的具体保存样式需要更明确一下;什么是Load Table?为什么要创建Load Table?
Pre-Test 模块
Pre-Test 界面介绍
每个按钮的使用焦老师在课堂上有非常详细地讲述,建议看完这个案例之后再返回去看焦老师对Pre-Test UI界面的介绍。
现在想要记录一下自己在观察True-Load Pre-test模块的一些浅显的猜测和思考。
Q: Pre-Test 模块主要为了什么目的而开发的?
为了确定应变片的位置。
Q: 这其中的关键点有哪些?Pre-Test 用了怎样的思路?
Pre-Test模块首先需要选中应变片的贴片区域(一般在1000个位置以下),进而通过稀疏得到一些离散位置(De-populate selection按钮), 然后载入这些位置在EFA文件中的单位载荷激励下应变张量,基于要贴片的数量,在条件数约束下(P-Scale针对相关性C矩阵,E-Scale针对Strain matrix),对离散的位置进行排序,越靠前的位置表示位置计算出来的条件数越低,测量效果将会越好。然后再Test similation 模块,通过单位载荷激励的应变数值5%的噪声来计算识别出来的载荷数值的敏感度。当条件数比较小的时候,输入的应变数值在其数值处加入5%的噪声时,识别出来的载荷Load也应当是在0~1之间的方波。如下图所示,说明使用Pre-Test模块找出来的贴片位置的条件数是让人满意的,输入的应变测量值有一定变化,输出的识别载荷也会有相应的变化,敏感度好且稳定。如果条件数不好,那么很容易就会失稳,不是规律的方波了。
由于从离散位置空间中进行寻优排序是最消耗时间的,所以后续如果想要增加和删减应变片的位置,如果不是点击全部重新寻找按钮,那么直接找到根据条件数约束已经排好序的位置直接替补和删除。
由于应变片的位置、角度还需要具备人为进行修改的功能,所以在找到的位置进行对应变片位置、旋转角度的显示、输入等是重要的,猜测这个角度数据将会用于将应变片的应变分解到两个正交的方向。
Q:什么是条件数?
并非是True-load所独有的,是线性代数中的一个术语,基本上是矩阵中最大值与最小值之比,用于度量方程线性稳定性的程度,最理想的情况是1,越大表示越不稳定,小于10代表非常稳定的系统,小于50也是可以接受的,大于等于100在工程上是不可接受的。对于一个模型,改善条件数有多种方法,贴加更多的应变片可以改善条件数。
条件数高,不稳定大致有以下几种原因:
- 应变片的数量不够,没有办法来捕捉传递路径从而生成唯一的应变特征
- Load case之间彼此本质上是相关的,一种情况彼此之间是相应的倍数关系,某些载荷会导致相同结构的变形,另一种情况是Load case之间相对的应变值差别非常大
载荷识别的原理
如下图所示,应变片的应变矩阵 ϵ × \epsilon \times ϵ× 相关矩阵 C C C等于彼此独立正交的Load矩阵。
相关矩阵 C C C怎么计算呢?根据单位载荷的应变值来进行计算。
在这个例子当中,最后Load assurance criteria (LAC 矩阵)应该和上图最后求出来的Load 矩阵有关系,越是正交,就说明计算出来的结果越好。