【ANSYS模态分析】Ansys Modal Analysis Tutorial Part3

已于 2024-08-16 18:27:40 修改
阅读量825 收藏 9
点赞数 26
分类专栏： MATLAB ANSYS 文章标签： matlab 矩阵
于 2024-08-16 18:22:25 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/RioiiRichard/article/details/141266769
版权

Ansys Modal Analysis Tutorial Part3

creation time: 2024/8/16

上一节，我们已经将 ANSYS 中获得的数据导出。

这一节，我们将想办法将 ANSYS 导出的数据导入 MATLAB，转化成 .mat 文件。

文章目录

读取 *PRINT 命令导出的数据

  1. 新建 RawTxt 和 MatData 文件夹，新建 cvtAnsysTxt2Mat.m 脚本：

    image-20240816163125767

  2. 将 ANSYS 导出的 .txt 文件复制进 RawTxt 文件夹：

    image-20240816163223190

  3. 修改并运行 cvtAnsysTxt2Mat.m 脚本，结果保存在 MatData 文件夹中：

    image-20240816163337298

    cvtAnsysTxt2Mat.m 脚本内容：

    %% 将 ANSYS 中使用 *PRINT 输出的数据转化为 mat 格式并保存
% author:   Richard
% version:  
% date:     2024-07-04 - 初次建立
%           2024-08-16 - 修改调整
% brief:    
% atention: 
%% cvtAnsysTxt2Mat.m ————————————————————%%
clc;
clear;
close all;


%% READ ————————————————————%%
MapForward = readAnsysTxt('RawTxt\MapForward.txt', 'VEC', 'R');
Phi = readAnsysTxt('RawTxt\Phi.txt', 'MAT', 'R');


%% SAVE ————————————————————%%
save MatData\MapForward.mat MapForward
save MatData\Phi.mat Phi;


