如何解多个方程?(管道设计)

2次查看(最近30天)

显示旧的评论

艾哈迈德穆罕默德Mansoor 2021年6月12日

编辑: dpb 2021年6月15日

答:接受 dpb

我需要从这些方程中计算出t。

其余的值都是给定的。但是这些方程是相互联系的。我只需要知道t的值是多少。

电脑= 1.324

规定的最小屈服强度，

= 448

v = 0.3

E = 2.1 * 10 ^ 5

D = 304.8 mm(12英寸)

管椭圆度, δ =我认为if条件足以根据D的大小提取出0.75或1.0的值

注意:忽略Dmin和Dmax

有人能帮我一下吗?谢谢你！

1评论
显示隐藏 None

dpb 2021年6月13日

写一个函数计算Pc机，然后使用 fsolve 求Pc - Pc_target =0时的临界t

登录评论。

登录回答这个问题。

接受的答案

dpb 2021年6月13日

1
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/854690-how-can-i-solve-multiple-equations-pipeline-design#answer_723940

编辑:dpb 2021年6月15日

                         函数电脑= pipelineExternalPC (d, t, sy)
                        
                         %用法:
                        
                         % Pc = pipelineExternalPC (d, t, sy)
                        
                         ％
                        
                         %返回Pc，管道外部坍塌压力(MPa)为管道
                        
                         直径和壁厚%;最小单位为英寸
                        
                         屈服强度sy (MPa)
                        
                         E = 207年e3;钢杨氏模量(MPa)
                        
                         ν= 0.30;%泊松比
                        
                         椭圆度= (0.75 + 0.25 * (d > = 20)) / 100;椭圆度修正系数
                        
                         D = 25.4 * D;t = 25.4 *;%转换尺寸为毫米
                        
                         托德= t / D;%厚度/直径比
                        
                         d =椭圆度/托德;直径与厚度之比
                        
                         g = @ (r d) (sqrt(一维* d) - d) * 12 + 1 + r * r) /√(1 + r * r * sqrt (1 + d * d) - d * d);
                        
                         Pe = 2 * E * tOD ^ 3 / (1.4 * (1-nu *ν)* (1-tOD) ^ 2);
                        
                         Py = 2 * sy *托托;
                        
                         r = Py /聚乙烯;
                        
                         电脑= g (r d) * Py * Pe /函数(Py、Pe);
                        
                         结束

现在写一个匿名函数来求解给定的约束条件…

                         > > fn = @ (t) pipelineExternalPC -1.324 (12, t, 448)
                        
                         fn =
                        
                         function_handle和值:
                        
                         @ (t) pipelineExternalPC -1.324 (12, t, 448)
                        
                         > > fsolve (fn, 0.25)
                        
                         方程解决了。
                        
                         fsolve完成是因为函数值的向量接近于零
                        
                         作为由功能公差的值来衡量,
                        
                         的通过梯度测量，问题似乎是规则的。
                        
                         <停止标准详细信息>
                        
                         ans =
                        
                         0.2576
                        
                         > > pipelineExternalPC (12, ans, 448)
                        
                         ans =
                        
                         1.3240
                        
                         >>

注:结果以英寸为单位，这是函数的输入尺寸单位。

2的评论
显示隐藏 1年长的评论

艾哈迈德穆罕默德Mansoor 2021年6月13日

非常感谢!!!!这工作。我对代码做了一些修改，出现了一些小错误。

这就是它现在的样子。

                             函数电脑= pipelineExternalPC (Di, E,ν,,t, y)
                            
                             %用法:
                            
                             % Pc = pipelineExternalPC (d, t, sy)
                            
                             ％
                            
                             %返回Pc，管道外部坍塌压力(MPa)为管道
                            
                             直径和壁厚%;最小单位为英寸
                            
                             屈服强度sy (MPa)
                            
                             椭圆度= (1. - -0.25 * (Di < 20));椭圆度修正系数
                            
                             托德= t /做;%厚度/直径比
                            
                             d =椭圆度/托德;直径与厚度之比
                            
                             g = @ (r d) ((sqrt (1 + d * d) - d) * sqrt (1 + r * r)) / (sqrt (1 + r * r * (sqrt (1 + d * d) - d) ^ 2));
                            
                             Pe = (2 * E * tOD ^ 3) / (1.4 * (1-nu *ν)* (1-tOD) ^ 2);
                            
                             Py = 2 * y *托托;
                            
                             r = Py /聚乙烯;
                            
                             电脑= g (r d) * Py * Pe /函数(Py、Pe);
                            
                             结束

现在我的调用函数是这样的

                             fn = @ (t) pipelineExternalPC (Di, E,ν,,t, y) -P_c;
                            
                             t_pc = fsolve (fn, 0.1);

谢谢你的帮助和支持。金宝app

dpb 2021年6月15日

椭圆度= (1. - -0.25 * (Di < 20));椭圆度修正系数

您确实捕获了尺寸的反向逻辑测试，对更大直径的管道使用1%;抢手货。

我觉得更好的写法是

椭圆度= (0.75 + 0.25 * (Di > = 20));椭圆度修正系数

增加1/4%的直径，而不是从较小的直径。似乎读得更好;对人类更有意义;电脑才不在乎呢!:)

但是，上面没有把百分比转换成绝对数字，这就产生了一个错误。

它被明确定义为比值

[Dmax-Dmin] / Dnominal——> [Dmax-Dmin] / {[Dmin +距离]/ 2}= = > 2 * [Dmax-Dmin] / [Dmin +距离)

这是绝对比率;如果用%表示，那就是那个数字的100倍。

注意:如果从实际管道测量中计算，参考文献中的公式缺少公称除数的二因子。

但是，通过删除函数中作为百分比输入的数的除数，引入了一个大了100的因数。

https://pipesak.com/2018/05/18/a-crash-course-part-3-pipeline-ovality/

登录评论。

如何解多个方程?(管道设计)

1评论
显示隐藏 None

接受的答案

2的评论
显示隐藏 1年长的评论

更多的答案(0)

另请参阅

类别

标签

下载188bet金宝搏

社区寻宝

如何解多个方程?(管道设计)

1评论 显示隐藏 None

接受的答案

2的评论 显示隐藏 1年长的评论

更多的答案(0)

另请参阅

类别

标签

下载188bet金宝搏

社区寻宝

1评论
显示隐藏 None

2的评论
显示隐藏 1年长的评论