一个结构可靠性分析实例

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

<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

谢谢楼主分享

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 一个结构可靠性分析实例

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。