版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 142.4
帖子: 39
在线: 6.3小时
虫号: 8623502
注册: 2018-04-25
专业: 岩土与基础工程

[求助] 用1stOpt解非线性方程组老是出错已有1人参与

用1stOpt解方程组，老是提示出错，大神们给指点一下
NewDivision;
Constant E=207*10^9,D=0.711,t=7.9*10^(-3),p=4*10^6,Es=0.035*10^9, vs=0.35, L0=40,q=4.35*10^3,v=28.7*10^(-3);
Parameter M0[1,],S0[1,],L[0.1,];
ConstStr d0=D-2*t;
      I=3.14*(D^4-d0^4)/64;
      s=0.25*3.1415*(D^2-d0^2);
      k=3.08/1.35*(Es*D^4/E/I)^(1/8)*Es/(1-vs^2);
      C11=-q*L0/(12*E*I);
      C12=-M0/(2*E*I);
      P=-k*v/(24*E*I);
      beta0=(k/(4*E*I))^0.25;
      phi=1+exp(L0*sqrt(S0/(E*I)));
      da=1-exp(L0*sqrt(S0/(E*I)));
      D1=(M0*S0+q*E*I)/(phi*S0^2);
      D2=(M0*S0+q*E*I)/(phi*S0^2)*exp(L0*sqrt(S0/(E*I)));
      C13=-2*C12-(6*C11+24*L*P)/(2*beta0^2)-3*C11*L^2-4*L^3*P-(12*P*L^2+6*C11*L+2*C12)/beta0;
      C14=L*(4*P*L^3+3*C11*L^2+2*C12)-C12*L^2-C11*L^3-L^4*P+((L*beta0+1)*(6*C11+24*L*P))/(2*beta0^3)+((2*L*beta0+1)*(12*P*L^2+6*C11*L+2*C12))/(2*beta0^2);
Function
  P*L^4+C11*L^3+C12*L^2+C13*L+C14=v;
  D1*sqrt(S0/(E*I))*da+q*L0/(2*S0)=abs(C13)；
  D1*exp(L0/2*sqrt(S0/(E*I)))+D2*exp(-L0/2*sqrt(S0/(E*I)))+q*L0^2/(8*S0)+abs(C14)-M0/S0-q*E*I/S0^2=sqrt(4*L0^2*S0/(3.14^2*E*s));

回复此楼

» 猜你喜欢

北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！已经有1人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收10人已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有108人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人已经有0人回复
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8 已经有0人回复
0702一志愿吉大B区求调剂有论文已经有4人回复
欢迎调剂到赣南师范大学智能制造与未来能源学院已经有28人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有4人回复
欢迎调剂到赣南师范大学智能制造与未来能源学院已经有21人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有0人回复

1楼 2019-10-21 15:26:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+10, 感谢热心指导 2019-10-23 08:58:46
尚玉杰: 金币+5, ★★★很有帮助 2019-10-23 14:48:46

ConstStr 里面除了最后一个以 ; 结尾，其他全部改为 ,
第二个方程结尾 ; 输入法不对。。所有标点都需要在英文输入法下输入
另外，1stOpt 不区分大写小，参数要重新定义。

赞一下(1人)

回复此楼

数值计算

2楼2019-10-21 15:48:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 142.4
帖子: 39
在线: 6.3小时
虫号: 8623502
注册: 2018-04-25
专业: 岩土与基础工程

引用回帖:

2楼: Originally posted by 独孤神宇 at 2019-10-21 15:48:26
ConstStr 里面除了最后一个以 ; 结尾，其他全部改为 ,
第二个方程结尾 ; 输入法不对。。所有标点都需要在英文输入法下输入
另外，1stOpt 不区分大写小，参数要重新定义。

