| 查看: 460 | 回复: 0 | ||
[求助]
可靠度模拟求助
|
|
最近在做可靠度分析，但是运行完了之后，在画柱状分析图的时候总是出现 “ERROR in EvalHistogramBars. Number of data points found is 0. Number data points expected is 2.” 请大家指教究竟是怎么回事。 命令流： Finish /clear *create,pdspavement,mac *set,ex1,7000e6 /prep7 et,1,solid185 !定义单元 mp,ex,1,ex1 mp,prxy,1,0.35 mp,dens,1,2100 l=4 !...........x方向宽度 w=5 !...........z方向长度 wh3=0.04 !...........y方向高度 Side_kuandu=18E-2 !荷载区域边长（18*20cm） Side_xingche=20E-2 blc4,0,0,l,wh3,w !*************荷载区域******* * b=Side_kuandu+0.3 !............荷载中心距离 d=Side_kuandu !............荷载边长 a=Side_xingche allsel,all wpave,0,h, w/2+5*a/2 vsbw,all wpave,0,h, w/2-5*a/2 vsbw,all WPROTA,0,0,90 wpave,l/2-b/2-d/2,h,0 vsbw,all wpave,l/2-b/2+d/2,h,0 vsbw,all wpave,l/2+b/2-d/2,h,0 vsbw,all wpave,l/2+b/2+d/2,h,0 vsbw,all wprota,0,0,-90 wpave,0,0,0 allsel,all vglue,all vsel,s,loc,y,0,wh3 VATT,1,,1,0 !*************划分网格********* lsel,s,length,,d lesize,all,,,4,,,,,1 lsel,s,length,,b-d LESIZE,all,,,2,,,,,1 lsel,s,length,,l/2-b/2-d/2 LESIZE,all,,,10,,,,,1 lsel,s,length,,w/2-5*a/2 LESIZE,all,,,14,,,,,1 lsel,s,length,,5*a LESIZE,all,,,20,,,,,1 lsel,s,length,,wh3 LESIZE,all,,,2,,,,,1 !............SMA-13 2份 MSHAPE,0,3d MSHKEY,1 allsel,all vmesh,all nsel,s,loc,x,0 !路面约束 nsel,a,loc,x,l d,all,ux,0 nsel,s,loc,z,0 nsel,a,loc,z,w d,all,uz,0 nsel,s,loc,y,0 d,all,all,0 Nsel,s,loc,x, l/2-b/2-d/2, l/2-b/2+d/2 Nsel,r,loc,y,wh3 Cm,movend1,node Nsel,s,loc,x, l/2+b/2-d/2, l/2+b/2+d/2 Nsel,r,loc,y,wh3 Cm,movend2,node po=0.7E6 !静荷载（Pa） elemdl=0.1 !每2个单元的长度 outres,all,all !输出所有荷载步的结果信息 NROPT,init !use the previously computed matrix(initial-stiffness) Sstif,off !关闭应力刚化效应 finish !****************施加荷载************* /solu ANTYPE,trans,new !分析类型为瞬态 Time,1 kbc,1 !加载为阶跃方式 !*****************添加荷载*****************左侧 Cmsel,s,movend1,node Nsel,r,loc,z, 20*elemdl,22*elemdl Esln,s,1 Esel,r,type,,2 sfe,all,1,pres,0, po !*****************添加荷载*****************右侧 Cmsel,s,movend2,node Nsel,r,loc,z, 20*elemdl,22*elemdl Esln,s,1 Esel,r,type,,2 sfe,all,1,pres,0, po allsel,all solve Save Finish /post1 Set,first Etable,s_y,s,y Etable,volu,volu Allsel,all Esel,r,mat,,1 Nsort,s,int,0,0 *get,max_int_SMA13,sort,0,max *get,max_int_node_SMA13,sort,0,imax Nsel,s,,,max_int_node_SMA13 Esln,s,0 *get,elem1,elem,0,num,max *get,elem2,elem,0,num,min *get,sig1,elem,elem1,etab,s_y *get,sig2,elem,elem2,etab,s_y Ssum *get,tvol,ssum,,item,volu Finish *end !结束宏定义 /input,pdspavement,mac !执行宏操作 /pds !进入可靠性分析模块 Pdanl,pdspavement,mac !指定分析文件 Pdvar,ex1,log1,7000e6,0.13 !服从对数正态分布 Pdvar,sig1,resp !输出变量 Pdvar,sig2,resp Pdvar,tvol,resp Pdmeth,mcs,dir !设置可靠度分析方法，采用monte carlo抽样法 Pddmcs,100,,all,,,,1234567 !设置分析循环次数，初始化种子值1234567 Pdexe,mcspavement !执行可靠度分析 Pdsens,mcspavement,tvol,both,rank,0.025 !创建敏感性图表 Pdshis,mcspavement,sig1,samp !sig1抽样过程显示 Pdhist,mcspavement,sig2 !sig2取值分布柱状图 Pdcdf,mcspavement,sig2 !绘图显示sig2分布函数 /pds Rsfit,rsmpavement,mcspavement,sig1,quad,none,0,fsr,0.95 !拟合响应面方程 Rsprnt,rsmpavement,sig1,no !查看响应面方程 |
回复此楼» 猜你喜欢
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
2026国自然函评费到账
已经有6人回复
-
科研孤儿太难了
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有12人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
40