CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    rhoSimpleFoam 求解可压缩流动时一个奇怪的现象

    OpenFOAM
    8
    53
    18323
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 宝
      宝丁 最后由 编辑

      一直尝试用OF求解叶栅这类的跨音速内流流动,原先用rhoCentralFoam和sonicFoam一直不收敛,然后尝试利用rhoPimpleFoam 计算,本来不抱什么希望,但是计算竟然收敛的挺好,但是计算的结果并不理想,和Fluent 计算得到的结果相比,其叶片前缘的局部高速区域明显很弱,这是不太符合一般的气动规律。
      原先认为是壁面yplus的影响,但是将yplus从(30,60)(标准壁面函数)调到了1(低雷诺数壁面函数),但是情况并没有太大改观。
      ebe960e3-e32d-4582-b4e6-8eeaa972d6b2-图片.png f2b5bb75-180c-4a8a-9d7d-0001a41594d7-图片.png

      感觉每走一步都会遇到很多困难,这应该是与我原先只是将OF 作为一个求解器使用的缘故吧。

      1 条回复 最后回复 回复 引用
      • 宝
        宝丁 最后由 编辑

        来流马赫数都是0.69

        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

          出数据图看看?这种云图看起来有时候并不太准。

          线上CFD课程开始报名:http://www.dyfluid.com/class.html

          CFD高性能服务器 http://dyfluid.com/servers.html

          1 条回复 最后回复 回复 引用
          • 宝
            宝丁 最后由 编辑

            013d07be-3988-41d6-9c99-52909d27e087-图片.png
            这是叶片表面等熵马赫数分布,可以明显看出叶片前缘那一块fluent和rhoPimpleFoam有很明显的差别

            李东岳 五好青年 2 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 @宝丁 最后由 编辑

              Fluent和OpenFOAM相符的时候,可以说明二者都对。
              但是不相符的时候,很难说哪个更好,最好跟实验值比较
              后面的趋势对的不错,可以发SCI了

              线上CFD课程开始报名:http://www.dyfluid.com/class.html

              CFD高性能服务器 http://dyfluid.com/servers.html

              宝 1 条回复 最后回复 回复 引用
              • 宝
                宝丁 @李东岳 最后由 编辑

                @东岳 这个离SCI还很遥远,应该是两个求解器之间的差别吧,实验和计算对起来也很麻烦,我再做两个算例:xiexie:

                1 条回复 最后回复 回复 引用
                • 李东岳
                  李东岳 管理员 最后由 编辑

                  这个后来有什么进展么 :xinxin3:

                  线上CFD课程开始报名:http://www.dyfluid.com/class.html

                  CFD高性能服务器 http://dyfluid.com/servers.html

                  宝 1 条回复 最后回复 回复 引用
                  • 宝
                    宝丁 @李东岳 最后由 编辑

                    @东岳 利用NACA0012,算了一把外流,速度和压力分布都很合理;
                    但是对与我要算的叶栅而言,不论如何调节差分格式,更改网格,计算结果都不让人满意,现在怀疑是算法的差异

                    F 1 条回复 最后回复 回复 引用
                    • F
                      Foamer24 @宝丁 最后由 编辑

                      @宝丁 您好!有改过湍流模型吗?不同模型结果也是有差别的

                      宝 1 条回复 最后回复 回复 引用
                      • 宝
                        宝丁 @Foamer24 最后由 编辑

                        @Foamer24 用过层流、S_A和kOmegeSST,结果并没有本质的区别,前缘加速不明显,低速区域影响范围过大

                        1 条回复 最后回复 回复 引用
                        • 宝
                          宝丁 最后由 编辑

                          我把我的求解文件放到这里,大家有兴趣的可以帮忙瞧一下。
                          链接:https://pan.baidu.com/s/15zFl76UsfD_XyWwUYEdyGQ
                          提取码:o3ae

                          1 条回复 最后回复 回复 引用
                          • 宝
                            宝丁 最后由 编辑

                            最近将离散格式全部改为中心差分格式(原先全部是二阶迎风),计算的结果还比较让人满意,但是这个差分格式还是不太让人满意,并不稳定,对于这个算例最终可以收敛,收敛的结果如下图所示。但是只要换一个叶片马上又发散了,现在真的是很疑惑,前缘的这类现象究竟要怎么调,或者需要动求解器内部的结构了吗?

                            divSchemes
                            {
                                default         none;
                            
                                div(phi,U)      Gauss limitedLinearV 1;//Gauss linearUpwind limited;
                                div(phi,k)      Gauss limitedLinear 1;
                                div(phi,omega)  Gauss limitedLinear 1;
                                div(phi,e)       Gauss limitedLinear 1;
                                div(phi,h)       Gauss limitedLinear 1;
                                div(phi,K)      Gauss limitedLinear 1;
                                div(phi,Ekp)     Gauss limitedLinear 1;
                                div(phiv,p)     Gauss limitedLinear 1;//Gauss linearUpwind limited;
                                div(phid,p)     Gauss limitedLinear 1;//Gauss linearUpwind limited;
                                div((phi|interpolate(rho)),p) Gauss limitedLinear 1;//Gauss linearUpwind limited;
                                div(((rho*nuEff)*dev2(T(grad(U)))))    Gauss linear;
                            }
                            

                            80cfc97e-89dc-45d9-8830-93539f96e8b4-图片.png

                            1 条回复 最后回复 回复 引用
                            • 李东岳
                              李东岳 管理员 最后由 编辑

                              我周末下载看看。

                              目前据我所知,你是对比OpenFOAM和Fluent的数据,同样都是模拟,为什么你更倾向让OpenFOAM的数据接近Fluent而不是反过来呢? :chigua2:

                              我觉得你最好跟实验比一下

                              线上CFD课程开始报名:http://www.dyfluid.com/class.html

                              CFD高性能服务器 http://dyfluid.com/servers.html

                              宝 1 条回复 最后回复 回复 引用
                              • 宝
                                宝丁 @李东岳 最后由 编辑

                                @东岳 首先非常感谢东岳老师
                                就我目前能找到的试验数据,Naca 65-810,fluent 计算和试验值相比较准一些(以下是该叶型计算值与试验值的对比),但是OpenFOAM 的计算结果在趋势上就差很多,所以个人倾向于fluent的计算结果。下面给出NACA65-810的计算网格以及fluent的计算结果与试验值的对比。
                                NACA65-810网格(y+~30-60) 链接:https://pan.baidu.com/s/1tZPl2NfspkOqjmoV4rSpvQ
                                提取码:jhg4

                                fd9d8d85-c8f5-4a4f-ba1b-3be81b40f720-图片.png

                                1 条回复 最后回复 回复 引用
                                • 宝
                                  宝丁 最后由 编辑

                                  最后一次自己顶一下这个帖子

                                  最近尝试了很多,从差分格式到求解和算法控制,一直都没有得到合理的结果,应该是求解器本身的问题了。但是自己时间限制,马上研三了,再做下去心里实在没底,准备去用自己ANSYS 那一系列的软件做事情去了。

                                  但是还是有些不甘心,毕竟了尝试了那么久。我的求解文件和一些网格都在上面,如果有人知道我的错误在哪里,请一定要告诉我,不论时间过去多久,我都愿意拿自己现在一个月的助研费作为酬劳(400,不多哈)。

                                  最后还是感谢东岳老师一直以来的帮助。

                                  1 条回复 最后回复 回复 引用
                                  • 李东岳
                                    李东岳 管理员 最后由 编辑

                                    @宝丁 我想起来了,周末的时候下载了你15楼这个文件,这个求解起来没问题啊?

                                    线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                    CFD高性能服务器 http://dyfluid.com/servers.html

                                    宝 2 条回复 最后回复 回复 引用
                                    • 宝
                                      宝丁 @李东岳 最后由 编辑

                                      @东岳 计我这里计算一开始没什么问题,但是一段时间以后就会发散啊

                                      1 条回复 最后回复 回复 引用
                                      • 宝
                                        宝丁 @李东岳 最后由 编辑

                                        @东岳 是这样,求解格式如果全部都是二阶迎风求解确实很稳定,但是和试验或者Fluent相差很大,所以想着改成中心差分,但是中心差分格式不稳定,第一个叶栅可以收敛,换一下就发散了,异常诡异。

                                        1 条回复 最后回复 回复 引用
                                        • 李东岳
                                          李东岳 管理员 最后由 编辑

                                          替代文字

                                          是说这个图么?实验值在哪里?

                                          线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                          CFD高性能服务器 http://dyfluid.com/servers.html

                                          宝 1 条回复 最后回复 回复 引用
                                          • 宝
                                            宝丁 @李东岳 最后由 编辑

                                            @东岳 是15楼的图,原来的叶栅试验结果和计算结果对不上,后来就放弃了。

                                            从文献找了NACA65-810的试验和计算结果(下图),发现fluent和试验对的非常好,但是rhoPimpleFoam在趋势上就差很多,计算结果就没有放在图上。

                                            f9d07bc1-d947-40ea-9554-cd82673a1258-1559548421462-fd9d8d85-c8f5-4a4f-ba1b-3be81b40f720-图片.png

                                            1 条回复 最后回复 回复 引用
                                            • 李东岳
                                              李东岳 管理员 最后由 编辑

                                              @宝丁

                                              1. 你需要把OpenFOAM的数据、Fluent数据、实验数据放在一起看区别

                                              2. 对不上的时候需要找原因。如何找对不上的原因需要从数据上分析,Cp对不上可能别的能对上。还需要看别的能不能对上,比如壁面剪切力,轴线湍流动能、进口速度等重要的变量。

                                              我近期要备课OKS,只能在7月后期抽出时间给你debug。你这个要写SCI么?还是就放在硕士论文中?我觉得你的研究很系统,如果你做出来,我愿意赞助 http://www.cfd-china.com/topic/1716/cfd验证算例支持计划

                                              线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                              CFD高性能服务器 http://dyfluid.com/servers.html

                                              宝 1 条回复 最后回复 回复 引用
                                              • 宝
                                                宝丁 @李东岳 最后由 编辑

                                                @东岳 嗯嗯 我这边会继续跟进的。

                                                原先想着毕设做这个的,但是没想到卡在这里。现在因为毕业的压力,所以只能用Fluent或者CFX先往前做事情。

                                                如果有新的进展,我也会发到这里跟您讨论的。

                                                1 条回复 最后回复 回复 引用
                                                • 李东岳
                                                  李东岳 管理员 最后由 编辑

                                                  因为你这个算例非常有意思我比较感兴趣,把你电话留下我联系你我们电话讨论

                                                  电话用英文留防止百度 13 six 45 seven 908 eight eight

                                                  线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                  CFD高性能服务器 http://dyfluid.com/servers.html

                                                  宝 1 条回复 最后回复 回复 引用
                                                  • 宝
                                                    宝丁 @李东岳 最后由 编辑

                                                    @东岳 好的 好的:ok:

                                                    13 six 36 seven 12 seven three two

                                                    微信和手机号同号

                                                    李东岳 1 条回复 最后回复 回复 引用
                                                    • 李东岳
                                                      李东岳 管理员 @宝丁 最后由 编辑

                                                      @宝丁 不好意思早上才想起来,你今天晚上5点打给我吧

                                                      线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                      CFD高性能服务器 http://dyfluid.com/servers.html

                                                      宝 1 条回复 最后回复 回复 引用
                                                      • 宝
                                                        宝丁 @李东岳 最后由 编辑

                                                        @东岳 嗯嗯,好的

                                                        1 条回复 最后回复 回复 引用
                                                        • 李东岳
                                                          李东岳 管理员 最后由 李东岳 编辑

                                                          过去半年了... 估计都已经毕业了吧。。。:zoule:

                                                          你觉得我这个结果看起来如何?我不太懂你们这个叶山关注的对象,你感兴趣的话我把算例文件发给你玩玩。用rhoSimpleFoam算的。

                                                          对了,这个结果应该捕获了你所说的

                                                          其叶片前缘的局部高速区域明显很弱,这是不太符合一般的气动规律。

                                                          但是我不太明白这个流体力学本质。能解释一下不?为什么叶片前缘存在局部高速?为什么这个东西前缘加速明显?

                                                          contour.jpeg

                                                          线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                          CFD高性能服务器 http://dyfluid.com/servers.html

                                                          宝 羽之下 2 条回复 最后回复 回复 引用
                                                          • 宝
                                                            宝丁 最后由 编辑

                                                            这个算的感觉很好啊,东岳大神能不能把文件发到我邮箱,baoding@mail.nwpu.edu.cn

                                                            还没有毕业,准备写打毕业论文呢,就是换了一个题目而已:baobao:

                                                            1 条回复 最后回复 回复 引用
                                                            • 宝
                                                              宝丁 @李东岳 最后由 编辑

                                                              @东岳 流体流至叶片前缘,受到阻滞而形成驻点,速度滞止,静压升高,因而会有一个顺压梯度;叶栅吸力面是一个凸曲面,流体流至叶栅前缘,存在一个局部压缩的现象,流体速度增加。在以上两者的共同作用下,叶栅前缘会出现局部高速。

                                                              叶栅流道是一个扩张型通道,起到减速增压的效果,因而整体而言速度会降低,同时静压升高,通道内会存在一个逆压梯度。在吸力面一侧,尤其是叶栅中后部,当流体惯性力不足以克服当地逆压梯度所产生的“阻力”时,会存在流动分离现象,在速度云图上就会出现一个低速区。

                                                              P 1 条回复 最后回复 回复 引用
                                                              • 李东岳
                                                                李东岳 管理员 最后由 编辑

                                                                专业!待我消化消化..

                                                                https://www.jianguoyun.com/p/DUWWPtQQ9s3ZBhiHg7cC

                                                                线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                                CFD高性能服务器 http://dyfluid.com/servers.html

                                                                宝 1 条回复 最后回复 回复 引用
                                                                • 宝
                                                                  宝丁 @李东岳 最后由 编辑

                                                                  @东岳 等我写完大论文再好好整理一下

                                                                  1 条回复 最后回复 回复 引用
                                                                  • P
                                                                    paopaoai11 @宝丁 最后由 编辑

                                                                    @宝丁 叶珊前缘吸力面局部加速跟流动滞止没太大关系,因为前缘滞止点压力最大也不会超过气流总压,所以由顺压梯度造成的加速最大也不会超过主流速度。吸力面近前缘处局部加速主要是由该处流线曲率决定的,曲率越大,即曲率半径越小,则局部流速越大(没有分离的情况下),该处流动法线方向的动能梯度与当地动能成正比,与曲率半径成反比。另外,该过程伴随着流速增大,密度减小,是局部膨胀过程,不是局部压缩过程。

                                                                    宝 我 4 条回复 最后回复 回复 引用
                                                                    • 宝
                                                                      宝丁 @paopaoai11 最后由 编辑

                                                                      @paopaoai11
                                                                      非常感谢您指出我的一些错误,受教了。我再结合您的回答说一下我的理解哈。
                                                                      第一点是前缘滞止所造成的顺压梯度只是一个流动加速的有利因素。
                                                                      第二点是流体在前缘加速,动压升高,静压降低,确实是对应一个膨胀过程,说压缩确实不合适。
                                                                      第三点是曲率对流动的影响,我当时的理解是受到前缘壁面的扰动,局部存在一个收缩的过程,也就对应曲率的影响。

                                                                      1 条回复 最后回复 回复 引用
                                                                      • 宝
                                                                        宝丁 @paopaoai11 最后由 编辑

                                                                        @paopaoai11
                                                                        另外,关于前缘流动的压缩,这个是指的流体虚拟流道的面积被压缩,因而流体流速增加,压强降低,密度增大。

                                                                        1 条回复 最后回复 回复 引用
                                                                        • 宝
                                                                          宝丁 @paopaoai11 最后由 编辑

                                                                          @paopaoai11 另外,关于前缘流动的压缩,这个是指的流体虚拟流道的面积被压缩,因而流体流速增加,压强降低,密度减小,(比体积增大:xinlei: )

                                                                          P 1 条回复 最后回复 回复 引用
                                                                          • 我
                                                                            我是河滩 @paopaoai11 最后由 编辑

                                                                            @paopaoai11 为什么曲率越大,流速越大?或者说为什么曲率越大,压力越小?

                                                                            动边界

                                                                            P 1 条回复 最后回复 回复 引用
                                                                            • P
                                                                              paopaoai11 @宝丁 最后由 编辑

                                                                              @宝丁 我明白你的意思,你是把叶片表面和主流间假想了一个流道,然后来解释这种现象。但这是不正确的,这也是很多同学常犯的理解上的错误。下面我举一个例子来说明一下。
                                                                              在均匀流场中放置一段曲线壁面,如下图所示,一种是ACB形式的上凸壁面,一种是ADB形式的下凹壁面,两者首尾端点位置相同,我们来考察在这两种情况下B点处的流动。试想一下,如果“假想流道”理论成立的话,由于B点位置一样,那么B点处速度应该相同才对,但事实上ACB所对应的B处的速度要远高于ADB所对应的B处的速度。所以,这种“假想流道”或“虚拟流道”理论是不对的。
                                                                              curvature.gif

                                                                              宝 1 条回复 最后回复 回复 引用
                                                                              • P
                                                                                paopaoai11 @我是河滩 最后由 编辑

                                                                                @我是河滩 给你提供一个思路来定量解释这个现象,
                                                                                考察一个二维无粘流动,建立流场中一点流动法线方向的压力梯度表达式(可通过动量方程获得),对于不可压流体,带入静压等于总压减去动压的关系式,即可获得该点法线方向上的速度或者动能梯度,该点的速度即可通过积分获得;对可压流体,则需要带入热焓表达式,同样可以消去压力梯度项,得到速度或者动能梯度,然后积分。

                                                                                1 条回复 最后回复 回复 引用
                                                                                • 宝
                                                                                  宝丁 @paopaoai11 最后由 编辑

                                                                                  @paopaoai11
                                                                                  首先,目前咱们共同的观点在于吸力面近前缘处局部加速主要是由该处流线曲率决定的,而不同的观点在于曲率导致前缘加速现象的机理是什么。
                                                                                  我的观点在于曲率的存在使得来流虚拟流道被压缩,因而流动流速增加。这里对虚拟流道加以解释:壁面的存在使得流场中产生了相应的扰动,不同壁面所对应的扰动传播的距离是不同的,因此将壁面至扰动传播的边界称之为虚拟流道,其是变化的,而非一成不变的。目前这个观点的问题在于我还没有定量的去描述这个虚拟流道,但就流动细节,这又是客观存在的。
                                                                                  但是凹凸曲面的虚拟边界相同吗?用凹凸曲面举例只是对为什么前缘加速只在吸力面一侧显著的复述罢了。注意,这并不是说压力面一侧没有。
                                                                                  其次,我不同意你的这一观点:
                                                                                  叶栅前缘吸力面局部加速跟流动滞止没太大关系,因为前缘滞止点压力最大也不会超过气流总压,所以由顺压梯度造成的加速最大也不会超过主流速度。
                                                                                  这个观点暗含了一个条件:流场中静压处处小于进口静压,但是这并不现实。在本问题中,叶栅前缘的最小静压可是要比进口静压小很多,因而并不一定能保证顺压梯度造成的加速最大也不会超过主流速度。
                                                                                  50585e00-8e2b-46e9-92df-c8406be643f0-图片.png
                                                                                  此外,关于前缘加速的原因我想通过圆柱扰流的算例加以说明:
                                                                                  下图给出了圆柱扰流瞬态静压与速度云图,流体首先在圆柱前缘滞止,静压升高,而在尾迹区,分离涡的存在使得该处为一个低压区,因此在圆柱上下表面存在一个较大的静压梯度,气流加速。
                                                                                  bbe0e601-504a-4eda-893b-da4b31f61cd7-图片.png
                                                                                  综上所述,曲率的存在使得来流虚拟流道被压缩,流体流速增加,其动力来源是静压差。

                                                                                  宝 1 条回复 最后回复 回复 引用
                                                                                  • 宝
                                                                                    宝丁 @宝丁 最后由 编辑

                                                                                    @宝丁 P*=115600:chitang:

                                                                                    1 条回复 最后回复 回复 引用
                                                                                    • 李东岳
                                                                                      李东岳 管理员 最后由 编辑

                                                                                      大佬们讨论的太精彩了,邀请发表在Journal of Fluid Mechanics: CFD-China issue

                                                                                      :chitang:

                                                                                      线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                                                      CFD高性能服务器 http://dyfluid.com/servers.html

                                                                                      宝 1 条回复 最后回复 回复 引用
                                                                                      • 宝
                                                                                        宝丁 @李东岳 最后由 编辑

                                                                                        @东岳 只求毕业啊,不说了,要去写论文了,现在在等一波计算结果出来,只求是最后一波了:xinlei:

                                                                                        1 条回复 最后回复 回复 引用
                                                                                        • 羽之下
                                                                                          羽之下 @李东岳 最后由 编辑

                                                                                          @东岳 查看了东岳老师贴出的case,关于“transonic no;”的设置有一些疑问,transonic的设置是否跟速度场超过340m/s有直接关系?结果显示速度场已经超过了340m/s,为何设置“transonic no;”。还有就是像“transonic”和“consistent”等选项的设置,如何查看其默认设置是否开启。

                                                                                          高速旋转机械/旋风分离器/气固分离

                                                                                          1 条回复 最后回复 回复 引用
                                                                                          • 李东岳
                                                                                            李东岳 管理员 最后由 编辑

                                                                                            transonic主要看你的压力是否具有非常强的传输性,高超那面肯定是了,近音速这面transonic收敛性不是特别好,我就没开

                                                                                            线上CFD课程开始报名:http://www.dyfluid.com/class.html

                                                                                            CFD高性能服务器 http://dyfluid.com/servers.html

                                                                                            1 条回复 最后回复 回复 引用
                                                                                            • 五好青年
                                                                                              五好青年 @宝丁 最后由 编辑

                                                                                              @宝丁 您好,我一直后处理的困惑。请教下这种翼型断面的Ma或者压力系数,类似这种图的后处理方法

                                                                                              请问这样图是怎么后处理的呢?

                                                                                              I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                                                                                              宝 1 条回复 最后回复 回复 引用
                                                                                              • 宝
                                                                                                宝丁 @五好青年 最后由 编辑

                                                                                                @五好青年 我一般是在cfd post输入马赫数求解公式算的

                                                                                                五好青年 1 条回复 最后回复 回复 引用
                                                                                                • 五好青年
                                                                                                  五好青年 @宝丁 最后由 编辑

                                                                                                  @宝丁 123123123.png

                                                                                                  您好,我知道后处理中可以得出压力系数(或者马赫数),但是应该是全流场的压力系数(或者马赫数)吧。

                                                                                                  我怎样才能得到围绕翼型一圈的压力系数(或者马赫数)呢?因为我做的风力机,需要截取50%半径位置处(此处也是二维翼型)的压力系数。

                                                                                                  我没办法选择翼型一圈呀,想问下你这个怎么选上的。

                                                                                                  感谢

                                                                                                  I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                                                                                                  宝 1 条回复 最后回复 回复 引用
                                                                                                  • 宝
                                                                                                    宝丁 @五好青年 最后由 编辑

                                                                                                    @五好青年 你是用什么软件做的计算和后处理

                                                                                                    五好青年 1 条回复 最后回复 回复 引用
                                                                                                    • 五好青年
                                                                                                      五好青年 @宝丁 最后由 编辑

                                                                                                      @宝丁 paraview和cfdpost和tecplot。
                                                                                                      你会哪个,就跟我讲下哪个怎么操作呗。感谢

                                                                                                      I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                                                                                                      宝 1 条回复 最后回复 回复 引用
                                                                                                      • 宝
                                                                                                        宝丁 @五好青年 最后由 编辑

                                                                                                        @五好青年 我用cfd post ,直接提取就好

                                                                                                        1 条回复 最后回复 回复 引用
                                                                                                        • First post
                                                                                                          Last post