谢谢，按照大神的建议修改了，能运行，但是出来的结果还是有问题，请问应该怎么调整啊
这是修改后的
NewDivision;
Constant E=207*10^9,D=0.711,t=7.9*10^(-3),p=4*10^6,Es=0.035*10^9, vs=0.35, L0=40,q=4.35*10^3,v=28.7*10^(-3);
Parameter M0[1,],S0[1,],L[0.1,];
ConstStr d0=D-2*t,
      I=3.14*(D^4-d0^4)/64,
      s=0.25*3.1415*(D^2-d0^2),
      k=3.08/1.35*(Es*D^4/E/I)^(1/8)*Es/(1-vs^2),
      C11=-q*L0/(12*E*I),
      C12=-M0/(2*E*I),
      P=-k*v/(24*E*I),
      beta0=(k/(4*E*I))^0.25,
      phi=1+Exp(L0*Sqrt(S0/(E*I))),
      da=1-Exp(L0*Sqrt(S0/(E*I))),
      D1=(M0*S0+q*E*I)/(phi*S0^2),
      D2=(M0*S0+q*E*I)/(phi*S0^2)*Exp(L0*Sqrt(S0/(E*I))),
      C13=-2*C12-(6*C11+24*L*P)/(2*beta0^2)-3*C11*L^2-4*L^3*P-(12*P*L^2+6*C11*L+2*C12)/beta0,
      C14=L*(4*P*L^3+3*C11*L^2+2*C12)-C12*L^2-C11*L^3-L^4*P+((L*beta0+1)*(6*C11+24*L*P))/(2*beta0^3)+((2*L*beta0+1)*(12*P*L^2+6*C11*L+2*C12))/(2*beta0^2);
Function
      P*L^4+C11*L^3+C12*L^2+C13*L+C14=v;
      D1*Sqrt(S0/(E*I))*da+q*L0/(2*S0)=Abs(C13);
      D1*Exp(L0/2*Sqrt(S0/(E*I)))+D2*Exp(-L0/2*Sqrt(S0/(E*I)))+q*L0^2/(8*S0)+Abs(C14)-M0/S0-q*E*I/S0^2=Sqrt(4*L0^2*S0/(3.14^2*E*s));

运行结果为

====== Results ======
Iterations: 94
Elapsed Time (Hr:Min:Sec:Msec): 00:00:01:319
Stop Reason: Convergence tolerance reached
Algorithms: Max Implementation Optimization (MIO1)
Function 1: ((-(3.08/1.35*(35000000*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.0079)^
4)/64))^(1/8)*35000000
/(1-0.35^2))*0.0287/(24*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))
)*l^4+((-4350*40/(12*207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))))*l^3+((-m0/(2*207000000000*(3.14*(0
.711^4-(0.711-2*0.0079)^4)
/64))))*l^2+((-2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))
)-(6*(-4350*40/(12*207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))+24*l*4000000)/(2*(((3.08/1.35*(3500
0000*0.711^4/207000000000
/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^2))/(4*207000
000000*(3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))^0.25)^2)-3*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0.7
11-2*0.0079)^4)/64)))*l^2
-4*l^3*4000000-(12*4000000*l^2+6*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0
.711-2*0.0079)^4)/64)))
*l+2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))))/(((3.08/1
.35*(35000000*0.711^4/207000000000
/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^2))/(4*207000
000000*(3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))^0.25)))*l+((l*(4*4000000*l^3+3*(-4350*40/(12*207000000000
*(3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))*l^2+2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079
)^4)/64))))-(-m0/(2*207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))*l^2-(-4350*40/(12*207000000000*(3.1
4*(0.711^4-(0.711-2*0.0079)^4)
/64)))*l^3-l^4*4000000+((l*(((3.08/1.35*(35000000*0.711^4/207000000000/(3.14*
(0.711^4-(0.711-2*0.0079)^4)
/64))^(1/8)*35000000/(1-0.35^2))/(4*207000000000*(3.14*(0.711^4-(0.711-2*0.00
79)^4)/64)))^0.25)+1)*(6
*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))+24*l*400
0000))/(2*(((3.08/1.35
*(35000000*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)
*35000000/(1-0.35^2))/(4
*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)^3)+((2*l*(((3.08
/1.35*(35000000*0.711^4
/207000000000/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^
2))/(4*207000000000*(3.14
*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)+1)*(12*4000000*l^2+6*(-4350*40/(12*
207000000000*(3.14*(0.711^4
-(0.711-2*0.0079)^4)/64)))*l+2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0
.0079)^4)/64)))))/(2*(((3.08
/1.35*(35000000*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^
(1/8)*35000000/(1-0.35^2))
/(4*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)^2)))-(0.0287)
= 13593990.66

2: (((m0*s0+4350*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))/((1+exp(40
*sqrt(s0/(207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))))*s0^2)))*sqrt(s0/(207000000000*((3
.14*(0.711^4-(0.711-2*0.0079)^4)
/64))))*((1-exp(40*sqrt(s0/(207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/6
4))))))+4350*40/(2*s0)
-(abs(((-2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))-(6*(
-4350*40/(12*207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))+24*l*4000000)/(2*(((3.08/1.35*(3500
0000*0.711^4/207000000000
/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^2))/(4*207000
000000*(3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))^0.25)^2)-3*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0.7
11-2*0.0079)^4)/64)))*l^2
-4*l^3*4000000-(12*4000000*l^2+6*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0
.711-2*0.0079)^4)/64)))
*l+2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))))/(((3.08/1
.35*(35000000*0.711^4/207000000000
/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^2))/(4*207000
000000*(3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))^0.25))))) = -17370977.29

