【有限元学习笔记】二维四边形单元位移动画云图可视化(有限元后处理)——MATLAB程序

已于 2024-08-16 22:20:41 修改
阅读量768 收藏 18
点赞数 16
分类专栏： 有限元学习笔记 MATLAB学习笔记 文章标签： 学习 笔记 matlab
于 2024-08-15 17:16:40 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_53350487/article/details/141225865
版权

目录

前言

  在有限元求解完成后，需要对结果进行可视化，比如位移云图、应力云图和反力云图等。在前面一篇文章中已经实现了位移云图可视化(静态)，本文章将继续使用MATLAB对二维四边形单元进行位移动画云图可视化。

准备数据

  第一步还是准备数据，节点和单元与前一篇文章相同，不同的是需要准备多个时刻的节点位移数据，下面提供了5个时刻的位移实例数据。

%% 准备数据
% 坐标
XYZ = [ 0 0;
        1 0;
        2 0;
        0 1;
        1 1;
        2 1;
        0 2;
        1 2;
        2 2;];

% 单元
EL = [  1 2	5 4;
        2 3	6 5;
        4 5	8 7;
        5 6	9 8;];

% 水平方向位移时刻t=0.00
DX0 = zeros(9,1);
% 水平方向位移时刻t=0.25
DX1 = [ 0; 0; 0; 0.0174337; 0.0151149; 0.0186553; 0.0318029; 0.0311055; 0.0330584;];
% 水平方向位移时刻t=0.50
DX2 = [ 0; 0; 0; 0.0348675; 0.0302299; 0.0373107; 0.0636058; 0.0622109; 0.0661167;];
% 水平方向位移时刻t=0.75
DX3 = [ 0; 0; 0; 0.0523012; 0.0453448; 0.055966; 0.0954087; 0.0933164; 0.0991751;];
% 水平方向位移时刻t=1.00
DX4 = [ 0; 0; 0; 0.0697349; 0.0604598; 0.0746214; 0.127212; 0.124422; 0.132233;];
% 所有时刻的水平位移
DX = [DX0, DX1, DX2, DX3, DX4];

% 竖直方向位移时刻t=0.00
DY0 = zeros(9,1);
% 竖直方向位移时刻t=0.25
DY1 = [ 0; 0; 0; -0.00104918; -0.0108555; -0.0209396; -0.00926313; -0.0219317; -0.0346498;];
% 竖直方向位移时刻t=0.50
DY2 = [ 0; 0; 0; -0.00209836; -0.021711; -0.0418792; -0.0185263; -0.0438634; -0.0692995;];
% 竖直方向位移时刻t=0.75
DY3 = [ 0; 0; 0; -0.00314755; -0.0325666; -0.0628188; -0.0277894; -0.0657951; -0.103949;];
% 竖直方向位移t=1
DY4 = [ 0; 0; 0; -0.00419673; -0.0434221; -0.0837584; -0.0370525; -0.0877267; -0.138599;];
% 所有时刻的竖直位移
DY = [DY0, DY1, DY2, DY3, DY4];

在这里插入图片描述
计算位移的模

% 位移的模
D = sqrt(DY.^2 + DX.^2);

  为了让整个动画云图的颜色栏采用同一个标准，还需要计算所有时刻位移的最大、最小值。同时，还需要计算坐标的最大、最小值，便于后面设置坐标轴范围。

% 位移模的最大最小值
Dmin = min(min(D));
Dmax = max(max(D));
Xmin = min(min(XYZ(:,1)));
Xmax = max(max(XYZ(:,1)));
Ymin = min(min(XYZ(:,2)));
Ymax = max(max(XYZ(:,2)));

  同样的，还有与Abaqus一致的配色方案。

% 颜色栏
Colormap = [0	0	1;
            0	0.364705882000000	1;
            0	0.725490196000000	1;
            0	1	0.909803922000000;
            0	1	0.545098039000000;
            0	1	0.180392157000000;
            0.180392157000000	1	0;
            0.545098039000000	1	0;
            0.909803922000000	1	0;
            1	0.725490196000000	0;
            1	0.364705882000000	0;
            1	0	0;];