%% Function：读取 ANSYS TXT
function mat = readAnsysTxt(file_name, MatorVec, RorC)
    %READANSYSTXT
    % date:     2024-08-16
    % brief:    读取 ANSYS TXT
    % atention: 
    % Input:
    %       file_name - 文件名
    %       MatorVec - 对象是矩阵还是向量，'MAT' \ 'VEC'
    %       RorC - 对象是实数还是复数，'R' \ 'C'
    % Output:
    %       mat - 以 mat 格式返回矩阵（或向量）
    % ---------

    % 读取数据文件
    txt = fileread(file_name);
    % 按行分割并获取数据的行数
    lines = strsplit(txt, '\n');
    line_num = length(lines);

    % 获取矩阵大小
    % 向量模式
    if strcmpi(MatorVec, 'vec')
        % 从后往前匹配 [row, col] 格式
        pattern1 = '\[(\d+)\]';
        for i = line_num:-1:1
            line = lines{i};
            if ~isempty(line)
                % 使用正则表达式匹配
                matches = regexp(line, pattern1, 'tokens');
                if ~isempty(matches)
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 获取矩阵大小
                    Mat_Size = [str2double(matches{end}), 1];
                    % 跳出循环
                    break;
                end
            end
        end
    % 矩阵模式
    elseif strcmpi(MatorVec, 'mat')
        % 从后往前匹配 [row, col] 格式
        pattern1 = '\[(\d+),\s*(\d+)\]';
        for i = line_num:-1:1
            line = lines{i};
            if ~isempty(line)
                % 使用正则表达式匹配
                matches = regexp(line, pattern1, 'tokens');
                if ~isempty(matches)
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 获取矩阵大小
                    Mat_Size = [str2double(matches{end}{1}), str2double(matches{end}{2})];
                    % 跳出循环
                    break;
                end
            end
        end
    else
        error('矩阵 \ 向量 选项设置错误')
    end
    % 提示矩阵大小
    disp(['文件 ', file_name, '，矩阵大小：']);
    disp(Mat_Size);
    
    % 创建空矩阵
    mat = zeros(Mat_Size);

    % 获取矩阵数据
    % 向量模式
    if strcmpi(MatorVec, 'vec')
        % 实数模式
        if strcmpi(RorC, 'r')
            % 匹配 [row]:  value 格式
            pattern1 = '\[(\d+)\]\s*:\s*(\d+)'; % value 为整数
            pattern2 = '\[(\d+)\]\s*:\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)'; % value 为科学计数法
            for i = 1:line_num
                line = lines{i};
                if ~isempty(line)
                    % 使用正则表达式匹配
                    matches = regexp(line, pattern1, 'tokens');
                    if isempty(matches)
                        matches = regexp(line, pattern2, 'tokens');
                        if isempty(matches)
                            continue;
                        end
                    end
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 赋值
                    for j = 1:length(matches)
                        idx = str2double(matches{j}(1));
                        value = str2double(matches{j}(2));
                        mat(idx) = value;
                    end
                end
            end
        % 复数模式
        elseif strcmpi(RorC, 'c')
            % 匹配 [row]:  value1 value2 格式
            pattern1 = '\[(\d+)\]\s*:\s*(\d+)\s*(\d+)'; % value 为整数
            pattern2 = '\[(\d+)\]\s*:\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)'; % value 为科学计数法
            for i = 1:line_num
                line = lines{i};
                if ~isempty(line)
                    % 使用正则表达式匹配
                    matches = regexp(line, pattern1, 'tokens');
                    if isempty(matches)
                        matches = regexp(line, pattern2, 'tokens');
                        if isempty(matches)
                            continue;
                        end
                    end
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 赋值
                    for j = 1:length(matches)
                        idx = str2double(matches{j}(1));
                        value1 = str2double(matches{j}(2));
                        value2 = str2double(matches{j}(3));
                        mat(idx) = value1 + 1i*value2;
                    end
                end
            end
        else
            error('实数 \ 复数 选项设置错误')
        end
    % 矩阵模式
    elseif strcmpi(MatorVec, 'mat')
        % 实数模式
        if strcmpi(RorC, 'r')
            % 匹配 [row,col]:  value 格式
            pattern1 = '\[\s*(\d+),\s*(\d+)\]\s*:\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)'; % value 为科学计数法
            pattern2 = '\[\s*(\d+),\s*(\d+)\]\s*:\s*(\d+)'; % value 为整数
            for i = 1:line_num
                line = lines{i};
                if ~isempty(line)
                    % 使用正则表达式匹配
                    matches = regexp(line, pattern1, 'tokens');
                    if isempty(matches)
                        matches = regexp(line, pattern2, 'tokens');
                        if isempty(matches)
                            continue;
                        end
                    end
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 赋值
                    for j = 1:length(matches)
                        idx1 = str2double(matches{j}(1));
                        idx2 = str2double(matches{j}(2));
                        value = str2double(matches{j}(3));
                        mat(idx1, idx2) = value;
                    end
                end
            end
        % 复数模式
        elseif strcmpi(RorC, 'c')
            % 匹配 [row,col]:  value1 value2 格式
            pattern1 = '\[\s*(\d+),\s*(\d+)\]\s*:\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)\s*([-+]?\d*\.\d+([eE][-+]?\d+)?)'; % value 为科学计数法
            pattern2 = '\[\s*(\d+),\s*(\d+)\]\s*:\s*(\d+)\s*(\d+)'; % value 为整数
            for i = 1:line_num
                line = lines{i};
                if ~isempty(line)
                    % 使用正则表达式匹配
                    matches = regexp(line, pattern1, 'tokens');
                    if isempty(matches)
                        matches = regexp(line, pattern2, 'tokens');
                        if isempty(matches)
                            continue;
                        end
                    end
                    % 提示匹配结果
                    disp(['文件 ', file_name, '第', num2str(i), '行匹配结果：'])
                    for j = 1:length(matches)
                        disp(matches{j})
                    end
                    % 赋值
                    for j = 1:length(matches)
                        idx1 = str2double(matches{j}(1));
                        idx2 = str2double(matches{j}(2));
                        value1 = str2double(matches{j}(3));
                        value2 = str2double(matches{j}(4));
                        mat(idx1, idx2) = value1 + 1i*value2;
                    end
                end
            end
        else
            error('实数 \ 复数 选项设置错误')
        end
    else
        error('矩阵 \ 向量 选项设置错误')
    end
    
end

