pds-pipe.mac

采用ANSYS软件的APDL语言编写的两个命令流文件

FINISH
/CLEAR
/FILNAM,PDS-OilPipe
/TITLE,PROBABILISTIC DESIGN OF AN UNDERGROUND OIL PIPELINE
/REPLOT

/inp,pds-pipe-loop,mac                                !执行初始化的分析过程

/PDS !进入PDS处理器
PDANL,pds-pipe-loop,mac

!定义随机输入变量
PDVAR,DIAMETER,GAUS,420,1.2                           !管道外径服从高斯分布,均值为420mm,标准差为1.2mm
PDVAR,THICKNESS,GAUS,10,0.2                           !管道壁厚服从高斯分布,均值为10mm,标准差为0.2mm
PDVAR,YOUNG,GAUS,2.1E5,0.063E5                        !弹性模量服从高斯分布,均值为2.1E5MPa,标准差为0.063E5MPa
PDVAR,TANG,GAUS,0.81E5,0.024E5                        !剪切模量服从高斯分布,均值为0.81E5MPa,标准差为0.024E5MPa
PDVAR,POSONG,GAUS,0.3,0.0075                          !泊松比服从高斯分布,均值为0.3,标准差为0.0075
PDVAR,DEPTH,TGAU,3.02,0.0459,0,THICKNESS              !腐蚀缺陷深度服从截断高斯分布,均值为3.02mm,标准差为0.0459mm
PDVAR,PRESSURE,GAUS,7.64,0.46                         !输送压力服从高斯分布,均值为7.64MPa,标准差为0.46MPa
PDVAR,YIELDSTRESS,GAUS,312,6.78                       !屈服强度服从高斯分布,均值为312MPa,标准差为6.78MPa
PDVAR,RESISTANCE,GAUS,450,9.56                        !抗拉强度服从高斯分布,均值为450MPa,标准差为9.56MPa

!绘制随机输入参数的分布函数图
/TITLE,PROBABILITY DENSITY FUNCTION AND CUMULATIVE DISTRIBUTION FUNCTION OF RESISTANCE
PDPLOT,RESISTANCE, , ,                                !抗拉强度的分布函数图

!定义随机输出变量
PDVAR,MAXSTRESS,RESP                                  !绝对值最大应力
PDVAR,LIMITSTATUSFUN,RESP                             !极限状态函数值

!设置最大循环数
MAXLOOP=1000

!选择概率设计方法
PDMETH,MCS,LHS                                        !选择拉丁超立方蒙特卡罗抽样法
PDLHS,MAXLOOP,1,RAND,,'ALL',,,,INIT                   !执行所有MAXLOOP个随机抽样点，初始化种子值为123457

!执行概率设计仿真循环计算，并定义结果标识字为LHSPIPE
PDEXE,LHSPIPE,SER,0,COPY,LHSPIPE

!概率设计结果后处理
!查看概率设计参数的抽样历史曲线
PDSHIS,LHSPIPE,MAXSTRESS,SAMP,1,                 
PDSHIS,LHSPIPE,MAXSTRESS,MEAN,1,
PDSHIS,LHSPIPE,LIMITSTATUSFUN,SAMP,1,
PDSHIS,LHSPIPE,LIMITSTATUSFUN,MEAN,1,

!查看概率参数值的柱状分布图
PDHIST,LHSPIPE,MAXSTRESS, ,REL
PDHIST,LHSPIPE,LIMITSTATUSFUN, ,REL

!查看概率参数值的累积分布函数图
PDCDF,LHSPIPE,MAXSTRESS,EMP,0.95,100,
PDCDF,LHSPIPE,LIMITSTATUSFUN,EMP,0.95,100,

!计算并列表显示概率响应参数在指定范围内出现的概率
PDPROB,LHSPIPE,MAXSTRESS,LT,450E6,'    ',0.95,
PDPROB,LHSPIPE,LIMITSTATUSFUN,LT,0,'    ',0.95,

!计算并列表显示概率响应参数在给定概率条件下的极限值
PDPINV,LHSPIPE,MAXSTRESS,0.9,'    ',0.95,
PDPINV,LHSPIPE,LIMITSTATUSFUN,0.9,'    ',0.95,

!计算并绘制概率参数的灵敏度图
PDSENS,LHSPIPE,MAXSTRESS,BOTH,RANK,0.025,
PDSENS,LHSPIPE,LIMITSTATUSFUN,BOTH,RANK,0.025, 

!计算并绘制概率参数的离散图
PDSCAT,LHSPIPE,DEPTH,MAXSTRESS,POLY,1,,          !绘制最大应力对腐蚀缺陷深度的散点图
PDSCAT,LHSPIPE,DEPTH,LIMITSTATUSFUN,POLY,1,,     !绘制极限状态函数对腐蚀缺陷深度的散点图

!计算并列表显示概率相关系数矩阵
PDCMAT,LHSPIPE,IO,'        ','        ',RANK,0.025,1

!自动生成概率设计报告
/TITLE,
PDROPT,0,0,2,3,2,2,2,5,0.95,
PDWRITE,Report,Liu,Xiaodong 
/REPLOT,RESIZE
SAVE
FINISH
/EXIT,ALL