绘图

  绘制动画云图只需要循环绘制每一时刻的位移云图即可，并且在每次绘制下一个云图前清除上一个的云图。下面这是在位移云图的基础上添加了循环的代码。位移云图的绘制可以查看上一篇文章——>点击转跳<——。

%% 绘图
f = figure('Name','位移云图');
hold on
width = 550;      %宽度，像素数
height = 400;     %高度
left = 200;       %距屏幕左下角水平距离
bottem = 200;     %距屏幕左下角垂直距离
f.Position = [left,bottem,width,height];

NT = size(D,2); % 位移云图动画帧数
for k = 1:NT

    XYZ1 = XYZ + [DX(:,k),DY(:,k)];    % 更新坐标

    patch('Faces',EL,'Vertices',XYZ1,'FaceVertexCData',D(:,k),'FaceColor','interp','EdgeColor','k');
    axis tight      % 坐标轴紧凑格式
    axis off        % 关闭坐标轴显示
    
    % 将图向左移动, 为绘制颜色栏做准备
    ax = gca;
    PositionNew = get(ax,'Position');
    PositionNew(1) = 0.08;      % 更改图形距离左边界的距离
    set(gca,'Position',PositionNew);
    
    % 绘制颜色栏
    Co = colorbar('eastoutside');   % 设置颜色栏位移图像外侧
    colormap(Colormap)  % 设置颜色栏
    set(gca, 'FontName','Times New Roman', 'FontSize',14)   % 设置字体字号
    Co.Label.String = 'U, Magnitude';  % 设置颜色栏标题
    
    % 设置颜色栏显示标签个数
    NTicks = 13;    % 标签个数
    [cmin, cmax] = caxis;
    Co.Ticks = linspace(cmin,cmax,NTicks)';
    
    % 设置颜色栏的刻度标签为科学计数法
    yt = Co.Ticks; % 获取刻度位置
    yt_labels = arrayfun(@(x) sprintf('%.3e', x), yt, 'UniformOutput', false); % 将刻度位置转换为科学计数法格式的字符串
    Co.YTickLabel = yt_labels; % 设置刻度标签

end

  下面是效果图，可以看到有几个问题：1.下一时刻的图绘制在了上一时刻的图上面；2.颜色栏的范围不统一；3.颜色栏的位置没有固定；4.原本最底下的边应该是固定的，但是现在产生了移动。
在这里插入图片描述

细节调整

清除上一时刻的绘图

  清除上一时刻的绘图只需要在循环开始时，使用cla清除绘图区即可，下面是部分代码。

NT = size(D,2); % 位移云图动画帧数
for k = 1:NT

    cla                 % 清除上一次的绘图

    XYZ1 = XYZ + [DX(:,k),DY(:,k)];    % 更新坐标

    patch('Faces',EL,'Vertices',XYZ1,'FaceVertexCData',D(:,k),'FaceColor','interp','EdgeColor','k');

在这里插入图片描述

统一颜色栏范围

  设置统一颜色栏范围只需要在绘制颜色栏后，使用caxis对颜色栏进行设置即可，下面是部分代码。代码中的Dmin和Dmax表示所有时刻位移的最大、最小值，在前面准备数据的时候已经提前计算了。

    % 绘制颜色栏
    Co = colorbar('eastoutside');   % 设置颜色栏位移图像外侧
    colormap(Colormap)  % 设置颜色栏
    set(gca, 'FontName','Times New Roman', 'FontSize',14)   % 设置字体字号
    Co.Label.String = 'U, Magnitude';  % 设置颜色栏标题
    
    caxis([Dmin,Dmax]); % 设置颜色栏的最大最小范围
    
    % 设置颜色栏显示标签个数
    NTicks = 13;    % 标签个数
    [cmin, cmax] = caxis;
    Co.Ticks = linspace(cmin,cmax,NTicks)';

在这里插入图片描述

