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

登录注册

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

24小时热门版块排行榜

返回列表

尚玉杰

新虫 (初入文坛)

应助: 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));

回复此楼

» 猜你喜欢

欢迎加入课题组已经有0人回复
散金币，求好运，祝面上顺利！已经有80人回复
物理学I论文润色/翻译怎么收费? 已经有265人回复
华东师范大学芯片设计徐珑真实水平咋样已经有1人回复
调剂已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有0人回复
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职已经有0人回复
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有2人回复
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。已经有0人回复
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！已经有1人回复

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的回帖

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

返回列表

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

[考研] 07化学303求调剂 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 材料与化工304求B区调剂 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 0854电子信息求调剂 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 0854人工智能方向招收调剂 +4	章小鱼567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 材料与化工考研调剂 +9	孅華 2026-03-22	9/450	2026-03-25 13:09 by cmz0325
[考研] 081700 调剂 267分 +10	迷人的哈哈 2026-03-23	10/500	2026-03-25 12:12 by userper
[考研] 0703化学调剂，求导师收 +7	天天好运来上岸� 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 材料专硕找调剂 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 070300化学求调剂 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 一志愿211 初试270分求调剂 +5	谷雨上岸 2026-03-23	6/300	2026-03-24 16:32 by laoshidan
[考研] 【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 种大树
[考研] 336化工调剂 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 350求调剂 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 求老师收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 一志愿070300浙大化学358分，求调剂！ +4	酥酥鱼.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 寻找调剂 +4	倔强芒? 2026-03-21	4/200	2026-03-22 16:14 by 木托莫露露
[考研] 初试 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 一志愿深大，0703化学，总分302，求调剂 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 学员8dgXkO
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng

24小时热门版块排行榜

尚玉杰

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

» 猜你喜欢

独孤神宇

【答案】应助回帖

尚玉杰

独孤神宇

【答案】应助回帖

尚玉杰

尚玉杰

独孤神宇

【答案】应助回帖

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