土木工程数值法结课作业姓名%%%学号100000000专业结构方向2015年10月具体操作图文解析如下一.Winkler地基上的方形板受中心集中力作用时的挠度1、进入preferences选中structuralGUI图形界面过滤对话框2、定义单元类型单元类型选择单元对话框3、定义实常数定义单元实常数对话框4、定义材料特性定义材料属性对话框5、建立单元模型地基板模型示意图6、划分网格生成有限元模型划分网格单元尺寸对话框7、网格划分后的有限元模型网格划分后的有限元模型8、求解选择分析类型GUI分析类型选择对话框施加约束采用点约束形式:固定四个顶点及各边中点施加位移约束对话框施加荷载施加压力荷载对话框求解求解完成对话框9、查看结果Z方向位移云图二.框架结构在自重作用下的变形及内力、结构在平面内水平振动的前5阶自振频率及相应的振型、结构在振动荷载(水平加速度为0.1*g*sin(2π*2*t))作用下的振动反应(10秒内)。1、自重作用下的变形及内力操作步骤如下:1)进入preferences选中structuralGUI图形界面过滤对话框2)定义单元类型3)定义实常数4)定义材料特性定义弹模和泊松比定义密度5)建立模型6、划分网格7)求解选择分析类型设置约束条件定义重力加速度求解完成8、查看结果结点位移云图内力图位移值如下所示:PRINTUNODALSOLUTIONPERNODE*****POST1NODALDEGREEOFFREEDOMLISTING*****LOADSTEP=1SUBSTEP=1TIME=1.0000LOADCASE=0THEFOLLOWINGDEGREEOFFREEDOMRESULTSAREINTHEGLOBALCOORDINATESYSTEMNODEUXUYUZUSUM10.00000.00000.00000.000020.17078E-02-0.28500E-010.00000.28551E-0130.42382E-01-0.53750E-020.00000.42721E-0140.13566-0.10500E-010.00000.1360750.22903-0.15375E-010.00000.2295560.27170-0.20000E-010.00000.2724370.21286-0.24375E-010.00000.214258-0.24698E-02-0.40500E-010.00000.40575E-0190.43133-0.39125E-010.00000.43310100.47235-0.37500E-010.00000.47383110.28524-0.35625E-010.00000.28746120.34666E-01-0.33500E-010.00000.48208E-0113-0.11473-0.31125E-010.00000.11888140.24698E-02-0.40500E-010.00000.40575E-0115-0.19758E-02-0.938310.00000.9383216-0.14819E-02-2.02550.00002.025517-0.98792E-03-3.01040.00003.010418-0.49396E-03-3.68470.00003.6847190.11121E-12-3.92330.00003.9233200.49396E-03-3.68470.00003.6847210.98792E-03-3.01040.00003.0104220.14819E-02-2.02550.00002.0255230.19758E-02-0.938310.00000.9383224-0.17078E-02-0.28500E-010.00000.28551E-0125-0.43133-0.39125E-010.00000.4331026-0.47235-0.37500E-010.00000.4738327-0.28524-0.35625E-010.00000.2874628-0.34666E-01-0.33500E-010.00000.48208E-01290.11473-0.31125E-010.00000.11888300.00000.00000.00000.000031-0.42382E-01-0.53750E-020.00000.42721E-0132-0.13566-0.10500E-010.00000.1360733-0.22903-0.15375E-010.00000.2295534-0.27170-0.20000E-010.00000.2724335-0.21286-0.24375E-010.00000.2142536-0.13662E-02-0.583560.00000.5835737-0.10247E-02-1.40420.00001.4042*****POST1NODALDEGREEOFFREEDOMLISTING*****LOADSTEP=1SUBSTEP=1TIME=1.0000LOADCASE=0THEFOLLOWINGDEGREEOFFREEDOMRESULTSAREINTHEGLOBALCOORDINATESYSTEMNODEUXUYUZUSUM38-0.68311E-03-2.19870.00002.198739-0.34155E-03-2.75870.00002.7587400.22519E-13-2.95920.00002.9592410.34155E-03-2.75870.00002.7587420.68311E-03-2.19870.00002.1987430.10247E-02-1.40420.00001.4042440.13662E-02-0.583560.00000.58357MAXIMUMABSOLUTEVALUESNODE1019019VALUE0.47235-3.92330.00003.9233结构在平面内水平振动的前5阶自振频率及相应的振型有限元模型的建立和前面所示一致在此不赘述按以下方法求解1、定义分析类型2、选择模态分析类型和提取数3、设置子空间法4、求解完成5、查看结果前5阶自振频率如下所示:*****INDEXOFDATASETSONRESULTSFILE*****SETTIME/FREQLOADSTEPSUBSTEPCUMULATIVE10.55991E-0111120.1932012230.2561313340.3130714450.64091155各阶振型如下所示结构在振动荷载作用下的振动反应有限元模型的建立和前面所示一致在此不赘述按以下方法求解1、定义分析类型2)瞬态分析方法3)定义参数4)输入命令5)求解完成6)查看结果三、预应力混凝土简支梁1)进入preferences选中structural2)定义单元类型3)定义实常数4)定义材料属性5)建立模型6)创建荷载施力点7)选择单元属性8)网格大小控制9)划分单元10)求解分析类型分析控制施加约束施加荷载结点压缩求解完成11)查看结果变形图云图四、桁架的优化设计1、进入preferences选中structural2、进行参数的初始值设定3、定义单元类型4、实常数定义点击ok,完成。5、定义材料特性6、建立单元模型7、求解施加约束:在1、3结点施加全约束,如下图所示施加荷载(在结点2处)求解完成8、后处理记录每个单元的体积(ElementTableDefineTable)求体积和(ElementTableSumofEachItem)SUMALLTHEACTIVEENTRIESINTHEELEMENTTABLETABLELABELTOTALEVOL0.241421E+07提取体积并赋予变量VTOT(ParametersGetScalarData)定义密度和重量参数(ParametersScalarParameters)定义轴向应力单元表(ElementTableDefineTable)提取轴向应力并赋予变量SIG1、SIG2(ParametersGetScalarData)点击Apply计算轴向应力(ParametersScalarParameters)输入sig1=abs(sig1)sig2=(sig2)9、生成文件(FileWriteDBLogFile)10、执行优化分析过程(DesignOptAnalysisFileAssign)定义优化设计变量(DesignOptDesignVariable)定义优化状态变量(DesignOptStateVariable)定义目标函数(DesignOptObjective)设置优化计算总体控制选项指定一阶优化方法执行优化分析