3: (((m0*s0+4350*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))/((1+exp(40
*sqrt(s0/(207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))))*s0^2)))*exp(40/2*sqrt(s0/(2070000
00000*((3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))))+(((m0*s0+4350*207000000000*(3.14*(0.711^4-(0.711-2*0.00
79)^4)/64))/((1+exp(40
*sqrt(s0/(207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))))*s0^2)*exp(4
0*sqrt(s0/(207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))))))*exp(-40/2*sqrt(s0/(207000000000*
((3.14*(0.711^4-(0.711
-2*0.0079)^4)/64)))))+4350*40^2/(8*s0)+abs(((l*(4*4000000*l^3+3*(-4350*40/(12
*207000000000*(3.14*(0.711^4
-(0.711-2*0.0079)^4)/64)))*l^2+2*(-m0/(2*207000000000*(3.14*(0.711^4-(0.711-2
*0.0079)^4)/64))))-(-m0
/(2*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))*l^2-(-4350*40/(12*2
07000000000*(3.14*(0.711^4
-(0.711-2*0.0079)^4)/64)))*l^3-l^4*4000000+((l*(((3.08/1.35*(35000000*0.711^4
/207000000000/(3.14*(0.711^4
-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/(1-0.35^2))/(4*207000000000*(3.14*(0
.711^4-(0.711-2*0.0079)^4)
/64)))^0.25)+1)*(6*(-4350*40/(12*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)
^4)/64)))+24*l*4000000))
/(2*(((3.08/1.35*(35000000*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.007
9)^4)/64))^(1/8)*35000000
/(1-0.35^2))/(4*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)^3
)+((2*l*(((3.08/1.35*(35000000
*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.0079)^4)/64))^(1/8)*35000000/
(1-0.35^2))/(4*207000000000
*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)+1)*(12*4000000*l^2+6*(-4350*4
0/(12*207000000000*(3.14
*(0.711^4-(0.711-2*0.0079)^4)/64)))*l+2*(-m0/(2*207000000000*(3.14*(0.711^4-(
0.711-2*0.0079)^4)/64)))))
/(2*(((3.08/1.35*(35000000*0.711^4/207000000000/(3.14*(0.711^4-(0.711-2*0.007
9)^4)/64))^(1/8)*35000000
/(1-0.35^2))/(4*207000000000*(3.14*(0.711^4-(0.711-2*0.0079)^4)/64)))^0.25)^2
))))-m0/s0-4350*207000000000
*((3.14*(0.711^4-(0.711-2*0.0079)^4)/64))/s0^2-(sqrt(4*40^2*s0/(3.14^2*207000
000000*((0.25*3.1415*(0.711^2
-(0.711-2*0.0079)^2)))))) = 11946658.43
Objective Function (Min.): 629270081574722
m0: 2.714223949163E15
s0: 60767551616.1643
l: 0.1
====== Finished ======

