C++编程有限元公式

适合于C++有限元编程初学者，利用C++编辑有限元力学公式，作用初学的参考，可以借鉴。

# include “stdio.h“
# define U（x） （（unsigned char）（x））
# define NLSTATE yyprevious=YYNEWLINE
# define BEGIN yybgin = yysvec + 1 +
# define INITIAL 0
# define YYLERR yysvec
# define YYSTATE （yyestate-yysvec-1）
# define YYOPTIM 1
# define YYLMAX 2048
# define output（c） （void）putc（cyyout）
#if defined（__cplusplus） || defined（__STDC__）

#ifdef __cplusplus
extern “C“ {
#endif
	int yylex（void）;
	int yyback（int * int）;
	int yyinput（void）;
	int yylook（void）;
	void yyoutput（int）;
	int yyracc（int）;
	int yyreject（void）;
	void yyunput（int）;
#ifndef yyless
	void yyless（long int）;
#endif
#ifndef yywrap
	int yywrap（void）;
#endif
#ifdef __cplusplus
}
#endif

#endif

# define input（） （（（yytchar=yysptr>yysbuf?U（*--yysptr）:getc（yyin））==10?（yylineno++yytchar）:yytchar）==EOF?0:yytchar）
# define unput（c） {yytchar= （c）;if（yytchar==‘\n‘）yylineno--;*yysptr++=yytchar;}
# define yymore（） （yymorfg=1）
# define ECHO （void）fprintf（yyout “%s“yytext）
# define REJECT { nstr = yyreject（）; goto yyfussy;}
int yyleng; extern char yytext[];
int yymorfg;
extern char *yysptr yysbuf[];
int yytchar;
FILE *yyin = {stdin} *yyout = {stdout};
extern int yylineno;
struct yysvf { 
	struct yywork *yystoff;
	struct yysvf *yyother;
	int *yystops;};
struct yysvf *yyestate;
extern struct yysvf yysvec[] *yybgin;

# line 2 “lexer.l“
#include 
#include 
#include 
#include 
#include 
#include 
#include 
int line_num = 1;
# define YYNEWLINE 10
yylex（void）{
int nstr; extern int yyprevious;
while（（nstr = yylook（）） >= 0）
yyfussy: switch（nstr）{
case 0:
if（yywrap（）） return（0）; break;
case 1:

# line 21 “lexer.l“
{ fprintf（stderr“ --- Reading Nodes\n“）; return NODE;    }
break;
case 2:

# line 23 “lexer.l“
{ return ZERO;    }
break;
case 3:

# line 25 “lexer.l“
{ return READMODE; }
break;
case 4:

# line 27 “lexer.l“
{ return PRESSURE; }
break;
case 5:

# line 29 “lexer.l“
{ fprintf（stderr“ --- Reading Control\n“）; return CONTROL;  }
break;
case 6:

# line 31 “lexer.l“
{ fprintf（stderr“ --- Reading Elements\n“）; return TOPOLOGY; }
break;
case 7:

# line 33 “lexer.l“
{fprintf（stderr“ --- Reading Attributes\n“）;return ATTRIBUTES; }
break;
case 8:

# line 35 “lexer.l“
{ fprintf（stderr“ --- Reading Element frames\n“）; return EframeS; }
break;
case 9:

# line 37 “lexer.l“
 { return IDIS;  }
break;
case 10:

# line 39 “lexer.l“
{ return IDIS6; }
break;
case 11:

# line 41 “lexer.l“
{ return IVEL; }
break;
case 12:

# line 43 “lexer.l“
{ return IACC; }
break;
case 13:

# line 45 “lexer.l“
 { fprintf（stderr“ --- Done Reading Input File\n“）; return END;  }
break;
case 14:

# line 47 “lexer.l“
 { return LOAD; }
break;
case 15:

# line 49 “lexer.l“
{ return LOAD; }
break;
case 16:

# line 51 “lexer.l“
  { return USE;  }
break;
case 17:

# line 53 “lexer.l“
{ return SENSORS; }
break;
case 18:

# line 55 “lexer.l“
{ return ACTUATORS; }
break;
case 19:

# line 57 “lexer.l“
{ fprintf（stderr“ --- Reading Materials\n“）; return MATER

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件         42  2013-05-13 13:07  FEMc++\basedir

     文件       2424  2013-05-13 13:08  FEMc++\Control.d\control.C

     文件      22608  2013-05-13 13:08  FEMc++\Control.d\control.so

     文件        707  2013-05-13 13:08  FEMc++\Control.d\ControlInterface.h

     文件        180  2013-05-13 13:08  FEMc++\Control.d\Makefile

     文件        122  2013-05-13 13:08  FEMc++\Control.d\so_locations

     文件      14130  2013-05-13 13:08  FEMc++\Driver.d\Domain.C

     文件      13314  2013-05-13 13:08  FEMc++\Driver.d\Domain.h

     文件      26262  2013-05-13 13:08  FEMc++\Driver.d\Dynam.C

     文件        897  2013-05-13 13:08  FEMc++\Driver.d\Dynam.h

     文件      10183  2013-05-13 13:08  FEMc++\Driver.d\DynamProbType.C

     文件        893  2013-05-13 13:08  FEMc++\Driver.d\DynamProbType.h

     文件       4860  2013-05-13 13:08  FEMc++\Driver.d\Eigen.C

     文件       4809  2013-05-13 13:08  FEMc++\Driver.d\EigenProbType.C

     文件        933  2013-05-13 13:08  FEMc++\Driver.d\EigenProbType.h

     文件      12650  2013-05-13 13:08  FEMc++\Driver.d\jacobi.f

     文件     323044  2013-05-13 13:08  FEMc++\Driver.d\libdriver.a

     文件       1686  2013-05-13 13:08  FEMc++\Driver.d\Makefile

     文件      10314  2013-05-13 13:08  FEMc++\Driver.d\OpMake.C

     文件      44414  2013-05-13 13:08  FEMc++\Driver.d\Static.C

     文件        602  2013-05-13 13:08  FEMc++\Driver.d\StaticProbType.C

     文件        615  2013-05-13 13:08  FEMc++\Driver.d\StaticProbType.h

     文件       1682  2013-05-13 13:07  FEMc++\Element.d\AddElem.C

     文件       1600  2013-05-13 13:07  FEMc++\Element.d\beamRotation.f

     文件       3333  2013-05-13 13:08  FEMc++\Element.d\Brick.d\br8mas.f

     文件       5039  2013-05-13 13:08  FEMc++\Element.d\Brick.d\brik8v.f

     文件        555  2013-05-13 13:08  FEMc++\Element.d\Brick.d\brkcmt.f

     文件       8888  2013-05-13 13:08  FEMc++\Element.d\Brick.d\EightNodeBrick.C

     文件        930  2013-05-13 13:08  FEMc++\Element.d\Brick.d\EightNodeBrick.h

     文件       3975  2013-05-13 13:08  FEMc++\Element.d\Brick.d\h8shpe.f

............此处省略253个文件信息