- 积分
- 9
- 注册时间
- 2006-3-24
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2012-11-2 16:50:11
|
显示全部楼层
来自 陕西西安
1.ANSYS通常使用的是蒙特卡罗法和响应平面法。
2.下面为一个使用蒙特卡罗法计算的可靠度实例,可根据各位的需要,作类似的修改即可。
问题描述:单边约束矩形板,一个角点受一垂直于板的集中载荷。
本例题将板长(威布尔分布)、板宽(均匀分布)和集中载荷(贝塔分布)定义为随机变量,计算最大应力不超过阈值ss(服从威布尔分布)的概率。
本例题需要两个文件夹aa.txt和bb.txt,将二者放置于工作目录,然后input from读入bb.txt,即可进行运行计算。计算结果可从
ansys main menu>prob design>prob results中相关选项读取。
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!将下列命令流保存为aa.txt!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
length=250
width=140
force=400
ss=200
/prep7
et,1,shell63
r,1,10
mp,ex,1,2.1e5
mp,nuxy,1,0.3
k,1
k,2,length
k,3,length,width
k,4,,width
a,1,2,3,4
esize,4
amesh,all
/solu
nsel,s,loc,x,0
d,all,all
nsel,s,loc,x,length
nsel,r,loc,y,width
f,all,fz,-force
allsel
solve
/post1
nsort,s,eqv,1,1,,0
*get,smax,sort,,max
subs=ss-smax
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!将下列命令流保存为bb.txt!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/inp,aa,txt
/pds
pdanl,aa,txt
pdvar,length,weib,2,250,100
pdvar,width,unif,139,141
pdvar,force,beta,1,3,,600
pdvar,ss,weib,2,280,150
pdvar,subs,resp
pdmeth,mcs,lhs
pdlhs,100,1,rand,,all,,,,init
pdexe,lhsrun
计算模型
部分计算结果示意图
注:上例蒙特卡罗法仅计算了100次,通常计算时所需要的计算次数为:100/失效概率,假如失效概率为1%,即可靠度是99%,那么需要的计算次数为100/0.01=10000次。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|