固定颜色栏位置

  在第一个循环时获取颜色栏的位置，然后在后面的每一个循环中将颜色栏设置成第一个颜色栏的位置，部分代码如下。其中，Co表示颜色栏的实例变量，前面绘制颜色栏的时候已经进行了赋值，具体见上一个代码块。

    % 让颜色栏保持同一个位置
    if k == 1
        colorbarPosition = get(Co,'Position');
    else
        Co.Position = colorbarPosition;
    end

在这里插入图片描述
  从上图中可以看到，虽然颜色栏的位置已经固定了，但是可以观察到颜色栏标签的颜色在不断加深。这是因为清除绘图区的时候没有把颜色栏标签清除，导致了颜色栏标签不断绘制叠加，致使颜色加深。

清除颜色栏标签叠加

  与清除绘图区类似，只需要在循环开始时将颜色栏清除即可，下面是部分代码。

    cla                 % 清除上一次的绘图
    colorbar('off')     % 关闭上一次的颜色栏

    XYZ1 = XYZ + [DX(:,k),DY(:,k)];    % 更新坐标

在这里插入图片描述

固定坐标轴范围

  上图中，原本最底下的边应该是固定的，但是现在产生了移动。这是因为每次绘制坐标轴范围都在改变，需要进行固定，在循环中添加下面这段代码即可。代码中的Xmin、Ymin、Xmax和Ymax表示坐标的最大、最小值，在前面准备数据的时候已经提前计算了，后面的1.1或者1.0是为了防止变形后云图超出坐标轴的最大范围而乘的系数。

    % 固定横纵坐标范围
    xlim([Xmin,Xmax*1.1])
    ylim([Ymin,Ymax*1.0])

在这里插入图片描述

保存动画云图

  将下面的代码放入循环的最后，即可保存动画云图，DelayTime表示一帧(一个时刻)的动画播放时间。

    % 保存动态图
    frame=getframe(gcf);
    imind=frame2im(frame);
    [imind,cm] = rgb2ind(imind,256);
    DelayTime = 0.50;   % 单帧播放时间
    if k == 1
         imwrite(imind,cm,'位移动画云图.gif','gif', 'Loopcount',inf,'DelayTime',DelayTime);
    else
         imwrite(imind,cm,'位移动画云图.gif','gif','WriteMode','append','DelayTime',DelayTime);
    end

MATLAB可视化结果
在这里插入图片描述
Abaqus可视化结果
在这里插入图片描述
  可以看到，MATLAB可视化结果与Abaqus可视化结果一致。

总代码

  位移动画云图可视化的总代码如下，后续还会更新MATLAB应力云图(高斯点外插、节点平均)的教程，制作不易，别忘了关注和点赞喔。

clc;clear;close all
set(0,'defaultfigurecolor','w');

%% 准备数据
% 坐标
XYZ = [ 0 0;
        1 0;
        2 0;
        0 1;
        1 1;
        2 1;
        0 2;
        1 2;
        2 2;];

% 单元
EL = [  1 2	5 4;
        2 3	6 5;
        4 5	8 7;
        5 6	9 8;];

% 水平方向位移时刻t=0.00
DX0 = zeros(9,1);
% 水平方向位移时刻t=0.25
DX1 = [ 0; 0; 0; 0.0174337; 0.0151149; 0.0186553; 0.0318029; 0.0311055; 0.0330584;];
% 水平方向位移时刻t=0.50
DX2 = [ 0; 0; 0; 0.0348675; 0.0302299; 0.0373107; 0.0636058; 0.0622109; 0.0661167;];
% 水平方向位移时刻t=0.75
DX3 = [ 0; 0; 0; 0.0523012; 0.0453448; 0.055966; 0.0954087; 0.0933164; 0.0991751;];
% 水平方向位移时刻t=1.00
DX4 = [ 0; 0; 0; 0.0697349; 0.0604598; 0.0746214; 0.127212; 0.124422; 0.132233;];
% 所有时刻的水平位移
DX = [DX0, DX1, DX2, DX3, DX4];

