- 博客(158)
- 收藏
- 关注
RSS订阅原创 ansys wokrbench热电制冷器分析
从计算结果看，制冷器顶部表面温度约为0.66摄氏度，与case2比较，相差约0.8摄氏度，说明实际工程应用中，考虑环境对流是非常必要的。Case1，用于计算热电制冷器制冷效率；从计算结果看，制冷器顶部表面温度约为零下0.15摄氏度，底部温度为50度，制冷效果较好；热冷气分析时，电场计算与制冷效率分析结果基本相同，这是因为此分析为理想分析，温度分布基本接近，热电效应逆反应产生的电流较小。5热电分析伴随温度梯度变化过程，因此，在材料属性定义时，应考虑随温度变化的各种参数，通过wb可快速与准确地创建。
2025-11-17 12:54:54
828
原创 ANSYS WB-零基础入门案例-方管悬臂梁静强度分析
然后在Y Component（即， Y方向）后面的输入框中输入-10000，即可完成端部1t重力的加载（此处，“-”代表受力方向为Y方向的负向，10000N由重力加速度10m/s^2，根据牛顿第二定律获得，即1000kg（1吨）* 10m/s^2 = 10000N ）Ansys workbench具有功能强大的单位系统（即，量纲），特别是对于复杂问题，单位系统的自由切换，为CAE工程师带来较大的便利，单位系统主要体现在求解与后处理部分，此处仅说明单位系统的重要性。
2025-11-17 12:53:56
1026
原创 ansys wokrbench温差发电器分析
Case4中P1=10.9W，P2= 1.0032W，发电器热效率为beta= P1/ P2=9.2%，约为case3中热效率10.8%的85%，由此可以看出，考虑热变化的材料属性对热效率计算有较大影响。但是随着铜带回路中电阻的增加，会对发热效率产生重要影响，因此，设计发电器时，应降低导线带来的电阻。5热电分析伴随温度梯度变化过程，因此，在材料属性定义时，应考虑随温度变化的各种参数，通过wb可快速与准确地创建。由于材料属性随着温度不同，而电压是因为热电效应产生的，与case3中对比，结果分布都有一定差异。
2025-11-11 11:45:00
887
原创 ANSYS WB-零基础入门案例-方管悬臂梁静强度分析
然后在Y Component（即， Y方向）后面的输入框中输入-10000，即可完成端部1t重力的加载（此处，“-”代表受力方向为Y方向的负向，10000N由重力加速度10m/s^2，根据牛顿第二定律获得，即1000kg（1吨）* 10m/s^2 = 10000N ）Ansys workbench具有功能强大的单位系统（即，量纲），特别是对于复杂问题，单位系统的自由切换，为CAE工程师带来较大的便利，单位系统主要体现在求解与后处理部分，此处仅说明单位系统的重要性。
2025-11-11 09:45:00
1107
原创 模态分析（频率分析）
比如此例中，1阶模态振动方向为Y方向，此时Y方向在前6阶模态的有效质量占比约为71%，是不满足计算要求的，也就是说，此例中模态计算是不准确的，采用当前的网格技术以及模态计算方法不能够进行后期的动力学计算（可通过提高网格精度以及提取更多的模态数提供有效质量占比）。动力学分析可分为时域分析与频域分析，模态分析是动力学频域分析的基础分析类型。其中，[M]=结构质量矩阵，[C]=结构阻尼矩阵，[K]=结构刚度矩阵，{F}=随时间变化的外力载荷函数，{u}=节点位移矢量，{}=节点速度矢量，{ü}=节点加速度矢量。
2025-11-10 10:14:55
627
原创 结构静力学分析
假设筷子是常规木材的，比较脆，基本上都能折断。假设一钢制悬臂梁结构（钢的基本材料属性假设为弹性模量=200GPa，泊松比=0.3，屈服强度235MPa，悬臂梁基本尺寸为300*50*16mm），其一端固定，另一端顶部放置70Kg的物体；如果有人偏胖，坐在椅子上，椅面会有明显下沉（谁坐谁知道...），此时，椅面的下沉量，可用位移矢量来表示。因此，一个结构静力学分析过程，就是在一种假设的情况下（工况），又假设结构在某种受力状态下，不考虑时间效应、不考虑惯性效应以及不考虑阻尼效应的一种理想情况下的结构分析过程。
2025-11-10 10:13:43
945
原创 结构计算的模型处理
通过“Tools” > “Analysis Tools” > “Interference Detection(Beta)”通过“Tools” > “Analysis Tools” > “Small Entiy Search”通过“Tools” > “Repair” > “Repair Hard Edges”- 选择一个元件，选择“Tools”,选择“Named Selections”通过“Tools” > “Repair” > “Repair Sliver”
2025-11-04 09:15:00
593
原创 接触非线性仿真关键参数设置
导致计算误差，尽量避免，工程经验：穿透量<1E-6KeyOpt, cid, 9, 4。- Workbench只支持1~4算法，算法5，须要APDL命令实现。- 球形区域内部，已发生接触，可发生滑动区域，定义为滑动接触。- 球形区域内部，已发生接触，且无滑动区域，定义为粘结接触。- 真实结构不存在此类接触方式，算法稳定性差，很少使用。- 法向接触刚度高，几何穿透小，求解精确，解收敛困难。- 可改变接触刚度(FKN)或穿透容差(FTOLN)- 球形区域内部，且没发生接触区域定义为近场接触。
2025-11-03 17:26:55
928
原创 结构瞬态动力学分析
我们期望在高频模型中使用可控的数值阻尼计算方法，因为使用有限元计算离散空间域的结果，在高频率的模式不太准确。在第一个载荷施加位移0.005mm，设置第一个载荷步的计算结束时间为0.001s，并且关闭时间积分，则初始速度为等于0.005mm/0.001s=5mm/s；ITS小到足够获取间隙“弹簧”频率；由于模态叠加法是基于模态计算得到的，因此在计算中如果存在非线性接触：摩擦，粗糙和无摩擦类型，程序都会转换为线性接触：绑定或不分离。在模态叠加法中，当施加高频的外部激励载荷时，结构的动力学响应是近似计算得到的。
2025-11-03 17:24:54
780
原创 UDF | 创建和使用UDF程序库
如果Interpreted UDF面板中的Display Assembly Listing选项被选中，当编译进行时，控制台窗口中将会显示汇编语言代码。对于每个面，面的质心（the face centroid）的坐标可由F_CENTROID 宏来获得。确认UDF 的case 文件和C语言源代码在工作路径下，如果不一样，则写全路径。UDF 中，begin_f_loop 被用来形成对边界区域上所有单元面的循环。begin_f_loop宏和F_PROFILE宏都是FLUENT提供的宏。
2025-10-28 09:30:00
918
原创 子模型技术与网格无关性
对于子模型切割边界上的所有结点，程序用粗糙模型网格中相应的单元确定自由度数值，然后这些数值用单元形状功能插值到切割边界上。也就是不断的改变网格的疏密，观察计算结果的变化，若其变化幅度在允许的范围之内，我们就可以说计算值已经与网格无关了。在有限元分析中往往出现这种情况，即对于用户关心的区域，如应力集中区域，网格太疏不能得到满意的结果，而对于这些区域之外的部分，网格密度已经足够了。切割边界就是子模型从整个较粗糙的模型分割开的边界。需要注意的是B3（全模型的几何）和C3（子模型的几何）必须基于统一的坐标系。
2025-10-27 10:15:47
684
原创 复合材料结构优化设计
所谓“最优设计”，指的是一种方案可以满足所有的设计要求，而且所需的支出（如重量，面积，体积，应力，费用等）最小。也就是说，最优设计方案就是一个最有效率的方案设计方案的任何方面都是可以优化的，比如说：尺寸（如厚度），形状（如过渡圆角的大小），支撑位置，制造费用，自然频率，材料特性等。设计变量（DVs）为自变量，优化结果的取得就是通过改变设计变量的数值来实现的。b和h都不可能为负值，因此其下限应为b，h>0，而且，h有上限hmax。它必须是设计变量的函数，也就是说，改变设计变量的数值将改变目标函数的数值。
2025-10-27 10:12:47
599
原创 复合材料压力容器热-结构耦合计算
因此，热通量是反映传热强度的指标，又称为热流密度，单位为w/m2，见图。用于3D 分析和2D平面应力及轴对称分析，完全绝热条件施加到表面上，可认为是零热流率加载，在热分析中，当不施加任何载荷时，它实际上是自然产生的边界条件。用户可以指定恒定温度，随时间变化的温度表格，随单空间变量变化的温度表格，也可以直接定义单变量温度函数，该函数可以关于时间或三个方向空间坐标的函数。密度，导热系数和比热,其中稳态传热计算，只需要知道材料的导热系数，并且用户可以定义恒定导热系数，也可以定义随时温度变化的导热系数；
2025-10-21 10:30:00
909
原创 复合材料界面失效和损伤计算方法
界面最大载荷对应的界面变形量与界面失效时对应的变形量比值；复合材料渐进损伤计算，需要在基础复合材料属性基础上，增加刚度衰减材料性能定义，ANSYS渐进损伤模拟支持失效准则包括：最大应力，最大应变，Puck，Hashin，LaRC03和LaRC04.虚拟裂纹闭合技术（VCCT）最初用于计算裂纹体的能量释放率，现在正被广泛应用于层合复合材料的界面裂纹扩展模拟中。该技术假定裂纹总是沿着预先定义的路径扩展，特别是在界面处。为建立裂纹扩展准则，用户必须定义裂纹开始和随后裂纹扩展的断裂准则。用界面单元离散裂纹路径。
2025-10-20 10:16:58
479
原创 复合材料压力容器建模方法
OSS中每一个单元的参考方向通过距离该单元中心最近的边的方向来确定。•可采用不同的rosettes来定义参考方向(0°fiber direction). 所有可用的rosettes在下页展示.径向坐标系的径向定义方向选择集（OSS）的参考方向。OSS中每一个单元的参考方向通过坐标系原点到该单元的中心径向向量来确定。OSS中的每一个单元的参考方向由单元中心点，按照右手定则确定切向向量来确定。采用球坐标定义参考方向，每一个单元的参考方向由该单元中心，绕球坐标Z轴的切向向量确定。
2025-10-20 10:15:36
951
原创 复合材料建模关键技术
Offset-定义CAD/锥面偏移量，即CAD模型/锥面和铺层的交界面可以定义偏移量，偏移的方向通过方向选择集的法向定义。位于体积包络在内的单元将被选择，体积包络根据CAD实体的尺寸或者CAD表面与捕获容差的组合来确定。-当Fabric和Stackups的材料处理设置成Global时，Global Drop-Off Material指定的材料属性将赋予错层单元。选择规则用于根据几何运算选择单元，以及定义结构的局部加强或错层，如图所示，选择规则与方向选择集和铺层组一起使用定义任意形状的铺层。
2025-10-14 09:30:00
1675
原创 复合材料失效准则
cft（纤维拉伸失效fiber tension failure），cfc（纤维压缩失效fiber compression failure），cmA（基体拉伸失效matrix tension failure），cmB（基体压缩失效，matrix compression failure），cmC（基体楔形失效matrix wedge shape failure）hf（纤维失效，fiber failure），hm（基体失效matrix failure），hd（分层失效，delamination failure）
2025-10-13 11:57:20
864
原创 结构计算常用单元介绍
Beam188用于分析细长或中等粗细的梁结构，该单元基于Timoshenko梁理论，包含剪切变形的影响，适用于线性分析及大转动、大应变非线性分析。Ansys Workbench中提供了丰富的单元类型，根据单元位移函数(形函数)的不同，可以分为一次单元和二次单元，这两种单元分别具有不同的特征。当需要修改单元类型时，如果是将二次单元修改为一次单元，比如186单元修改为185单元，那么可以直接在Mesh详细信息面板修改。ansys WB支持的单元类型有：点单元、 线单元 、平面单元、壳体单元、体单元。
2025-10-10 13:15:00
1122
原创 单向/双向流-热耦合
流体会影响热量的传递，同时温度会影响流体的流动，两者会相互影响造成最后才会达到平衡，这种情况我们叫双向流-热耦合，例如电机的自然冷却，电机产生的温度造成空气密度变化及重力作用下，空气会运动起来，运动的空气会带走从电机散出来的热量，可以看成双向流-热耦合。流体只起传递热作用，而热量不影响流体的流动情况，这种情况我们叫单向流-热耦合，例如电机的强制水冷（压力或者水流较高时），可以看成单向流-热耦合。– 用于给定自由流的可压缩流动状态，给定自由流的马赫数和静压，静温等。
2025-10-09 16:50:04
660
原创 ABAQUS复合材料建模及基本分析流程
案例1：创建开孔矩形复合材料常规壳层合板，层合板一端固定，另一端施加拉伸载荷，对模型进行分析，查看每层单方向的应力，对比云图和加载时的铺层额方向，理解铺层方向与lamina材料的概念。Inp同文件夹下创建.bat文件，内容为Call abaqus job=job_convshell cpus=4 int可批量提交计算。利用Abaqus中composite layup快速创建复合材料模型。创建边界条件：固定左边界edge平动自由度。创建载荷件：右边界edge施加拉力。划分网格及选择单元类型。
2025-10-09 16:48:36
663
原创 Fluent流体热分析
产生接触热阻的主要原因是，任何外表上看来接触良好的两物体，直接接触的实际面积只是交界面的一部分（见图）,其余部分都是缝隙。而热边界层的形成与流动边界层的形成类似，只不过形成温度梯度的范围一般比形成速度梯度的范围要小，因而大多数情况下，流动边界层的厚度大于传热边界层的厚度。一般来说，增大流速对传热有利，因为流速越大，对流换热系数增大，传递的热量就越多，反之亦然。主要是实体，其内部会进行传热，两接触的实体的传热类似于一实体内部，能量从高流向低的位置在流体热计算中，仅需要设置实体材料的相关属性就能考虑实体传热。
2025-09-23 15:15:00
1065
原创 Abaqus轮廓积分法与裂纹网格
• 在该区域之外（即，离断裂过程区足够“远”），只要塑性/损伤区域“足够小”，LEFM是可以准确描述断裂状态的，这就是小规模屈服理论（ small scale yielding ，SSY）• 为了能够预测粘结强度。• K用于预测裂纹尖端附近的应力状态（“应力强度”），当这种应力状态（即K）达到临界值时，一个小裂纹就会导致开裂（裂纹扩展），此时材料就会失效。• 但是，如果裂纹尖端附近的非弹性变形区域足够小，以至于在这个区域之外有一个有限区域，那么在这个区域外，LEFM渐近解是精确的。
2025-09-22 10:18:22
920
原创 结构热分析[基本理论/热接触分析/载荷与边界条件/结果与后处理]
用户可以为纯罚函数和增广拉格朗日方程定义一个有限热接触传导（TCC）。– 在在细节窗口 ，为每个接触域指指定定TCC输入值– 如果已知接触热阻，那么它的相反数除以接触面积就可得到TCC值。
2025-09-22 10:16:23
808
原创 APDL在WB中的应用
ANYS提供两种工作方式，GUI图形用户界面（Graphical User Interface又称图形用户接口）操作和命令流。在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成同GUI方式一样甚至GUI不能完成的的操作。
2025-09-17 09:45:00
950
原创 结构网格划分[网格划分原理、整体网格控制、局部网格控制、网格质量评估]
全局网格设置通常用于整体网格划分的部署，包括网格尺寸函数、inflation、平滑度、模型简化、参数输入、激活等。设置合适的全局网格参数可以减小后面具体网格参数的设置工作量，下面以结构分析为例对其展开描述。
2025-09-16 09:39:04
1161
原创 Fluent | 光滑法
平滑方法移动内部节点以吸收运动/变形边 界的运动● 节点数量及其连接不会改变可以用：● tri，tet，quad，hex，prism和poly（单元分割）网格（即对2D和3D）内部区域有两种不同的平滑方法：● 基于弹簧的平滑● 基于扩散的例子：● 变形膜● 摆动结构。
2025-09-10 12:15:00
1040
原创 Fluent | 动网格区域
o 通用程序选择边界/单元区域ü 在区域名称下（Zone names ）-点击绘图可以看到它选择正确的动作类型设置参数ü 与选定的运动类型相关联点击创建点击创建-验证动区域如果在条件上做了任何更改-点击创建要删除一个或多个动区域，-点击删除（全部删除）o 用户定义动网格区域条件必须与UDF结合使用o 允许定义每个节点的位置运动/变形区域中用户定义的运动o 定义用户定义区域选择区域选择用户定义设置运动属性参数-运动属性（UDF / Profile / ...）
2025-09-09 11:39:28
1492
原创 Fluent | 动网格介绍
运动/变形网格（MDM），又叫动网格。一种Fluent模型，允许在模拟过程中改变CFD域的形状，因此在每个时间步中都会添加/移除/拉伸一些单元以适应运动。Ø从初始网格（读取或导入） 从初始网格（读取或导入）Ø求解器处理网在许多工程应用中，零件会以简单或复杂的运动（空间和时间）移动或变形。MDMü 更好地代表和考 第五级 虑系统的真实行为● 更好地代表和考虑系统的真实行为● 更好地预测零件之间的相互作用● 给予的影响力- 流体流动中的零件运动，或- 零件运动中流体流动，
2025-09-09 11:38:00
712
原创 基于pairs公式的裂纹扩展寿命计算
利用尖缺口并带有预制疲劳裂纹的标准试样，如中心裂纹拉伸试样( CCT 试样 -Center Crack Tension)或者紧凑拉伸试样( CT 试样-Compact Tension), 在给定载荷条件下进行恒幅疲劳实验，记录裂纹扩展过程中的裂纹尺寸 a和循环次数 N。（1）由于ΔK=(1-R)Kmax, 当Kmax趋近于Kc，则式中右端分母趋近于零，裂纹扩展速率da/dN趋近于无穷大，裂纹失稳扩展而发生断裂，该公式可以更好的考虑裂纹扩展第III阶段的寿命。a∼N 曲线的斜率，就是裂纹扩展速率da/dN。
2025-09-09 11:36:42
958
原创 常用结构断裂参数的计算
定义沿裂纹面的路径，应以裂纹顶端作为路径的第一点，对于半个裂纹模型而言，沿裂纹面需有两个附加点，对于整个模型，则应包括两个裂纹面共需四个附加点，两个点沿一个裂纹面，其他两个点沿另一个裂纹面，下面图形给出了二维模型的情况。由上面的原理可知，对于非对称模型的平面问题应力强度因子，需要提前5个位置的点的位移，对于对称问题，只需要提取3个数据点，同理在工程中还有一种应力法来求解应力强度因子，但是该方法的精度要比位移法的精度差。VCCT基于裂纹扩展需要的能量等于裂纹闭合需要的能量这样的假设。
2025-09-03 09:00:00
1479
原创 裂纹的建模方法
裂纹尖端第一圈单元的数量（NTHET）和裂纹尖端单元中间节点位置（KCTIP）。单元尺寸的推荐值同2D的情况。对于曲线形式的裂纹前缘，沿裂纹前缘的单元尺寸取决于局部曲率的大小。目前工程中，常见的裂纹类型包括平面裂纹，三维贯穿裂纹，三维表面裂纹和三维埋藏裂纹等形式，其中ANSYS除了三维埋藏裂纹外，其他裂纹模型都可以在软件中直接建立，三维埋藏裂纹模型，需要借助APDL语言编程，才可以建立。裂纹尖端（或裂纹前缘）周围的单元应该是二次单元，在1/4点（at the quarter points）处有中间节点。
2025-09-02 16:29:23
801
原创 转子动力学 | 模态分析
在垂直于自旋轴的平面中，节点的轨道是椭圆形它由三个特征定义：局部坐标系（x，y，z）中的半轴A，B和相位ψ，其中x是旋转轴角度φ是节点相对于主轴A的初始位置。• 目标是估计0到4000 rpm之间的前六个固有频率和模态形状，创建Campbell图，观察由于模态对的分离而看到的陀螺效应（前旋和后旋）–由此产生的陀螺矩阵[g]对垂直于旋转轴Z的平面上的自由度进行了配对。–由此产生的陀螺矩阵[g]将是斜对称的。轴中心椭圆旋转轨迹（本例显示了后旋（BW） ，由于轴承支承是各向异性的，因此会导致椭圆振荡）
2025-08-26 10:00:00
1232
原创 转子动力学 | 后处理
模态的旋转方向计算为属于旋转部件的节点的旋转方向的平均值。 例如，对于在垂直于旋转速度轴的平面上的壳单元中进行网格划分的旋转结构，薄圆盘的旋转效应不会由PRCAMP或PLCAMP命令绘制或打印。例如，如果频率在旋转速度范围内显示出显著的变化，则必须确保已经执行了足够的模态分析，以准确地表示这些变化。发出OUTRES，ALL，NONE，然后发出OUTRES，ITEM，FREQ，CNAME，其中ITEM=NSOL，FREQ=ALL，CNAME是基于节点的组件的名称。在这个例子中，所有的频率都是稳定的。
2025-08-25 17:45:14
1068
原创 转子动力学 | 求解
在多负载步模态分析（例如坎贝尔分析）中使用QRDAMP特征求解器时，可以激活QRDOPT命令上的REUSEKEY，以重用第一步对称特征求解器中的特征量。但是，请注意，如果存在可变轴承（具有表格特性的COMBI214），考虑到每个负载步骤中刚度的变化，此选项可能会导致不正确的结果。虽然QRDAMP特征求解器的计算效率比阻尼特征求解器高，但它仅限于阻尼（粘性、材料等）不重要的情况。当分析中包含旋转阻尼（CORIOLIS命令中的ROTMAST=ON）并且结构的旋转部件使用固体单元时，建议使用阻尼特征求解器。
2025-08-25 17:44:20
776
原创 转子动力学 | 载荷和约束
例如，如果您的模型是线单元模型（每个节点6个自由度），因为所有节点都位于旋转速度轴上，所以使用IC 输入直接强制节点旋转速度（OMGX, OMGY或OMGZ）。如果您的模型是实体单元模型（每个节点有3个自由度），请使用ICROTATE 输入指定旋转速度和旋转速度轴，以生成等效的节点平移速度。在使用基于组件旋转（CMOMEGA）的分析中定义不平衡力时，请确保在连接到属于旋转组件的元素的节点上施加力。在瞬态分析中，使用表阵列参数定义旋转力，以指定每个方向、每个时间步的力的振幅。–在中心节点施加储存反作用力。
2025-08-25 17:43:18
575
原创 转子动力学的建模和分析
在使用过程中，您可以使用DMPRAT，ALPHAD，BETAD，MP，ALPD和MP，BETD指定阻尼。不支持其他类型的阻尼，例如单元，阻尼和瑞利阻尼（DMPRAT，ALPHAD，BETAD，MP，ALPD，MP，BETD，TB，SDAMP ,,, ALPD和TB，SDAMP,,, BETD）。在使用过程中，您可以使用DMPRAT，ALPHAD，BETAD，MP，ALPD和MP，BETD指定阻尼。如果在生成和使用过程中都指定了真实旋转速度，则由于陀螺矩阵的内部缩放，所得到的陀螺效应将基于旋转速度的平方。
2025-08-22 13:15:00
1206
原创 概述 | 转子动力学分析
在较高的旋转速度下，例如在燃气涡轮发动机中，必须始终如一地表示旋转部件的惯性效应，以便准确地预测转子的行为。这种方法的一个主要限制是质量和惯性的位置和分布的不精确近似，以及由此导致的零件本身内力和应力计算的不准确性。陀螺矩阵[G]取决于旋转速度（或者如果结构的零件具有不同的旋转速度），并且是旋转动力学分析的主要贡献者。• 航空发动机，涡轮机，航天器等应用中的机械故障会造成巨大的维修费用，更重要的可能会危及人的生命。[B]是旋转阻尼矩阵，它也取决于旋转速度，改变结构的表观刚度并且可以产生不稳定的运动。
2025-08-21 15:59:22
1056
原创 干货 | 结构热分析
它们之间的净热量传递可以用斯蒂芬—波尔兹曼方程来计算：，式中为热流率，为辐射率（黑度），为斯蒂芬－波尔兹曼常数，约为5.67×10-8W/m2.K4，A1为辐射面1的面积，为由辐射面1到辐射面2的形状系数，为辐射面1的绝对温度，为辐射面2的绝对温度。• 热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。用于 3D 分析和 2D平面应力及轴对称分析，完全绝热条件施加到表面上， 可认为是零热流率加载，在热分析中，当不施加任何载荷时，它实际上是自然产生的边界条件。
2025-08-21 15:57:42
1551
原创 Fluent | 热分析
● 流动传热是指流体在流动过程中与外界发生热量交换的现象。● 流体力学与传热虽然分属于两个不同的领域和研究方向，它们之间有着密切的内在联系。
2025-08-12 11:15:00
1478
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人