这句话就是刚度矩阵KLL主元比过大,有可能是矩阵奇异,或者网格质量太差。一般要先检查是否存在约束不足的问题。或者你算的模型,存在弯曲刚度局部连接处相差太大的原因,要是这个原因,你可以试试下面方法:
在NXNastran程序文件中找到nast7.rcf(其路径为NASTRAN\conf\nast7.rcf),先去掉其只读的属性,然后用记事本打开,在里面加入一行PARAM,BAILOUT,-1在进行计算!
小知识:
“PARAM,BAILOUT,-1”这个参数一般算是debug用的,用这个参数算出来的结果常常是不可信的,
刚度矩阵奇异了,你可以检查下结构约束是否完全拿掉,另外也可以用AUTOSPC这个参数去拿掉些GridSingularity,如果刚度矩阵有mechanismsingularity,在f06中会打印出来一些表格,表明哪个自由度导致了这个奇异的发生,Nastran的一般做法是把刚度矩阵对角化处理,如果本身刚度矩阵奇异的话,那么理论上,对角化以后在对角线上会出现0值,但由于roundoff的误差,会是一个很小的数,把刚度矩阵原来对角线上项的值除以这个对角化以后在对角线上的值,会变得很大的一个数,这个就算是MAXRATIO,如果这个比值很大的话,就说明刚度矩阵有奇异项,有时调高参数MAXRATIO的值也能让求解继续下去,这个值一般默认是1.0E07,一般调到1.0E10也也算是安全的(比如在考虑接触的时候)
用函数加载器加载,刚试了下,我是用柱坐标系Y(θ)变化加载的力,也可以用其他参数time,x,y。
先定义函数,然后点file〉保存(随便写个英文函数名字),关闭函数加载器。然后在点utility menu>functions>read form file,选刚保存的函数文件。设定table parameter name为加载函数名字(如volt),坐标系也换成变量(y)所需要的坐标系。OK。再去solution中加载,选existing table。
PS:楼上那个加载不是正弦加载,而是每次加载值用计算出来的力的正弦值再次加载。比如,t=1,得出力加载上去,求解;t=2,得出力再加载上去,求解。属于运用正弦变力加载。而lz应该是加载正弦载荷(载荷成正弦分布)。
首先在input 界面里把荷载类型设置成动荷载,在calculations界面里,calculation type选项里选择dynamic analysis,在multipliers选项里点击红色波形曲线的图标后,输入正弦波的振幅和频率即可欢迎分享,转载请注明来源:内存溢出
评论列表(0条)