赞一下

回复此楼

3楼2019-10-22 10:15:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

1stOpt 不区分大写小，参数要重新定义，不要用大小写区分。

赞一下

回复此楼

数值计算

4楼2019-10-22 10:22:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 142.4
帖子: 39
在线: 6.3小时
虫号: 8623502
注册: 2018-04-25
专业: 岩土与基础工程

引用回帖:

4楼: Originally posted by 独孤神宇 at 2019-10-22 10:22:48
1stOpt 不区分大写小，参数要重新定义，不要用大小写区分。

不同参数使用了不同字母表示，大小写的使用是为了跟自身定义相符

发自小木虫IOS客户端

赞一下

回复此楼

5楼2019-10-22 10:50:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 142.4
帖子: 39
在线: 6.3小时
虫号: 8623502
注册: 2018-04-25
专业: 岩土与基础工程

我在修改试试

发自小木虫IOS客户端

回复此楼

6楼2019-10-22 10:52:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

引用回帖:

5楼: Originally posted by 尚玉杰 at 2019-10-22 10:50:27
不同参数使用了不同字母表示，大小写的使用是为了跟自身定义相符
...

不是说了吗，软件本身不区分大小写。。这么定义肯定有问题的。。

赞一下

回复此楼

数值计算

7楼2019-10-22 10:52:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主尚玉杰的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 考研调剂 +5	Amber00 2026-03-31	5/250	2026-03-31 22:43 by zhengzhi8909
[考研] 调剂申请 +7	张张张张zy 2026-03-31	7/350	2026-03-31 22:25 by piaoyi415
[考研] 352分-085602-一志愿985 +6	海纳百川Ly 2026-03-29	6/300	2026-03-31 21:06 by yuq
[考研] 284求调剂 +9	小熊～～ 2026-03-31	9/450	2026-03-31 18:22 by 253863592
[考研] 求调剂推荐材料 304 +18	荷包蛋hyj 2026-03-26	18/900	2026-03-31 18:08 by 544594351
[考研] 0856 335分 \| 封装or激光加工老师看看我！！ +9	cccchenso 2026-03-29	9/450	2026-03-31 16:37 by lishahe
[考研] 311求调剂 +9	蓝月亮亮 2026-03-30	9/450	2026-03-31 16:32 by yedezhan
[考研] 化学工程085602 305分求调剂 +28	RichLi_ 2026-03-25	36/1800	2026-03-31 14:56 by JourneyLucky
[考博] 材料专业申博 +5	杜雨婷dyt 2026-03-29	5/250	2026-03-31 11:19 by oooqiao
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +6	星沉uu 2026-03-26	7/350	2026-03-31 10:19 by GdShizy
[考研] 0817化工学硕调剂 +7	努力上岸中！ 2026-03-31	7/350	2026-03-31 09:58 by nalakaiqi
[考研] 332求调剂 +14	032500 2026-03-25	14/700	2026-03-30 20:03 by 无际的草原
[考研] 293求调剂 +3	末未mm 2026-03-30	5/250	2026-03-30 17:23 by 王保杰33
[考研] 071010 323 分求调剂 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 求调剂 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐儿
[硕博家园] 招收生物学/细胞生物学调剂 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 数一英一271专硕（085401）求调剂，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木虫tim
[考研] 286求调剂 +12	PolarBear11 2026-03-26	12/600	2026-03-28 12:14 by zllcz
[考研] 348求调剂 +4	小懒虫不懒了 2026-03-27	5/250	2026-03-27 12:47 by 果果妈咪
[考研] 求调剂 +3	李李不服输 2026-03-25	3/150	2026-03-25 13:03 by cmz0325

24小时热门版块排行榜

尚玉杰

[求助] 用1stOpt解非线性方程组老是出错 已有1人参与

» 猜你喜欢

独孤神宇

【答案】应助回帖

尚玉杰

独孤神宇

【答案】应助回帖

尚玉杰

尚玉杰

独孤神宇

【答案】应助回帖

[求助] 用1stOpt解非线性方程组老是出错已有1人参与