第二节求解器基础
- 格式:ppt
- 大小:1.40 MB
- 文档页数:18
结构力学求解器重要知识点看结构力学书了解:自振频率、几何构造(怎样判断多余约束数与自由度)。
1、输入数值时(如刚度等),可用科学记数法输入:如1.25e52、自己学习时,可打开“入门向导.inp”文件,其中包含两个例题及创建该例题的详细的操作步骤。
可按其注释说明进行同样的操作。
3、变量定义,如L=6,H=L/2,可方便后续的数值的输入4、删除命令后,可按ctrl+z撤消5、荷载条件中,均布荷载沿杆轴通用,沿水平则是均布荷载水平分布(不能在竖直杆上分布),沿竖直是均布荷载竖直分布(不能在水平杆上分布)。
6、0.4*0.6**3/12=(0.4*0.63)/127、自振频率参数,10,1,0.00000005?8、修改命令:调用对话框以修改当前光标所在的命令行;9、一下命令在编辑菜单中:复原:使用该命令,可以取消上一次所做的编辑,如果无编辑操作,则该命令无法使用;删除:从文档中删除选定的文本,但不将删除的文本放到剪贴板中;查找:使用该命令,可以在活动文档中搜索指定的文字;查找下一个:不需要打开“查找”对话框即可重复前一个搜索动作,若选择该命令前未使用“查找”命令,则该命令无法使用;全选:使用该命令,可以选定整个文档;10、以下命令在命令菜单中:尺寸线:打开尺寸线命令对话框,输入尺寸线定义相关的参数(可在“标注值”输入单位,如:2m);文本:可在观览器上输入文字或公式;全文翻译:将所有命令的关键词转换成选中的语言(简体中文、英文)。
11、观览器命令:单步显示:根据设计文档,逐行(命令行)显示当前问题的结构图形;停止单步显示:在单步显示按钮按下后,位于工具栏的最右方的该按钮将由灰变亮,可用于中段单步显示;连续显示:根据设计文档,一次性地将当前问题地整个结构图画出;暂停显示:使用该命令,可以暂时中断结构图形的显示;定制结点码、单元码:命令-定制编码-标注-定制结点码、单元码杆件内力图划分段数:可以确定绘制杆件内力图、变形图时的数据点数,数据点的采集首先保证取到图形的关键点。
1.变步长(Variable—Step)求解器可以选择的变步长求解器有:ode45,ode23,ode113,odel5s,ode23s和discret.缺省情况下,具有状态的系统用的是ode45;没有状态的系统用的是discrete.1)ode45基于显式Runge—Kutta(4,5)公式,Dormand—Prince对.它是—个单步求解器(solver)。
也就是说它在计算y(tn)时,仅仅利用前一步的计算结果y(tn-1).对于大多数问题.在第一次仿真时、可用ode45试一下.2)ode23是基于显式Runge—Kutta(2,3).Bogackt和Shampine对.对于宽误差容限和存在轻微刚性的系统、它比ode45更有效一些.ode23也是单步求解器.3)odell3是变阶Adams-Bashforth—Moulton PECE求解器.在误差容限比较严时,它比ode45更有效.odell3是一个多步求解器,即为了计算当前的结果y(tn),不仅要知道前一步结果y(tn-1),还要知道前几步的结果y(tn-2),y(tn-3),…;4)odel5s是基于数值微分公式(NDFs)的变阶求解器.它与后向微分公式BDFs(也叫Gear方法)有联系.但比它更有效.ode15s是一个多步求解器,如果认为一个问题是刚性的,或者在用ode45s时仿真失败或不够有效时,可以试试odel5s。
odel5s是基于一到五阶的NDF 公式的求解器.尽管公式的阶数越高结果越精确,但稳定性会差一些.如果模型是刚性的,并且要求有比较好的稳定性,应将最大的阶数减小到2.选择odel5s求解器时,对话框中会显示这一参数.可以用ode23求解器代替。
del5s,ode23是定步长、低阶求解器.5)ode23s是基于一个2阶改进的Rosenbrock公式.因为它是一个单步求解器,所以对于宽误差容限,它比odel5s更有效.对于一些用odel5s不是很有效的刚性问题,可以用它解决.6)ode23t是使用“自由”内插式梯形规则来实现的.如果问题是适度刚性,而且需要没有数字阻尼的结果,可采用该求解器.7)ode23tb是使用TR—BDF2来实现的,即基于隐式Runge—Kutta公式,其第一级是梯形规则步长和第二级是二阶反向微分公式.两级计算使用相同的迭代矩阵.与ode23s相似,对于宽误差容限,它比odtl5s更有效.8)discrete(变步长)是simulink在检测到模型中没有连续状态时所选择的一种求解器.====================================================================== ===================2.定步长(Flxed—Step)求解器可以选择的定步长求解器有:ode5,ode4,ode3,ode2,ode1和discrete.1)ode5是ode45的一个定步长版本,基于Dormand—Prince公式.2)ode4是RK4,基于四阶Runge—Kutta公式.3) ode3是ode23的定步长版本,基于Bogacki-Sbampine公式.4) ode2是Heun方法,也叫作改进Euler公式.5) odel是Euler方法.6) discrete(定步长)是不执行积分的定步长求解器.它适用于没有状态的模型,以及对过零点检测和误差控制不重要的模型.====================================================================== ===================3.诊断页(Diagnostics)可以通过选择Simulation Parameters对话框的Diagnostics标签来指明在仿真期间遇到一些事件或者条件时希望执行的动作.对于每一事件类型,可以选择是否需要提示消息,是警告消息还是错误消息.警告消息不会终止仿真,而错误消息则会中止仿真的运行.(1)一致性检查一致性检查是一个调试工具.用它可以验证Simulink的0DE求解器所做的某些假设.它的主要用途是确保s函数遵循Simulink内建模块所遵循的规则.因为一致性检查会导致性能的大幅度下阵(高达40%),所以一般应将它设为关的状态.使用一致性检查可以验证s函数,并有助于确定导致意外仿真结果的原因.为了执行高效的积分运算,Simulink保存一些时间步的结果,并提供给下一时间步使用.例如,某一时间步结束的导数通常可以放下一时间步开始时再使用.求解器利用这一点可以防止多余的导数运算.一致性检查的另一个目的是保证当模块被以一个给定的t(时间)值调用时.它产生一常量输出.这对于刚性求解器(ode23s和odel5s)非常重要,因为当计算Jacobi行列式时.模块的输出函数可能会被以相同的t值调用多次.如果选择了一致性检查,Simulink置新计算某些值,并将它们与保存在内存中的值进行比较,如果这些值有不相同的,将会产生一致性错误.Simulink比较下列量的计算值:1)输出;2)过零点3)导数;4)状态.(2)关闭过零点检测可以关闭一个仿真的过零点检测.对于一个有过零点的模型,关闭过零点检测会加快仿真的速度,但是可能影响仿真结果的精度.这一选项关闭那些本来就有过零点检测的模块的过零点检测.它不能关闭Hir crossing模块的过零点检测.(3)关闭优化I/O存储选择该选项,将导致Simulink为每个模块约I/()值分配单独的缓存,而不是重新利用援存.这样可以充分增加大模型仿真所需内存的数量.只有需要调试模型时才选择该选项.在下列情况下,应当关闭缓存再利用;1)调试一个C-MEX S-函数;2)使用浮点scope或display模块来察看调试模型中的信号.如果缓存再利用打开,并且试图使用浮点scope或display模块来显示缓存已被再利用的信号,将会打开一个错误对话框.(4)放松逻辑类型检验选择该选项,可使要求逻辑类型输入的模块接受双精度类型输入.这样可保证与Simulink 3版本之前的模型的兼容性.====================================================================== ===================4.提高仿真性能和精度仿值性能相精度由多种因素决定,包括模型的设计和仿真参数的选择.求解器使用它们的缺省参数值可以使大多数模型的仿真比较精确有效,然而,对于一些模型如果调整求解器相仿真参数将会产生更好的结果.而且,如果对模型的性能比较熟悉,并且将这些信息提供给求解器,得到的仿真效果将会提高。
结构力学求解器重要知识点结构力学求解器重要知识点看结构力学书了解:自振频率、几何构造(怎样判断多余约束数与自由度)。
1、输入数值时(如刚度等),可用科学记数法输入:如1.25e52、自己学习时,可打开“入门向导.inp”文件,其中包含两个例题及创建该例题的详细的操作步骤。
可按其注释说明进行同样的操作。
3、变量定义,如L=6,H=L/2,可方便后续的数值的输入4、删除命令后,可按ctrl+z撤消5、荷载条件中,均布荷载沿杆轴通用,沿水平则是均布荷载水平分布(不能在竖直杆上分布),沿竖直是均布荷载竖直分布(不能在水平杆上分布)。
6、0.4*0.6**3/12=(0.4*0.63)/127、自振频率参数,10,1,0.00000005?8、修改命令:调用对话框以修改当前光标所在的命令行;9、一下命令在编辑菜单中:复原:使用该命令,可以取消上一次所做的编辑,如果无编辑操作,则该命令无法使用;删除:从文档中删除选定的文本,但不将删除的文本放到剪贴板中;查找:使用该命令,可以在活动文档中搜索指定的文字;查找下一个:不需要打开“查找”对话框即可重复前一个搜索动作,若选择该命令前未使用“查找”命令,则该命令无法使用;全选:使用该命令,可以选定整个文档;10、以下命令在命令菜单中:尺寸线:打开尺寸线命令对话框,输入尺寸线定义相关的参数(可在“标注值”输入单位,如:2m);文本:可在观览器上输入文字或公式;全文翻译:将所有命令的关键词转换成选中的语言(简体中文、英文)。
11、观览器命令:单步显示:根据设计文档,逐行(命令行)显示当前问题的结构图形;停止单步显示:在单步显示按钮按下后,位于工具栏的最右方的该按钮将由灰变亮,可用于中段单步显示;连续显示:根据设计文档,一次性地将当前问题地整个结构图画出;暂停显示:使用该命令,可以暂时中断结构图形的显示;定制结点码、单元码:命令-定制编码-标注-定制结点码、单元码杆件内力图划分段数:可以确定绘制杆件内力图、变形图时的数据点数,数据点的采集首先保证取到图形的关键点。
求解器的使用FLUENT提供了三种不同的求解器Segregated,coupled implicit,coupled explicit(显式格式主要用于激波等波动解的捕捉问题)传统上,分离解法(segregated)主要用于不可压缩以及适度压缩的流动中。
相反,耦合算法是为高速可压流体设计的。
默认情况下,fluent使用分离求解器。
对于高速可压流体,与很强的体积力高度耦合的流动,或者是在非常精确的网格上求解流动情况,可以考虑使用耦合隐式算法代替。
对于需要使用耦合隐式算法(coupled implicit)的case,如果电脑没有足够的内存,可以使用分离解法(segregated)或者耦合显式算法(coupled explicit)代替。
(显示算法节约内存,但是需要更多的计算步数达到收敛。
)选择离散格式1.一阶迎风格式v.s. 二阶迎风格式当流动与网格匹配(校准)时,一阶迎风格式是可以接受的。
对于三角形和四面体网格,由于流动不会与网格匹配,通常使用二阶离散格式会得到更准确的结果。
对于四边形/六面体网格,使用二阶离散格式会取得更好的结果,尤其是复杂的流动情况。
对于大多数情况,可以在计算初始,使用二阶的离散格式。
然而在一些情况下,可以开始使用一阶的离散格式然后在一些计算之后转变为二阶格式。
例如,如果正在运行一个高马赫数的流动计算,这个的初始解与期望的解相差很大,最终,如果二阶离散格式很难收敛,应该尝试使用一阶离散格式。
2.Quick格式v.s. Upwind(Quick格式适用于网络结构,流动方向与网格一致,对于非结构网格推荐使用二阶迎风)对于在四边形或者六面体网格中的旋转或者回旋流,Quick离散格相比于二阶离散格式可以提供更准确的结果。
对于存在震动的可压缩流动(网格为四边形,六面体或者混合网格),推荐对所有的变量使用Quick离散格式,包括密度。
3.中心差分格式v.s. 迎风格式当使用LES湍流模型时,是可以使用中心差分格式的,并且只有当网格间距足够好,以至于局部的Peclet数的大小小于1时才可以使用。
之杨若古兰创作结构力学求解器(使用指南)结构力学求解器(SM Solver of Windows)是一个关于结构力学分析计算的计算机软件,其功能包含求解平面杆件结构(体系)的几何构成、静定和超静定结构的内力、位移,影响线、自在振动的自振频率和振型,和弹性波动等结构力学课程中所涉及的绝大部分成绩.对几何可变体系可作静态或动态显示机构模态;能绘制结构内力图和位移图;能静态或动态显示结构自在振动的各阶振型和弹性波动分析的失稳模态;能绘制结构的影响线图.该软件的版本为V1.5.清华大学土木系研制.高教出版社发行.一.运转环境Windows 98/NT. 8M内存. 2M硬盘空间.二.装机与运转将软件光盘置入光驱,在Windows环境下运转光盘上的SMsetup.exe,然后按提示操纵即可完成装机.装机完成后,桌面上将出现一个名为"求解器"的图标.双击桌面上的"求解器"图标,再单击软件的封面,即可使用该求解器.三.输入数椐先对结构的结点及单元进行编码,然后按以下诸项输入数椐:1.结点定义N,Nn,x,yNn---结点编码;x---结点的x 坐标;y---结点的y 坐标.结构全体坐标系为xoy,普通取结构左下支座结点为坐标原点(0,0).2.结点生成(即成批输入结点坐标)NGEN,Ngen,Nincr,N1,N2,N12incr,Dx,DYNgen---结点生成的次数;Nincr---每次生成的结点码增量;N1、N2---基础结点范围;N12incr---基础结点的编码增量;Dx,DY---生成结点的x ,y坐标增量.3.单元定义E,N1,N2[,DOF11,DOF12,DOF13,DOF21,DOF22,DOF23]N1,N2---单元两端的结点码;以下连接方式:1为连接,0为不连接;DOF11,DOF12---分别为单元在杆端1处的x、y方向自在度的连接方式,缺省值=1;DOF13---单元在杆端1处的转角方向自在度的连接方式,缺省值=0;DOF21,DOF22---分别为单元在杆端2处的x、y方向自在度的连接方式,缺省值=1;DOF23---单元在杆端2处的转角方向自在度的连接方式,缺省值=0.4.单元生成(即成批输入单元两端的连接方式)EGEN,Ngen,E1,E2,NincrEgen---生成次数;E1,E2---基础单元范围;Nincr---生成中单元两端点对应的结点码增量.5.支座束缚定义NSUPT,Sn,Stype,Sdir,[,Sdisx,Sdisy,SdisR]Sn---支座的结点码;Stype---支座类型码;Sdir---支座方向,以图示方向为零,绕结点逆时针扭转为正;Sdisx---x方向的支坐位移,缺省值=0;Sdisy---y方向的支坐位移,缺省值=0;SdisR---转角方向的支坐位移,缺省值=0.以上(1)~(6)为支座类型码.6.单元材料性质ECHAR,ElemStart,ElemEnd,EA,EI,mElemStart---单元起始码;ElemEnd---单元终止码;EA,EI---分别为单元的抗拉和抗弯刚度;m---单元的均布质量(kg/m).7.结点荷载NLOAD,Ln,Ltype,Lsize[,Ldir]Ln---荷载感化的结点码;Ltype---荷载类型;Ltype=1(-1),集中荷载,指向(叛变)结点;Ltype=2(-2),逆时(顺时)针方向的集中力矩;Lsize---荷载大小(kn,kn-m);Ldir---荷载方向(度),仅当Ltype=1或-1时入,缺省值=0 . 说明:竖向集中力,感化在结点上方时,取=-90 ,反之,取=90 ;水平集中力,感化在结点左方时,取=0 ,反之,取=180 .8.单元荷载ELOAD,Ln,Ltype,Lsize1[,Lpos1[Lpos2[,Ldir]]]Ln---荷载感化的单元码;Ltype---荷载类型;Ltype=1(-1),集中荷载,指向(叛变)单元;Ltype=2(2),逆时(顺时)针方向的集中力矩;Ltype=3(-3),均布荷载,指向(叛变)单元;Lsize1---荷载大小;Lpos1---荷载起点至单元杆端1的距离与单元杆长的比值,缺省值=0;Lpos2---荷载起点至单元杆端1的距离与单元杆长的比值,缺省值=1;(仅对均布荷载输入Lpos2)Ldir---荷载方向(度),仅当Ltype=1,3或-1,-3时输入,缺省值=0.(注:按局部坐标系定义,其角度以逆时针方向为正)9.频率计算参数FREQ,Nfreq,FreqStart,TolNfreq---欲求的频率数目;FreqStart---频率起始阶数;Tol---精度误差限,如0.0005.10.影响线参数IL,LoadDOF,En,pos,FdofLoadDOF---单位荷载的方向(全体坐标系):1为水平,2为竖直,3为转角;En---单元码;pos---单元上截面地位:距杆端1的距离与杆长的比值;Fdof---欲求影响线的内力自在度(局部坐标系),1为轴力,2为剪力,3为弯矩.说明:1.计算结构的内力和位移时,仅输入1(或及2),3(或及4),5,6,7,8项;2.当单元的抗拉刚度(EA)或抗弯刚度(EI)为无量大时,则分别填-1;3.当斜杆单元感化沿水平线的均布荷载时,需按合力相等的准绳,变换成沿杆轴线分布的均布荷载输入,荷载类型码仍为3(见例5).四.上机操纵步调1.双击桌面上的"求解器"图标,再单击"求解器"的封面进入使用形态;2.键入数椐文件名(如TITLE,XXXX),逐行输入数椐(也可用命令方式输入);3.将数椐文件存盘单击桌面上方的"文件",在文件菜单中点"保管"或"另存为",键入文件名,点"确定";4.再单击"文件",在文件菜单中点"退出";5.见提示"?此命令将结束本次SM Solver!"点"取消"或'确定",从头进入SM Solver;6.单击"文件",在文件菜单中点"打开";7.点所要运转的数椐文件名,并单击"确定";8.单击桌面下方的"观览器"图标,(桌面上显示结构计算简图的外形),并单击"最大化"按钮,将图形放大;9.单击桌面上方的"标注",在"标注"菜单中点所要显示的参数;(如无误,则进行下一步,若有误,则进行点窜)10.单击"观览器"图标,点桌面上方的"求解";11.在"求解"菜单中,点所要计算的内容(如内力计算、位移计算等),即可显示计算结果(如各杆杆端的内力或位移,对照结构的单元编号或结点编号浏览);12.单击所要显示的内力类型(轴力、剪力、弯矩)及显示对象(如"结构"或"单元");13.单击"观览器"图标,则显示出内力图或位移图;14.反复单击"观览器"图标,即可选定和显示分歧的内力图;15.逐层单击题目栏右侧的"关闭"按钮,当显示:"此命令将结束本次SM Solver的运转"或提示"结力求解器!Overflow"时,则点"确定"退出.五.计算例题例1 求图示刚架的内力.各杆的EA=3.12X10 KN,EI=4.16X10 KN-M.TITLE,AAA-1N,1,0,0N,2,4,0N,3,0,4N,4,4,4E,1,3,1,1,1,1,1,0E,3,4,1,1,0,1,1,1E,4,2,1,1,1,1,1,1NSUPT,1,6,0,0,0,0NSUPT,2,6,0,0,0,0ECHAR,1,3,3.12E+06,4.16E+04NLOAD,3,1,30,0ELOAD,1,3,20,0,1,90例2 求图示组合结构的内力.设各杆的EA = EI= 1.TITLE,AAA-2N,1,0,0 NSUPT,1,2,-90,0,0N,2,2,0 NSUPT,5,1,0,0N,3,4,0 ECHAR,1,9,1,1N,4,6,0 ELOAD,1,3,1,0,1,90N,5,8,0 ELOAD,2,3,1,0,1,90N,6,2,-2 ELOAD,3,3,1,0,1,90N,7,6,-2 ELOAD,4,3,1,0,1,90E,1,2,1,1,0,1,1,1E,2,3,1,1,1,1,1,0E,3,4,1,1,0,1,1,1E,4,5,1,1,1,1,1,0E,6,7,1,1,0,1,1,0E,6,2,1,1,0,1,1,0E,7,4,1,1,0,1,1,0E,6,1,1,1,0,1,1,0E,7,5,1,1,0,1,1,0例3.求图示桁架各杆的轴力. TITLE,,AAA-3N,1,1,0NGEN,4,1,1,1,1,1,0NGEN,1,5,1,5,1,0,1E,1,2,1,1,0,1,1,0EGEN,3,1,1,1EGEN,1,1,4,5E,6,1,1,1,0,1,1,0EGEN,4,9,9,1E,1,7,1,1,0,1,1,0E,1,8,1,1,0,1,1,0E,8,5,1,1,0,1,1,0E,5,9,1,1,0,1,1,0NSUPT,1,2,-90,0,0NSUPT,5,1,0,0NLOAD,8,1,1,-90NLOAD,9,-1,2,-90(注:此题系静定结构,其内力与材料性质有关,故可不输入ECHAR项)例4.求图示桁架的轴力.提示:支座束缚和结点荷载信息为NSUPT,1,2,-90,0,0NSUPT,2,1,0,0NLOAD,7,1,8,-90NLOAD,8,1,4,180NLOAD,5,1,4,180例5.求图示三铰刚架的内力.提示:支座束缚及斜杆单元的荷载信息为NSUPT,1,2,-90,0,0NSUPT,2,2,0,0,0ELOAD,2,3,9.48682596,0,1,71.565注:将沿水平线均布荷载(q)变换成沿杆轴线的均布荷载(q ) 即q =qcos =10X6 / 40 = 9.48682596例6.求图示刚架的内力.设EI=1.TITLE,AAA-6N,1,0,0N,2,4,0N,3,8,0N,4,4,-4E,1,2,1,1,1,1,1,1E,2,3,1,1,1,1,0,1E,4,2,1,1,0,1,1,1NSUPT,1,6,-90,0,0,0NSUPT,3,5,0,0,0NSUPT,4,4,90,0,0ECHAR,1,1,-1,1ECHAR,2,2,-1,2ECHAR,3,3,-1,1ELOAD,1,3,30,0,1,90ELOAD,2,1,50,0.5,90(注:取EA= ,填-1)例7.求图示梁的内力和位移.EI=5X10 KN-M .TITLE,AAA-7N,1,0,0N,2,6,0N,3,7.5,0E,1,2,1,1,0,1,1,1E,2,3,1,1,1,0,0,0NSUPT,1,3,0,0NSUPT,2,1,0,0ECHAR,1 2 -1,5E+04ELOAD,2,3,16,0,1,90ELOAD,2,1,20,1,90例8.求图示铰接排架的内力.EI = 1, EI = 6(设横梁的EI=1,柱子的EA= )TITLE,AAA-8N,1,0,0 NSUPT,1,6,0,0,0,0N,2,6,0 NSUPT,2,6,0,0,0,0N,3,16,0 NSUPT,3,6,0,0,0,0N,4,0,6 ECHAR,1,1,-1,1N,5,6,6 ECHAR,2,4,-1,6N,6,6,7 ECHAR,5,6,-1,1N,7,6,10 ECHAR,7,8,-1,1N,8,16,10 ELOAD,6,1,20,1/3,90N,9,16,7E,1,4,1,1,1,1,1,0E,2,5,1,1,1,1,1,1E,3,9,1,1,1,1,1,1E,5,6,1,1,1,1,1,1E,6,7,1,1,1,1,1,0E,9,8,1,1,1,1,1,0E,4,5,1,1,0,1,1,0E,7,8,1,1,0,1,1,0例9.计算图示两层刚架的自振频率和主振型横梁的均布质量m = m = 15X10 kg/m 柱子的抗弯刚度EI =1X10 kn.m设EA =TITLE,AAA-9N,1,0,0N,2,4,0N,3,0,3N,4,4,3N,5,0,6N,6,4,6E,1,3,1,1,1,1,1,1E,3,5,1,1,1,1,1,1E,2,4,1,1,1,1,1,1E,4,6,1,1,1,1,1,1E,3,4,1,1,1,1,1,1E,5,6,1,1,1,1,1,1NSUPT,1,6,0,0,0,0NSUPT,2,6,0,0,0,0ECHAR,1,4,-1,1E+08,1E-08ECHAR,5,6,-1,-1,1.5E+04FREQ,2,1,0.0005(注:柱子的质量不克不及填0,可填一个很小的数,如10 ) 例10.对图示两跨四层框架结构,分别计算竖向荷载和水平荷载感化下的内力.各杆的EA、EI值见下表:框架梁柱计算参数表截面弹性模量惯性矩EA EI构件A=bXh(m ) E(kn/m ) I(m ) (kn) (kn-m )底层0.25xX0.5 3X10 0.521X10 0.375X10 1.563X10梁其它层0.25X0.5 2.8X10 0.521X10 0.350X10 1.459X10 底边柱0.4X0.4 3X10 0.213X10 0.480X10 0.639X10 层中柱0.45X0.45 3X10 0.342X10 0.608X10 1.026X10柱其边柱0.4X0.4 2.8X10 0.213X10 0.448X10 0.596X10它层中柱0.45X0.45 2.8X10 0.342X10 0.567X10 0.958X10TITLE,AAA-10N,1,0,0 水平荷载感化N,2,5.0,0 NLOAD,4,1,8.05,0N,3,10.0,0 NLOAD,7,1,11.17,0NGEN,1,3,1,3,1,0,4.5 NLOAD,10,1,15.20,0NGEN,3,3,4,6,1,0,3 NLOAD,13,1,19.10,0E,1,4,1,1,1,1,1,1 竖向荷载感化EGEN,2,1,1,1 ELOAD,13,3,19.30,0,1,9EGEN,3,1,3,3 ELOAD,14,3,19.30,0,1,90E,4,5,1,1,1,1,1,1 ELOAD,15,3,19.30,0,1,90E,5,6,1,1,1,1,1,1 ELOAD,16,3,19.30,0,1,90E,7,8,1,1,1,1,1,1 ELOAD,17,3,19.30,0,1,90EGEN,2,15,15,3 ELOAD,18,3,19.30,0,1,90E,8,9,1,1,1,1,1,1 ELOAD,19,3,19.50,0,1,90EGEN,2,18,18,3 ELOAD,20,3,19.50,0,1,90NSUPT,1,6,0,0,0,0 NLOAD,4,1,53.79,-90NSUPT,2,6,0,0,0,0 NLOAD,7,1,53.79,-90NSUPT,3,6,0,0,0,0 NLOAD,10,1,53.79,-90ECHAR,1,1,4.8E+06,6.39E+04 NLOAD,6,1,53.79,-90 ECHAR,2,2,6.08E+06,10.26E+04 NLOAD,9,1,53.79,-90 ECHAR,3,3,4.8E+06,6.39E+04 NLOAD,12,1,53.79,-90 ECHAR,4,4,4.48E+06,5.96E+04 NLOAD,5,1,71.97,-90 ECHAR,6,6,4.48E+06,5.96E+04 NLOAD,8,1,71.97,-90ECHAR,7,7,4.48E+06,5.96E+04 NLOAD,11,1,71.97,-90 ECHAR,9,9,4,48E+06,5.96E+04 NLOAD,13,1,44.08,-90 ECHAR,10,10,4.48E+06,5.96E+04 NLOAD,15,1,44.08,-90 ECHAR,12,12,4.48E+06,5.96E+04 NLOAD,14,1,50.86,-90 ECHAR,5,5,5.67E+06,9.58E+04ECHAR,8,8,5.67E+06,9.58E+04ECHAR,11,11,5.67E+06,9.58E+04ECHAR,13,14,3.75E+06,15.63E+04ECHAR,15,20,3.50E+06,14.59E+04。
求解器使用手册该求解器(HTSolver)配普通高等教育“十一五”国家级规划教材赵镇南编《传热学》(第二版)使用。
版权属于作者所有。
这是一套完全自主研制的传热学习题解题软件,在国内是第一个,也是唯一的一个。
它具有强大的计算功能和与传热学教材紧密配合的众多辅助功能,其用法符合一般的思维习惯,因而简单易学。
尤其在处理比较复杂的复合传热问题,必须做迭代计算的问题(包括正常问题的逆向求解),以及本教材习题中大量需要进行变参数分析或参数优选的问题时,它的优点显得格外突出。
它的基本功能和特点可以归纳为:⏹建立了很多基本类型问题的数学模型-计算模板,可供直接套用。
⏹用户可以自主扩展、增添新的模型,从而拓宽软件的应用范围。
⏹对迭代计算的问题,既可以手动迭代(即看到迭代过程中每一步得出的数值),也可以全自动迭代,一键得出迭代结果。
⏹对肋效率、角系数、温差修正系数和黑体辐射函数等常用的图表,均实现了数字化,免除了查图表的麻烦,同时自变量取值范围加宽,计算精度也大幅提高。
⏹提供常用介质的物性子程序,免除了手工插值计算的麻烦。
⏹具备传热学常用数学辅助计算工具,如求解超越方程和高次方程,线性代数方程组,数值积分等。
HTSolver的主界面如下图所示,所有计算工作均围绕计算表进行。
计算表本质上是一个行数可以任意扩展的计算程序。
一般按照与手工计算相同或相近的思路逐步写出计算过程即可。
对于需要迭代求解的问题,可以用人工方式一步步地进行,也可以设置一个循环回路,迭代过程将自动进行。
迭代初值原则上没有限制,不违反物理规律即可。
计算表的构造设计如下:1.左第1栏为行号(编号)。
2.第2栏为说明文字,对计算过程无实质性影响,只为本人或者其他人日后容易看明白,可以省略。
3.第3栏为“目标量”,即对每个变量(常数或计算式)赋予一个名字。
注意该变量名是区分大小写的,即认为大写A与小写a不是同一个量。
4.第4栏“计算式”中,写入给定常量或者所采用的计算公式。
求解器介绍求解器Krylov、Sparse、CG_ILU、Mumps之间的区别如下:Krylov 求解器:利用阻抗Z(ω) 与频率f(ω)的结构形式Z(ω) :二阶函数f(ω) = f0+ωf1:线性函数直接求解4个频率,Krylov 子空间重构,显著减少CPU计算时间,对内存需求不变,精度没有影响,对模态提取与多孔介质不适用。
Krylov求解器需要与子求解器MUMPS或SPARSE连用,适用多频率点计算。
例如计算范围10-1000Hz,分析步1Hz和2Hz时时间几乎一致,但用其它求解器时间会存在倍数关系。
ACTRAN_12.2_Users_Guide_Vol1(6.2.7) 局限:不能用于DIFFUSE_FIELD和TURBULENT分析;不能用于管道模态声源分析;不能用于来自外部求解器载荷如BC_MESH和Aero-acoustic sources;SPARSE and MUMPS求解器:直接求解,快速,高效,利用稀疏矩阵的特点,大规模问题的求解。
SPARSE是Actran默认求解器,但MUMPS比SPARSE性能好大多数情况下取代它。
MUMPS (Multifrontal Massively Parallel Solver),需要内存少,MUMPS适用所有应用场合,特别适用对集群并行求解很有优势。
The RHS_BLOCK_SIZE option is a common option for all the direct solvers like SPARSE and MUMPS (not supported with iterative solvers like CG_ILU and KRYLOV), and allows to handle multiple right-hand-sides by blocks during the backtransformation step. This allows to reduce the memory used during the factorization and the back-transformation when computing a high number of loadcases.The consequence is an increase of the time required by the back-transformation.The default value of this parameter is 0, which corresponds to no partitioning of the right hand side.The MEMORY_RELAXATION corresponds to a percentageincrease in the estimated working space (default value 0.2). Decreasing this parameter allows consuming less memory, but a minimal overestimation should always be provided. This parameter should be increased when extra memory usage is encountered (for instance due to additional numerical pivoting).The ORDERING parameters selects the ordering method used inside the solver. Four ordering methods are available: CG_ILU求解器:迭代求解器,使用频率不多,优势对内存要求更低,计算时间也少,但是迭代不容易收敛。