ANSYS程序变截面梁与耦合自由度的应用CAE产品创新数字化(PLM)619.doc

ansys程序变截面梁与耦合自由度的应用_cae_产品创新数字化(plm)一、前言变截面梁，即两端截面不同的梁。在工程中通常将梁受力吃紧的一端采用较大的截面，另一端采用较小的截面，以实现等强度设计，节省原材料。这种变截面梁生活中普遍存在，远至古代娶亲用的花轿轿杆、农村大车车辕，直至当今建筑钢结构中各种采光大蓬横梁（包括弯曲横梁）、各种连杆摇臂结构等。特别是悬臂梁几乎处处可见，就连输电用的水泥电杆也做成根粗尖细的。由于这种变截面梁在工程中普遍存在，在ansys程序结构分析中，专门设置了变截面梁单元的功能，使用起来特别方便。具体操作过程在下面的实例分析中详细介绍。自由度耦合即构件连接处两个节点的自由度（包括移动自由度和转动自由度）变化是一致的，主节点如何变化，从节点随着同样变化。自由度耦合在静力分析时常用在连接件上，特别用在具有转动的连接件上。例如：汽车挂挡手柄连接端的球铰，各种销钉与耳环的连接，各种转盘与转轴的连接等。自由度耦合的概念与自由度释放的概念正好相反。在ansys程序中没有自由度释放功能，只有自由度耦合功能，但用自由度耦合功能完全可以达到自由度释放的目的。这种功能具体操作过程，也在下面的实例分析中详细介绍。二、雨蓬计算分析该雨蓬结构是由11根变截面工字钢梁、1根等截面工字钢梁、3根圆管钢梁、1根槽钢梁和6根拉杆钢组成。其结构如下：图1雨蓬结构1原始数据（1）坐标数据关键点号xyz100020063003045000 （第一拉杆上端点）40-10000 （变截面梁z轴方向）5010000 （第一排工字钢z轴方向）60-10006300（第六排槽钢z轴方向）（2）材料数据主钢梁变截面工字钢：大端30015086小端15015086第一横梁工字钢：15010086第二、第三、第四横梁圆管钢：1024 ri=47 r0=51第五横梁槽钢：1608086拉杆：1024a=1231.5 =0.0001(3) 载荷由玻璃均布载荷计算而得:中间节点: f=5560n边节点: f=2780n角节点: f=1690n(4) 单元梁单元: beam188杆单元: link8(5) 边界条件主梁固定于墙上:位移和转角全约束；拉杆与主梁连接: 位移耦合,转角自由；拉杆与墙连接: 位移约束, 转角自由。2 计算方法有限元几何非线性方法3操作命令（1）前处理a. 选择单元类型选择单元类型菜单路径：【preprocessor】【element type】【add/edit/delete】执行选择单元类型菜单路径命令，弹出定义单元类型对话框，单击add按钮，弹出单元类型库对话框，然后选择“beam 2node 188”梁单元，单击ok按钮。用同样方法选择“link spar 8”杆单元，单击ok按钮。b. 设置单元实常数设置单元实常数菜单路径：【preprocessor】【real constants】【add/edit/delete】执行设置单元实常数菜单路径命令，弹出对话框后，单击add按钮，选择“type 2 link8”单元，并单击ok按钮，弹出定义实常数对话框，在“id”中输入“1”，在“area”中输入“1231.5”，单击ok按钮。c. 设置单元截面数据变截面工字钢梁截面设置：大端截面设置：执行菜单路径命令【preprocessor】【sections】【beam】【common sections】，弹出定义截面特性对话框后，在“id”中设置梁单元截面代号“1”，在“name”中设置梁单元截面名称“m1”，在sub-type中选择工字钢截面图形。在“w1”中输入“150”，在“w2”中输入“150”，在“w3”中输入“300”，在“t1”中输入“8”，在“t2”中输入“8”，在“t3”中输入“6”，单击ok按钮。小端截面设置：执行菜单路径命令【preprocessor】【sections】【beam】【common sections】，弹出定义截面特性对话框后，在“id”中设置梁单元截面代号“2”，在“name”中设置梁单元截面名称“m2”，在sub-type中选择工字钢截面图形。在“w1”中输入“150”，在“w2”中输入“150”，在“w3”中输入“150”，在“t1”中输入“8”，在“t2”中输入“8”，在“t3”中输入“6”，单击ok按钮。尖劈梁截面设置：执行菜单路径命令【preprocessor】【sections】【beam】【taper section】【by xyz location】，弹出定义变截面特性对话框后，在“id”中设置变截面梁单元截面代号“3”，在“name”中设置变截面梁单元截面名称“m3”。在“beginning section id”中选择变截面梁起始截面代号“1 m1”，在“xyz location of beginning sect”中输入变截面梁起始截面坐标“0，0，0”；在“ending section id”中选择变截面梁结束截面代号“2 m2”，在“xyz location of ending section”中输入变截面梁结束截面坐标“0，0，6300”，单击ok按钮。普通工字钢梁截面设置：执行菜单路径命令【preprocessor】【sections】【beam】【common sections】，弹出定义截面特性对话框后，在“id”中设置梁单元截面代号“4”，在“name”中设置梁单元截面名称“m4”，在sub-type中选择工字钢截面图形。在“w1”中输入“100”，在“w2”中输入“100”，在“w3”中输入“150”，在“t1”中输入“8”，在“t2”中输入“8”，在“t3”中输入“6”，单击ok按钮。圆管梁截面设置：设置梁单元截面数据菜单路径：【preprocessor】【sections】【beam】【common】【sections】执行设置梁单元截面数据菜单路径命令，弹出对话框后，在“id”中设置梁单元截面代号“5”，在“sub-type”中选择圆管钢截面图形。在“ri”中输入“47”，在“r0”中输入“51”，单击ok按钮。槽钢梁截面设置：执行菜单路径命令【preprocessor】【sections】【beam】【common sections】，弹出定义截面特性对话框后，在“id”中设置梁单元截面代号“6”，在sub-type中选择槽钢截面图形。在“w1”中输入“80”，在“w2”中输入“80”，在“w3”中输入“160”，在“t1”中输入“8”，在“t2”中输入“8”，在“t3”中输入“6”，单击ok按钮。d. 设置材料性能数据设置材料性能数据菜单路径：【preprocessor】【material props】【material models】执行设置材料性能数据菜单路径命令，弹出定义材料性能数据对话框，选择右侧材料有关变量依次双击“structural”“linaer”“elastic”“isotropic” (也可以依次双击“favorites”“linaer static”“linaer isotropic”) 后，在出现的对话框中，输入弹性模量“ex=2.06e5”和泊松比“prxy=0.3”数值，输入密度数值“7.85e-9”，单击ok按钮，并退出对话框。e. 建立几何模型（a）生成关键点菜单路径：【preprocessor】【modeling】【create】【keypoints】【on working plane】执行生成关键点菜单路径命令，在弹出的工作平面对话框中，依照上述坐标次序，输入各关键点坐标值：“x，y，z”，最后单击工作平面对话框中的ok按钮。轴侧显示菜单路径：【utility menu】【plotctrls】【pan zoom rotate】【iso】(b) 生成直线菜单路径：【preprocessor】【modeling】【create】【lines】【lines】【straight line】执行生成直线菜单路径命令，弹出生成直线对话框，采用鼠标在图形窗口依次点击前、后两个关键点，生成一条直线，单击对话框中的ok按钮。(c) 将直线分成21等份菜单路径：【preprocessor】【operate】【booleans】【divide】【line into n lns】，弹出对话框后，选择直线，单击apply按钮，输入划分数“21”，单击ok按钮。（d）直线相加菜单路径：【preprocessor】【modeling】【operate】【booleans】【add】【lines】执行直线相加菜单路径命令，弹出对话框后，选择除第一段直线外的5段直线，单击apply按钮，接下来再次选择5段直线，单击apply按钮，直到共生成5段直线为止。(e) 5段直线赋置属性菜单路径：【preprocessor】【meshing】【meshing attribute】【picked lines】执行赋属性菜单路径命令，弹出对话框后，单击对话框中的pick all按钮，弹出赋属性对话框，选择属性：材料号“mat”为“1”，实常数号“real”为默认值，单元类型号“type”为“1 beam188”，截面号“sect”为“3”，并激活“pick orientation keypoints”，单击apply按钮，用鼠标在主梁下选择确定梁单元z轴方向的“k”点，单击ok按钮。(f) 5段直线设置划分尺寸菜单路径：【preprocessor】【meshing】【size ctrls】【manualsize】【lines】【picked lines】执行设置划分尺寸菜单路径命令，采用鼠标在图形窗口拾取全部直线，单击apply按钮，在弹出的对话框“ndiv”中输入尺寸数值“1”，单击ok按钮。(g) 5段直线划分网格采用分网工具，其菜单路径：【preprocessor】【meshing】【meshtool】执行分网工具菜单路径命令，弹出分网工具对话框。先在分网工具窗口中确定网格实体：“lines”，单击mesh按钮后，单击pick all按钮。单击工具栏中的save-db按钮，保存实体网格文件。(h) 复制直线菜单路径：【preprocessor】【modeling】【copy】【lines】执行复制直线菜单路径命令，弹出对话框后，单击pick all按钮，在弹出对话框“number of copies-”中输入“11”，在“dx”中输入“2250”，单击ok按钮(目的在于生成11根变截面梁)。(i) 生成横线菜单路径：【preprocessor】【modeling】【create】【lines】【lines】【straight line】执行生成直线菜单路径命令，弹出生成直线对话框，采用鼠标在图形窗口依次点击前、后两个关键点，生成所有横直线，单击对话框中的ok按钮。（j）给各横直线赋置属性菜单路径：【preprocessor】【meshing】【meshing attribute】【picked lines】执行赋属性菜单路径命令，弹出对话框后，采用鼠标在图形窗口拾取第一横直线，单击apply按钮，弹出赋属性对话框，选择属性：材料号“mat”为“1”，实常数号“real”为默认值，单元类型号“type”为“1 beam188”，截面号“sect”为“4”，并激活“pick orientation keypoints”，单击apply按钮，用鼠标在第一横直线起点上选择确定梁单元z轴方向的“k”点，单击ok按钮。第二、第三、第四横直线为圆管钢，截面号为“5”，梁单元z轴方向的“k”点选择第五横直线上任意一点即可；第五横直线为槽钢，截面号为“6”，梁单元z轴方向的“k”点选择第五横直线起点下一点，单击ok按钮。(k) 各横直线划分尺寸设置为“1”、网格划分同前。(l) 合并菜单路径：【preprocessor】【numbering ctrls】【merge items】执行合并菜单路径命令，弹出对话框后，在“label”中选择“all”，单击ok按钮。(m) 复制关键点菜单路径：【preprocessor】【modeling】【copy】【keypoints】执行复制关键点菜单路径命令，弹出对话框后，用鼠标在图形窗口选择拉杆上端点，单击apply按钮，在弹出对话“number of copies”中输入“6”，在“dx”中输入“4500”，单击ok按钮。执行复制关键点菜单路径命令，弹出对话框后，用鼠标在图形窗口选择拉杆下端点之间的关键点，复制距离为“dx=2250”，复制后使各拉杆下端点位置形成两个关键点号。(n) 生成直线菜单路径：【preprocessor】【modeling】【create】【lines】【lines】【straight line】执行生成直线菜单路径命令，弹出对话框后，用鼠标在图形窗口选择拉杆上、下端点，生成拉杆线。注意：在选择下端点时，要选择后生成的关键点号。(o ) 拉杆线的属性赋置、划分尺寸设置和网格划分同前。( p) 自由度耦合设置该题6根拉杆下端与变截面梁的连接为铰链连接，即每个连接节点3个方向的转动自由度应释放为自由。在ansys程序中，自由度耦合的概念与其他程序中自由度释放的概念正相反。对于该题而言，拉杆下端与变截面梁铰链连接，应设置连接点处为两个节点，一个节点为拉杆的，另一个节点为变截面梁的，并使他们3个方向的位移自由度各自耦合，3个方向的转动自由度不作任何处理。自由度耦合设置菜单路径：【preprocessor】【conpling/ceqn】【coupling dofs】执行自由度耦合设置菜单路径命令，弹出对话框后，在图形窗口拾取耦合节点，单击ok按钮，再次在该位置拾取耦合点，单击apply按钮，弹出自由度耦合对话框，在“nset”中输入耦合集号“”，在“lab”中选择耦合自由度方向“”，单击apply按钮，直到完成所有6对节点位移自由度的耦合（6对节点共建18个耦合集），单击ok按钮。单击工具栏中的save-db按钮，保存实体网格文件。（2）加载求解a. 施加位移约束施加位移约束菜单路径：【solution】【define loads】【apply】【structural】【displacement】【on nodes】执行施加位移约束菜单路径命令，弹出施加位移约束选择节点对话框，采用鼠标在图形窗口拾取主梁与墙连接点，单击apply按钮，弹出施加位移约束选择分量窗口，选择节点位移分量：“ux=uy=uz=0，rotx= roty=rotz=0”，单击apply按钮。再拾取拉杆与墙连接点，选择节点位移分量：“ux=uy=uz=0”，单击ok按钮。b. 施加载荷选择节点菜单路径【solution】【define loads】【apply】【structural】【force/moment】【on nodes】执行选择节点菜单路径命令，弹出对话框后，在图形窗口拾取雨蓬中间各节点，单击apply按钮，弹出的对话框后选择“fy”，在value框中输入“-5560”，单击apply按钮；拾取雨蓬边缘各节点，单击apply按钮，弹出的对话框后选择“fy”，在value框中输入“-2780”，单击apply按钮；拾取雨蓬各角节点，单击apply按钮，弹出对话框后选择“fy”，在value框中输入“-1690”，单击ok按钮。c. 施加重力加速度（不考虑结构的重量可以不加）选择菜单路径：【solution】【define loads】【apply】【structural】【inertia】【gravity】【global】，弹出拾取对话框后, 在acelx中输入“”（水平方向），在acely中输入“”（垂直方向），单击ok按钮。d. 选择求解方法选择菜单路径：【main menu】【solution】【analysis type】【soln control】，弹出对话框后，在“analysis options”中选择“large displacemant static”，单击ok按钮。e. 用当前载荷步解题选择菜单路径：【main menu】【solution】【current ls】，弹出状态窗口和求解窗口，退出状态窗口，单击求解窗口中的ok按钮，进行解题。当出现此句英语：“solution is done !”则结果正确，单击close按钮。（3）后处理a. 显示变形选择菜单路径：【general postproc】【read results】【first set】【plot results】【contour plot】【nodal solu】，在弹出的对话框中，单击“dof solution”项，选择变形分量（例如：ycomponent of displacement），单击ok按钮，出现带标尺的变形图。b. 显示当量应力选择菜单路径：【general postproc】【read results】【first set】【plot results】【nodal solu】，在弹出的对话框中，单击“stress”项，选择“von mises stress”项，单击ok按钮。再选择菜单路径：【utility menu】【plotctrls】【style】【size and shape】，在弹出的对话框中，激活“display of elemet”项，单击ok按钮，出现带标尺的应力图。c. 显示弯曲应力定义单元表格菜单路径：【general postproc】【element table】【define table】执行定义单元表格菜单路径命令，弹出对话框以后单击add按钮，在滚动框中选择“by sequence num”，在右框中选择“smisc”，在右下框中“smisc，”后面输入应力的序列号“34”（beam4梁单元弯曲应力序列号为“ls,4”，各种梁单元应力的序列号是不同的，见在线帮助），单击ok按钮。绘制表格（或列表）菜单路径：【general postproc】【element table】【plot table】（【list elem table】）执行绘制表格（或列表）菜单路径命令，弹出对话框以后，单击ok按钮，绘制表内应力图形（或列出表格数据）。d. 支反力显示列表方式显示支反力菜单路径：【general postproc】【list results】【reaction solu】执行列表方式显示支反力菜单路径命令，弹出用列表方式显示支反力窗口，选择“all items”项，单击ok按钮。e. 将图形存入文件首先选择菜单路径：【utility menu】【plotctris】【style】【backgroud】，单击“display picture backgruod”（改变背景颜色，如果背景颜色已经改变则不要执行此步）。再选择菜单路径：【utility menu】【plotctris】【hard copy】【to file】，弹出的对话框后，在“save to：”中定义用户文件名，单击ok按钮。f. 输出结构总体积和总重量菜单路径【general postproc】【element table】【define table】【general postproc】【element table】【sum of each item】【utility menu】【paramenters】【get scalar data】【utility menu】【paramenters】【scalar paramenter】执行菜单路径命令1，弹出对话框以后，单击add, 弹出单元表对话框，在“lab”中输入用户定义名“evol”，在左下列选择“geometry”，在右列选择“elem volume volu”，单击ok按钮。执行菜单路径命令2，弹出对话框以后，单击ok按钮，得到单元总体积值。执行菜单路径命令3，弹出对话框以后，在左列选择“results data”，在右列选择“elem table sums”，单击ok按钮，弹出对话框以后输入用户定义名“vtot”，单击ok按钮。执行菜单路径命令4，弹出对话框以后，在“items”中选择“vtot=”，再在“selection”中乘以材料比重数值“例如：7.85e-6”（vtot=*7.85e-6），单击accept按钮，在“items”中得到单元总重量值（单位：vtot=千克重）。该雨蓬结构需用钢材吨。4 计算结果（1）位移分布云图：图2 雨蓬变形分布云图（2）应力分布云图：图3 雨蓬当量应力分布云图图4 雨蓬的梁弯曲应力分布云图5 结果分析及结论由图2图4可知，在当前载荷的作用下，雨蓬的最大变形在前缘为5.2mm；最大当量应力在根部为37.9mpa（最大弯曲应力37.2 mpa）。应力和变形都很小，可见用钢材浪费严重，结构设计很有改进的余地。三、总结1雨蓬计算分析例题前处理操作顺序（1）生成关键点（6个）；（2）生成一条直线；（3）分直线为21等份（具体尺寸数据需要）；（4）除第一段线外，每5段线相加，共形成5段线；（5）5段线赋置属性（注意z轴方向）；（6）设置划分尺寸为1；（7）划分单元；（8）储存单元；（9）复制5段线为11条直线（225011）；（10）生成所有横线；（11）横线分别赋置属性（注意z轴方向）；（12）所有横线设置划分尺寸为1；（13）所有横线划分单元；（14）所有合并（在此步一定要合并）；（15）储存所有单元；（16）复制拉杆上端节点（45006）；（17）利用拉杆下端点之间的节点复制拉杆下端节点，形成一个位置两个节点号；（18）利用新生成的拉杆下端节点和拉杆上端节点生成拉杆线；（19）拉杆线赋置属性；（20）拉杆线设置划分尺寸为1；（21）拉杆线划分单元；（22）储存所有单元；（23）从第一个拉杆下端节点位置开始，做每对节点的位移耦合，每做一个建立一个耦合集，共建立6对节点的18个位移耦合集。（24）施加其余边界条件和集中载荷；（25）使用几何非线性进行计算。2雨蓬分析注意下面几点（1）形成变截面工字钢梁时，与建立一般梁一样，建立两个梁截面，再建立尖劈梁。不要利用偏置功能，目的在于使尖劈梁轴水平放置。（2）beam188梁单元需要定义z轴方向，各梁放置位置不同，z轴的方向也不同。（3）拉杆预应变设置为0.0001，即给一个比较小的预紧力就可以了。（4）拉杆与主梁连接为铰链形式，采用位移自由度耦合，即每两个节点的ux、uy和uz耦合，转角自由。（5）拉杆与墙连接也为铰链形式，采用位移自由度约束，即ux=uy=uz=0，转角自由。（6）当量应力为合成应力，永远是正值；弯曲应力是有方向的，不过应取绝对值进行分析。（7）显示弯曲应力