读取 HBMAT 命令导出的数据

  1. 如果有 hb_to_msm.m 函数文件，可以使用以下函数：

    function matrix = fctn_readHBMAT(file_name)
    % brief: 读取ANSYS数据
    % author: Richard
    % time: 2024/7/4
    % 读取数据文件
    matrix = hb_to_msm(file_name);
    % 将稀疏矩阵转换为double类型
    matrix = full(matrix);
    % 利用下三角阵生成对称矩阵
    matrix = matrix + matrix' - diag(diag(matrix));
end

  2. 或者使用以下函数（可能存在错误）：

    function [matrix, right_vector] = fctn_readHBMAT(file_name)
    % brief: 读取ANSYS数据
    % author: Richard
    % time: 2024/7/4
    % 读取数据文件
    dats = fileread(file_name);
    % 按行分割
    lines = strsplit(dats, '\n');
    % 获取数据的行数
    num_of_rows = length(lines);
    
    % 第二行包含索引信息，使用正则表达式提取
    pattern = '[-+]?\d*\.?\d+([eE][-+]?\d+)?';
    numbers = regexp(lines{2}, pattern, 'match');
    rows_of_data = str2double(numbers{1});              % 总行数，不包括文件头
    rows_of_matrixColPointer = str2double(numbers{2});  % 矩阵列指针的总行数
    rows_of_matrixRowIndex = str2double(numbers{3});    % 矩阵行索引的总行数
    rows_of_matrixElement = str2double(numbers{4});     % 矩阵元素数值总行数
    rows_of_rightVector = str2double(numbers{5});       % 右边项总行数
    
    % 第三行包含矩阵信息，使用正则表达式提取
    pattern = '[-+]?\d*\.?\d+([eE][-+]?\d+)?';
    numbers = regexp(lines{3}, pattern, 'match');
    rows_of_matrix = str2double(numbers{1});            % 矩阵行数
    cols_of_matrix = str2double(numbers{2});            % 矩阵列数
    num_of_matrixRowIndex = str2double(numbers{3});     % 矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）
    
    
    % Step1：重建列指针
    % 从第 idx1_start = num_of_rows - rows_of_data 开始，往后是列指针索引
    idx1_start = num_of_rows - rows_of_data;    % 列指针索引起始行
    % 新建索引向量
    idx1_vector = zeros(rows_of_matrixRowIndex,1);  % 零初始化列指针索引
    col_cnt = 1;    % 列索引计数器
    % 读取索引
    for i = idx1_start+1 : idx1_start + rows_of_matrixColPointer-1
        idx1_vector(str2double(lines{i-1}):str2double(lines{i})-1) = col_cnt;     % 将 x 个元素赋值为当前的列索引
        col_cnt = col_cnt+1;    % 列索引计数器 + 1
    end
    
    % Step2：重建行索引
    % 从第 idx2_start = idx1_start + rows_of_matrixColPointer 开始，往后是行索引
    idx2_start = idx1_start + rows_of_matrixColPointer;
    % 新建索引向量
    idx2_vector = zeros(rows_of_matrixRowIndex,1);
    % 读取索引
    for i = idx2_start : idx2_start + rows_of_matrixRowIndex-1
        idx2_vector(i - idx2_start + 1) = str2double(lines{i});
    end

    % Step3：重建数据
    % 从第 idx3_start = idx2_start + rows_of_matrixRowIndex 开始，往后是数据
    idx3_start = idx2_start + rows_of_matrixRowIndex;
    % 新建索引向量
    idx3_vector = zeros(rows_of_matrixElement,1);
    % 读取索引
    for i = idx3_start : idx3_start + rows_of_matrixElement-1
        % 使用正则表达式提取数值
        pattern = '[-+]?\d*\.?\d+([eE][-+]?\d+)?';
        numbers = regexp(lines{i}, pattern, 'match');
        % 数据采用了科学计数法表示
        if length(numbers) == 1
            idx3_vector(i - idx3_start + 1) = str2double(numbers{1});
        elseif length(numbers) == 2
            idx3_vector(i - idx3_start + 1) = str2double(numbers{1}) * 10^(str2double(numbers{2}));
        else
            disp('Wrong.')
        end
    end
    
    % Step4：重建右端向量
    if rows_of_rightVector ~= 0
        % 从第 idx4_start = idx3_start + rows_of_matrixElement 开始，往后是右端项
        idx4_start = idx3_start + rows_of_matrixElement;
        % 新建右端项向量
        idx4_vector = zeros(rows_of_rightVector,1);
        % 读取索引
        for i = idx4_start : idx4_start + rows_of_rightVector-1
            % 使用正则表达式提取数值
            pattern = '[-+]?\d*\.?\d+([eE][-+]?\d+)?';
            numbers = regexp(lines{i}, pattern, 'match');
            % 数据采用了科学计数法表示
            if length(numbers) == 1
                idx4_vector(i - idx4_start + 1) = str2double(numbers{1});
            elseif length(numbers) == 2
                idx4_vector(i - idx4_start + 1) = str2double(numbers{1}) * 10^(str2double(numbers{2}));
            else
                disp('Wrong.')
            end
        end
    end
    % 复制右端项
    right_vector = idx4_vector;
    % 新建稀疏矩阵
    matrix = sparse ( idx2_vector, idx1_vector, idx3_vector, rows_of_matrix, cols_of_matrix);
    % 将稀疏矩阵转换为 Double 类型
    matrix = full(matrix);
    % 利用下三角阵生成对称矩阵
    matrix = matrix + matrix' - diag(diag(matrix));
