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

登录注册

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

24小时热门版块排行榜石溪大学接受考研调剂申请>

>论坛更新日志 (7230)
>休闲灌水 (594)
>虫友互识 (505)
>硕博家园 (486)
>文献求助 (467)
>导师招生 (457)
>论文投稿 (300)
>基金申请 (247)
>考研 (243)
>考博 (219)
>招聘信息布告栏 (194)
>博后之家 (182)
>公派出国 (119)
>教师之家 (114)
>有机交流 (111)
>论文道贺祈福 (77)

【调剂】北京石油化工学院2024年16个专业接受调剂

返回列表

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 137.4
帖子: 38
在线: 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));

回复此楼

» 猜你喜欢

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

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

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29537.8
散金: 802
红花: 120
沙发: 1
帖子: 5338
在线: 850.2小时
虫号: 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 (幼儿园)
金币: 137.4
帖子: 38
在线: 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的回帖

独孤神宇

版主 (知名作家)

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

【答案】应助回帖

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

赞一下

回复此楼

数值计算

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

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

尚玉杰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 137.4
帖子: 38
在线: 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 (幼儿园)
金币: 137.4
帖子: 38
在线: 6.3小时
虫号: 8623502
注册: 2018-04-25
专业: 岩土与基础工程

我在修改试试

发自小木虫IOS客户端

回复此楼

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

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

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29537.8
散金: 802
红花: 120
沙发: 1
帖子: 5338
在线: 850.2小时
虫号: 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的回帖

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

返回列表

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

[找工作] 普通院校药学硕士，做合成的，感觉找不到工作 +8	pom戴墨镜 2024-04-24	15/750	2024-04-25 20:50 by pom戴墨镜
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +12	彩虹初见 2024-04-24	12/600	2024-04-25 19:21 by 雪径踏青
[教师之家] 刚博士毕业，入职普通一本高校 +18	保持肌肤 2024-04-20	19/950	2024-04-25 18:30 by lyfbangong
[考博] 取博导收留 5+3	zzb777888 2024-04-20	8/400	2024-04-25 17:40 by 那年的旧城
[教师之家] 博士论文被抄袭 +23	和尚敲小木鱼 2024-04-22	39/1950	2024-04-25 15:48 by 和尚敲小木鱼
[博后之家] 南京大学-广州大学联合招聘博士后欢迎广大优秀人才！！！ +4	黑魔变身啾 2024-04-20	12/600	2024-04-25 11:18 by dodonaomi
[基金申请] "颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人 +3	zju2000 2024-04-24	3/150	2024-04-25 10:00 by shenyeqing
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函 +33	bjdxyxy 2024-04-19	53/2650	2024-04-25 08:05 by 大雁高飞
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+4	yifeng11 2024-04-20	10/500	2024-04-25 02:14 by yifeng11
[考研] 博士生招生！！ +5	2458827884 2024-04-19	5/250	2024-04-24 22:28 by 不懂的风清
[考博] 24年申博化学/材料一作6篇sci +9	wangyp123 2024-04-23	11/550	2024-04-24 19:01 by bangbangbiu
[基金申请] 化学结清有情发出来了 +3	starboy7286 2024-04-20	3/150	2024-04-24 15:52 by 嘿，黑贝
[有机交流] 紧急求助，有谁用过三甲基硅醇钾TMSOK吗？ 1000+3	genius2008 2024-04-19	10/500	2024-04-24 12:37 by genius2008
[教师之家] 大家访学都是怎么找的啊？ +3	luokereng 2024-04-22	3/150	2024-04-24 11:40 by xuechenli
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4	我属驴核动力驴 2024-04-22	5/250	2024-04-24 08:47 by 晓目崇
[教师之家] 找工作过程中碰到好多人问南京某新换校长211的政策，规定如下 +15	ZHONGWU_U 2024-04-21	17/850	2024-04-24 05:12 by juta256
[教师之家] 发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？ +11	河西夜郎 2024-04-19	11/550	2024-04-23 15:01 by mylifecg
[考博] 研二光催化6月底4篇2区 +7	wjtab 2024-04-22	11/550	2024-04-23 06:59 by byron2012
[论文投稿] 研究光催化的，好中的三四区 20+3	sl.0117 2024-04-20	3/150	2024-04-22 09:53 by bobvan
[硕博家园] 材料工程，二本硕士，值得读吗 +3	简风23 2024-04-20	6/300	2024-04-21 12:19 by 简风23

24小时热门版块排行榜 石溪大学接受考研调剂申请>

尚玉杰

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

» 猜你喜欢

独孤神宇

【答案】应助回帖

尚玉杰

独孤神宇

【答案】应助回帖

尚玉杰

尚玉杰

独孤神宇

【答案】应助回帖

24小时热门版块排行榜石溪大学接受考研调剂申请>