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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (2567)
>考研 (293)
>导师招生 (255)
>虫友互识 (59)
>休闲灌水 (47)
>基金申请 (31)
>考博 (22)
>文献求助 (21)
>博后之家 (18)
>教师之家 (17)
>论文投稿 (17)
>论文道贺祈福 (16)
>硕博家园 (13)
>找工作 (11)
>公派出国 (9)
>生物科学 (5)

北京石油化工学院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));

回复此楼

» 猜你喜欢

调剂已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有261人回复
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职已经有0人回复
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有2人回复
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。已经有0人回复
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！已经有1人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收10人已经有0人回复
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人已经有0人回复

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

已阅回复此楼关注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的回帖

查看全部 7 个回答

独孤神宇

版主 (知名作家)

应助: 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的回帖

查看全部 7 个回答

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

[考研] 求调剂推荐材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 289求调剂 +8	新时代材料 2026-03-27	8/400	2026-03-28 04:09 by fmesaito
[考研] 材料292调剂 +12	橘颂思美人 2026-03-23	12/600	2026-03-27 15:44 by caszguilin
[考研] 283求调剂（080500） +4	A child 2026-03-27	4/200	2026-03-27 15:34 by XPU李庆
[考研] 08开头275求调剂 +4	拉谁不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 333求调剂 +3	question挽风 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚的貓
[考研] 求调剂 +3	刘柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 求调剂 +6	林之夕 2026-03-24	6/300	2026-03-27 08:38 by hypershenger
[考研] 333求调剂 +6	wfh030413@ 2026-03-23	6/300	2026-03-26 22:45 by 学员8dgXkO
[考研] 321求调剂 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 297求调剂 +6	田洪有 2026-03-26	6/300	2026-03-26 15:55 by 不吃魚的貓
[考研] 一志愿南京邮电大学 288分材料考研求调剂 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7	niko- 2026-03-22	7/350	2026-03-25 20:14 by qingfeng258
[考研] 080500求调剂 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 305分求调剂（食品工程） +5	Sxy112 2026-03-21	7/350	2026-03-24 12:27 by 544594351
[考研] 一志愿吉大化学322求调剂 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 环境学硕288求调剂 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 热情沙漠
[考研] 276求调剂 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在风里
[考研] 285求调剂 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 求调剂 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike

24小时热门版块排行榜

尚玉杰

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

» 猜你喜欢

尚玉杰

独孤神宇

【答案】应助回帖

尚玉杰

独孤神宇

【答案】应助回帖

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