Ansys讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

gghhjj 发表于 2005-10-17 07:45

一个结构可靠性分析实例

<DIV >问题描述：桁架结构如图所示，各节点坐标为：1（0，0，0），2（10，0，0），3（20，0，0），4（10，-10，0）（单位：cm）。1、2、3节点完全固定，节点4受集中力作用，FX＝20000N，FY＝20000N。杨氏模量EX＝2.1e5N/cm2。单元1的横截面积A1~N（10，0.5），单元2的横截面积A2服从三角分布，Xmin=10，Xmlv=11，Xmax=12；单元3的横截面积A3服从均匀分布，Xmin=9，Xmax=11；A1、A3的相关系数为0.2。试分析杆件横截面积的不确定性对结构体积、杆单元轴向应力的影响。 本题采用LINK1单元，直接法构建有限元模型求解，提取结构体积、各单元的轴向力作为可靠性分析输出结果变量，选择各杆横截面积作为输入变量，采用蒙特卡罗法中的直接法进行分析。以下是求解的命令流： *CREATE, pds3bar !生成宏pds3bar作为可靠性性分析文件 *SET, A1, 10 !初始化设计变量 *SET, A2, 10 *SET, A3, 10 /PREP7 ET, 1, 1 EX, 1, 2.1E5 R, 1, A1 R, 2, A2 R, 3, A3 N, 1 N, 2, 10 N, 3, 20 N, 4, 10, -10 !生成有限元模型 REAL, 1 E, 1, 4 REAL, 2 E, 2, 4 REAL, 3 E, 3, 4 FINISH /SOLU D, 1, ALL,,, 3 F, 4, FX, 20000 F, 4, FY, -20000 SOLVE FINISH /POST1 SET, 1 ETABLE, VOLU, VOLU !将单元的体积放入表VOLU中 ETABEL, AXST, LS, 1 !将单元轴向应力放入表AXST中 *GET, SIG1, ELEM, 1, ETAB, AXST !SIG1＝单元1的轴向应力 *GET, SIG2, ELEM, 2, ETAB, AXST *GET, SIG3, ELEM, 3, ETAB, AXST SSUM !将单元表格内数据求和 *GET,TVOL, SSUM,, ITEM, VOLU !提取结构的总体积 FINISH *end !结束宏定义 *use, pds3bar !执行宏操作 /PDS !进入可靠性分析模块 pdanl, pds3bar !指定分析文件 !定义输入变量 PDVAR, A1, GAUS, 10,0.5 PDVAR, A2, tria, 10,11,12 PDVAR, A3, unif, 9, 11 PDCOR, A1, A3, 0.2 !定义输出变量 PDVAR, SIG1, resp PDVAR, SIG2, resp PDVAR, SIG3, resp PDVAR, TVOL, resp PDMETH, MCS, DIR !指定分析方法 PDDMCS, 100, NONE, ALL,,,, 123457 !设定分析循环次数 PDEXE, mcs3bar !执行分析 PDSENS, MCS3BAR, TVOL　　　！TVOL灵敏度显示 PDSHIS, MCS3BAR, SIG1, SMAP　　！SIG1抽样过程显示 PDHIST, MCS3BAR, SIG2　　　　　　！绘制SIG2取值分布柱状图 PDCDF, MCS3BAR, SIG3　　　　　　　！绘制显示SIG3分布函数</DIV>

jhmgirl 发表于 2008-4-29 13:40

谢谢分享！正在学习，论坛里关于可靠性的好少啊！不知有没有关于振动可靠性的例子？

aircraft 发表于 2008-4-29 15:31

大家有做疲劳分析的吗?
一起讨论讨论,呵呵.

zhoushengren 发表于 2013-5-23 11:07

谢谢楼主分享

页: [1]

声振论坛's Archiver

一个结构可靠性分析实例