% 竖直方向位移时刻t=0.00
DY0 = zeros(9,1);
% 竖直方向位移时刻t=0.25
DY1 = [ 0; 0; 0; -0.00104918; -0.0108555; -0.0209396; -0.00926313; -0.0219317; -0.0346498;];
% 竖直方向位移时刻t=0.50
DY2 = [ 0; 0; 0; -0.00209836; -0.021711; -0.0418792; -0.0185263; -0.0438634; -0.0692995;];
% 竖直方向位移时刻t=0.75
DY3 = [ 0; 0; 0; -0.00314755; -0.0325666; -0.0628188; -0.0277894; -0.0657951; -0.103949;];
% 竖直方向位移t=1
DY4 = [ 0; 0; 0; -0.00419673; -0.0434221; -0.0837584; -0.0370525; -0.0877267; -0.138599;];
% 所有时刻的竖直位移
DY = [DY0, DY1, DY2, DY3, DY4];

% 位移的模
D = sqrt(DY.^2 + DX.^2);

% 位移模的最大最小值
Dmin = min(min(D));
Dmax = max(max(D));
Xmin = min(min(XYZ(:,1)));
Xmax = max(max(XYZ(:,1)));
Ymin = min(min(XYZ(:,2)));
Ymax = max(max(XYZ(:,2)));

% 颜色栏
Colormap = [0	0	1;
            0	0.364705882000000	1;
            0	0.725490196000000	1;
            0	1	0.909803922000000;
            0	1	0.545098039000000;
            0	1	0.180392157000000;
            0.180392157000000	1	0;
            0.545098039000000	1	0;
            0.909803922000000	1	0;
            1	0.725490196000000	0;
            1	0.364705882000000	0;
            1	0	0;];

%% 绘图
f = figure('Name','位移云图');
hold on
width = 550;      %宽度，像素数
height = 400;     %高度
left = 200;       %距屏幕左下角水平距离
bottem = 200;     %距屏幕左下角垂直距离
f.Position = [left,bottem,width,height];

NT = size(D,2); % 位移云图动画帧数
for k = 1:NT

    cla                 % 清除上一次的绘图
    colorbar('off')     % 关闭上一次的颜色栏

    XYZ1 = XYZ + [DX(:,k),DY(:,k)];    % 更新坐标

    patch('Faces',EL,'Vertices',XYZ1,'FaceVertexCData',D(:,k),'FaceColor','interp','EdgeColor','k');
    axis tight      % 坐标轴紧凑格式
    axis off        % 关闭坐标轴显示

    % 固定横纵坐标范围
    xlim([Xmin,Xmax*1.1])
    ylim([Ymin,Ymax*1.0])
    
    % 将图向左移动, 为绘制颜色栏做准备
    ax = gca;
    PositionNew = get(ax,'Position');
    PositionNew(1) = 0.08;      % 更改图形距离左边界的距离
    set(gca,'Position',PositionNew);
    
    % 绘制颜色栏
    Co = colorbar('eastoutside');   % 设置颜色栏位移图像外侧
    colormap(Colormap)  % 设置颜色栏
    set(gca, 'FontName','Times New Roman', 'FontSize',14)   % 设置字体字号
    Co.Label.String = 'U, Magnitude';  % 设置颜色栏标题

    caxis([Dmin,Dmax]); % 设置颜色栏的最大最小范围
    
    % 设置颜色栏显示标签个数
    NTicks = 13;    % 标签个数
    [cmin, cmax] = caxis;
    Co.Ticks = linspace(cmin,cmax,NTicks)';
    
    % 设置颜色栏的刻度标签为科学计数法
    yt = Co.Ticks; % 获取刻度位置
    yt_labels = arrayfun(@(x) sprintf('%.3e', x), yt, 'UniformOutput', false); % 将刻度位置转换为科学计数法格式的字符串
    Co.YTickLabel = yt_labels; % 设置刻度标签

    % 让颜色栏保持同一个位置
    if k == 1
        colorbarPosition = get(Co,'Position');
    else
        Co.Position = colorbarPosition;
    end

    % 保存动态图
    frame=getframe(gcf);
    imind=frame2im(frame);
    [imind,cm] = rgb2ind(imind,256);
    DelayTime = 0.50;   % 单帧播放时间
    if k == 1
         imwrite(imind,cm,'位移动画云图.gif','gif', 'Loopcount',inf,'DelayTime',DelayTime);
    else
         imwrite(imind,cm,'位移动画云图.gif','gif','WriteMode','append','DelayTime',DelayTime);
    end