end

  3. HBMAT 导出的质量、刚度、阻尼矩阵应与 *PRINT 导出的矩阵一致。

YunchengLi
关注
专栏目录
ansys 命令流学习
weixin_50061663的博客
06-10 5242
ansys 基础命令流
【ANSYS模态分析】Ansys Modal Analysis Tutorial Part1
RioiiRichard的博客
08-16 1426
柔性航天器控制涉及到对卫星模型质量矩阵、刚度矩阵、模态振型矩阵的有限元分析，整个ANSYS的学习过程十分痛苦，阅读了书籍和在线资料，花了很多时间但进展缓慢。本专栏文章将记录作者利用 ANSYS Workbench 进行模态分析的过程。在这一节中，我们将利用 ANSYS Workbench 新建一片简单的太阳能帆板，并对其做模态分析。
【ANSYS模态分析】Ansys Modal Analysis Tutorial Part2
RioiiRichard的博客
08-16 1050
上这一节中，我们利用 ANSYS Workbench 新建一片简单的太阳能帆板，并对其做模态分析。 这一节，我们将想办法将 ANSYS 中获得的数据进行导出。由于是边学 ANSYS 边记录，这一节可能存在错误或者问题，请大家批评指正。
基于ANSYS的无约束梁的模态分析与实验结果比较
ikhui7的博客
03-04 2759
通过借助ANSYS建立无约束梁有限元模型进行模态分析，计算前3阶固有频率，并与实验结果进行比较，结果表明，ANSYS能够在较少的划分单元数的情况下，较精确地与实验结果吻合。
ansys中模态扩展是什么意思_ansys模态分析步骤
weixin_39661353的博客
12-19 1255
模态分析步骤第1步：载入模型Plot>Volumes第2步：指定分析标题并设置分析范畴1设置标题等UtilityMenu>File>ChangeTitleUtilityMenu>File>ChangeJobnameUtilityMenu>File>ChangeDirectory2选取菜单途径MainMenu>Preference,单击...
ANSYS workbench的模态分析基本原理和步骤
qq_42384612的博客
04-21 2万+
模态分析是用来确定结构的振动特性（固有频率和振型）的技术，是其他动力学分析的基础。 作用：使结构避免共振，或以特定的频率进行振动，通过模态分析，可以认识到结构对不同动力载荷的响应特性。 模态分析重点关注：固有频率和振型。
ANSYS_APDL——实例001-模态分析
alsj123456的博客
09-16 4615
学习APDL无从下手，怎么办，就从实例出发，一点一点了解 /clear /prep7 et,1,solid186 mp,ex,1,2e11 !mp后面都为材料参数 mp,prxy,1,0.3 mp,dens,1,7800 block,0,0.01, 0,0.01, 0,0.1 !创建体，基于工作平面，x1，x2，y1，y2，z1，z2 lesize,1,,,3 !将编号为1的线划分为3份， lesize,2,,,3 !将编号为2的线划分3份 lesize,9,,,15 !将编号为9的线划分15份 m
ANSYS模态分析
热门推荐
qq_47470595的博客
02-12 3万+
目录 一、ANSYS求解模态分析步骤  建模  施加载荷和求解 扩展模态（这是区别于其他分析的一个关键） 查看结果 二、模态分析实例（均匀直杆的子空间法模态分析） 1.问题描述 2.有限元分析 1）进入ansys经典界面 2）设置单元类型：Solid Brick 20 node 186 3）设置材料属性 4）建模 5）划分网格单元 6）加载：模态分析...
ANSYS_APDL——实例002-模态分析
alsj123456的博客
09-16 2718
/clear /prep7 et,1,link1 mp,ex,1,2e11 !材料属性 mp,prxy,1,0.3 mp,dens,1,7800 R,1,1e-6 !定义元素实常数 k,1,0,0,0 !定义1的坐标为000 k,2,1,0,0 !定义2的坐标为100 lstr,1,2 !定义点1到点2的直线 lesize,1,,,50 !线1，中间50结 lmesh,1 !沿着线生成节点 finish !============== /solu dk,1,ux !定义约束 dk,1,uy.
ANSYS 有限元分析 轴承座的模态分析 分析过程及常见问题解决
m0_68737530的博客
02-04 740
1.确定分析类型为 static structural。8.固定轴承座进行分析，得到约束状态下的结果如图7-12所示。3.建立几何模型，模型为一个part多个body的形式。7.得到自由状态下轴承座的态分析结果如图1-6所示。2.工程材料采用Q235。5.设置边界条件和载荷。
ANSYS简单框架结构地震分析算例.rar_ansys_ansys 地震_ansys框架_ansys模态分析_地震时程ansys
07-15
在ANSYS中，可以使用`/MODAL`命令启动模态分析，通过计算出的固有频率，工程师可以判断结构是否容易在地震波的某个特定频率下共振，从而采取必要的设计措施。 **时程分析**是地震工程中更为实际的分析方法，它考虑...
ansys支架静力、屈曲及模态分析
06-18
《ANSYS支架静力、屈曲及模态分析详解》 在工程领域，ANSYS是一款广泛应用的有限元分析软件，能够对结构、流体、热传递等问题进行精确模拟。本篇将深入探讨如何利用ANSYS进行支架的静力分析、屈曲分析以及模态分析...
Ansys模态分析专章讲解.rar
04-23
4. 模态求解器设置：在Ansys Mechanical中选择“Modal”分析类型，设置求解参数，如最大模态数量、频率范围等。 5. 运行分析：启动求解器，Ansys将计算出结构的所有固有频率和对应的振型。 6. 后处理：通过图形...
Desktop.rar_ANSYS导出模态_APDL_APDL模态_NGY_ansys 模态振型 命令流
07-15
3. **运行求解器**：使用`/SOLU`命令启动求解器，然后执行`/MODAL,DYNAMIC`命令开始模态分析。 4. **结果提取**：使用APDL命令如`/POST1`进入后处理模式，通过`MODE,SPEC,1,NGY`命令获取前NGY个模态的频率和振型。 ...
IWO-Kmeans聚类 | MATLAB实现IWO-Kmeans侵入性杂草优化K均值聚类算法
关注并私信文章链接，获取对应文章源码和数据，机器学习之心的博客。
10-14 97
IWO-Kmeans聚类 | MATLAB实现IWO-Kmeans侵入性杂草优化K均值聚类算法
记录搜罗到的Matlab 对散点进行椭圆拟合
learner_jj的博客
10-11 257
散点的椭圆拟合
MATLAB(Octave)混电动力能耗评估
最新发布
跨学科知识视角展现
10-15 1089
1. 处理电动和混动汽车能耗的后向和前向算法模型(simulink)，以及图形函数、后处理函数等实现。 2. 构建储能元数据信息：电池标称特性、电池标识符等以及静止、恒定电流和恒定电压等特征阶段。 3. 使用电流脉冲或要识别的等效电路模型类型配置阻抗分析。 4. 依据恒定电压阶段相关或不相关的恒定电流阶段等，分析电池的特征量。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值