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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (2922)
>虫友互识 (344)
>导师招生 (161)
>考博 (148)
>文献求助 (130)
>论文投稿 (79)
>基金申请 (74)
>硕博家园 (66)
>博后之家 (53)
>找工作 (53)
>仿真模拟 (49)
>休闲灌水 (48)
>招聘信息布告栏 (37)
>公派出国 (34)
>教师之家 (33)
>考研 (33)

返回列表

尚玉杰

新虫 (初入文坛)

应助: 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年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有134人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复
大豆异黄酮分离已经有0人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有3人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复

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	Peterchao 2026-05-18	8/400	2026-05-24 11:58 by 晓晓爱翠翠
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3	bjvtcliu 2026-05-24	5/250	2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3	bjvtcliu 2026-05-24	5/250	2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！ +4	zju2000 2026-05-21	5/250	2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4	ZXW0202 2026-05-22	9/450	2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	hvkbtfonbv 2026-05-23	3/150	2026-05-24 08:01 by 9ps9vgkqva
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	pmo95bazuy 2026-05-23	7/350	2026-05-24 06:35 by fpo5ljpv91
[基金申请] 揭秘青基评审内幕：几个A才能顺利中标 +3	国自然国社科中 2026-05-23	4/200	2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5	chemBioBro 2026-05-22	7/350	2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3	焦晓明 2026-05-21	3/150	2026-05-23 11:26 by mlc840311
[文学芳草园] 献血感触 +7	呀呀好傻 2026-05-19	13/650	2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +14	1234567wang 2026-05-17	16/800	2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4	无名者登山 2026-05-20	5/250	2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7	无名者登山 2026-05-18	11/550	2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9	fummck 2026-05-20	10/500	2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15	popular289 2026-05-18	26/1300	2026-05-21 10:35 by 西葫芦炒鸡蛋
[有机交流] 反应很差，大量原料没有反应 5+3	Mr.Zot 2026-05-19	8/400	2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3	铁达火车 2026-05-18	5/250	2026-05-20 09:33 by tfang
[考博] 博士申请 +5	星…… 2026-05-18	6/300	2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3	Love_MH 2026-05-17	3/150	2026-05-18 02:22 by 竹林孤影

24小时热门版块排行榜

尚玉杰

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

» 猜你喜欢

独孤神宇

【答案】应助回帖

尚玉杰

独孤神宇

【答案】应助回帖

尚玉杰

尚玉杰

独孤神宇

【答案】应助回帖

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