发布网友 发布时间:2023-09-06 22:31
共1个回答
热心网友 时间:2023-09-11 13:23
#include "stdlib.h"#include "alloc.h"#include "math.h"#include "stdio.h"#define pi 3.1415926fstruct elems{int nedge;float x[4],y[4];float s11[4],s22[4],s12[4];};int get_s(float norg[],float ang,float pt[],float *sx,float *sy, float *txy,struct elems *elem_stress);void mnq(float center[],float radius,float orient_ang,float thick, float result[],struct elems *elem_stress);void convertxytor(float sx,float sy,float t,float ang,float *sr,float *tr);float getangle(float pt[],float p1[],float p2[]);void xytox1y1(float old[],float neworg[],float ang);int is_in_box(float pt[],float p[][2],int num);void gaussjordan(float coefary[][4],float constary[],int numcol,float solcoef[]);float secu_k(float m,float n,float h,float rr2f,float rl);int nelement=130,nnode,nsum_node;