多场耦合问题的程序生成.pptVIP

应力场描述方程已知变形和温度分布后，代入几何方程求应变，再通过本构方程求出热应力几何方程本构方程第29页，共42页，星期日，2025年，2月5日dispsx,sy,sxycoorx,ycoefu,v,tnshap%1%2gaus%3mass%11.0matepepvalfafxfy1.0e9;0.3;1.0e-5;0.0;0.0;$c6fact=pe/(1.+pv)/(1.-pv)vectfssxsysxyvectfffsxfsyfsxyvectemexeyvectfmfsxfsyvectfefttfttmatrde221pvpv1应力场ETC.VDE文件第30页，共42页，星期日，2025年，2月5日第1页，共42页，星期日，2025年，2月5日学习目的1.掌握GCN文件的填写格式；2.掌握用GCN功能生成程序的一般过程；3.了解非线性问题线性化的过程。第2页，共42页，星期日，2025年，2月5日GCN文件的结构aNFE文件名耦合场名bNFE文件名耦合场名cNFE文件名耦合场名............空一行命令流第3页，共42页，星期日，2025年，2月5日库中NFE文件命名基本原则NFE算法文件是FEPG系统用来生成单元计算程序(即E元件程序)和后处理计算程序(即U元件程序)描叙算法，系统提供大量供选择的算法。算法是与方程的类型有关的（方程类型可参看关于《PDE文件填写》的powerpoint文档），库中NFE文件命名基本原则：ell——椭圆par——抛物wave——波动n*——非线性第4页，共42页，星期日，2025年，2月5日ELL.NFE求解线性椭圆型方程的算法程序。PAR.NFE求解线性抛物型方程的算法程序，时间离散采用Crank-Nicolson格式。PARB.NFE求解线性抛物型方程的算法程序，时间离散采用向后差分格式。WAVE.NFE求解线性双曲型方程的算法程序，时间离散采用Wave速度法格式NEWMARK.NFE求解线性双曲型方程的算法程序，时间离散采用Newmark格式。NELL.NFE求解非线性椭圆型方程的算法程序。NPAR.NFE求解非线性抛物型方程的算法程序,时间离散采用Crank-Nicolson格式。第5页，共42页，星期日，2025年，2月5日NPARB.NFE求解非线性抛物型方程的算法程序，时间离散采用向后差分格式。NWAVE.NFE求解非线性双曲型方程的算法程序，时间离散采用Wave速度法格式。NNW.NFE求解非线性双曲型方程的算法程序，时间离散采用Newmark格式。STR.NFE已知位移场求应力场的算法程序，采用最小二乘法。也是求解已知场量的梯度的算法程序。第6页，共42页，星期日，2025年，2月5日命令流命令流段给出耦合场的计算流程。该段由两种命令组成，一种是DOS命令；另一种是本系统的专有命令。较常用的专有命令可参看《FEPG中级教程》关于GCN文件填写的部分。第7页，共42页，星期日，2025年，2月5日线性稳态问题命令流段startcsolvcSTART元件程序E元件程序SOLV求解器U元件程序命令流段线性稳态问题流程第8页，共42页，星期日，2025年，2月5日线性动态问题命令流段startcifexiststopdelstop:1bftsolvcifnotexiststopgoto1命令流段线性动态问题流程START元件程序BFT元件程序E元件程序SOLV求解器U元件程序时间计算是否结束否是第9页，共42页，星期日，2025年，2月5日非线性稳态问题命令流段startcifexistenddelend:1solvcifnotexistendgoto1命令流段非线性稳态问题流程START元件程序E元件程序SOLV求解器U元件程序迭代是否结束否第10页，共42页，星期日，2025年，2月5日非线性动态问题命令流段startcifexiststopdelstop:1bftifexistenddelend:2solvcifnotexistendgoto2ifnotexiststopgoto1