end

总结

  这只是一个基础的示例，实际中还会有更具体的、更细致的要求，这就需要再做额外调整；另外本人也仍在学习中，这只是个人的学习笔记，可能还有一些不足之处，欢迎指正。

Infww
关注
  • 16
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
Matlab二维三角形单元有限元程序
qq_42909159的博客
07-21 4580
Matlab二维平面应力三角形单元有限元程序 问题：带孔薄板的平面应力问题 带孔薄板参数：h = 3 in，a = 0.5 in，w = 6 in，σ = 1000 psi，E = 10E6 psi。 求解平板孔的应力情况与弹性力学解对比。 求解代码 求解代码分为输入，单刚计算，总刚计算，求解，输出。 输入代码 %读取input文件按顺序读取节点坐标，单元信息，受力节点信息，固定节点 fid = fopen('input.txt','rt'); %打开文件 NodeCoordinate = fscanf(
matlab流场可视化后处理
热门推荐
hyhhyh21的博客
03-12 3万+
matlab流场可视化后处理1流体中标量的可视化1.1 云图1.2 切片图绘制1.3 三维等值面图绘制2流体中矢量的可视化2.1 箭头图或速度图2.2 流线图2.4 带节点的流线图2.5流管图和流带图2.6 圆锥体图2.7 粒子动画图3 参考链接
【有限元学习笔记】二维四边形单元位移云图可视化(有限元后处理)——MATLAB程序
最新发布
qq_53350487的博客
08-14 721
在有限元求解完成后，需要对结果进行可视化，比如位移云图、应力云图和反力云图等。本文章使用MATLAB对二维四边形单元进行位移云图的可视化。
【VTK+有限元后处理】可视化结果云图
MSTIFIY的博客
11-18 3718
成员函数分别用于读取从ProCAST软件中导出的inp文件（存储有限元的几何拓扑数据）和ntl文件（存储有限元的节点属性数据，如温度、应力等）。类是VTK中的非结构化网格类，可用于有限元分析、计算几何和几何建模这类领域。为了能一次性读取有限元模型整体数据，我们可以将其保存为vtk格式的文件。这个类的作用是，方便我们读取并构建用于在VTK中显示的。为了读取不同格式的有限元计算结果文件，我们先写一个。类来管理有限元的几何拓扑和属性信息。下一次只需读取对应的vtk文件即可。
matlab三维图、二维云图的可视化绘图
fantasty9413的博客
07-16 3万+
关于matlab三维图、二维云图的可视化。 测试和对比了surf、surf+view、surfl、pcolo和contourf等五种绘图可视化的方案，并且还对比了shading interp、shading flat和shading faceted三种类型的可视化效果。
Matlab有限元结果后处理 | 不规则云图绘制(二维/三维)|【Matlab源码+视频教程】| 三维云图 | 二维云图
u010542847的博客
03-14 2175
本程序通过matlab实现有限元计算结果的后处理显示，具体通过两个案例来介绍：案例一、实现一个不规则形状的二维云图的绘制，其中还包含三角形网格剖分方法，然后通过节点解来绘制云图，具体后处理结果如图1所示；案例二、实现三维云图（四面体单元）的绘制，由于这个案例给出单元节点组成，所以无需划分网格，直接通过节点解来绘制云图，如图2所示。
平面四边形单元有限元实现——有限元实践笔记（6）
weixin_43891316的博客
04-10 9790
引言 以上一篇博客平面三角形单元有限元实现中讲过的平面应力问题为例，利用四边形单元及其等参单元实现平面应力问题的有限元求解。 本文由以下几个部分组成： 等单元概念 数值积分方法 平面四边形单元的刚度矩阵和方程的推导 平面结构的有限元分析实例及matlab代码 1 等参元概念 在有限元分析中，实际划分出来的网格很多是不规则形状的单元且形状各不相同，造成位移插值函数也不相同，求解计算效率较低。通过将形状不规整的实际单元和形状规整的标准单元之间可以建立坐标映射关系，映射后的单元称为参数单元，如果坐标变换和试函
【VTK+有限元后处理】实时剖切视图
MSTIFIY的博客
11-19 1572
VTK+有限元后处理：云图实时剖切
二维三角形元有限元的MATLAB实现
weixin_55143442的博客
10-19 1433
三角形元
MATLAB有限元应用-四边形八节点梁受力弯曲
corn1949的博客
01-30 2146
MATLAB有限元应用-四边形八节点梁受力弯曲
static2d.rar_matlab二维云图_二维云图_二维有限元_应力 云图_应力云图
07-14
有限元二维静力计算程序，采用四节点线性单元，可以显示变形和应力云图。
FEM-Matlab.zip_matlab 二维fem_三角形单元_三角形有限元_有限元_有限元三角形
07-14
学习和理解这个MATLAB程序，对于掌握二维有限元方法及其在MATLAB中的实现具有重要意义。你可以通过运行这个程序，并逐步调试和修改，来深入理解有限元方法的工作原理和MATLAB编程技巧。此外，该程序还可以作为进一步...
static3d20.rar_20节点_matlab应力云图_matlab结构云图_云图_有限元 3维
07-15
三维有限元静力问题分析程序，用Matlab编写，采用20节点3维单元，可绘制应力云图。
有限元-六面体单元.zip_FEM_matlab 六面体_六面体单元_有限元_等参元
07-14
有限元六面体单元MATLAB代码，有限元，FEM，六面体单元，等参元
VTK-vtkStructuredGrid学习
weixin_69505365的博客
08-14 346
可以理解为单元格顺序沿着坐标轴排列，但是每个单元格可以不一样。维度理解为X/Y/Z方向的点数，如草图所示分别是4，3，1。此时只提供了12个点，只能构造一层，5.点图、线图、面图如下，对应行的代码改一下出现不同效果图。7.增加Z方向维度，点也要增加，点序保持一致。6.点可以不在同一个平面，随机改变点Z数值。2.获取这些点信息，点顺序要严格按照图示。直接赋值第一层的12个点，改一下Z值。所以Z方向维度是1，不能是其他数据。1.使用CAD随意绘制个网格草图。4.将上面的网格显示出来。
fme处理空间数据入门v0
weixin_43229258的博客
08-13 716
本教程是最最最基本的fme处理空间数据内容，给新同事介绍我以前是怎么干活用的，基本啥有用的东西。本教程主要说如何在 FME Workbench中创建模板文件（一般我们把fme写的工作空间就叫模板），教程里面主要是利用GIS数据。按照分步说明操作，1.了解如何从读取数据；2.添加转换器以在空间上过滤数据；3.将过滤后的数据写入 Esri Shapefile;4.如何运行模板文件。
VBA学习(26)：引用单元格/单元格区域的方法
nianfen的博客
08-13 456
实际上也可以将这些属性组合在一起，例如要选择数据底部的下一个新行，可以使用Range(“A1”).End(xlDown).Offset(1,0)，首先转到A1，向下到数据的底部，然后将结果向下偏移一行，以选择下一个空白行。引用连续单元格区域中最底部的单元格（即该单元格下方的单元格为空）。引用目标单元格的所有相关单元格，包括相关单元格的相关单元格。引用目标单元格的所有从属单元格，包括从属单元格的从属单元格。在VBA中，可以通过多种不同的方式来引用工作表中的单元格/单元格区域。引用目标单元格的从属单元格。
【Python】python中抽象类学习笔记
qq_42761751的博客
08-12 403
抽象类是一个不能被实例化的类，通常作为其他类的基类抽象方法是一个没有具体实现的方法一个抽象类可以有或者没有抽象方法, 但是一般来说定义抽象类的时候总会在类中定义抽象方法Python并没有直接支持抽象类，但是提供了一个模块（abc）来允许定义抽象类通过继承abc.ABC类定义一个抽象类pass上述代码中类虽然继承了ABC但是没有定义抽象方法，仍然可以被实例化。"""使用@abstractmethod来定义抽象方法"""passmain()代码报错：不能实例化一个带有抽象方法的抽象类。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Infww

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值